Preparing an unsupervised massive analysis of SPHERE high contrast data with the PACO algorithm

We aim at searching for exoplanets on the whole ESO/VLT-SPHERE archive with improved and unsupervised data analysis algorithm that could allow to detect massive giant planets at 5 au. To prepare, test and optimize our approach, we gathered a sample of twenty four solar-type stars observed with SPHERE using angular and spectral differential imaging modes. We use PACO, a new generation algorithm recently developed, that has been shown to outperform classical methods. We also improve the SPHERE pre-reduction pipeline, and optimize the outputs of PACO to enhance the detection performance. We develop custom built spectral prior libraries to optimize the detection capability of the ASDI mode for both IRDIS and IFS. Compared to previous works conducted with more classical algorithms than PACO, the contrast limits we derived are more reliable and significantly better, especially at short angular separations where a gain by a factor ten is obtained between 0.2 and 0.5 arcsec. Under good observing conditions, planets down to 5 MJup, orbiting at 5 au could be detected around stars within 60 parsec. We identified two exoplanet candidates that require follow-up to test for common proper motion. In this work, we demonstrated on a small sample the benefits of PACO in terms of achievable contrast and of control of the confidence levels. Besides, we have developed custom tools to take full benefits of this algorithm and to quantity the total error budget on the estimated astrometry and photometry. This work paves the way towards an end-to-end, homogeneous, and unsupervised massive re-reduction of archival direct imaging surveys in the quest of new exoJupiters.Comment: Accepted for publication in A&

An imaged 15Mjup companion within a hierarchical quadruple system

Since 2019, the direct imaging B-star Exoplanet Abundance Study (BEAST) at SPHERE@VLT has been scanning the surroundings of young B-type stars in order to ascertain the ultimate frontiers of giant planet formation. Recently, the $17^{+3}_{-4}$ Myr HIP 81208 was found to host a close-in (~50 au) brown dwarf and a wider (~230 au) late M star around the central 2.6Msun primary. Alongside the continuation of the survey, we are undertaking a complete reanalysis of archival data aimed at improving detection performances so as to uncover additional low-mass companions. We present here a new reduction of the observations of HIP 81208 using PACO ASDI, a recent and powerful algorithm dedicated to processing high-contrast imaging datasets, as well as more classical algorithms and a dedicated PSF-subtraction approach. The combination of different techniques allowed for a reliable extraction of astrometric and photometric parameters. A previously undetected source was recovered at a short separation from the C component of the system. Proper motion analysis provided robust evidence for the gravitational bond of the object to HIP 81208 C. Orbiting C at a distance of ~20 au, this 15Mjup brown dwarf becomes the fourth object of the hierarchical HIP 81208 system. Among the several BEAST stars which are being found to host substellar companions, HIP 81208 stands out as a particularly striking system. As the first stellar binary system with substellar companions around each component ever found by direct imaging, it yields exquisite opportunities for thorough formation and dynamical follow-up studies.Comment: 12 pages, 9 figures, 5 tables. Accepted for publication as a Letter in Astronomy and Astrophysics, section 1. Letters to the Edito

A scaled-up planetary system around a supernova progenitor

Stars and planetary system

A scaled-up planetary system around a supernova progenitor

Context. Virtually all known exoplanets reside around stars with M M⊙ either due to the rapid evaporation of the protostellar disks or to selection effects impeding detections around more massive stellar hosts. Aims. To clarify if this dearth of planets is real or a selection effect, we launched the planet-hunting B-star Exoplanet Abundance STudy (BEAST) survey targeting B stars (M > 2.4 M⊙) in the young (5−20 Myr) Scorpius-Centaurus association by means of the high-contrast spectro-imager SPHERE at the Very Large Telescope. Methods. In this paper we present the analysis of high-contrast images of the massive (M ~ 9 M⊙) star μ2 Sco obtained within BEAST. We carefully examined the properties of this star, combining data from Gaia and from the literature, and used state-of-the-art algorithms for the reduction and analysis of our observations. Results. Based on kinematic information, we found that μ2 Sco is a member of a small group which we label Eastern Lower Scorpius within the Scorpius-Centaurus association. We were thus able to constrain its distance, refining in turn the precision on stellar parameters. Around this star we identify a robustly detected substellar companion (14.4 ± 0.8 MJ)at a projected separation of 290 ± 10 au, and a probable second similar object (18.5 ± 1.5 MJ) at 21 ± 1 au. The planet-to-star mass ratios of these objects are similar to that of Jupiter to the Sun, and the flux they receive from the star is similar to those of Jupiter and Mercury, respectively. Conclusions. The robust and the probable companions of μ2 Sco are naturally added to the giant 10.9 MJ planet recently discovered by BEAST around the binary b Cen system. While these objects are slightly more massive than the deuterium burning limit, their properties are similar to those of giant planets around less massive stars and they are better reproduced by assuming that they formed under a planet-like, rather than a star-like scenario. Irrespective of the (needed) confirmation of the inner companion, μ2 Sco is the first star that would end its life as a supernova that hosts such a system. The tentative high frequency of BEAST discoveries is unexpected, and it shows that systems with giant planets or small-mass brown dwarfs can form around B stars. When putting this finding in the context of core accretion and gravitational instability formation scenarios, we conclude that the current modeling of both mechanisms is not able to produce this kind of companion. The completion of BEAST will pave the way for the first time to an extension of these models to intermediate and massive stars

