Search

Direct imaging of exoplanets is a challenging task that involves distinguishing faint planetary signals from the overpowering glare of their host stars, often obscured by time-varying stellar noise known as "speckles". The predominant algorithms for speckle noise subtraction employ principal-based point spread function (PSF) fitting techniques to discern planetary signals from stellar speckle noise. We introduce torchKLIP, a benchmark package developed within the machine learning (ML) framework PyTorch. This work enables ML techniques to utilize extensive PSF libraries to enhance direct imaging post-processing. Such advancements promise to improve the post-processing of high-contrast images from leading-edge astronomical instruments like the James Webb Space Telescope and extreme adaptive optics systems., Comment: 8 pages, 6 figures, SPIE Optics and Photonics 2024

Real-time control (RTC) is pivotal for any Adaptive Optics (AO) system, including high-contrast imaging of exoplanets and circumstellar environments. It is the brain of the AO system, and what wavefront sensing and control (WFS\&C) techniques need to work with to achieve unprecedented image quality and contrast, ultimately advancing our understanding of exoplanetary systems in the context of high contrast imaging (HCI). Developing WFS\&C algorithms first happens in simulation or a lab before deployment on-sky. The transition to on-sky testing is often challenging due to the different RTCs used. Sharing common RTC standards across labs and telescope instruments would considerably simplify this process. A data architecture based on the interprocess communication method known as shared memory is ideally suited for this purpose. The CACAO package, an example of RTC based on shared memory, was initially developed for the Subaru-SCExAO instrument and now deployed on several benches and instruments. This proceeding discusses the challenges, requirements, implementation strategies, and performance evaluations associated with integrating a shared memory-based RTC. The Santa Cruz Extreme AO Laboratory (SEAL) bench is a platform for WFS\&C development for large ground-based segmented telescopes. Currently, SEAL offers the user a non-real-time version of CACAO, a shared-memory based RTC package initially developed for the Subaru-SCExAO instrument, and now deployed on several benches and instruments. We show here the example of the SEAL RTC upgrade as a precursor to both RTC upgrade at the 3-m Shane telescopes at Lick Observatory (Shane-AO) and a future development platform for the Keck II AO. This paper is aimed at specialists in AO, astronomers, and WFS\&C scientists seeking a deeper introduction to the world of RTCs., Comment: SPIE Astronomical Telescope and Intrumentation 2024

Continuous wavefront sensing benefits space observatories in on-orbit optical performance maintenance. To measure the phase of a wavefront, phase retrieval is an attractive technique as it uses multiple point spread function (PSF) images that are acquired by the telescope itself without extra metrology systems nor complicated calibration. The focus diverse phase retrieval utilizes PSFs from predetermined defocused positions to enhance the dynamic range of the algorithm. We describe an updated visible light active optics testbed with the addition of a linear motorized focus stage. The performance of the phase retrieval algorithm in broadband is tested under various cases. While broadband pass filters have advantages in higher signal-to-noise ratio (SNR), the performance of phase retrieval can be restricted due to blurred image caused by diffraction and increased computing cost. We used multiple bandpass filters (10 nm, 88 nm, and 150 nm) and investigated effects of bandwidth on the accuracy and required image acquisition conditions such as SNR, reaching accuracies below 20 nm RMS wavefront error at the widest bandwidth. We also investigated the dynamic range of the phase retrieval algorithm depending on the bandwidth and required amount of defocus to expand dynamic range. Finally, we simulated the continuous wavefront sensing and correction loop with a range of statistically generated representative telescope disturbance time series to test for edge cases.

When a telescope doesn't reach a reasonable point spread function on the detector or detectable wavefront quality after initial assembly, a coarse phase alignment on-sky is crucial. Before utilizing a closed loop adaptive optics system, the observatory needs a strategy to actively align the telescope sufficiently for fine wavefront sensing. This paper presents a method of early-stage alignment using a stochastic parallel-gradient-descent (SPGD) algorithm which performs random perturbations to the optics of a three mirror anastigmat telescope design. The SPGD algorithm will drive the telescope until the wavefront error is below the acceptable range of the fine adaptive optics system to hand the telescope over. The focused spot size over the field of view is adopted as a feed parameter to the SPGD algorithm and wavefront peak-to-valley error values are monitored to directly compare our mechanical capabilities to our alignment goal of diffraction limited imaging and fine wavefront sensing., Comment: 7 pages, 4 figures, presented at SPIE Optics + Photonics 2024

