K-Stacker: an algorithm to hack the orbital parameters of planets hidden in high-contrast imaging. First applications to VLT/SPHERE multi-epoch observations

Context. Recent high-contrast imaging surveys, using the Spectro- Polarimetic High contrast imager for Exoplanets REsearch (SPHERE) or the Gemini Planet Imager in search of planets in young, nearby systems, have shown evidence of a small number of giant planets at relatively large separation beyond 10-30 au, where those surveys are the most sensitive. Access to smaller physical separations between 5 and 30 au is the next step for future planet imagers on 10 m telescopes and the next generation of extremely large telescopes in order to bridge the gap with indirect techniques such as radial velocity, transit, and soon astrometry with Gaia. In addition to new technologies and instruments, the development of innovative observing strategies combined with optimized data processing tools is participating in the improvement of detection capabilities at very close angular separation. In that context, we recently proposed a new algorithm, Keplerian-Stacker, which combines multiple observations acquired at different epochs and takes into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. We showed that this algorithm is able to find planets in time series of simulated images of the SPHERE InfraRed Dual-band Imager and Spectrograph (IRDIS) even when a planet remains undetected at one epoch. Aims: Our goal is to test and validate the K-Stacker algorithm performances on real SPHERE datasets to demonstrate the resilience of this algorithm to instrumental speckles and the gain offered in terms of true detection. This will motivate future dedicated multi-epoch observation campaigns of well- chosen, young, nearby systems and very nearby stars carefully selected to search for planets in emitted and reflected light, respectively, to open a new path concerning the observing strategy used with current and future planet imagers. Methods: To test K-Stacker, we injected fake planets and scanned the low signal-to-noise ratio (S/N) regime in a series of raw observations obtained by the SPHERE/IRDIS instrument in the course of the SPHERE High-contrast ImagiNg survey for Exoplanets. We also considered the cases of two specific targets intensively monitored during this campaign: β Pictoris and HD 95086. For each target and epoch, the data were reduced using standard angular differential imaging processing techniques and then recombined with K-Stacker to recover the fake planetary signals. In addition, the known exoplanets β Pictoris b and HD 95086 b previously identified at lower S/N in single epochs have also been recovered by K-Stacker. Results: We show that K-Stacker achieves a high success rate of ≈100% when the S/N of the planet in the stacked image reaches ≈9. The improvement of the S/N is given as the square root of the total exposure time contained in the data being combined. At S/N Based on observations collected at the European Southern Observatory under programs: 095.C-0298, 096.C-0241, 097.C-0865, 198.C-0209, 099.C-0127

In-depth direct imaging and spectroscopic characterization of the young Solar System analog HD 95086

Context. HD 95086 is a young nearby Solar System analog hosting a giant exoplanet orbiting at 57 au from the star between an inner and outer debris belt. The existence of additional planets has been suggested as the mechanism that maintains the broad cavity between the two belts. Aims. We present a dedicated monitoring of HD 95086 with the VLT/SPHERE instrument to refine the orbital and atmospheric properties of HD 95086 b, and to search for additional planets in this system. Methods. SPHERE observations, spread over ten epochs from 2015 to 2019 and including five new datasets, were used. Combined with archival observations, from VLT/NaCo (2012-2013) and Gemini/GPI (2013-2016), the extended set of astrometric measurements allowed us to refine the orbital properties of HD 95086 b. We also investigated the spectral properties and the presence of a circumplanetary disk around HD 95086 b by using the special fitting tool exploring the diversity of several atmospheric models. In addition, we improved our detection limits in order to search for a putative planet c via the K-Stacker algorithm. Results. We extracted for the first time the JH low-resolution spectrum of HD 95086 b by stacking the six best epochs, and confirm its very red spectral energy distribution. Combined with additional datasets from GPI and NaCo, our analysis indicates that this very red color can be explained by the presence of a circumplanetary disk around planet b, with a range of high-temperature solutions (1400–1600 K) and significant extinction (Av ≳ 10 mag), or by a super-solar metallicity atmosphere with lower temperatures (800–300 K), and small to medium amount of extinction (Av ≲ 10 mag). We do not find any robust candidates for planet c, but give updated constraints on its potential mass and location

In-depth direct imaging and spectroscopic characterization of the young Solar System analog HD 95086

Context. HD 95086 is a young nearby Solar System analog hosting a giant exoplanet orbiting at 57 au from the star between an inner and outer debris belt. The existence of additional planets has been suggested as the mechanism that maintains the broad cavity between the two belts. Aims. We present a dedicated monitoring of HD 95086 with the VLT/SPHERE instrument to refine the orbital and atmospheric properties of HD 95086 b, and to search for additional planets in this system. Methods. SPHERE observations, spread over ten epochs from 2015 to 2019 and including five new datasets, were used. Combined with archival observations, from VLT/NaCo (2012-2013) and Gemini/GPI (2013-2016), the extended set of astrometric measurements allowed us to refine the orbital properties of HD 95086 b. We also investigated the spectral properties and the presence of a circumplanetary disk around HD 95086 b by using the special fitting tool exploring the diversity of several atmospheric models. In addition, we improved our detection limits in order to search for a putative planet c via the K-Stacker algorithm. Results. We extracted for the first time the JH low-resolution spectrum of HD 95086 b by stacking the six best epochs, and confirm its very red spectral energy distribution. Combined with additional datasets from GPI and NaCo, our analysis indicates that this very red color can be explained by the presence of a circumplanetary disk around planet b, with a range of high-temperature solutions (1400–1600 K) and significant extinction (Av ≳ 10 mag), or by a super-solar metallicity atmosphere with lower temperatures (800–300 K), and small to medium amount of extinction (Av ≲ 10 mag). We do not find any robust candidates for planet c, but give updated constraints on its potential mass and location.ISSN:0004-6361ISSN:1432-074