SPOTS: The Search for Planets Orbiting Two Stars: II. First constraints on the frequency of sub-stellar companions on wide circumbinary orbits

A large number of direct imaging surveys for exoplanets have been performed in recent years, yielding the first directly imaged planets and providing constraints on the prevalence and distribution of wide planetary systems. However, like most of the radial velocity ones, these surveys generally focus on single stars, hence binaries and higher-order multiples have not been studied to the same level of scrutiny. This motivated the SPOTS (Search for Planets Orbiting Two Stars) survey, which is an ongoing direct imaging study of a large sample of close binaries, started with VLT/NACO and now continuing with VLT/SPHERE. To complement this survey, we have identified the close binary targets in 24 published direct imaging surveys. Here we present our statistical analysis of this combined body of data. We analysed a sample of 117 tight binary systems, using a combined Monte Carlo and Bayesian approach to derive the expected values of the frequency of companions, for different values of the companion's semi-major axis. Our analysis suggest that the frequency of sub-stellar companions in wide orbit is moderately low ($\lesssim$13% with a best value of 6% at 95% confidence level) and not significantly different between single stars and tight binaries. One implication of this result is that the very high frequency of circumbinary planets in wide orbits around post-common envelope binaries, implied by eclipse timing (up to 90% according to Zorotovic & Schreiber 2013), can not be uniquely due to planets formed before the common-envelope phase (first generation planets), supporting instead the second generation planet formation or a non-Keplerian origin of the timing variations.Comment: 21 pages, 3 figure

ALICE Data Release: A Revaluation of HST-NICMOS Coronagraphic Images

The Hubble Space Telescope NICMOS instrument was used from 1997 to 2008 to perform coronagraphic observations of about 400 targets. Most of them were part of surveys looking for substellar companions or resolved circumstellar disks to young nearby stars, making the NICMOS coronagraphic archive a valuable database for exoplanets and disks studies. As part of the Archival Legacy Investigations of Circumstellar Environments program, we have consistently reprocessed a large fraction of the NICMOS coronagrahic archive using advanced starlight subtraction methods. We present here the high-level science products of these re-analyzed data, which we delivered back to the community through the Mikulski Archive for Space Telescopes: doi:10.17909/T9W89V. We also present the second version of the HCI-FITS format (for High-Contrast Imaging FITS format), which we developed as a standard format for data exchange of imaging reduced science products. These re-analyzed products are openly available for population statistics studies, characterization of specific targets, or detected point-source identification

Spectral cube extraction for the VLT/SPHERE IFS: Open-source pipeline with full forward modeling and improved sensitivity

We present a new open-source data-reduction pipeline to reconstruct spectral data cubes from raw SPHERE integral-field spectrograph (IFS) data. The pipeline is written in Python and based on the pipeline that was developed for the CHARIS IFS. It introduces several improvements to SPHERE data analysis that ultimately produce significant improvements in postprocessing sensitivity. We first used new data to measure SPHERE lenslet point spread functions (PSFs) at the four laser calibration wavelengths. These lenslet PSFs enabled us to forward-model SPHERE data, to extract spectra using a least-squares fit, and to remove spectral crosstalk using the measured lenslet PSFs. Our approach also reduces the number of required interpolations, both spectral and spatial, and can preserve the original hexagonal lenslet geometry in the SPHERE IFS. In the case of least-squares extraction, no interpolation of the data is performed. We demonstrate this new pipeline on the directly imaged exoplanet 51 Eri b and on observations of the hot white dwarf companion to HD 2133. The extracted spectrum of HD 2133B matches theoretical models, demonstrating spectrophotometric calibration that is good to a few percent. Postprocessing on two 51 Eri b data sets demonstrates a median improvement in sensitivity of 80% and 30% for the 2015 and 2017 data, respectively, compared to the use of cubes reconstructed by the SPHERE Data Center. The largest improvements are seen for poorer observing conditions. The new SPHERE pipeline takes less than three minutes to produce a data cube on a modern laptop, making it practical to reprocess all SPHERE IFS data.Comment: 17 pages, 11 figures. Software available at: https://github.com/PrincetonUniversity/charis-de

The First Scattered Light Image of the Debris Disk around the Sco-Cen target HD 129590

