E. Sturm, R.Davies, J. Alves, Y. Clénet, J. Kotilainen, A. Monna, H. Nicklas, J.-U. Pott, E. Tolstoy, B. Vulcani, J. Achren, S. Annadevara, H. Anwand-Heerwart, C. Arcidiacono, S. Barboza, L. Barl, P. Baudoz, R. Bender, N. Bezawada, F. Biondi, et al (134) MICADO is a first light instrument for the Extremely Large Telescope (ELT), set to start operating later this decade. It will provide diffraction limited imaging, astrometry, high contrast imaging, and long slit spectroscopy at near-infrared wavelengths. During the initial phase operations, adaptive optics (AO) correction will be provided by its own natural guide star wavefront sensor. In its final configuration, that AO system will be retained and complemented by the laser guide star multi-conjugate adaptive optics module MORFEO (formerly known as MAORY). Among many other things, MICADO will study exoplanets, distant galaxies and stars, and investigate black holes, such as Sagittarius A* at the centre of the Milky Way. After their final design phase, most components of MICADO have moved on to the manufacturing and assembly phase. Here we summarize the final design of the instrument and provide an overview about its current manufacturing status and the timeline. Some lessons learned from the final design review process will be presented in order to help future instrumentation projects to cope with the challenges arising from the substantial differences between projects for 8-10m class telescopes (e.g. ESO-VLT) and the next generation Extremely Large Telescopes (e.g. ESO-ELT). Finally, the expected performance will be discussed in the context of the current landscape of astronomical observatories and instruments. For instance, MICADO will have similar sensitivity as the James Webb Space Telescope (JWST), but with six times the spatial resolution.
Matteo Simioni, Daniel Jodlbauer, Carmelo Arcidiacono, Andrea Grazian, Marco Gullieuszik, Elisa Portaluri, Benedetta Vulcani, Roland Wagner, Anita Zanella, Johanna Hartke, Tapio Helin, Hanindyo Kuncarayakti, Elena Masciadri, Fernando Pedichini, Roberto Piazzesi, Alessio Turchi, Piero Vaccari The highest scientific return, for adaptive optics (AO) observations, is achieved with a reliable reconstruction of the PSF. This is especially true for MICADO@ELT. In this presentation, we will focus on extending the MICADO PSF reconstruction (PSF-R) method to the off-axis case. Specifically, a novel approach based on temporal-based tomography of AO telemetry data has been recently implemented. Results from the PSF-R of both simulated and real data show that, at half isoplanatic angle distances, a precision of about 10-15% is achievable in both Strehl ratio and full-width at half maximum, paving the way to extend the MICADO PSF-R tool also to the multi-conjugated AO case.
Time-delay error is a significant error source in adaptive optics (AO) systems. It arises from the latency between sensing the wavefront and applying the correction. Predictive control algorithms reduce the time-delay error, providing significant performance gains, especially for high-contrast imaging. However, the predictive controller's performance depends on factors such as the WFS type, the measurement noise, the AO system's geometry, and the atmospheric conditions. This work studies the limits of prediction under different imaging conditions through spatiotemporal Gaussian process models. The method provides a predictive reconstructor that is optimal in the least-squares sense, conditioned on the fixed times series of WFS data and our knowledge of the atmosphere. We demonstrate that knowledge is power in predictive AO control. With an SHS-based extreme AO instrument, perfect knowledge of Frozen Flow evolution (wind and Cn2 profile) leads to a reduction of the residual wavefront phase variance up to a factor of 3.5 compared to a non-predictive approach. If there is uncertainty in the profile or evolution models, the gain is more modest. Still, assuming that only effective wind speed is available (without direction) led to reductions in variance by a factor of 2.3. We also study the value of data for predictive filters by computing the experimental utility for different scenarios to answer questions such as: How many past data frames should the prediction filter consider, and is it always most advantageous to use the most recent data? We show that within the scenarios considered, more data consistently increases prediction accuracy. Further, we demonstrate that given a computational limitation on how many past frames we can use, an optimized selection of $n$ past frames leads to a 10-15% additional improvement in RMS over using the n latest consecutive frames of data.
Direct imaging of Earth-like exoplanets is one of the most prominent scientific drivers of the next generation of ground-based telescopes. Typically, Earth-like exoplanets are located at small angular separations from their host stars, making their detection difficult. Consequently, the adaptive optics (AO) system's control algorithm must be carefully designed to distinguish the exoplanet from the residual light produced by the host star. A new promising avenue of research to improve AO control builds on data-driven control methods such as Reinforcement Learning (RL). RL is an active branch of the machine learning research field, where control of a system is learned through interaction with the environment. Thus, RL can be seen as an automated approach to AO control, where its usage is entirely a turnkey operation. In particular, model-based reinforcement learning (MBRL) has been shown to cope with both temporal and misregistration errors. Similarly, it has been demonstrated to adapt to non-linear wavefront sensing while being efficient in training and execution. In this work, we implement and adapt an RL method called Policy Optimization for AO (PO4AO) to the GHOST test bench at ESO headquarters, where we demonstrate a strong performance of the method in a laboratory environment. Our implementation allows the training to be performed parallel to inference, which is crucial for on-sky operation. In particular, we study the predictive and self-calibrating aspects of the method. The new implementation on GHOST running PyTorch introduces only around 700 microseconds in addition to hardware, pipeline, and Python interface latency. We open-source well-documented code for the implementation and specify the requirements for the RTC pipeline. We also discuss the important hyperparameters of the method, the source of the latency, and the possible paths for a lower latency implementation.
Andrea Grazian, Matteo Simioni, Carmelo Arcidiacono, Jani Achren, Yann Clenet, Yixian Cao, Richard Davies, Marco Gullieuszik, Tapio Helin, Daniel Jodlbauer, Hanindyo Kuncarayakti, Miska Le Louarn, Seppo Mattila, Fernando Pedichini, Roberto Piazzesi, Elisa Portaluri, Aleksi Salo, Gijs Verdoes Kleijn, Benedetta Vulcani, Roland Wagner, et al (3) MICADO is a workhorse instrument for the ESO ELT, allowing first light capability for diffraction limited imaging and long-slit spectroscopy at near-infrared wavelengths. The PSF Reconstruction (PSF-R) Team of MICADO is currently implementing, for the first time within all ESO telescopes, a software service devoted to the blind reconstruction of the PSF. This tool will work independently of the science data, using adaptive optics telemetry data, both for Single Conjugate (SCAO) and Multi-Conjugate Adaptive Optics (MCAO) allowed by the MORFEO module. The PSF-R service will support the state-of-the-art post-processing scientific analysis of the MICADO imaging and spectroscopic data. We provide here an update of the status of the PSF-R service tool of MICADO, after successfully fulfilling the Final Design Review phase, and discuss recent results obtained on simulated and real data gathered on instruments similar to MICADO.
Matteo Simioni, Carmelo Arcidiacono, Roland Wagner, Andrea Grazian, Marco Gullieuszik, Elisa Portaluri, Benedetta Vulcani, Anita Zanella, Guido Agapito, Richard Davies, Tapio Helin, Fernando Pedichini, Roberto Piazzesi, Enrico Pinna, Ronny Ramlau, Fabio Rossi, Aleksi Salo Current state-of-the-art adaptive optics (AO) provides ground-based, diffraction-limited observations with high Strehl ratios (SR). However, a detailed knowledge of the point spread function (PSF) is required to fully exploit the scientific potential of these data. This is even more crucial for the next generation AO instruments that will equip 30-meter class telescopes, as the characterization of the PSF will be mandatory to fulfill the planned scientific requirements. For this reason, there is a growing interest in developing tools that accurately reconstruct the observed PSF of AO systems, the so-called PSF reconstruction. In this context, a PSF-R service is a planned deliverable for the MICADO@ELT instrument and our group is in charge of its development. In the case of MICADO, a blind PSF-R approach is being pursued to have the widest applicability to science cases. This means that the PSF is reconstructed without extracting information from the science data, relying only on telemetry and calibrations. While our PSF-R algorithm is currently being developed, its implementation is mature enough to test performances with actual observations. In this presentation we will discuss the reliability of our reconstructed PSFs and the uncertainties introduced in the measurements of scientific quantities for bright, on-axis observations taken with the SOUL+LUCI instrument of the LBT. This is the first application of our algorithm to real data. It demonstrates its readiness level and paves the way to further testing. Our PSF-R algorithm is able to reconstruct the SR and full-width at half maximum of the observed PSFs with errors smaller than 2% and 4.5%, respectively. We carried out the scientific evaluation of the obtained reconstructed PSFs thanks to a dedicated set of simulated observations of an ideal science case.
Matteo Simioni, Carmelo Arcidiacono, Roland Wagner, Andrea Grazian, Marco Gullieuszik, Elisa Portaluri, Benedetta Vulcani, Anita Zanella, Guido Agapito, Richard Davies, Tapio Helin, Fernando Pedichini, Roberto Piazzesi, Enrico Pinna, Ronny Ramlau, Fabio Rossi, Aleksi Salo This paper presents the status of an ongoing project aimed at developing a PSF reconstruction software for adaptive optics (AO) observations. In particular, we test for the first time the implementation of pyramid wave-front sensor data on our algorithms. As a first step in assessing its reliability, we applied the software to bright, on-axis, point-like sources using two independent sets of observations, acquired with the single-conjugated AO upgrade for the Large Binocular Telescope. Using only telemetry data, we reconstructed the PSF by carefully calibrating the instrument response. The accuracy of the results has been first evaluated using the classical metric: specifically, the reconstructed PSFs differ from the observed ones by less than 2% in Strehl ratio and 4.5% in full-width at half maximum. Moreover, the recovered encircled energy associated with the PSF core is accurate at 4% level in the worst case. The accuracy of the reconstructed PSFs has then been evaluated by considering an idealized scientific test-case consisting in the measurements of the morphological parameters of a compact galaxy. In the future, our project will include the analysis of anisoplanatism, low SNR regimes, and the application to multi-conjugated AO observations.
J. Nousiainen, C. Rajani, M. Kasper, T. Helin, S. Y. Haffert, C. Vérinaud, J. R. Males, K. Van Gorkom, L. M. Close, J. D. Long, A. D. Hedglen, O. Guyon, L. Schatz, M. Kautz, J. Lumbres, A. Rodack, J.M. Knight, K. Miller The direct imaging of potentially habitable Exoplanets is one prime science case for the next generation of high contrast imaging instruments on ground-based extremely large telescopes. To reach this demanding science goal, the instruments are equipped with eXtreme Adaptive Optics (XAO) systems which will control thousands of actuators at a framerate of kilohertz to several kilohertz. Most of the habitable exoplanets are located at small angular separations from their host stars, where the current XAO systems' control laws leave strong residuals.Current AO control strategies like static matrix-based wavefront reconstruction and integrator control suffer from temporal delay error and are sensitive to mis-registration, i.e., to dynamic variations of the control system geometry. We aim to produce control methods that cope with these limitations, provide a significantly improved AO correction and, therefore, reduce the residual flux in the coronagraphic point spread function. We extend previous work in Reinforcement Learning for AO. The improved method, called PO4AO, learns a dynamics model and optimizes a control neural network, called a policy. We introduce the method and study it through numerical simulations of XAO with Pyramid wavefront sensing for the 8-m and 40-m telescope aperture cases. We further implemented PO4AO and carried out experiments in a laboratory environment using MagAO-X at the Steward laboratory. PO4AO provides the desired performance by improving the coronagraphic contrast in numerical simulations by factors 3-5 within the control region of DM and Pyramid WFS, in simulation and in the laboratory. The presented method is also quick to train, i.e., on timescales of typically 5-10 seconds, and the inference time is sufficiently small (< ms) to be used in real-time control for XAO with currently available hardware even for extremely large telescopes.
Reinforcement Learning (RL) presents a new approach for controlling Adaptive Optics (AO) systems for Astronomy. It promises to effectively cope with some aspects often hampering AO performance such as temporal delay or calibration errors. We formulate the AO control loop as a model-based RL problem (MBRL) and apply it in numerical simulations to a simple Shack-Hartmann Sensor (SHS) based AO system with 24 resolution elements across the aperture. The simulations show that MBRL controlled AO predicts the temporal evolution of turbulence and adjusts to mis-registration between deformable mirror and SHS which is a typical calibration issue in AO. The method learns continuously on timescales of some seconds and is therefore capable of automatically adjusting to changing conditions.
Markus Kasper, Nelly Cerpa Urra, Prashant Pathak, Markus Bonse, Jalo Nousiainen, Byron Engler, Cédric Taïssir Heritier, Jens Kammerer, Serban Leveratto, Chang Rajani, Paul Bristow, Miska Le Louarn, Pierre-Yves Madec, Stefan Ströbele, Christophe Verinaud, Adrian Glauser, Sascha P. Quanz, Tapio Helin, Christoph Keller, Frans Snik, et al (5) The Planetary Camera and Spectrograph (PCS) for the Extremely Large Telescope (ELT) will be dedicated to detecting and characterising nearby exoplanets with sizes from sub-Neptune to Earth-size in the neighbourhood of the Sun. This goal is achieved by a combination of eXtreme Adaptive Optics (XAO), coronagraphy and spectroscopy. PCS will allow us not only to take images, but also to look for biosignatures such as molecular oxygen in the exoplanets' atmospheres. This article describes the PCS primary science goals, the instrument concept and the research and development activities that will be carried out over the coming years.
Matteo Simioni, Carmelo Arcidiacono, Andrea Grazian, Yann Clenet, Richard Davies, Marco Gullieuszik, Gijs Verdoes Kleijn, Fernando Pedichini, Roland Wagner, Ronny Ramlau, Werner W. Zeilinger, Fabrice Vidal, Benedetta Vulcani, Roberto Ragazzoni, Arnaud Sevin, Bernardo Salasnich, Andrea Baruffolo, Lorenzo Busoni, Simone Esposito, Éric Gendron, et al (8) The point spread function reconstruction (PSF-R) capability is a deliverable of the MICADO@ESO-ELT project. The PSF-R team works on the implementation of the instrument software devoted to reconstruct the point spread function (PSF), independently of the science data, using adaptive optics (AO) telemetry data, both for Single Conjugate (SCAO) and Multi-Conjugate Adaptive Optics (MCAO) mode of the MICADO camera and spectrograph. The PSF-R application will provide reconstructed PSFs through an archive querying system to restore the telemetry data synchronous to each science frame that MICADO will generate. Eventually, the PSF-R software will produce the output according to user specifications. The PSF-R service will support the state-of-the-art scientific analysis of the MICADO imaging and spectroscopic data.
Adaptive optics (AO) is a technology in modern ground-based optical telescopes to compensate the wavefront distortions caused by atmospheric turbulence. One method that allows to retrieve information about the atmosphere from telescope data is so-called SLODAR, where the atmospheric turbulence profile is estimated based on correlation data of Shack--Hartmann wavefront measurements. This approach relies on a layered Kolmogorov turbulence model. In this article, we propose a novel extension of the SLODAR concept by including a general non-Kolmogorov turbulence layer close to the ground with an unknown power spectral density. We prove that the joint estimation problem of the turbulence profile above ground simultaneously with the unknown power spectral density at the ground is ill-posed and propose three numerical reconstruction methods. We demonstrate by numerical simulations that our methods lead to substantial improvements in the turbulence profile reconstruction, compared to standard SLODAR-type approach. Also, our methods can accurately locate local perturbations in non-Kolmogorov power spectral densities.
The next generation ground-based telescopes rely heavily on adaptive optics for overcoming the limitation of atmospheric turbulence. In the future adaptive optics modalities, like multi-conjugate adaptive optics (MCAO), atmospheric tomography is the major mathematical and computational challenge. In this severely ill-posed problem a fast and stable reconstruction algorithm is needed that can take into account many real-life phenomena of telescope imaging. We introduce a novel reconstruction method for the atmospheric tomography problem and demonstrate its performance and flexibility in the context of MCAO. Our method is based on using locality properties of compactly supported wavelets, both in the spatial and frequency domain. The reconstruction in the atmospheric tomography problem is obtained by solving the Bayesian MAP estimator with a conjugate gradient based algorithm. An accelerated algorithm with preconditioning is also introduced. Numerical performance is demonstrated on the official end-to-end simulation tool OCTOPUS of European Southern Observatory.