Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

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.

Context. HD72946 is a bright and nearby solar-type star hosting a low-mass companion at long period (P~16 yr) detected with the radial velocities (RV) method. The companion has a minimum mass of 60.4+/-2.2 MJ and might be a brown dwarf. Its expected semi-major axis of ~243 mas makes it a suitable target for further characterization with high-contrast imaging, in particular to measure its inclination, mass, and spectrum and thus definitely establish its substellar nature. Aims. We aim to further characterize the orbit, atmosphere, and physical nature of HD72946B. Methods. We present high-contrast imaging data in the near-infrared with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. We also use proper motion measurements of the star from Hipparcos and Gaia. Results. The SPHERE data reveal a point source with a contrast of ~9 mag at a projected separation of ~235 mas. No other point sources are detected in the field of view. By jointly fitting the RV, imaging, and proper motion data, we constrain all the orbital parameters of HD72946B and assess a dynamical mass of 72.4+/-1.6 MJ and a semi-major axis of 6.45$^{+0.08}_{-0.07}$ au. Empirical comparison of its IFS spectrum to template dwarfs indicates a spectral type of L5.0+/-1.5. The J-H3 color is close to the expectations of the DUSTY models and suggests a cloudy atmosphere. Comparison with atmospheric models of the spectrophotometry suggests an effective temperature of ~1700 K. The bolometric luminosity (log(L/LS)=-4.11+/-0.10 dex) and dynamical mass of HD72946B are more compatible with evolutionary models for an age range of ~0.9-3 Gyr. The formation mechanism of the companion is currently unclear as it appears slightly away from the bulk of model predictions. HD72946B is currently the closest benchmark brown dwarf companion to a solar-type star with imaging, RV, and proper motion measurements.