We present the first scattered light image of the debris disk around HD 129590, a ~1.3 M$_\odot$ G1V member of the Scorpius Centaurus association with age ~10-16 Myr. The debris disk is imaged with the high contrast imaging instrument SPHERE at the Very Large Telescope, and is revealed by both the IRDIS and IFS subsytems, operating in the H and YJ bands respectively. The disk has a high infrared luminosity of $L_{\textrm{IR}}/L_{\textrm{star}}$~5$\times$10$^{-3}$, and has been resolved in other studies using ALMA. We detect a nearly edge on ring, with evidence of an inner clearing. We fit the debris disk using a model characterized by a single bright ring, with radius ~60-70 AU, in broad agreement with previous analysis of the target SED. The disk is vertically thin, and has an inclination angle of ~75$^\circ$. Along with other previously imaged edge-on disks in the Sco-Cen association such as HD 110058, HD 115600, and HD 111520, this disk image will allow of the structure and morphology of very young debris disks, shortly after the epoch of planet formation has ceased.Comment: 8 pages, 4 figures. Accepted for publication in ApJ

Early Results from VLT-SPHERE: Long-Slit Spectroscopy of 2MASS 0122-2439B, a Young Companion Near the Deuterium Burning Limit

We present 0.95-1.80 $\mu$m spectroscopy of the $\sim$12-27 $M_{\rm Jup}$ companion orbiting the faint ($R$$\sim$13.6), young ($\sim$120 Myr) M-dwarf 2MASS J01225093--2439505 ("2M0122--2439 B") at 1.5 arcsecond separation (50 AU). Our coronagraphic long-slit spectroscopy was obtained with the new high contrast imaging platform VLT-SPHERE during Science Verification. The unique long-slit capability of SPHERE enables spectral resolution an order of magnitude higher than other extreme AO exoplanet imaging instruments. With a low mass, cool temperature, and very red colors, 2M0122-2439 B occupies a particularly important region of the substellar color-magnitude diagram by bridging the warm directly imaged hot planets with late-M/early-L spectral types (e.g. $\beta$ Pic b and ROXs 42Bb) and the cooler, dusty objects near the L/T transition (e.g. HR 8799bcde and 2MASS 1207b). We fit BT-Settl atmospheric models to our $R$$\approx$350 spectrum and find $T_{\rm eff}$=1600$\pm$100 K and $\log(g)$=4.5$\pm$0.5 dex. Visual analysis of our 2M0122-2439 B spectrum suggests a spectral type L3-L4, and we resolve shallow $J$-band alkali lines, confirming its low gravity and youth. Specifically, we use the Allers & Liu (2013) spectral indices to quantitatively measure the strength of the FeH, VO, KI, spectral features, as well as the overall $H$-band shape. Using these indices, along with the visual spectral type analysis, we classify 2M0122-2439 B as an intermediate gravity (INT-G) object with spectral type L3.7$\pm$1.0.Comment: Accepted to ApJ Letters, 8 pages, 4 figures, some minor typographical issues were fixe

Constraining the presence of giant planets in two-belt debris disk systems with VLT/SPHERE direct imaging and dynamical arguments

Giant, wide-separation planets often lie in the gap between multiple, distinct rings of circumstellar debris: this is the case for the HR 8799 and HD 95086 systems, and even the solar system where the Asteroid and Kuiper belts enclose the four gas and ice giants. In the case that a debris disk, inferred from an infrared excess in the SED, is best modelled as two distinct temperatures, we infer the presence of two spatially separated rings of debris. Giant planets may well exist between these two belts of debris, and indeed could be responsible for the formation of the gap between these belts. We observe 24 such two-belt systems using the VLT/SPHERE high contrast imager, and interpret our results under the assumption that the gap is indeed formed by one or more giant planets. A theoretical minimum mass for each planet can then be calculated, based on the predicted dynamical timescales to clear debris. The typical dynamical lower limit is ˜0.2MJ in this work, and in some cases exceeds 1MJ. Direct imaging data, meanwhile, is typically sensitive to planets down to ˜3.6MJ at 1", and 1.7MJ in the best case. Together, these two limits tightly constrain the possible planetary systems present around each target, many of which will be detectable with the next generation of high-contrast imagers