As wavefront quality demands tighten on space systems for applications such as astronomy and laser communication, mounting small optics such that the wavefront is undisturbed, positioning is adjustable and the design is producible, while surviving harsh space environments, is a continuing challenge. We designed multiple candidate flexure mounts to support small optics (up to 50 mm diameter, and over 100 grams) to survive the qualification and acceptance tests of small spacecraft and units as defined in ISO 19683 and a mounting structure which is adjustable in decenter [+/-0.5mm], tip/tilt +/-0.5deg, and piston [+/-0.25mm]. We will present design details along with measurements showing less than approximately lambda/10 wavefront contribution from the optic bonding process, along with thermal and multi-axis vibration test data showing the mounted optics survived the acceptance testing loads and are suitable for operation in a wide range of harsh environments.

Epsilon Eridani ($\epsilon$ Eri) is one of the first debris disk systems detected by the Infrared Astronomical Satellite (IRAS). However, the system has thus far eluded detection in scattered light with no components having been directly imaged. Its similarity to a relatively young Solar System combined with its proximity makes it an excellent candidate to further our understanding of planetary system evolution. We present a set of coronagraphic images taken using the Space Telescope Imaging Spectrograph (STIS) coronagraph on the Hubble space telescope at a small inner working angle to detect a predicted warm inner debris disk inside 1". We used three different post-processing approaches; Non-negative Matrix Factorization (NMF), Karhunen-Lo`eve Image Processing (KLIP), and Classical reference differential imaging (RDI), to best optimize reference star subtraction, and find that NMF performed the best overall while KLIP produced the absolute best contrast inside 1". We present limits on scattered light from warm dust, with constraints on surface brightness at 6 mJy/as$^2$ at our inner working angle of 0.6". We also place a constraint of 0.5 mJy/as$^2$ outside 1", which gives us an upper limit on the brightness for outer disks and substellar companions. Finally, we calculated an upper limit on the dust albedo at $\omega <$ 0.487., Comment: 13+2 pages, 7+2 figures; Accepted for publication in the Astronomical Journal

Future space observatories will likely have segmented primaries, causing diffraction effects that reduce coronagraph performance. Reflective binary pupil apodizer masks can mitigate these, with the metamaterial black silicon (BSi) showing promise as a strong absorber. To bring contrast ratios to the $10^-{10}$ level as needed to observe Earth-like exoplanets, feature sizes on these BSi masks will need to be less than $5$ microns when paired with MEMS (micro-electromechanical systems) deformable mirrors. As scalar diffraction cannot reliably model this feature size, we developed a Finite-Difference Time-Domain (FDTD) model of BSi masks using Meep software. We characterize the FDTD-derived polarization-dependent bidirectional reflectance distribution function of BSi and discuss the model's shortcomings., Comment: 8 pages, 10 figures, submitted to SPIE Astronomical Telescopes and Instrumentation (AS24)

The Space Coronagraph Optical Bench (SCoOB) is a high-contrast imaging testbed built to demonstrate starlight suppression techniques at visible wavelengths in a space-like vacuum environment. The testbed is designed to achieve ${<}10^{-8}$ contrast from $3-10\lambda/D$ in a one-sided dark hole using a liquid crystal vector vortex waveplate and a 952-actuator Kilo-C deformable mirror (DM) from Boston Micromachines (BMC). We have recently expanded the testbed to include a field stop for mitigation of stray/scattered light, a precision-fabricated pinhole in the source simulator, a Minus K passive vibration isolation table for jitter reduction, and a low-noise vacuum-compatible CMOS sensor. We report the latest contrast performance achieved using implicit electric field conjugation (iEFC) at a vacuum of ${\sim}10^{-6}$ Torr and over a range of bandpasses with central wavelengths from 500 to 650nm and bandwidths (BW) from $\ll 1\%$ to 15\%. Our jitter in vacuum is $<3\times10^{-3} \lambda/D$, and the best contrast performance to-date in a half-sided D-shaped dark hole is $2.2\times10^{-9}$ in a $\ll 1 \%$ BW, $4\times10^{-9}$ in a 2\% BW, and $2.5\times10^{-8}$ in a 15\% BW., Comment: 16 pages, 12 figures, SPIE Astronomical Telescopes and Instrumentation 2024