The SPHERE infrared survey for exoplanets (SHINE). III. The demographics of young giant exoplanets below 300 au with SPHERE

The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses between 1 and 75 MJup and semimajor axes between 5 and 300 au. We adopt detection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using the BEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections in the sample, we use a MCMC tool to compare our observations to two different types of models. The first is a parametric model based on observational constraints, and the second type are numerical models that combine advanced core accretion and gravitational instability planet population synthesis. Using the parametric model, we show that the frequencies of systems with at least one substellar companion are $23.0_{-9.7}^{+13.5}\%$, $5.8_{-2.8}^{+4.7}\%$, and $12.6_{-7.1}^{+12.9}\%$ for BA, FGK, and M stars, respectively. We also demonstrate that a planet-like formation pathway probably dominates the mass range from 1-75 MJup for companions around BA stars, while for M dwarfs, brown dwarf binaries dominate detections. In contrast, a combination of binary star-like and planet-like formation is required to best fit the observations for FGK stars. Using our population model and restricting our sample to FGK stars, we derive a frequency of $5.7_{-2.8}^{+3.8}\%$, consistent with predictions from the parametric model. More generally, the frequency values that we derive are in excellent agreement with values obtained in previous studies.Comment: 24 pages, 14 figures, 3 tables. Accepted for publication in A&

Unveiling the β Pictoris system, coupling high contrast imaging, interferometric, and radial velocity data

Context. The nearby and young β Pictoris system hosts a well resolved disk, a directly imaged massive giant planet orbiting at ≃9 au, as well as an inner planet orbiting at ≃2.7 au, which was recently detected through radial velocity (RV). As such, it offers several unique opportunities for detailed studies of planetary system formation and early evolution. Aims: We aim to further constrain the orbital and physical properties of β Pictoris b and c using a combination of high contrast imaging, long base-line interferometry, and RV data. We also predict the closest approaches or the transit times of both planets, and we constrain the presence of additional planets in the system. Methods: We obtained six additional epochs of SPHERE data, six additional epochs of GRAVITY data, and five additional epochs of RV data. We combined these various types of data in a single Markov-chain Monte Carlo analysis to constrain the orbital parameters and masses of the two planets simultaneously. The analysis takes into account the gravitational influence of both planets on the star and hence their relative astrometry. Secondly, we used the RV and high contrast imaging data to derive the probabilities of presence of additional planets throughout the disk, and we tested the impact of absolute astrometry. Results: The orbital properties of both planets are constrained with a semi-major axis of 9.8 ± 0.4 au and 2.7 ± 0.02 au for b and c, respectively, and eccentricities of 0.09 ± 0.1 and 0.27 ± 0.07, assuming the HIPPARCOS distance. We note that despite these low fitting error bars, the eccentricity of β Pictoris c might still be over-estimated. If no prior is provided on the mass of β Pictoris b, we obtain a very low value that is inconsistent with what is derived from brightness-mass models. When we set an evolutionary model motivated prior to the mass of β Pictoris b, we find a solution in the 10-11 M[SUB]Jup[/SUB] range. Conversely, β Pictoris c's mass is well constrained, at 7.8 ± 0.4 M[SUB]Jup[/SUB], assuming both planets are on coplanar orbits. These values depend on the assumptions on the distance of the β Pictoris system. The absolute astrometry HIPPARCOS-Gaia data are consistent with the solutions presented here at the 2σ level, but these solutions are fully driven by the relative astrometry plus RV data. Finally, we derive unprecedented limits on the presence of additional planets in the disk. We can now exclude the presence of planets that are more massive than about 2.5 M[SUB]Jup[/SUB] closer than 3 au, and more massive than 3.5 M[SUB]Jup[/SUB] between 3 and 7.5 au. Beyond 7.5 au, we exclude the presence of planets that are more massive than 1-2 M[SUB]Jup[/SUB]. Conclusions: Combining relative astrometry and RVs allows one to precisely constrain the orbital parameters of both planets and to give lower limits to potential additional planets throughout the disk. The mass of β Pictoris c is also well constrained, while additional RV data with appropriate observing strategies are required to properly constrain the mass of β Pictoris b.Peer reviewe