ALICE Data Release: A Revaluation of HST-NICMOS Coronagraphic Images

The Hubble Space Telescope NICMOS instrument was used from 1997 to 2008 to perform coronagraphic observations of about 400 targets. Most of them were part of surveys looking for substellar companions or resolved circumstellar disks to young nearby stars, making the NICMOS coronagraphic archive a valuable database for exoplanets and disks studies. As part of the Archival Legacy Investigations of Circumstellar Environments program, we have consistently reprocessed a large fraction of the NICMOS coronagrahic archive using advanced starlight subtraction methods. We present here the high-level science products of these re-analyzed data, which we delivered back to the community through the Mikulski Archive for Space Telescopes: doi:10.17909/T9W89V. We also present the second version of the HCI-FITS format (for High-Contrast Imaging FITS format), which we developed as a standard format for data exchange of imaging reduced science products. These re-analyzed products are openly available for population statistics studies, characterization of specific targets, or detected point-source identification

Impacts of high-contrast image processing on atmospheric retrievals

Many post-processing algorithms have been developed in order to better separate the signal of a companion from the bright light of the host star, but the effect of such algorithms on the shape of exoplanet spectra extracted from integral field spectrograph data is poorly understood. The resulting spectra are affected by noise that is correlated in wavelength space due to both optical and data processing effects. Within the framework of Bayesian atmospheric retrievals, we aim to understand how these correlations and other systematic effects impact the inferred physical parameters. We consider three algorithms (KLIP, PynPoint and ANDROMEDA), optimizing the choice of algorithmic parameters using a series of injection tests into archival SPHERE and GPI data of the HR 8799 system. The wavelength-dependent covariance matrix is calculated to provide a measure of instrumental and algorithmic systematics. We perform atmospheric retrievals using petitRADTRANS on optimally extracted spectra to measure how these data processing systematics influence the retrieved parameter distributions. The choice of data processing algorithm and parameters significantly impact the accuracy of retrieval results, with the mean posterior parameter bias ranging from 1 to 3 $\sigma$ from the true input parameters. Including the full covariance matrix in the likelihood improves the accuracy of inferred parameters, and cannot be accounted for using ad hoc scaling parameters in the retrieval framework. Using the Bayesian information criterion and other statistical measures as a heuristic goodness-of-fit metrics, the retrievals including the full covariance matrix are favoured when compared to using only the diagonal elements.Comment: 22 pages, 13 figures, accepted to Astronomy & Astrophysic

Exoplanet imaging with ELTs: exploring a second-stage AO with a Zernike wavefront sensor on the ESO/GHOST testbed

We propose to explore a cascade extreme Adaptive optics (ExAO) approach with a second stage based on a Zernike wavefront sensor (ZWFS) for exoplanet imaging and spectroscopy. Most exoplanet imagers currently use a single-stage ExAO to correct for the effects of atmospheric turbulence and produce high-Strehl images of observed stars in the near-infrared. While such systems enable the observation of warm gaseous companions around nearby stars, adding a second-stage AO enables to push the wavefront correction further and possibly observe colder or smaller planets. This approach is currently investigated in different exoplanet imagers (VLT/SPHERE, Mag-AOX, Subaru/SCExAO) by considering a Pyramid wavefront sensor (PWFS) in the second arm to measure the residual atmospheric turbulence left from the first stage. Since these aberrations are expected to be very small (a few tens of nm in the near-infrared domain), we propose to investigate an alternative approach based on the ZWFS. This sensor is a promising concept with a small capture range to estimate residual wavefront errors thanks to its large sensitivity, simple phase reconstruction and easiness of implementation. In this contribution, we perform preliminary tests on the GHOST testbed at ESO to validate this approach experimentally. Additional experiments with petalling effects are also showed, giving promising wavefront correction results. Finally, we briefly discuss a first comparison between PWFS-based and ZWFS-based second-stage AO to draw preliminary conclusions on the interests of both schemes for exoplanet imaging and spectroscopy with the upgrade of the current exoplanet imagers and the envisioned ExAO instruments for ELTs.Comment: 17 pages, 10 figures, pre-print of the proceeding of the AO4ELT7 conference held in June 2023 in Avignon, Franc