Context. The temperate Earth-mass planet Proxima b is the closest exoplanet to Earth and represents what may be our best ever opportunity to search for life outside the Solar System. Aims. We aim at ...directly detecting Proxima b and characterizing its atmosphere by spatially resolving the planet and obtaining high-resolution reflected-light spectra. Methods. We propose to develop a coupling interface between the SPHERE high-contrast imager and the new ESPRESSO spectrograph, both installed at ESO VLT. The angular separation of 37 mas between Proxima b and its host star requires the use of visible wavelengths to spatially resolve the planet on a 8.2-m telescope. At an estimated planet-to-star contrast of ~10-7 in reflected light, Proxima b is extremely challenging to detect with SPHERE alone. However, the combination of a ~103–104 contrast enhancement from SPHERE to the high spectral resolution of ESPRESSO can reveal the planetary spectral features and disentangle them from the stellar ones. Results. We find that significant but realistic upgrades to SPHERE and ESPRESSO would enable a 5σ detection of the planet and yield a measurement of its true mass and albedo in 20–40 nights of telescope time, assuming an Earth-like atmospheric composition. Moreover, it will be possible to probe the O2 bands at 627, 686 and 760 nm, the water vapour band at 717 nm, and the methane band at 715 nm. In particular, a 3.6σ detection of O2 could be made in about 60 nights of telescope time. Those would need to be spread over three years considering optimal observability conditions for the planet. Conclusions. The very existence of Proxima b and the SPHERE-ESPRESSO synergy represent a unique opportunity to detect biosignatures on an exoplanet in the near future. It is also a crucial pathfinder experiment for the development of extremely large telescopes and their instruments, in particular the E-ELT and its high-resolution visible and near-IR spectrograph.
Context. The direct detection of exoplanets with high-contrast imaging requires advanced data processing methods to disentangle potential planetary signals from bright quasi-static speckles. Among ...them, angular differential imaging (ADI) permits potential planetary signals with a known rotation rate to be separated from instrumental speckles that are either statics or slowly variable. The method presented in this paper, called ANDROMEDA for ANgular Differential OptiMal Exoplanet Detection Algorithm, is based on a maximum likelihood approach to ADI and is used to estimate the position and the flux of any point source present in the field of view. Aims. In order to optimize and experimentally validate this previously proposed method, we applied ANDROMEDA to real VLT/NaCo data. In addition to its pure detection capability, we investigated the possibility of defining simple and efficient criteria for automatic point source extraction able to support the processing of large surveys. Methods. To assess the performance of the method, we applied ANDROMEDA on VLT/NaCo data of TYC-8979-1683-1 which is surrounded by numerous bright stars and on which we added synthetic planets of known position and flux in the field. In order to accommodate the real data properties, it was necessary to develop additional pre-processing and post-processing steps to the initially proposed algorithm. We then investigated its skill in the challenging case of a well-known target, β Pictoris, whose companion is close to the detection limit and we compared our results to those obtained by another method based on principal component analysis (PCA). Results. Application on VLT/NaCo data demonstrates the ability of ANDROMEDA to automatically detect and characterize point sources present in the image field. We end up with a robust method bringing consistent results with a sensitivity similar to the recently published algorithms, with only two parameters to be fine tuned. Moreover, the companion flux estimates are not biased by the algorithm parameters and do not require a posteriori corrections. Conclusions. ANDROMEDA is an attractive alternative to current standard image processing methods that can be readily applied to on-sky data.
Context. 51 Eridani b is an exoplanet around a young (20 Myr) nearby (29.4 pc) F0-type star, which was recently discovered by direct imaging. It is one of the closest direct imaging planets in ...angular and physical separation (~0.5′′, ~13 au) and is well suited for spectroscopic analysis using integral field spectrographs. Aims. We aim to refine the atmospheric properties of the known giant planet and to constrain the architecture of the system further by searching for additional companions. Methods. We used the extreme adaptive optics instrument SPHERE at the Very Large Telescope (VLT) to obtain simultaneous dual-band imaging with IRDIS and integral field spectra with IFS, extending the spectral coverage of the planet to the complete Y- to H-band range and providing additional photometry in the K12-bands (2.11, 2.25 μm). The object is compared to other known cool and peculiar dwarfs. The posterior probability distributions for parameters of cloudy and clear atmospheric models are explored using MCMC. We verified our methods by determining atmospheric parameters for the two benchmark brown dwarfs Gl 570D and HD 3651B. We used archival VLT-NACO (L′) Sparse Aperture Masking data to probe the innermost region for additional companions. Results. We present the first spectrophotometric measurements in the Y and K bands for the planet and revise its J-band flux to values 40% fainter than previous measurements. Cloudy models with uniform cloud coverage provide a good match to the data. We derive the temperature, radius, surface gravity, metallicity, and cloud sedimentation parameter fsed. We find that the atmosphere is highly super-solar (Fe/H = 1.0 ± 0.1 dex), and the low \hbox{${f_{\rm sed} = 1.26^{+0.36}_{-0.29}}$}fsed=1.26-0.29+0.36 value is indicative of a vertically extended, optically thick cloud cover with small sized particles. The model radius and surface gravity estimates suggest higher planetary masses of \hbox{${M_\mathrm{gravity} = 9.1^{+4.9}_{-3.3} \, {M}_\mathrm{J}}$}Mgravity=9.1-3.3+4.9 MJ. The evolutionary model only provides a lower mass limit of > 2 MJ (for pure hot-start). The cold-start model cannot explain the luminosity of the planet. The SPHERE and NACO/SAM detection limits probe the 51 Eri system at solar system scales and exclude brown-dwarf companions more massive than 20 MJ beyond separations of ~2.5 au and giant planets more massive than 2 MJ beyond 9 au.
In recent years, there has been intensive research into the direct detection of exoplanets. Data obtained in the future with high-contrast imaging instruments, optimized for the direct detection of ...giant planets, may be strongly limited by speckle noise. Specific observing strategies and data analysis methods, such as angular and spectral differential imaging, are required to attenuate the noise level and possibly to detect the flux of faint planets. Even though these methods are very efficient at suppressing the speckles, the photometry of faint planets is dominated by the speckle residuals. The determination of the effective temperature and surface gravity of the detected planets from photometric measurements in different bands is then limited by the photometric error on the planet flux. In this paper, we investigate this photometric error and the consequences on the determination of the physical parameters of the detected planets. We perform detailed end-to-end simulation with the caos-based software package for spectro-polarimetric high-contrast exoplanet research (SPHERE) to obtain realistic data representing typical observing sequences in the Y, J, H and Ks bands with a high-contrast imager. The simulated data are used to measure the photometric accuracy as a function of contrast for planets detected with angular and spectral+angular differential methods. We apply this empirical accuracy to study the characterization capabilities of a high-contrast differential imager. We show that the expected photometric performances will allow the detection and characterization of exoplanets down to a Jupiter mass at angular separations of 1.0 and 0.2 arcsec, respectively, around high-mass and low-mass stars with two observations in different filter pairs. We also show that the determination of the physical parameters of the planets from photometric measurements in different filter pairs is essentially limited by the error on the determination of the surface gravity.
Here, we show that the approximately 10-million-year-old β Pictoris system hosts a massive giant planet, β Pictoris b, located 8 to 15 astronomical units from the star. This result confirms that gas ...giant planets form rapidly within disks and validates the use of disk structures as fingerprints of embedded planets. Among the few planets already imaged, β Pictoris b is the closest to its parent star. Its short period could allow for recording of the full orbit within 17 years.
Understanding the formation and evolution of giant planets (≥1 MJup) at wide orbital separation (≥5 AU) is one of the goals of direct imaging. Over the past 15 yr, many surveys have placed strong ...constraints on the occurrence rate of wide-orbit giants, mostly based on non-detections, but very few have tried to make a direct link with planet formation theories. In the present work, we combine the results of our previously published VLT/NaCo large program with the results of 12 past imaging surveys to constitute a statistical sample of 199 FGK stars within 100 pc, including three stars with sub-stellar companions. Using Monte Carlo simulations and assuming linear flat distributions for the mass and semi-major axis of planets, we estimate the sub-stellar companion frequency to be within 0.75–5.70% at the 68% confidence level (CL) within 20–300 AU and 0.5–75 MJup, which is compatible with previously published results. We also compare our results with the predictions of state-of-the-art population synthesis models based on the gravitational instability (GI) formation scenario with and without scattering. We estimate that in both the scattered and non-scattered populations, we would be able to detect more than 30% of companions in the 1–75 MJup range (95% CL). With the threesub-stellar detections in our sample, we estimate the fraction of stars that host a planetary system formed by GI to be within 1.0–8.6% (95% CL). We also conclude that even though GI is not common, it predicts a mass distribution of wide-orbit massive companions that is much closer to what is observed than what the core accretion scenario predicts. Finally, we associate the present paper with the release of the Direct Imaging Virtual Archive (DIVA), a public database that aims at gathering the results of past, present, and future direct imaging surveys.
Context. Direct imaging of circumstellar disks requires high-contrast and high-resolution techniques. The angular differential imaging (ADI) technique is one of them, initially developed for ...point-like sources but now increasingly applied to extended objects such as disks. This new field of application raises many questions because the disk images reduced with ADI depend strongly on the amplitude of field rotation and the ADI data reduction strategy. Both of them directly affect the disk observable properties. Aims. Our aim is to characterize the applicability and biases of some ADI data reduction strategies for different disk morphologies. A particular emphasis is placed on parameters mostly used for disks such as their surface brightness distribution, their width if the disk is a ring, and local features such as gaps or asymmetries. We first present a general method for predicting and quantifying those biases. In a second step we illustrate them for some widely used ADI algorithms applied to typical debris disk morphologies: inclined rings with various inner/outer slopes and width. Last, our aim is also to propose improvements of classical ADI to limit the biases on extended objects. Methods. Simulated fake disks seen under various observing conditions were used to reduce ADI data and quantify the resulting biases. These conclusions are complemented by previous results from NaCo L’ real-disk images of HR 4796A. Results. As expected, ADI induces flux losses on disks. This makes this technique appropriate only for low- to medium-inclination disks. A theoretical criterion is derived to predict the amount of flux loss for a given disk morphology, and quantitative estimates of the biases are given in some specific configurations. These biases alter the disk observable properties, such as the slopes of the disk surface brightness or the radial/azimuthal extent of the disk. Additionally, this work demonstrates that ADI can very easily create artificial features without involving astrophysical processes. For example, a particularly striking feature appears for a ring when the amplitude of field rotation is too small. The two ring ansae are surrounded by two flux-depleted regions, which makes them appear as bright blobs. This observation does not require any astrophysical process such as dust blown by radiation pressure, as previously proposed in H-band images of HR 4796A. Conclusions. The ADI techniques behave as spatial filtering algorithms and can bias disk observables. Therefore, the filtering process needs to be properly calibrated when deriving disk parameters from processed images.
Context. Understanding the diversity of planets requires studying the morphology and physical conditions in the protoplanetary disks in which they form. Aims. We aim to study the structure of the ~10 ...Myr old protoplanetary disk HD 100453, to detect features that can trace disk evolution and to understand the mechanisms that drive these features. Methods. We observed HD 100453 in polarized scattered light with VLT/SPHERE at optical (0.6 μm, 0.8 μm) and near-infrared (1.2 μm) wavelengths, reaching an angular resolution of ~0.02′′, and an inner working angle of ~0.09′′. Results. We spatially resolve the disk around HD 100453, and detect polarized scattered light up to ~0.42′′ (~48 au). We detect a cavity, a rim with azimuthal brightness variations at an inclination of ~38° with respect to our line of sight, 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 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°. 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. Conclusions. 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. However, the origin of such a misalignment in HD 100453, and of the spirals, is still unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit. Observations over a few years will allow us to measure the spiral pattern speed, and determine if the shadows are fixed or moving, which may constrain their origin.
Context. The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced ...stage of disk clearing and may harbor direct signatures of disk evolution. Aims. We aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. Methods. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with VLT/SPHERE at 1.04 μm to resolve scattered light at a smaller inner working angle (0.093′′) and a higher angular resolution (0.027′′) than previously achieved. Results. We observe polarized scattered light within 0.53′′ (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r ~ 0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Conclusions. Our observations reveal the complex morphology of the disk MWC 758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.
Context. Debris disks are observed around 10 to 20% of FGK main-sequence stars as infrared excess emission. They are important signposts for the presence of colliding planetesimals and therefore ...provide important information about the evolution of planetary systems. Direct imaging of such disks reveals their geometric structure and constrains their dust-particle properties. Aims. We present observations of the known edge-on debris disk around HIP 79977 (HD 146897) taken with the ZIMPOL differential polarimeter of the SPHERE instrument. We measure the observed polarization signal and investigate the diagnostic potential of such data with model simulations. Methods. SPHERE-ZIMPOL polarimetric data of the 15 Myr-old F star HIP 79977 (Upper Sco, 123 pc) were taken in the Very Broad Band (VBB) filter (λc = 735 nm, Δλ = 290 nm) with a spatial resolution of about 25 mas. Imaging polarimetry efficiently suppresses the residual speckle noise from the AO system and provides a differential signal with relatively small systematic measuring uncertainties. We measure the polarization flux along and perpendicular to the disk spine of the highly inclined disk for projected separations between 0.2′′ (25 AU) and 1.6′′ (200 AU). We perform model calculations for the polarized flux of an optically thin debris disk which are used to determine or constrain the disk parameters of HIP 79977. Results. We measure a polarized flux contrast ratio for the disk of (Fpol)disk/F∗ = (5.5 ± 0.9) × 10-4 in the VBB filter. The surface brightness of the polarized flux reaches a maximum of SBmax = 16.2 mag arcsec-2 at a separation of 0.2′′–0.5′′ along the disk spine with a maximum surface brightness contrast of 7.64 mag arcsec-2. The polarized flux has a minimum near the star <0.2′′ because no or only little polarization is produced by forward or backward scattering in the disk section lying in front of or behind the star. The width of the disk perpendicular to the spine shows a systematic increase in FWHM from 0.1′′ (12 AU) to 0.3′′−0.5′′, when going from a separation of 0.2′′ to >1′′. This can be explained by a radial blow-out of small grains. The data are modelled as a circular dust belt with a well defined disk inclination i = 85( ± 1.5)° and a radius between r0 = 60 and 90 AU. The radial density dependence is described by (r/r0)α with a steep (positive) power law index α = 5 inside r0 and a more shallow (negative) index α = −2.5 outside r0. The scattering asymmetry factor lies between g = 0.2 and 0.6 (forward scattering) adopting a scattering-angle dependence for the fractional polarization such as that for Rayleigh scattering. Conclusions. Polarimetric imaging with SPHERE-ZIMPOL of the edge-on debris disk around HIP 79977 provides accurate profiles for the polarized flux. Our data are qualitatively very similar to the case of AU Mic and they confirm that edge-on debris disks have a polarization minimum at a position near the star and a maximum near the projected separation of the main debris belt. The comparison of the polarized flux contrast ratio (Fpol)disk/F∗ with the fractional infrared excess provides strong constraints on the scattering albedo of the dust.