Spectral Differential Imaging Research Articles

REXPACO ASDI: joint unmixing and deconvolution of the circumstellar environment by angular and spectral differential imaging

ABSTRACT Angular and spectral differential imaging is an observational technique of choice to investigate the immediate vicinity of stars. By leveraging the relative angular motion and spectral scaling between on-axis and off-axis sources, post-processing techniques can separate residual star light from light emitted by surrounding objects such as circumstellar discs or point-like objects. This paper introduces a new algorithm that jointly unmixes these components and deconvolves disc images. The proposed algorithm is based on a statistical model of the residual star light, accounting for its spatial and spectral correlations. These correlations are crucial yet remain inadequately modelled by existing reconstruction algorithms. We employ dedicated shrinkage techniques to estimate the large number of parameters of our correlation model in a data-driven fashion. We show that the resulting separable model of the spatial and spectral covariances captures very accurately the star light, enabling its efficient suppression. We apply our method to data sets from the Very Large Telescope/Spectro-Polarimetry High-contrast Exoplanet REsearch instrument and compare its performance with standard algorithms (median subtraction, PCA, PACO). We demonstrate that considering the multiple correlations within the data significantly improves reconstruction quality, resulting in better preservation of both disc morphology and photometry. With its unique joint spectral modelling, the proposed algorithm can reconstruct discs with circular symmetry (e.g. rings, spirals) at intensities one million times fainter than the star, without needing additional reference data sets free from off-axis objects.

Preparation for an unsupervised massive analysis of SPHERE high-contrast data with PACO

Context. Despite tremendous progress in the detection and characterization of extrasolar planetary systems in the last 25 yr, we have not pinpointed any Solar System analogues. In particular, Jupiter-like planets (either mature or old) are barely detectable beyond 5 au with indirect techniques and they are still out of the reach of direct imaging techniques. Aims. Our study is aimed at a search for exoplanets throughout the whole ESO/VLT-SPHERE archive with an improved and unsupervised data analysis algorithm that could allow us to detect massive giant planets at 5 au. To prepare, test, and optimize our approach, we gathered a sample of 24 solar-type stars observed with SPHERE using angular and spectral differential imaging modes. Methods. We used PACO, a recently developed new-generation algorithm that has been shown to outperform classical methods. We also improved the SPHERE pre-reduction pipeline and optimized the outputs of PACO to enhance the detection performance. We developed custom-built spectral prior libraries to optimize the detection capability of the ASDI mode for both IRDIS and IFS. Results. Compared to previous works conducted with more classical algorithms, the contrast limits we derived with PACO are more reliable and significantly improved, especially at short angular separations, where a gain by a factor ten has been 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 will require a follow-up to test for a common proper motion. Conclusions. In this work, we use a small sample to demonstrate the benefits of PACO in terms of achievable contrast and of control of the confidence levels. In addition, we have developed custom tools to take full advantage 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 for new exo-Jupiters.

Spectrally dispersed kernel phase interferometry with SCExAO/CHARIS: proof of concept and calibration strategies

Kernel phase interferometry (KPI) is a data processing technique that allows for the detection of asymmetries (such as companions or disks) in high-Strehl images, close to and within the classical diffraction limit. We show that KPI can successfully be applied to hyperspectral image cubes generated from integral field spectrographs (IFSs). We demonstrate this technique of spectrally-dispersed kernel phase by recovering a known binary with the SCExAO/CHARIS IFS in high-resolution K-band mode. We also explore a spectral differential imaging (SDI) calibration strategy that takes advantage of the information available in images from multiple wavelength bins. Such calibrations have the potential to mitigate high-order, residual systematic kernel phase errors, which currently limit the achievable contrast of KPI. The SDI calibration presented here is applicable to searches for line emission or sharp absorption features, and is a promising avenue toward achieving photon-noise-limited kernel phase observations. The high angular resolution and spectral coverage provided by dispersed kernel phase offers novel opportunities for science observations which would have been challenging to achieve otherwise.

AF Lep b: The lowest-mass planet detected by coupling astrometric and direct imaging data

