Kernel Phase and Coronagraphy with Automatic Differentiation

Abstract

The accumulation of aberrations along the optical path in a telescope produces distortions and speckles in the resulting images, limiting the performance of cameras at high angular resolution. It is important to achieve the highest possible sensitivity to faint sources such as planets, using both hardware and data analysis software. While analytic methods are efficient, real systems are better-modelled numerically, but such models with many parameters can be hard to understand, optimize and apply. Automatic differentiation software developed for machine learning now makes calculating derivatives with respect to aberrations straightforward for arbitrary optical systems. We apply this powerful new tool to enhance high-angular-resolution astronomical imaging. Self-calibrating observables such as the 'closure phase' or 'bispectrum' have been widely used in optical and radio astronomy to mitigate optical aberrations and achieve high-fidelity imagery. Kernel phases are a generalization of closure phases in the limit of small phase errors. Using automatic differentiation, we reproduce existing kernel phase theory within this framework and demonstrate an extension to the Lyot coronagraph, finding self-calibrating combinations of speckles which are resistant to phase noise, but only in the very high-wavefront-quality regime. As an illustrative example, we reanalyze Palomar adaptive optics observations of the binary alpha Ophiuchi, finding consistency between the new pipeline and the existing standard. We present a new Python package 'morphine' that incorporates these ideas, with an interface similar to the popular package poppy, for optical simulation with automatic differentiation. These methods may be useful for designing improved astronomical optical systems by gradient descent.

Accepted ApJ

References in corpus (14)

• Array Programming with NumPy
• Polarimetry with the Gemini Planet Imager: Methods, Performance at First Light, and the Circumstellar Ring around HR 4796A
• Fast computation of Lyot-style coronagraph propagation
• Using Automatic Differentiation as a General Framework for Ptychographic Reconstruction
• Beyond the Kepler/K2 bright limit: variability in the seven brightest members of the Pleiades
• A PSF-based Approach to TESS High quality data Of Stellar clusters (PATHOS) -- I. Search for exoplanets and variable stars in the field of 47 Tuc
• High precision astrometry with a diffractive pupil telescope
• Comparing Non-Redundant Masking and Filled-Aperture Kernel Phase for Exoplanet Detection and Characterization
• Kernel phase imaging with VLT/NACO: high-contrast detection of new candidate low-mass stellar companions at the diffraction limit
• Astrophotonics: molding the flow of light in astronomical instruments
• Recovering saturated images for high dynamic Kernel-Phase analysis Application to the determination of dynamical masses for the system Gl 494AB
• Laboratory Demonstration of Spatial Linear Dark Field Control For Imaging Extrasolar Planets in Reflected Light
• Kernel Phase and Kernel Amplitude in Fizeau Imaging
• Wavefront sensing from the image domain with the Oxford-SWIFT integral field spectrograph