Context. PDS 70 is a young (5.4 Myr), nearby (~113 pc) star hosting a known transition disk with a large gap. Recent observations with SPHERE and NACO in the near-infrared (NIR) allowed us to detect a planetary mass companion, PDS70b, within the disk cavity. Moreover, observations in H_alpha with MagAO and MUSE revealed emission associated to PDS70b and to another new companion candidate, PDS70c, at a larger separation from the star. Aims. Our aim is to confirm the discovery of the second planet PDS70c using SPHERE at VLT, to further characterize its physical properties, and search for additional point sources in this young planetary system. Methods. We re-analyzed archival SPHERE NIR observations and obtained new data in Y, J, H and K spectral bands for a total of four different epochs. The data were reduced using the data reduction and handling pipeline and the SPHERE data center. We then applied custom routines (e.g. ANDROMEDA and PACO) to subtract the starlight. Results. We re-detect both PDS 70 b and c and confirm that PDS70c is gravitationally bound to the star. We estimate this second planet to be less massive than 5 M Jup and with a T_eff around 900 K. Also, it has a low gravity with log g between 3.0 and 3.5 dex. In addition, a third object has been identified at short separation (~0.12") from the star and gravitationally bound to the star. Its spectrum is however very blue, so that we are probably seeing stellar light reflected by dust and our analysis seems to demonstrate that it is a feature of the inner disk. We, however, cannot completely exclude the possibility that it is a planetary mass object enshrouded by a dust envelope. In this latter case, its mass should be of the order of few tens of M_Earth. Moreover, we propose a possible structure for the planetary system based on our data that, however, cannot be stable on a long timescale.

The current sparse wavelength range coverage of exoplanet direct imaging observations, and the fact that models are defined using a finite wavelength range, lead both to uncertainties on effective temperature determination.We study these effects using black-bodies and atmospheric models and we detail how to infer this parameter. Through highlighting the key wavelength coverage that allows for a more accurate representation of the effective temperature, our analysis can be used to mitigate or manage extra uncertainties being added in the analysis from the models. We find that the wavelength range coverage will soon no longer be a problem. An effective temperature computed by integrating the spectroscopic observations of the James Webb Space Telescope (JWST) will give uncertainties similar to, or better than, the current state-of-the-art, which is to fit models to data. Accurately calculating the effective temperature will help to improve current modelling approaches. Obtaining an independent and precise estimation of this crucial parameter will help the benchmarking process to identify the best practice to model exoplanet atmospheres.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

High contrast imaging enables the determination of orbital parameters for substellar companions (planets, brown dwarfs) from the observed relative astrometry and the estimation of model and age-dependent masses from their observed magnitudes or spectra. Combining astrometric positions with radial velocity gives direct constraints on the orbit and on the dynamical masses of companions. A brown dwarf was discovered with the VLT/SPHERE instrument in 2017, which orbits at $\sim$ 11 au around HD 206893. Its mass was estimated between 12 and 50 $M_{Jup}$ from evolutionary models and its photometry. However, given the significant uncertainty on the age of the system and the peculiar spectrophotometric properties of the companion, this mass is not well constrained. We aim at constraining the orbit and dynamical mass of HD 206893 B. We combined radial velocity data obtained with HARPS spectra and astrometric data obtained with the high contrast imaging VLT/SPHERE and VLT/NaCo instruments, with a time baseline less than three years. We then combined those data with astrometry data obtained by Hipparcos and Gaia with a time baseline of 24 years. We used a MCMC approach to estimate the orbital parameters and dynamical mass of the brown dwarf from those data. We infer a period between 21 and 33{\deg} and an inclination in the range 20-41{\deg} from pole-on from HD 206893 B relative astrometry. The RV data show a significant RV drift over 1.6 yrs. We show that HD 206893 B cannot be the source of this observed RV drift as it would lead to a dynamical mass inconsistent with its photometry and spectra and with Hipparcos and Gaia data. An additional inner (semimajor axis in the range 1.4-2.6 au) and massive ($\sim$ 15 $M_{Jup}$) companion is needed to explain the RV drift, which is compatible with the available astrometric data of the star, as well as with the VLT/SPHERE and VLT/NaCo nondetection.

Jets are rarely associated with pre-main-sequence intermediate-mass stars. Optical and near-IR observations of jet-driving sources are often hindered by the presence of a natal envelope. Jets around partly embedded sources are a useful diagnostic to constrain the geometry of the concealed protoplanetary disk. In fact, the jet-driving mechanisms are affected by both spatial anisotropies and episodic variations at the (sub-)au scale from the star. We obtained a rich set of high-contrast VLT/SPHERE observations from 0.6 micron to 2.2 micron of the young intermediate-mass star RY Tau. Given the proximity to the Sun of this source, our images have the highest spatial resolution ever obtained for an atomic jet. Optical observations in polarized light show no sign of the protoplanetary disk detected by ALMA. Instead, we observed a diffuse signal resembling a remnant envelope with an outflow cavity. The jet is detected in four spectral lines. The jet appears to be wiggling and its radial width increasing with the distance is complementary to the shape of the outflow cavity suggesting a strong jet/envelope interaction. Through the estimated tangential velocity, we revealed a possible connection between the launching time of the jet sub-structures and the stellar activity of RY Tau. RY Tau is at an intermediate stage toward the dispersal of the natal envelope. This source shows episodic increases of mass accretion/ejection similarly to other known intermediate-mass stars. The amount of observed jet wiggle is consistent with the presence of a precessing disk warp or misaligned inner disk that would be induced by an unseen planetary/sub-stellar companion at sub-/few-au scales. The high disk mass of RY Tau and of two other jet-driving intermediate-mass stars, HD163296 and MWC480, suggests that massive, full disks are more efficient at launching prominent jets.

Current observational data of exoplanets are providing increasing detail of their 3D atmospheric structures. As characterisation efforts expand in scope, the need to develop consistent 3D radiative-transfer methods becomes more pertinent as the complex atmospheric properties of exoplanets are required to be modelled together consistently. We aim to compare the transmission and emission spectra results of a 3D Monte Carlo Radiative Transfer (MCRT) model to contemporary radiative-transfer suites. We perform several benchmarking tests of a MCRT code, Cloudy Monte Carlo Radiative Transfer (CMCRT), to transmission and emission spectra model output. We add flexibility to the model through the use of k-distribution tables as input opacities. We present a hybrid MCRT and ray tracing methodology for the calculation of transmission spectra with a multiple scattering component. CMCRT compares well to the transmission spectra benchmarks at the 10s of ppm level. Emission spectra benchmarks are consistent to within 10% of the 1D models. We suggest that differences in the benchmark results are likely caused by geometric effects between plane-parallel and spherical models. In a practical application, we post-process a cloudy 3DHD 189733b GCM model and compare to available observational data. Our results suggest the core methodology and algorithms of CMCRT produce consistent results to contemporary radiative transfer suites. 3D MCRT methods are highly suitable for detailed post-processing of cloudy and non-cloudy 1D and 3D exoplanet atmosphere simulations in instances where atmospheric inhomogeneities, significant limb effects/geometry or multiple scattering components are important considerations.

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 Mid-Infrared instrument (MIRI) on board the James Webb Space Telescope will perform the first ever characterization of young giant exoplanets observed by direct imaging in the 5-28 microns spectral range. This wavelength range is key for both determining the bolometric luminosity of the cool known exoplanets and for accessing the strongest ammonia bands. In conjunction with shorter wavelength observations, MIRI will enable a more accurate characterization of the exoplanetary atmospheric properties. Here we consider a subsample of the currently known exoplanets detected by direct imaging and we discuss their detectability with MIRI, either using the coronagraphic or the spectroscopic modes. By using the Exo-REM atmosphere model we calculate the mid-infrared emission spectra of fourteen exoplanets, and we simulate MIRI coronagraphic or spectroscopic observations. Specifically we analyze four coronagraphic observational setups, which depend on (i) the target-star and reference-star offset (0, 3, 14 mas), (ii) the wave-front-error (130, 204 nm rms), (iii) the telescope jitter amplitude (1.6, 7 mas). We then determine the signal-to-noise and integration time values for the coronagraphic targets whose planet-to-star contrasts range from 3.9 to 10.1 mag. We conclude that all the MIRI targets should be observable with different degrees of difficulty, which depends on the final in-flight instrument performances. Furthermore, we test for detection of ammonia in the atmosphere of the coolest targets. Finally, we present the case of HR 8799 b to discuss what MIRI observations can bring to the knowledge of a planetary atmosphere, either alone or in combination with shorter wavelength observations.

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.

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

The study of high contrast imaged brown dwarfs and exoplanets depends strongly on evolutionary models. To estimate the mass of a directly imaged substellar object, its extracted photometry or spectrum is used and adjusted with model spectra together with the estimated age of the system. These models still need to be properly tested and constrained. HD 4747B is a brown dwarf close to the H burning mass limit, orbiting a nearby, solar-type star and has been observed with the radial velocity method over almost two decades now. Its companion was also recently detected by direct imaging, allowing a complete study of this particular object. We aim to fully characterize HD 4747B by combining a well constrained dynamical mass and a study of its observed spectral features in order to test evolutionary models for substellar objects and characterize its atmosphere. We combine the radial velocity measurements of HIRES and CORALIE taken over two decades and high contrast imaging of several epochs from NACO, NIRC2 and SPHERE to obtain a dynamical mass. From the SPHERE data we obtain a low resolution spectrum of the companion from Y to H band, as well as two narrow band-width photometric measurements in the K band. A study of the primary star allows in addition to constrain the age of the system as well as its distance. Thanks to the new SPHERE epoch and NACO archival data combined with previous imaging data and high precision radial velocity measurements, we have been able to derive a well constrained orbit. We derive a dynamical mass of mB=70.0$\pm$1.6 MJup which is higher than a previous study, but in better agreement with the models. By comparing the object with known brown dwarfs spectra, we derive a spectral type of L9 and an effective temperature of 1350$\pm$50 K. With a retrieval analysis we constrain the oxygen and carbon abundances and compare them with the ones from the HR 8799 planets.

HR2562 is an F5V star located at ~33 pc from the Sun hosting a substellar companion that was discovered using the GPI instrument. The main objective of the present paper is to provide an extensive characterisation of the substellar companion, by deriving its fundamental properties. We observed HR 2562 with the near-infrared branch (IFS and IRDIS) of SPHERE at the VLT. During our observations IFS was operating in the YJ band, while IRDIS was observing with the $H$ broad-band filter. The data were reduced with the dedicated SPHERE GTO pipeline, which is custom-designed for this instrument. On the reduced images, we then applied the post-processing procedures that are specifically prepared to subtract the speckle noise. The companion is clearly detected in both IRDIS and IFS datasets. We obtained photometry in three different spectral bands. The comparison with template spectra allowed us to derive a spectral type of T2-T3 for the companion. Using both evolutionary and atmospheric models we inferred the main physical parameters of the companion obtaining a mass of 32+/-14$~MJup, T_{eff}=1100+/-200 K and log(g)=4.75+/-0.41.

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 developed a simple, physical and self-consistent cloud model for brown dwarfs and young giant exoplanets. We compared different parametrisations for the cloud particle size, by either fixing particle radii, or fixing the mixing efficiency (parameter fsed) or estimating particle radii from simple microphysics. The cloud scheme with simple microphysics appears as the best parametrisation by successfully reproducing the observed photometry and spectra of brown dwarfs and young giant exoplanets. In particular, it reproduces the L-T transition, due to the condensation of silicate and iron clouds below the visible/near-IR photosphere. It also reproduces the reddening observed for low-gravity objects, due to an increase of cloud optical depth for low gravity. In addition, we found that the cloud greenhouse effect shifts chemical equilibriums, increasing the abundances of species stable at high temperature. This effect should significantly contribute to the strong variation of methane abundance at the L-T transition and to the methane depletion observed on young exoplanets. Finally, we predict the existence of a continuum of brown dwarfs and exoplanets for absolute J magnitude=15-18 and J-K color=0-3, due to the evolution of the L-T transition with gravity. This self-consistent model therefore provides a general framework to understand the effects of clouds and appears well-suited for atmospheric retrievals.

Given the forthcoming launch of the James Webb Space Telescope (JWST) which will allow observ- ing exoplanet atmospheres with unprecedented signal-over-noise ratio, spectral coverage and spatial resolution, the uncertainties in the atmosphere modelling used to interpret the data need to be as- sessed. As the first step, we compare three independent 1D radiative-convective models: ATMO, Exo-REM and petitCODE. We identify differences in physical and chemical processes taken into ac- count thanks to a benchmark protocol we developed. We study the impact of these differences on the analysis of observable spectra. We show the importance of selecting carefully relevant molecular linelists to compute the atmospheric opacity. Indeed, differences between spectra calculated with Hitran and ExoMol exceed the expected uncertainties of future JWST observations. We also show the limitation in the precision of the models due to uncertainties on alkali and molecule lineshape, which induce spectral effects also larger than the expected JWST uncertainties. We compare two chemical models, Exo-REM and Venot Chemical Code, which do not lead to significant differences in the emission or transmission spectra. We discuss the observational consequences of using equilibrium or out-of- equilibrium chemistry and the major impact of phosphine, detectable with the JWST.Each of the models has benefited from the benchmarking activity and has been updated. The protocol developed in this paper and the online results can constitute a test case for other models.

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.

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.

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

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

We aim to interpret future photometric and spectral measurements from these instruments, in terms of physical parameters of the planets, with an atmospheric model using a minimal number of assumptions and parameters. We developed Exoplanet Radiative-convective Equilibrium Model (Exo-REM) to analyze the photometric and spectro- scopic data of directly imaged planets. The input parameters are a planet's surface gravity (g), effective temperature (Teff ), and elemental composition. The model predicts the equilibrium temperature profile and mixing ratio profiles of the most important gases. Opacity sources include the H2-He collision-induced absorption and molecular lines from eight compounds (including CH4 updated with the Exomol line list). Absorption by iron and silicate cloud particles is added above the expected condensation levels with a fixed scale height and a given optical depth at some reference wavelength. Scattering was not included at this stage. We applied Exo-REM to photometric and spectral observations of the planet beta Pictoris b obtained in a series of near-IR filters. We derived Teff = 1550 +- 150 K, log(g) = 3.5 +- 1, and radius R = 1.76 +- 0.24 RJup (2-{\sigma} error bars from photometric measurements). These values are comparable to those found in the literature, although with more conservative error bars, consistent with the model accuracy. We were able to reproduce, within error bars, the J- and H-band spectra of beta Pictoris b. We finally investigated the precision to which the above parameter

The intermediate-mass star Beta Pictoris is known to be surrounded by a structured edge-on debris disk within which a gas giant planet was discovered orbiting at 8-10 AU. The physical properties of Beta Pic b were previously inferred from broad and narrow-band 0.9-4.8 microns photometry. We used commissioning data of the Gemini Planet Imager (GPI) to obtain new astrometry and a low-resolution (R=35-39) J-band (1.12-1.35 microns) spectrum of the planet. We find that the planet has passed the quadrature. We constrain its semi-major axis to $\leq$ 10 AU (90 % prob.) with a peak at 8.9+0.4-0.6 AU. The joint fit of the planet astrometry and the most recent radial velocity measurements of the star yields a planet's dynamical mass $\leq$ 20 MJup (greater than 96 % prob.). The extracted spectrum of Beta Pic b is similar to those of young L1-1.5+1 dwarfs. We use the spectral type estimate to revise the planet luminosity to log(L/Lsun)=-3.90+-0.07. The 0.9-4.8 microns photometry and spectrum are reproduced for Teff=1650+-150 K and a log g lower than 4.7 dex by 12 grids of PHOENIX-based and LESIA atmospheric models. If we adopt the most recent system age estimate (21+-4 Myr), the bolometric luminosity and the constraints on the dynamical mass of Beta Pic b are only reproduced by warm- and hot-start tracks with initial entropies Si greater than 10.5 kB/baryon. Such initial conditions may result from an inefficient accretion shock and/or a planetesimal density at formation higher than in the classical core accretion model. Considering a younger age for the system or a conservative formation time for Beta Pic b does not change these conclusions.

HD 95086 is an intermediate-mass debris-disk-bearing star. VLT/NaCo $3.8 \mu m$ observations revealed it hosts a $5\pm2 \mathrm{M}_{Jup}$ companion (HD 95086 b) at $\simeq 56$ AU. Follow-up observations at 1.66 and 2.18 $\mu m$ yielded a null detection, suggesting extremely red colors for the planet and the need for deeper direct-imaging data. In this Letter, we report H- ($1.7 \mu m$) and $\mathrm{K}_1$- ($2.05 \mu m$) band detections of HD 95086 b from Gemini Planet Imager (GPI) commissioning observations taken by the GPI team. The planet position in both spectral channels is consistent with the NaCo measurements and we confirm it to be comoving. Our photometry yields colors of H-L'= $3.6\pm 1.0$ mag and K$_1$-L'=$2.4\pm 0.7$ mag, consistent with previously reported 5-$\sigma$ upper limits in H and Ks. The photometry of HD 95086 b best matches that of 2M 1207 b and HR 8799 cde. Comparing its spectral energy distribution with the BT-SETTL and LESIA planet atmospheric models yields T$_{\mathrm{eff}}\sim$600-1500 K and log g$\sim$2.1-4.5. Hot-start evolutionary models yield M=$5\pm2$ M$_{Jup}$. Warm-start models reproduce the combined absolute fluxes of the object for M=4-14 M$_{Jup}$ for a wide range of plausible initial conditions (S$_{init}$=8-13 k$_{B}$/baryon). The color-magnitude diagram location of HD 95086 b and its estimated T$_{\mathrm{eff}}$ and log g suggest that the planet is a peculiar L-T transition object with an enhanced amount of photospheric dust.

We present a radiative-convective equilibrium model for young giant exoplanets. Model predictions are compared with the existing photometric measurements of planet beta Pictoris b in the J, H, Ks, L', NB 4.05, M' bands . This model will be used to interpret future photometric and spectroscopic observations of exoplanets with SPHERE, mounted at the VLT with a first light expected mid-2014.