Aims. Using the direct-imaging technique, we searched for low-mass companions around the star AF Lep, which presents a significant proper-motion anomaly (PMa) signal obtained from the comparison of HIPPARCOS and Gaia eDR3 catalogs. Methods. We observed AF Lep in two epochs with VLT/SPHERE using its subsystems IFS and IRDIS in the near-infrared, covering wavelengths ranging from the Y to the K spectral bands (between 0.95 and 2.3 μm). We then reduced the data using the high-contrast imaging techniques angular differential imaging (ADI) and spectral differential imaging in order to be able to retrieve the signal from low-mass companions of the star. Results. A faint companion was retrieved at a separation of ~0.335″ from the star and with a position angle of ~70.5° in the first epoch and with a similar position in the second epoch. This corresponds to a projected separation of ~9 au. The extracted photometry allowed us to estimate a mass for the companion of between 2 and 5 MJup. This mass is in good agreement with astrometric measurements of the dynamic mass of the companion, which give 5.2–5.5 MJup. This is the first companion with a mass well below the deuterium burning limit that was discovered by coupling direct imaging with PMa measurements. Orbital fitting done using the orvara tool allowed us to further confirm the companion mass and to define its main orbital parameters.

Post-processing CHARIS integral field spectrograph data with pyklip

Abstract We present the pyKLIP-CHARIS post-processing pipeline, a Python library that reduces high contrast imaging data for the CHARIS integral field spectrograph used with the SCExAO project on the Subaru Telescope. The pipeline is a part of the pyklip package, a Python library dedicated to the reduction of direct imaging data of exoplanets, brown dwarfs, and discs. For PSF subtraction, the pyKLIP-CHARIS post-processing pipeline relies on the core algorithms implemented in pyklip but uses image registration and calibrations that are unique to CHARIS. We describe the pipeline procedures, calibration results, and capabilities in processing imaging data acquired via the angular differential imaging and spectral differential imaging observing techniques. We showcase its performance on extracting spectra of injected synthetic point sources as well as compare the extracted spectra from real data sets on HD 33632 and HR 8799 to results in the literature. The pipeline is a python-based complement to the SCExAO project supported, widely used (and currently IDL-based) CHARIS data post-processing pipeline (CHARIS DPP) and provides an additional approach to reducing CHARIS data and extracting calibrated planet spectra.

Searching for Hα-emitting sources in the gaps of five transitional disks

Context. (Pre-)Transitional disks show gaps and cavities that can be related to ongoing planet formation. According to theory, young embedded planets can accrete material from the circumplanetary and circumstellar disks and can be detected using accretion tracers, such as the Hα emission line. Aims. We aim to detect accreting protoplanets within the cavities of five (pre-)transitional disks through adaptive-optics(AO)-assisted spectral angular differential imaging in the optical regime. Methods. We performed simultaneous AO observations in the Hα line and the adjacent continuum using the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) with the Zurich Imaging Polarimeter (ZIMPOL) at the Very Large Telescope (VLT). We combined spectral and angular differential imaging techniques to increase the contrast in the innermost regions close to the star and search for the signature of young accreting protoplanets. Results. The reduced images show no clear Hα point source around any of the targets. We report the presence of faint Hα emission around TW Hya and HD163296: while the former is most probably an artifact related to a spike, the nature of the latter remains unclear. The spectral and angular differential images yield contrasts of 6–8 magnitudes at ~100 mas from the central stars, except in the case of LkCa15, with values of ~3 mag. We used the contrast curves to estimate average upper limits to the Hα line luminosity of LHα ~ 5 × 10-6 L⊙ at separations ≥200 mas for TW Hya, RXJ1615, and T Cha, while for HD163296 and LkCa15 we derive values of ~3 × 10−5 L⊙. We estimated upper limits to the accretion luminosity of potential protoplanets, obtaining that planetary models provide an average value of Lacc ~ 10−4 L⊙ at 200 mas, which is about two orders of magnitude higher than the Lacc estimated from the extrapolation of the LHα - Lacc stellar relationship. Conclusions. When considering all the objects observed with SPHERE/ZIMPOL in the Hα line, 5 in this work and 13 from the literature, we can explain the lack of protoplanet detections by a combination of factors, such as a majority of low-mass, low-accreting planets; potential episodic accretion; significant extinction from the circumstellar and circumplanetary disks; and the fact that the contrast is less favorable at separations of smaller than 100 mas, where giant planets are more likely to form.