The SPHERE Infrared Survey for Exoplanets (SHINE): II. Observations, Data Reduction and Analysis, Detection Performances, and Initial Results

Context. In recent decades, direct imaging has confirmed the existence of substellar companions (exoplanets or brown dwarfs) on wide orbits (>10 au) around their host stars. In striving to understand their formation and evolution mechanisms, in 2015 we initiated the SPHERE infrared survey for exoplanets (SHINE), a systematic direct imaging survey of young, nearby stars that is targeted at exploring their demographics. Aims. We aim to detect and characterize the population of giant planets and brown dwarfs beyond the snow line around young, nearby stars. Combined with the survey completeness, our observations offer the opportunity to constrain the statistical properties (occurrence, mass and orbital distributions, dependency on the stellar mass) of these young giant planets. Methods. In this study, we present the observing and data analysis strategy, the ranking process of the detected candidates, and the survey performances for a subsample of 150 stars that are representative of the full SHINE sample. Observations were conducted in a homogeneous way between February 2015 and February 2017 with the dedicated ground-based VLT/SPHERE instrument equipped with the IFS integral field spectrograph and the IRDIS dual-band imager, covering a spectral range between 0.9 and 2.3 μm. We used coronographic, angular, and spectral differential imaging techniques to achieve the best detection performances for this study, down to the planetary mass regime. Results. We processed, in a uniform manner, more than 300 SHINE observations and datasets to assess the survey typical sensitivity as a function of the host star and of the observing conditions. The median detection performance reached 5σ-contrasts of 13 mag at 200 mas and 14.2 mag at 800 mas with the IFS (YJ and YJH bands), and of 11.8 mag at 200 mas, 13.1 mag at 800 mas, and 15.8 mag at 3 as with IRDIS in H band, delivering one of the deepest sensitivity surveys thus far for young, nearby stars. A total of sixteen substellar companions were imaged in this first part of SHINE: seven brown dwarf companions and ten planetary-mass companions.These include two new discoveries, HIP 65426 b and HIP 64892 B, but not the planets around PDS70 that had not been originally selected for the SHINE core sample. A total of 1483 candidates were detected, mainly in the large field of view that characterizes IRDIS. The color-magnitude diagrams, low-resolution spectrum (when available with IFS), and follow-up observations enabled us to identify the nature (background contaminant or comoving companion) of about 86% of our subsample. The remaining cases are often connected to crowded-field follow-up observations that were missing. Finally, even though SHINE was not initially designed for disk searches, we imaged twelve circumstellar disks, including three new detections around the HIP 73145, HIP 86598, and HD 106906 systems. Conclusions. Nowadays, direct imaging provides a unique opportunity to probe the outer part of exoplanetary systems beyond 10 au to explore planetary architectures, as highlighted by the discoveries of: one new exoplanet, one new brown dwarf companion, and three new debris disks during this early phase of SHINE. It also offers the opportunity to explore and revisit the physical and orbital properties of these young, giant planets and brown dwarf companions (relative position, photometry, and low-resolution spectrum in near-infrared, predicted masses, and contrast in order to search for additional companions). Finally, these results highlight the importance of finalizing the SHINE systematic observation of about 500 young, nearby stars for a full exploration of their outer part to explore the demographics of young giant planets beyond 10 au and to identify the most interesting systems for the next generation of high-contrast imagers on very large and extremely large telescopes. © M. Langlois et al. 2021.SPHERE is an instrument designed and built by a consortium consisting of IPAG (Grenoble, France), MPIA (Heidelberg, Germany), LAM (Marseille, France), LESIA (Paris, France), Laboratoire Lagrange (Nice, France), INAF – Osservatorio di Padova (Italy), Observatoire de Genève (Switzerland), ETH Zürich (Switzerland), NOVA (Netherlands), ONERA (France) and ASTRON (Netherlands) in collaboration with ESO. SPHERE was funded by ESO, with additional contributions from CNRS (France), MPIA (Germany), INAF (Italy), FINES (Switzerland) and NOVA (Netherlands). SPHERE also received funding from the European CommissionSixth and Seventh Framework Programmes as part of the Optical Infrared Coordination Network for Astronomy (OPTICON) under grant number RII3-Ct-2004-001566 for FP6 (2004-2008), grant number 226604 for FP7 (2009-2012) and grant number 312430 for FP7 (2013-2016). This paper is based on observations collected at the European Southern Observatory under ESO programmes 198.C-0209, 097.C-0865, 095.C-0298, 095.C-0309,096.C-0241. This work has made use of the SPHERE Data Centre, jointly operated by OSUG/IPAG (Grenoble), PYTHEAS/LAM/CeSAM (Marseille), OCA/Lagrange (Nice), Observatoire de Paris/LESIA (Paris), and Observatoire de Lyon (OSUL/CRAL). This work is supported by the French National Research Agency in the framework of the Investissements d’Avenir program (ANR-15-IDEX-02), through the funding of the “Origin of Life” project of the Univ. Grenoble-Alpes. This work is jointly supported by the French National Programms (PNP and PNPS) and by the Action Spécifique Haute Résolution Angulaire (ASHRA) of CNRS/INSU co-funded by CNES. We also thank the anonymous referee for her/his careful reading of the manuscript as well as her/his insightful comments and suggestions. AV acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 757561). A.-M.L. acknowledges funding from French National Research Agency (GIPSE project). C.P. acknowledges financial support from Fondecyt (grant 3190691) and financial support from the ICM (Iniciativa Científica Milenio) via the Núcleo Milenio de Formación Planetaria grant, from the Universidad de Valparaíso. T.H. acknowledges support from the European Research Council under the Horizon 2020 Framework Program via the ERC Advanced Grant Origins 832428