High-contrast imaging in the next decade aims to image exoplanets at smaller angular separations and deeper contrasts than ever before. A problem that has recently garnered attention for telescopes equipped with high-contrast coronagraphs is polarization aberration arising from the optics. These aberrations manifest as low-order aberrations of different magnitudes for orthogonal polarization states and spread light into the dark hole of the coronagraph that cannot be fully corrected. The origin of polarization aberrations has been modeled at the telescope level. However, we don't fully understand how polarization aberrations arise at the instrument level. To directly measure this effect, we construct a dual-rotating-retarder polarimeter around the SCoOB high-contrast imaging testbed to measure its Mueller matrix. With this matrix, we directly characterize the diattenuation, retardance, and depolarization of the instrument as a function of position in the exit pupil. We measure the polarization aberrations in the Lyot plane to understand how polarization couples into high-contrast imaging residuals., Comment: 16 pages, 13 figures

High-order wavefront sensing and control (HOWFSC) is key to create a dark hole region within the coronagraphic image plane where high contrasts are achieved. The Roman Coronagraph is expected to perform its HOWFSC with a ground-in-the-loop scheme due to the computational complexity of the Electric Field Conjugation (EFC) algorithm. This scheme provides the flexibility to alter the HOWFSC algorithm for given science objectives. The baseline HOWFSC scheme involves running EFC while observing a bright star such as {\zeta} Puppis to create the initial dark hole followed by a slew to the science target. The new implicit EFC (iEFC) algorithm removes the optical diffraction model from the controller, making the final contrast independent of model accuracy. While previously demonstrated with a single DM, iEFC is extended to two deformable mirror systems in order to create annular dark holes. The algorithm is then applied to the Wide-Field-of-View Shaped Pupil Coronagraph (SPC-WFOV) mode designed for the Roman Space Telescope using end-to-end physical optics models. Initial monochromatic simulations demonstrate the efficacy of iEFC as well as the optimal choice of modes for the SPC-WFOV instrument. Further simulations with a 3.6% wavefront control bandpass and a broader 10% bandpass then demonstrate that iEFC can be used in broadband scenarios to achieve contrasts below 1E-8 with Roman. Finally, an EMCCD model is implemented to estimate calibration times and predict the controller's performance. Here, 1E-8 contrasts are achieved with a calibration time of about 6.8 hours assuming the reference star is {\zeta} Puppis. The results here indicate that iEFC can be a valid HOWFSC method that can mitigate the risk of model errors associated with space-borne coronagraphs, but to maximize iEFC performance, lengthy calibration times will be required to mitigate the noise accumulated during calibration., Comment: 20 pages, 15 figures

Paraxial diffraction modeling based on the Fourier transform has seen widespread implementation for simulating the response of a diffraction-limited optical system. For systems where the paraxial assumption is not sufficient, a class of algorithms has been developed that employs hybrid propagation physics to compute the propagation of an elementary beamlet along geometric ray paths. These "beamlet decomposition" algorithms include the well-known Gaussian Beamlet Decomposition (GBD) algorithm, of which several variants have been created. To increase the computational efficiency of the GBD algorithm, we derive an alternative expression of the technique that utilizes the analytical propagation of beamlets to tilted planes. We then use this accelerated algorithm to conduct a parameter-space search to find the optimal combination of free parameters in GBD to construct the analytical Airy function. The experiment is conducted on a consumer-grade CPU, and a high-performance GPU, where the new algorithm is 34 times faster than the previously published algorithm on CPUs, and 67,513 times faster on GPUs., Comment: 19 pages, 10 figures, Accepted to Optics Express on April 18 2024

Directly observing exoplanets with coronagraphs is impeded by the presence of speckles from aberrations in the optical path, which can be mitigated in hardware with wavefront control as well as in post-processing. This work explores using an instrument model in post-processing to separate astrophysical signals from residual aberrations in coronagraphic data. The effect of wavefront error (WFE) on the coronagraphic intensity consists of a linear contribution and a quadratic contribution. When either of the terms is much larger than the other, the instrument response can be approximated by a transfer matrix mapping WFE to detector plane intensity. From this transfer matrix, a useful projection onto instrumental modes that removes the dominant error modes can be derived. We apply this projection to synthetically generated Roman Space Telescope hybrid Lyot coronagraph (HLC) data to extract "robust observables," which can be used instead of raw data for applications such as detection testing. The projection improves planet flux ratio detection limits by about 28% in the linear regime and by over a factor of 2 in the quadratic regime, illustrating that robust observables can increase sensitivity to astrophysical signals and improve the scientific yield from coronagraphic data. While this approach does not require additional information such as observations of reference stars or modulations of a deformable mirror, it can and should be combined with these other techniques, acting as a model-informed prior in an overall post-processing strategy., Comment: 19 pages, 10 figures