Orbital and dynamical analysis of the system around HR 8799

Context. HR 8799 is a young planetary system composed of four planets and a double debris belt. Being the first multi-planetary system discovered with the direct imaging technique, it has been observed extensively since 1998. This wide baseline of astrometric measurements, counting over 50 observations in 20 years, permits a detailed orbital and dynamical analysis of the system. Aims. To explore the orbital parameters of the planets, their dynamical history, and the planet-to-disk interaction, we made follow-up observations of the system during the VLT/SPHERE guaranteed time observation program. We obtained 21 observations, most of them in favorable conditions. In addition, we observed HR 8799 with the instrument LUCI at the Large Binocular Telescope (LBT). Methods. All the observations were reduced with state-of-the-art algorithms implemented to apply the spectral and angular differential imaging method. We re-reduced the SPHERE data obtained during the commissioning of the instrument and in three open-time programs to have homogeneous astrometry. The precise position of the four planets with respect to the host star was calculated by exploiting the fake negative companions method. We obtained an astrometric precision of the order of 6 mas in the worst case and 1 mas in the best case. To improve the orbital fitting, we also took into account all of the astrometric data available in the literature. From the photometric measurements obtained in different wavelengths, we estimated the masses of the planets following the evolutionary models. Results. We obtained updated parameters for the orbits with the assumption of coplanarity, relatively small eccentricities, and periods very close to the 2:1 resonance. We also refined the dynamical mass of each planet and the parallax of the system (24.49 ± 0.07 mas), which overlap with the recent Gaia eDR3/DR3 estimate. Hydrodynamical simulations suggest that inward migration of the planets caused by the interaction with the disk might be responsible for the planets being locked in resonance. We also conducted detailed N-body simulations indicating possible positions of a putative fifth planet with a mass below the present detection limits of ≃3 MJup.

Carina High-contrast Imaging Project for massive Stars (CHIPS)

