The ELT is a 39m large, ground-based optical and near- to mid-infrared telescope under construction in the Chilean Atacama desert. Operation is planned to start around the middle of the next decade. All first light instruments will come with wavefront sensing devices that allow control of the ELT's intrinsic M4 and M5 wavefront correction units, thus building an adaptive optics (AO) system. To take advantage of the ELT's optical performance, full diffraction-limited operation is required and only a high performance AO system can deliver this. Further technically challenging requirements for the AO come from the exoplanet research field, where the task to resolve the very small angular separations between host star and planet, has also to take into account the high-contrast ratio between the two objects. We present in detail the results of our simulations and their impact on high-contrast imaging in order to find the optimal wavefront sensing device for the METIS instrument. METIS is the mid-infrared imager and spectrograph for the ELT with specialised high-contrast, coronagraphic imaging capabilities, whose performance strongly depends on the AO residual wavefront errors. We examined the sky and target sample coverage of a generic wavefront sensor in two spectral regimes, visible and near-infrared, to pre-select the spectral range for the more detailed wavefront sensor type analysis. We find that the near-infrared regime is the most suitable for METIS. We then analysed the performance of Shack-Hartmann and pyramid wavefront sensors under realistic conditions at the ELT, did a balancing with our scientific requirements, and concluded that a pyramid wavefront sensor is the best choice for METIS. For this choice we additionally examined the impact of non-common path aberrations, of vibrations, and the long-term stability of the SCAO system including high-contrast imaging performance.

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.