In this paper, we address the inverse problem of fast, stable, and high-quality wavefront reconstruction from pyramid wavefront sensor data for Adaptive Optics systems on Extremely Large Telescopes. For solving the indicated problem we apply well-known iterative mathematical algorithms, namely conjugate gradient, steepest descent, Landweber, Landweber-Kaczmarz and steepest descent-Kaczmarz iteration based on theoretical studies of the pyramid wavefront sensor. We compare the performance (in terms of correction quality and speed) of these algorithms in end-to-end numerical simulations of a closed adaptive loop. The comparison is performed in the context of a high-order SCAO system for METIS, one of the first-light instruments currently under design for the Extremely Large Telescope. We show that, though being iterative, the analyzed algorithms, when applied in the studied context, can be implemented in a very efficient manner, which reduces the related computational effort significantly. We demonstrate that the suggested analytically developed approaches involving iterative algorithms provide comparable quality to standard matrix-vector-multiplication methods while being computationally cheaper.

We consider the mathematical background of the wavefront sensor type that is widely used in Adaptive Optics systems for astronomy, microscopy, and ophthalmology. The theoretical analysis of the pyramid sensor forward operators presented in this paper is aimed at a subsequent development of fast and stable algorithms for wavefront reconstruction from data of this sensor type. In our analysis we allow the sensor to be utilized in both the modulated and non-modulated fashion. We derive detailed mathematical models for the pyramid sensor and the physically simpler roof wavefront sensor as well as their various approximations. Additionally, we calculate adjoint operators which build preliminaries for the application of several iterative mathematical approaches for solving inverse problems such as gradient based algorithms, Landweber iteration or Kaczmarz methods.

Accurate and robust wavefront reconstruction methods for pyramid wavefront sensors are in high demand as these sensors are planned to be part of many instruments currently under development for ground based telescopes. The pyramid sensor relates the incoming wavefront and its measurements in a non-linear way. Nevertheless, almost all existing reconstruction algorithms are based on a linearization of the model. The assumption of a linear pyramid sensor response is justifiable in closed loop AO when the measured phase information is small but may not be reasonable in reality due to unpreventable errors depending on the system such as non common path aberrations. In order to solve the non-linear inverse problem of wavefront reconstruction from pyramid sensor data we introduce two new methods based on the non-linear Landweber and Landweber-Kaczmarz iteration. Using these algorithms we experience high-quality wavefront estimation especially for the non-modulated sensor by still keeping the numerical effort feasible for large-scale AO systems.

MICADO will enable the ELT to perform diffraction limited near-infrared observations at first light. The instrument's capabilities focus on imaging (including astrometric and high contrast) as well as single object spectroscopy. This contribution looks at how requirements from the observing modes have driven the instrument design and functionality. Using examples from specific science cases, and making use of the data simulation tool, an outline is presented of what we can expect the instrument to achieve.