Context.Most massive stars belong to multiple systems, yet the formation process leading to such high multiplicity remains insufficiently understood. To help constrain the different formation scenarios that exist, insights into the low-mass end of the companion mass function of such stars is crucial. However, this is a challenging endeavour as (sub-)solar mass companions at angular separations (ρ) below 1″ (corresponding to 1000–3000 au in nearby young open clusters and OB associations) are difficult to detect due to the large brightness contrast with the central star.Aims.With the Carina High-contrast Imaging Project of massive Stars (CHIPS), we aim to obtain statistically significant constraints on the presence and properties of low-mass companions around massive stars in a previously unreachable observing window (Δmag ≳ 10 atρ ≲ 1″). In the second paper of the series, we focus on the Trumpler 14 cluster, which harbours some of the youngest and most massive O-type stars in the Milky Way.Methods.We obtained VLT-SPHERE observations of seven O-type objects in Trumpler 14 using IRDIFS_EXT mode. These provide us with a 12″ × 12″ field of view (approximately ((3 × 3)×104au) centred on each O star and allow us to search for companions at separations larger than 0″​​.15 (approx. 360 au) and down to magnitude contrast > 10 mag in the near-infrared. We used angular and spectral differential imaging along with Point Spread Function (PSF) fitting to detect sources and measure their flux relative to that of the central object. We then used grids of ATLAS9 and PHOENIX Local Thermodynamic Equilibrium (LTE) atmosphere models combined with (pre-)main-sequence evolutionary tracks to estimate the mass of the detected candidate companions.Results.We detected 211 sources with near-infrared magnitude contrast in the range of 2–12. Given the large surface number density of stars in Trumpler 14, one cannot reliably distinguish between cluster members and genuine companions for most of the detected sources. The closest companion, at only 0″​​.26, is characterised as a 1.4M⊙star with an age of 0.6 Myr, in excellent agreement with previous age estimates for Tr 14. The mass function peaks at about 0.4M⊙and presents a dearth of stars in the 0.5–0.8M⊙mass range compared to previous estimates of the initial mass function in Tr 14. While statistically significant, part of these differences may result from contamination of theK-band fluxes by circumstellar material.Conclusions.SPHERE is clearly suitable to probe the low-mass end of the mass function in the vicinity of massive stars. Follow-up SPHERE observations to obtain the full Y to K spectral energy distribution would allow for better constraints on the masses of the detected sources, and to confirm (or invalidate) the curious mass function that we derived for low-mass stars in the vicinity of the O-type objects in Trumpler 14.

Signs of late infall and possible planet formation around DR Tau using VLT/SPHERE and LBTI/LMIRCam

Context. Protoplanetary disks around young stars often contain substructures like rings, gaps, and spirals that could be caused by interactions between the disk and forming planets. Aims. We aim to study the young (1–3 Myr) star DR Tau in the near-infrared and characterize its disk, which was previously resolved through submillimeter interferometry with ALMA, and to search for possible substellar companions embedded into it. Methods. We observed DR Tau with VLT/SPHERE both in polarized light (H broad band) and total intensity (in Y, J, H, and K spectral bands). We also performed L′ band observations with LBTI/LMIRCam on the Large Binocular Telescope (LBT). We applied differential imaging techniques to analyze both the polarized data, using dual beam polarization imaging, and the total intensity data, using angular and spectral differential imaging. Results. We found two previously undetected spirals extending north-east and south of the star, respectively. We further detected an arc-like structure north of the star. Finally a bright, compact and elongated structure was detected at a separation of 303 ± 10 mas and a position angle 21.2 ± 3.7 degrees, just at the root of the north-east spiral arm. Since this feature is visible both in polarized light and total intensity and has a blue spectrum, itis likely caused by stellar light scattered by dust. Conclusions. The two spiral arms are at different separations from the star, have very different pitch angles, and are separated by an apparent discontinuity, suggesting they might have a different origin. The very open southern spiral arm might be caused by infalling material from late encounters with cloudlets into the formation environment of the star itself. The compact feature could be caused by interaction with a planet in formation still embedded in its dust envelope and it could be responsible for launching the north–east spiral. We estimate a mass of the putative embedded object of the order of few MJup.

Keck/OSIRIS Paβ High-contrast Imaging and Updated Constraints on PDS 70b

Abstract We present a high-contrast imaging search for Paβ line emission from protoplanets in the PDS 70 system with Keck/OSIRIS integral field spectroscopy. We applied the high-resolution spectral differential imaging technique to the OSIRIS J-band data but did not detect the Paβ line at the level predicted using the parameters of Hashimoto et al. (2020). This lack of Paβ emission suggests the MUSE-based study may have overestimated the line width of Hα. We compared our Paβ detection limits with the previous Hα flux and Hβ limits and estimated A V to be ∼0.9 and 2.0 for PDS 70 b and c, respectively. In particular, PDS 70 b’s A V is much smaller than implied by high-contrast near-infrared studies, which suggests the infrared-continuum photosphere and the hydrogen-emitting regions exist at different heights above the forming planet.

Spectral and angular differential imaging with SPHERE/IFS

Context. Direct imaging of exoplanets is a challenging task that requires state-of-the-art instrumentation and advanced image-processing techniques. Differential imaging techniques have proven useful for the detection of exoplanet companions around stars. Angular differential imaging (ADI) and spectral differential imaging (SDI) are commonly used for direct detection and characterisation of young, Jovian exoplanets in datasets obtained with the SPHERE/IFS instrument. Aims. We compare the performance of ADI, SDI, and three combinations of ADI and SDI to find which technique achieves the highest signal-to-noise ratio (S/N), and we analyse their performance as functions of integration time, field rotation, and wavelength range. Methods. We analyse SPHERE/IFS observations of three known exoplanets, namely β Pictoris b, 51 Eridani b, and HR 8799 e, with five differential imaging techniques. We split the datasets into subsets to vary each parameter before the data are processed with each technique. The differential imaging techniques are applied using principal component analysis (PCA). Results. The tests show that a combination of SDI and ADI consistently achieves better results than ADI alone, and using SDI and ADI simultaneously (combined differential imaging; CODI) achieved the best results. The integration time test shows that targets with a separation larger than 0.24 arcsec observed with an integration time of more than 103 s were photon-noise limited. Field rotation shows a strong correlation with S/N for field rotations up to 1 full width at half maximum (FWHM), after which no significant increase in S/N with field rotation is observed. Wavelength range variation shows a general increase in S/N for broader wavelength ranges, but no clear correlation is seen. Conclusions. Spectral information is essential to boost S/N compared to regular ADI. Our results suggest that CODI should be the preferred processing technique to search for new exoplanets with SPHERE/IFS. To optimise direct-imaging observations, the field rotation should exceed 1 FWHM to detect exoplanets at small separations.

The SPHERE infrared survey for exoplanets (SHINE)

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 (YJandYJHbands), and of 11.8 mag at 200 mas, 13.1 mag at 800 mas, and 15.8 mag at 3 as with IRDIS inHband, 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.

A MUSE view of the asymmetric jet from HD 163296

Context. Jets and outflows are thought to play important roles in regulating star formation and disk evolution. An important question is how the jets are launched. HD 163296 is a well-studied Herbig Ae star that hosts proto-planet candidates, a protoplanetary disk, a protostellar jet, and a molecular outflow, which makes it an excellent laboratory for studying jets. Aims. We aim to characterize the jet at the inner regions and check if there are large differences with the features at large separations. A secondary objective is to demonstrate the performance of Multi Unit Spectroscopic Explorer (MUSE) in high-contrast imaging of extended line emission. Methods. MUSE in the narrow field mode (NFM) can provide observations at optical wavelengths with high spatial (∼75 mas) and medium spectral (R ∼ 2500) resolution. With the high-resolution spectral differential imaging technique, we can characterize the kinematic structures and physical conditions of jets down to 100 mas. Results. We detect multiple atomic lines in two new knots, B3 and A4, at distances of < 4″ from the host star with MUSE. The derived Ṁjet/Ṁacc is about 0.08 and 0.06 for knots B3 and A4, respectively. The observed [Ca II]/[S II] ratios indicate that there is no sign of dust grains at distances of < 4″. Assuming the A4 knot traced the streamline, we can estimate a jet radius at the origin by fitting the half width half maximum of the jet, which sets an upper limit of 2.2 au on the size of the launching region. Although MUSE has the ability to detect the velocity shifts caused by high- and low-velocity components, we found no significant evidence of velocity decrease transverse to the jet direction in our 500 s MUSE observation. Conclusions. Our work demonstrates the capability of using MUSE NFM observations for the detailed study of stellar jets in the optical down to 100 mas. The derived Ṁjet/Ṁacc, no dust grain, and jet radius at the star support the magneto-centrifugal models as a launching mechanism for the jet.

Improved Contrast in Images of Exoplanets Using Direct Signal-to-noise Ratio Optimization

Abstract Direct imaging of exoplanets is usually limited by quasi-static speckles. These uncorrected aberrations in a star’s point-spread function (PSF) obscure faint companions and limit the sensitivity of high-contrast imaging instruments. Most current approaches to processing differential imaging sequences like angular differential imaging and spectral differential imaging produce a self-calibrating data set that is combined using a linear least-squares solution to minimize the noise. Due to temporal and chromatic evolution of a telescope’s PSF, the best correlated reference images are usually the most contaminated by the planet, leading to self-subtraction and reducing the planet throughput. In this paper, we present an algorithm that directly optimizes the nonlinear equation for planet signal-to-noise ratio (S/N). This new algorithm does not require us to reject adjacent reference images and optimally balances noise reduction with self-subtraction. We then show how this algorithm can be applied to multiple images simultaneously for a further reduction in correlated noise, directly maximizing the S/N of the final combined image. Finally, we demonstrate the technique on an illustrative sequence of HR8799 using the new Julia-based Signal to Noise Analysis Pipeline. We show that S/N optimization can provide up to a 5× improvement in contrast close to the star. Applicable to both new and archival data, this technique will allow for the detection of fainter, lower mass, and closer-in companions, or achieve the same sensitivity with less telescope time.

Constraining PDS 70b’s Formation Mechanism with Multi-hydrogen-emission Observations

Abstract We present our Keck/OSIRIS observations of the Paβ emission line (1.282 μm) to investigate accretion mechanisms of PDS 70 planetary system. Our spectral differential imaging reduction to remove the stellar PSF resulted in null detection of Paβ at the locations of PDS 70b and c. The 5σ detection limit of Paβ compared with the theoretical model of Aoyama & Ikoma (2019) indicates the gas velocity onto PDS 70b is smaller than 70 km s−1, which suggests MUSE-based Hα studies overestimated the gas velocity and the mass of PDS 70b.

Searching for proto-planets with MUSE

Context. Protoplanetary disks contain structures such as gaps, rings, and spirals, which are thought to be produced by the interaction between the disk and embedded protoplanets. However, only a few planet candidates are found orbiting within protoplanetary disks, and most of them are being challenged as having been confused with disk features. Aims. The VLT/MUSE discovery of PDS 70 c demonstrated a powerful way of searching for still-forming protoplanets by targeting accretion signatures with medium-resolution integral field spectroscopy. We aim to discover more proto-planetary candidates with MUSE, with a secondary aim of improving the high-resolution spectral differential imaging (HRSDI) technique by analyzing the instrumental residuals of MUSE. Methods. We analyzed MUSE observations of five young stars with various apparent brightnesses and spectral types. We applied the HRSDI technique to perform high-contrast imaging. The detection limits were estimated using fake planet injections. Results. With a 30 min integration time, MUSE can reach 5σ detection limits in apparent Hα line flux down to 10−14 and 10−15 erg s−1 cm−2 at 0.075′′ and 0.25′′, respectively. In addition to PDS 70 b and c, we did not detect any clear accretion signatures in PDS 70, J1850-3147, and V1094 Sco down to 0.1′′. MUSE avoids the small sample statistics problem by measuring the noise characteristics in the spatial direction at multiple wavelengths. We detected two asymmetric atomic jets in HD 163296 with a very high spatial resolution (down to 8 au) and medium spectral resolution (R ~ 2500). Conclusions. The HRSDI technique when applied to MUSE data allows us to reach the photon noise limit at small separations (i.e., <0.5′′). With the combination of high-contrast imaging and medium spectral resolution, MUSE can achieve fainter detection limits in apparent line flux than SPHERE/ZIMPOL by a factor of ~5. MUSE has some instrumental issues that limit the contrast that appear in cases with strong point sources, which can be either a spatial point source due to high Strehl observations or a spectral point source due to a high line-to-continuum ratio. We modified the HRSDI technique to better handle the instrumental artifacts and improve the detection limits. To avoid the instrumental effects altogether, we suggest faint young stars with relatively low Hα line-to-continuum ratio to be the most suitable targets for MUSE to search for potential protoplanets.

High-contrast Hα imaging with Subaru/SCExAO + VAMPIRES

We present the current status of Hα high-contrast imaging observations with Subaru/Subaru Coronagraphic Extreme Adaptive Optics + VAMPIRES. Our adaptive optics correction at optical wavelengths in combination with (double) spectral differential imaging (SDI) and angular differential imaging (ADI) was capable of resolving a ring feature around omi Cet and detect the Hα counterpart of jet around RY Tau. We tested the post-processing by changing the order of ADI and SDI and both of the contrast limits achieved ∼10 − 3 to 5 × 10 − 4 at 0.3″, which is comparable to other Hα high-contrast imaging instruments in the Southern Hemisphere such as very large telescope (VLT)/spectro-polarimetric high-contrast exoplanet research, VLT/MUSE, and Magellan AO. Current wavefront sensing and adaptive optics correction at optical wavelengths empirically depend on airmass, and Subaru/VAMPIRES provide great opportunities for Hα high-contrast imaging for Northern Hemisphere targets.

Carina High-contrast Imaging Project for massive Stars (CHIPS)

Context.Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations (ρ) smaller than 1″ (approx. 1000–3000 au, considering the typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are needed, however, for helping, for example, to distinguish among the various scenarios that describe the formation of massive stars.Aims.With the aim of obtaining a statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys (Δmag ≲ 5 atρ ≲ 1″), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) coronagraphic instrument on the Very Large Telescope (VLT). In this, the first paper of the series, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ Car.Methods.We obtained VLT-SPHERE snapshot observations in the IRDIFS_EXT mode, which combines the IFS and IRDIS sub-systems and simultaneously provides us four-dimensional (4D) data cubes in two different fields-of-view: 1.73″ × 1.73″ for IFS (39 spectral channels across theYJHbands) and 12″ × 12″ for IRDIS (two spectral channels across theKband). Angular- and spectral-differential imaging techniques as well as PSF-fitting were applied to detect and measure the relative flux of the companions in each spectral channel. The latter were then flux-calibrated using theoretical SED models of the central object and compared to a grid of ATLAS9 atmosphere model and (pre-)main-sequence evolutionary tracks, providing a first estimate of the physical properties of the detected companions.Results.Detection limits of 9 mag atρ > 200 mas for IFS, and as faint as 13 mag atρ > 1.​″8 for IRDIS (corresponding to sub-solar masses for potential companions), can be reached in snapshot observations of only a few minutes integration times, allowing us to detect 19 sources around the QZ Car system. All but two are reported here for the first time. With near-IR magnitude contrasts in the range of 4 to 7.5 mag, the three brightest sources (Ab, Ad, and E) are most likely to be physically bound. They have masses in the range of 2 to 12M⊙and are potentially co-eval with QZ Car central system. The remaining sources have flux contrast of 1.5 × 105to 9.5 × 106(ΔK ≈ 11 to 13 mag). Their presence can be explained by the local source density and they are, thus, likely to be chance alignments. If they were members of the Carina nebula, they would be sub-solar-mass pre-main sequence stars.Conclusions.Based on this proof of concept, we show that the VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. It paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

VLT/SPHERE survey for exoplanets around young early-type stars, including systems with multi-belt architectures

Context. Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant submillimeter near-infrared images of young protoplanetary and cool dusty debris disks around main-sequence stars that usually show substantial spatial structures. With two decades of direct imaging of exoplanets already studied, it is striking to note that a majority of recent discoveries of imaged giant planets have been obtained around young early-type stars hosting a circumstellar disk. Aims. Our aim was to create a direct imaging program designed to maximize our chances of giant planet discovery and target 22 young early-type stars. About half of them show indications of multi-belt architectures. Methods. Using the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images, we conducted a systematic search in the close environment of these young, dusty, and early-type stars. We used a combination of angular and spectral differential imaging to reach the best detection performances down to the planetary mass regime. Results. We confirm that companions detected around HIP 34276, HIP 101800, and HIP 117452 are stationary background sources and binary companions. The companion candidates around HIP 8832, HIP 16095, and HIP 95619 are determined as background contaminations. Regarding the stars for which we infer the presence of debris belts, a theoretical minimum mass for planets required to clear the debris gaps can be calculated. The dynamical mass limit is at least 0.1 MJ and can exceed 1 MJ. Direct imaging data is typically sensitive to planets down to ~3.6 MJ at 1′′, and 1.7 MJ in the best case. These two limits tightly constrain the possible planetary systems present around each target. These systems will be probably detectable with the next generation of planet imagers.

PACO ASDI: an algorithm for exoplanet detection and characterization in direct imaging with integral field spectrographs

Context.Exoplanet detection and characterization by direct imaging both rely on sophisticated instruments (adaptive optics and coronagraph) and adequate data processing methods. Angular and spectral differential imaging (ASDI) combines observations at different times and a range of wavelengths in order to separate the residual signal from the host star and the signal of interest corresponding to off-axis sources.Aims.Very high contrast detection is only possible with an accurate modeling of those two components, in particular of the background due to stellar leakages of the host star masked out by the coronagraph. Beyond the detection of point-like sources in the field of view, it is also essential to characterize the detection in terms of statistical significance and astrometry and to estimate the source spectrum.Methods.We extend our recent methodPACO, based on local learning of patch covariances, in order to capture the spectral and temporal fluctuations of background structures. From this statistical modeling, we build a detection algorithm and a spectrum estimation method:PACO ASDI. The modeling of spectral correlations proves useful both in reducing detection artifacts and obtaining accurate statistical guarantees (detection thresholds and photometry confidence intervals).Results.An analysis of several ASDI datasets from the VLT/SPHERE-IFS instrument shows thatPACO ASDIproduces very clean detection maps, for which setting a detection threshold is statistically reliable. Compared to other algorithms used routinely to exploit the scientific results of SPHERE-IFS, sensitivity is improved and many false detections can be avoided. Spectrally smoothed spectra are also produced byPACO ASDI. The analysis of datasets with injected fake planets validates the recovered spectra and the computed confidence intervals.Conclusions.PACO ASDIis a high-contrast processing algorithm accounting for the spatio-spectral correlations of the data to produce statistically-grounded detection maps and reliable spectral estimations. Point source detections, photometric and astrometric characterizations are fully automatized.