Robustness to bad frames in angular differential imaging: a local weighting approach

International audienc

REXPACO: An algorithm for high contrast reconstruction of the circumstellar environment by angular differential imaging

International audienceContext: Direct imaging is a method of choice for probing the close environment of young stars. Even with the coupling of adaptive optics and coronagraphy, the direct detection of off-axis sources such as circumstellar disks and exoplanets remains challenging due to the required high contrast and small angular resolution. Angular differential imaging (ADI) is an observational technique that introduces an angular diversity to help disentangle the signal of off-axis sources from the residual signal of the star in a post-processing step.Aims: While various detection algorithms have been proposed in the last decade to process ADI sequences and reach high contrast for the detection of point-like sources, very few methods are available to reconstruct meaningful images of extended features such as circumstellar disks. The purpose of this paper is to describe a new post-processing algorithm dedicated to the reconstruction of the spatial distribution of light (total intensity) received from off-axis sources, in particular from circumstellar disks.Methods: Built on the recent PACO algorithm dedicated to the detection of point-like sources, the proposed method is based on the local learning of patch covariances capturing the spatial fluctuations of the stellar leakages. From this statistical modeling, we develop a regularized image reconstruction algorithm (REXPACO) following an inverse problems approach based on a forward image formation model of the off-axis sources in the ADI sequences.Results: Injections of fake circumstellar disks in ADI sequences from the VLT/SPHERE-IRDIS instrument show that both the morphology and the photometry of the disks are better preserved by REXPACO compared to standard post-processing methods such as cADI. In particular, the modeling of the spatial covariances proves useful in reducing typical ADI artifacts and in better disentangling the signal of these sources from the residual stellar contamination. The application to stars hosting circumstellar disks with various morphologies confirms the ability of REXPACO to produce images of the light distribution with reduced artifacts. Finally, we show how REXPACO can be combined with PACO to disentangle the signal of circumstellar disks from the signal of candidate point-like sources.Conclusions .REXPACO is a novel post-processing algorithm for reconstructing images of the circumstellar environment from high contrast ADI sequences. It produces numerically deblurred images and exploits the spatial covariances of the stellar leakages and of the noise to efficiently eliminate this nuisance term. The processing is fully unsupervised, all tuning parameters being directly estimated from the data themselves

Exoplanet imaging data challenge: benchmarking the various image processing methods for exoplanet detection

International audienceThe Exoplanet Imaging Data Challenge is a community-wide effort meant to offer a platform for a fair and common comparison of image processing methods designed for exoplanet direct detection. For this purpose, it gathers on a dedicated repository (Zenodo), data from several high-contrast ground-based instruments worldwide in which we injected synthetic planetary signals. The data challenge is hosted on the CodaLab competition platform, where participants can upload their results. The specifications of the data challenge are published on our website https://exoplanet-imaging-challenge.github.io/. The first phase, launched on the 1st of September 2019 and closed on the 1st of October 2020, consisted in detecting point sources in two types of common data-set in the field of high-contrast imaging: data taken in pupil-tracking mode at one wavelength (subchallenge 1, also referred to as ADI) and multispectral data taken in pupil-tracking mode (subchallenge 2, also referred to as ADI+mSDI). In this paper, we describe the approach, organisational lessons-learnt and current limitations of the data challenge, as well as preliminary results of the participants' submissions for this first phase. In the future, we plan to provide permanent access to the standard library of data sets and metrics, in order to guide the validation and support the publications of innovative image processing algorithms dedicated to high-contrast imaging of planetary systems