skip to main content


Title: Direct Exoplanet Detection using Convolutional Image Reconstruction (ConStruct): A New Algorithm for Post-processing High-contrast Images
Abstract

We present a novel machine-learning approach for detecting faint point sources in high-contrast adaptive optics (AO) imaging data sets. The most widely used algorithms for primary subtraction aim to decouple bright stellar speckle noise from planetary signatures by subtracting an approximation of the temporally evolving stellar noise from each frame in an imaging sequence. Our approach aims to improve the stellar noise approximation and increase the planet detection sensitivity by leveraging deep learning in a novel direct imaging post-processing algorithm. We show that a convolutional autoencoder neural network, trained on an extensive reference library of real imaging sequences, accurately reconstructs the stellar speckle noise at the location of a potential planet signal. This tool is used in a post-processing algorithm we call Direct Exoplanet Detection with Convolutional Image Reconstruction, orConStruct. The reliability and sensitivity ofConStructare assessed using real Keck/NIRC2 angular differential imaging data sets. Of the 30 unique point sources we examine,ConStructyields a higher signal-to-noise ratio than traditional principal component analysis-based processing for 67% of the cases and improves the relative contrast by up to a factor of 2.6. This work demonstrates the value and potential of deep learning to take advantage of a diverse reference library of point-spread function realizations to improve direct imaging post-processing.ConStructand its future improvements may be particularly useful as tools for post-processing high-contrast images from JWST and extreme AO instruments, both for the current generation and those being designed for the upcoming 30 m class telescopes.

 
more » « less
NSF-PAR ID:
10489710
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astronomical Journal
Volume:
167
Issue:
3
ISSN:
0004-6256
Format(s):
Medium: X Size: Article No. 92
Size(s):
["Article No. 92"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    High-contrast imaging has afforded astronomers the opportunity to study light directly emitted by adolescent (tens of megayears) and “proto” (<10 Myr) planets still undergoing formation. Direct detection of these planets is enabled by empirical point-spread function (PSF) modeling and removal algorithms. The computational intensity of such algorithms, as well as their multiplicity of tunable input parameters, has led to the prevalence of ad hoc optimization approaches to high-contrast imaging results. In this work, we present a new, systematic approach to optimization vetted using data of the high-contrast stellar companion HD 142527 B from the Magellan Adaptive Optics Giant Accreting Protoplanet Survey (GAPlanetS). More specifically, we present a grid search technique designed to explore three influential parameters of the PSF subtraction algorithmpyKLIP: annuli, movement, and KL modes. We consider multiple metrics for postprocessed image quality in order to optimally recover at Hα(656 nm) synthetic planets injected into contemporaneous continuum (643 nm) images. These metrics include peak (single-pixel) signal-to-noise ratio (S/N), average (multipixel average) S/N, 5σcontrast, and false-positive fraction. We apply continuum-optimized KLIP reduction parameters to six Hαdirect detections of the low-mass stellar companion HD 142527 B and recover the companion at a range of separations. Relative to a single-informed, nonoptimized set of KLIP parameters applied to all data sets uniformly, our multimetric grid search optimization led to improvements in companion S/N of up to 1.2σ, with an average improvement of 0.6σ. Since many direct imaging detections lie close to the canonical 5σthreshold, even such modest improvements may result in higher yields in future imaging surveys.

     
    more » « less
  2. Abstract

    TESS has proven to be a powerful resource for finding planets, including those that orbit the most prevalent stars in our galaxy: M dwarfs. Identification of stellar companions (both bound and unbound) has become a standard component of the transiting planet confirmation process in order to assess the level of light-curve dilution and the possibility of the target being a false positive. Studies of stellar companions have also enabled investigations into stellar multiplicity in planet-hosting systems, which has wide-ranging implications for both exoplanet detection and characterization, as well as for the formation and evolution of planetary systems. Speckle and AO imaging are some of the most efficient and effective tools for revealing close-in stellar companions; we therefore present observations of 58 M-dwarf TOIs obtained using a suite of speckle imagers at the 3.5 m WIYN telescope, the 4.3 m Lowell Discovery Telescope, and the 8.1 m Gemini North and South telescopes. These observations, as well as near-infrared adaptive optics images obtained for a subset (14) of these TOIs, revealed only two close-in stellar companions. Upon surveying the literature, and cross-matching our sample with Gaia, SUPERWIDE, and the catalog from El-Badry et al., we reveal an additional 15 widely separated common proper motion companions. We also evaluate the potential for undetected close-in companions. Taking into consideration the sensitivity of the observations, our findings suggest that the orbital period distribution of stellar companions to planet-hosting M dwarfs is shifted to longer periods compared to the expected distribution for field M dwarfs.

     
    more » « less
  3. One of the top priorities in observational astronomy is the direct imaging and characterization of extrasolar planets (exoplanets) and planetary systems. Direct images of rocky exoplanets are of particular interest in the search for life beyond the Earth, but they tend to be rather challenging targets since they are orders-of-magnitude dimmer than their host stars and are separated by small angular distances that are comparable to the classicalλ<#comment/>/Ddiffraction limit, even for the coming generation of 30 m class telescopes. Current and planned efforts for ground-based direct imaging of exoplanets combine high-order adaptive optics (AO) with a stellar coronagraph observing at wavelengths ranging from the visible to the mid-IR. The primary barrier to achieving high contrast with current direct imaging methods is quasi-static speckles, caused largely by non-common path aberrations (NCPAs) in the coronagraph optical train. Recent work has demonstrated that millisecond imaging, which effectively “freezes” the atmosphere’s turbulent phase screens, should allow the wavefront sensor (WFS) telemetry to be used as a probe of the optical system to measure NCPAs. Starting with a realistic model of a telescope with an AO system and a stellar coronagraph, this paper provides simulations of several closely related regression models that take advantage of millisecond telemetry from the WFS and coronagraph’s science camera. The simplest regression model, called the naïve estimator, does not treat the noise and other sources of information loss in the WFS. Despite its flaws, in one of the simulations presented herein, the naïve estimator provides a useful estimate of an NCPA of∼<#comment/>0.5radian RMS (≈<#comment/>λ<#comment/>/13), with an accuracy of∼<#comment/>0.06radian RMS in 1 min of simulated sky time on a magnitude 8 star. Thebias-corrected estimatorgeneralizes the regression model to account for the noise and information loss in the WFS. A simulation of the bias-corrected estimator with 4 min of sky time included an NCPA of∼<#comment/>0.05radian RMS (≈<#comment/>λ<#comment/>/130) and an extended exoplanet scene. The joint regression of the bias-corrected estimator simultaneously achieved an NCPA estimate with an accuracy of∼<#comment/>5×<#comment/>10−<#comment/>3radian RMS and an estimate of the exoplanet scene that was free of the self-subtraction artifacts typically associated with differential imaging. The5σ<#comment/>contrast achieved by imaging of the exoplanet scene was∼<#comment/>1.7×<#comment/>10−<#comment/>4at a distance of3λ<#comment/>/Dfrom the star and∼<#comment/>2.1×<#comment/>10−<#comment/>5at10λ<#comment/>/D. These contrast values are comparable to the very best on-sky results obtained from multi-wavelength observations that employ both angular differential imaging (ADI) and spectral differential imaging (SDI). This comparable performance is despite the fact that our simulations are quasi-monochromatic, which makes SDI impossible, nor do they have diurnal field rotation, which makes ADI impossible. The error covariance matrix of the joint regression shows substantial correlations in the exoplanet and NCPA estimation errors, indicating that exoplanet intensity and NCPA need to be estimated self-consistently to achieve high contrast.

     
    more » « less
  4. The success of ground-based, high contrast imaging for the detection of exoplanets in part depends on the ability to differentiate between quasi-static speckles caused by aberrations not corrected by adaptive optics (AO) systems, known as non-common path aberrations (NCPAs), and the planet intensity signal. Frazin (ApJ, 2013) introduced a post-processing algorithm demonstrating that simultaneous millisecond exposures in the science camera and wavefront sensor (WFS) can be used with a statistical inference procedure to determine both the series expanded NCPA coefficients and the planetary signal. We demonstrate, via simulation, that using this algorithm in a closed-loop AO system, real-time estimation and correction of the quasi-static NCPA is possible without separate deformable mirror (DM) probes. Thus the use of this technique allows for the removal of the quasi-static speckles that can be mistaken for planetary signals without the need for new optical hardware, improving the efficiency of ground-based exoplanet detection. In our simulations, we explore the behavior of the Frazin Algorithm (FA) and the dependence of its convergence to an accurate estimate on factors such as Strehl ratio, NCPA strength, and number of algorithm search basis functions. We then apply this knowledge to simulate running the algorithm in real-time in a nearly ideal setting. We then discuss adaptations that can be made to the algorithm to improve its real-time performance, and show their efficacy in simulation. A final simulation tests the technique’s resilience against imperfect knowledge of the AO residual phase, motivating an analysis of the feasibility of using this technique in a real closed-loop Extreme AO system such as SCExAO or MagAO-X, in terms of computational complexity and the accuracy of the estimated quasi-static NCPA correction. 
    more » « less
  5. Abstract

    We introduce a new framework for point-spread function subtraction based on the spatiotemporal variation of speckle noise in high-contrast imaging data where the sampling timescale is faster than the speckle evolution timescale. One way that space–time covariance arises in the pupil is as atmospheric layers translate across the telescope aperture and create small, time-varying perturbations in the phase of the incoming wavefront. The propagation of this field to the focal plane preserves some of that space–time covariance. To utilize this covariance, our new approach uses a Karhunen–Loève transform on an image sequence, as opposed to a set of single reference images as in previous applications of Karhunen–Loève Image Processing (KLIP) for high-contrast imaging. With the recent development of photon-counting detectors, such as microwave kinetic inductance detectors, this technique now has the potential to improve contrast when used as a post-processing step. Preliminary testing on simulated data shows this technique can improve contrast by at least 10%–20% from the original image, with significant potential for further improvement. For certain choices of parameters, this algorithm may provide larger contrast gains than spatial-only KLIP.

     
    more » « less