Here we describe a simple, efficient, and most importantly fully operational point-spread-function(PSF)-reconstruction approach for laser-assisted ground layer adaptive optics (GLAO) in the frame of the Multi Unit Spectroscopic Explorer (MUSE) Wide Field Mode. Based on clear astrophysical requirements derived by the MUSE team and using the functionality of the current ESO Adaptive Optics Facility we aim to develop an operational PSF-reconstruction (PSFR) algorithm and test it both in simulations and using on-sky data. The PSFR approach is based on a Fourier description of the GLAO correction to which the specific instrumental effects of MUSE Wide Field Mode (pixel size, internal aberrations, etc.) have been added. It was first thoroughly validated with full end-to-end simulations. Sensitivity to the main atmospheric and AO system parameters was analysed and the code was re-optimised to account for the sensitivity found. Finally, the optimised algorithm was tested and commissioned using more than one year of on-sky MUSE data. We demonstrate with an on-sky data analysis that our algorithm meets all the requirements imposed by the MUSE scientists, namely an accuracy better than a few percent on the critical PSF parameters including full width at half maximum and global PSF shape through the kurtosis parameter of a Moffat function. The PSFR algorithm is publicly available and is used routinely to assess the MUSE image quality for each observation. It can be included in any post-processing activity which requires knowledge of the PSF.

Point spread function (PSF) reconstruction (PSF-R) is a well established technique to determine reliably and accurately the PSF from Adaptive Optics (AO) control loop data. We have successfully applied this technique to improve the precision on photometry and astrometry to observation of NGC6121 obtained with SPHERE/ZIMPOL as it will be presented in a forthcoming letter. Firstly, we present the methodology we followed to reconstruct the PSF combining pupil-plane and focal-plane measurements using using our PSF-R method PRIME (Beltramo-Martin et al. 2019), with upgrade of both the model and best-fitting steps compared to previous papers. Secondly, we highlight that PRIME allows to maintain the PSF fitting residual below 0.2% over 2 hours of observation and using only 30 s of AO telemetry, which may have important consequences for telemetry storage for PSF-R purpose on future 30-40 m class telescopes. Finally, we deploy PRIME in a more realistic regime using faint stars so as to identify the precision needed on the initial guess parameters to ensure the convergence towards the optimal solution.

Precise photometric and astrometric measurements on astronomical images require an accurate knowledge of the Point Spread Function (PSF). When the PSF cannot be modelled directly from the image, PSF-reconstruction techniques become the only viable solution. So far, however, their performance on real observations has rarely been quantified. Aims. In this Letter, we test the performance of a novel hybrid technique, called PRIME, on Adaptive Optics-assisted SPHERE/ZIMPOL observations of the Galactic globular cluster NGC6121. Methods. PRIME couples PSF-reconstruction techniques, based on control-loop data and direct image fitting performed on the only bright point-like source available in the field of view of the ZIMPOL exposures, with the aim of building the PSF model. Results. By exploiting this model, the magnitudes and positions of the stars in the field can be measured with an unprecedented precision, which surpasses that obtained by more standard methods by at least a factor of four for on-axis stars and by up to a factor of two on fainter, off-axis stars. Conclusions. Our results demonstrate the power of PRIME in recovering precise magnitudes and positions when the information directly coming from astronomical images is limited to only a few point-like sources and, thus, paving the way for a proper analysis of future Extremely Large Telescope observations of sparse stellar fields or individual extragalactic objects.

Context. Adaptive optics (AO) systems greatly increase the resolution of large telescopes, but produce complex point spread function (PSF) shapes, varying in time and across the field of view. This PSF must be accurately known since it provides crucial information about optical systems for design, characterisation, diagnostics and image post processing. Aims. We develop here a model of the AO long exposure PSF, adapted to various seeing conditions and any AO system. This model is made to match accurately both the core of the PSF and its turbulent halo. Methods. The PSF model we develop is based on a parsimonious parameterization of the phase power spectral density with only five parameters to describe circularly symmetric PSFs and seven parameters for asymmetrical ones. Moreover, one of the parameters is directly the Fried parameter r0 of the turbulence s strength. This physical parameter is an asset in the PSF model since it can be correlated with external measurements of the r0, such as phase slopes from the AO real time computer (RTC) or site seeing monitoring. Results. We fit our model against endtoend simulated PSFs using OOMAO tool, and against on sky PSFs from the SPHERE ZIMPOL imager and the MUSE integral field spectrometer working in AO narrowfield mode. Our model matches the shape of the AO PSF both in the core and the halo, with a sub 1 percent relative error for simulated and experimental data. We also show that we retrieve the r0 parameter with subcentimeter precision on simulated data. For ZIMPOL data, we show a correlation of 97 percent between our r0 estimation and the RTC estimation. Finally, MUSE allows us to test the spectral dependency of the fitted r0 parameter. It follows the theoretical $\lambda^{6/5}$ evolution with a standard deviation of 0.3 cm. Evolution of other PSF parameters, such as residual phase variance or aliasing, is also discussed.

Atmospheric profiling is a requirement for controlling wide-field Adaptive Optics (AO) instruments, analyzing the AO performance with respect to the observing conditions and predicting the Point Spread Function (PSF) spatial variations. We present PEPITO, a new concept for profiling atmospheric turbulence from {\em post~facto} tip-tilt (TT) corrected short-exposure images. PEPITO utilizes the anisokinetism effect in the images between several stars separated from a reference star, and then produces the profile estimation using a model-fitting methodology, by fitting to the long exposure TT-corrected PSF. PEPITO has a high sensitivity to both $C_n^2(h)$ and $L_0(h)$ by relying on the full telescope aperture and a large field of view. It then obtains a high vertical resolution (1\,m-400\,m) configurable by the camera pixel scale, taking advantage of fast statistical convergence (of order of tens of seconds). With only a short exposure-capable large format detector and a numerical complexity independent of the telescope diameter, PEPITO perfectly suits accurate profiling for night optical turbulence site characterization or adaptive optics instruments operations. We demonstrate, in simulation, that the $C_n^2(h)$ and $L_0(h)$ can be estimated to better than 1\% accuracy, from fitted PSFs of magnitude V=11 on a D=0.5\,m telescope with a 10 arcmin field of view.

In order to enhance accuracy of astrophysical estimates obtained on Adaptive-optics (AO) images, such as photometry and astrometry, we investigate a new concept to constrain the Point Spread Function (PSF) model called PSF Reconstruction and Identification for Multi-sources characterization Enhancement (PRIME), that handles jointly the science image and the AO control loop data. We present in this paper the concept of PRIME and validate it on Keck II telescope NIRC2 images. We show that by calibrating the PSF model over the scientific image, PSF reconstruction achieves 1\% and 3 mas of accuracy on respectively the Strehl-ratio and the PSF full width at half maximum. We show on NIRC2 binary images that PRIME is sufficiently robust to noise to retain photometry and astrometry precision below 0.005 mag and 100$\mu$as on a $m_H=$ 14 mag object. Finally, we also validate that PRIME performs a PSF calibration on the triple system Gl569BAB which provides a separation of 66.73$\pm 1.02$ and a differential photometry of 0.538$\pm 0.048$, compared to the reference values obtained with the extracted PSF which are 66.76 mas $\pm$ 0.94 and 0.532 mag $\pm$ 0.041.

Advanced adaptive optics (AO) instruments on ground-based telescopes require accurate knowledge of the atmospheric turbulence strength as a function of altitude. This information assists point spread function reconstruction, AO temporal control techniques and is required by wide-field AO systems to optimize the reconstruction of an observed wavefront. The variability of the atmosphere makes it important to have a measure of the optical turbulence profile in real time. This measurement can be performed by fitting an analytically generated covariance matrix to the cross-covariance of Shack-Hartmann wavefront sensor (SHWFS) centroids. In this study we explore the benefits of reducing cross-covariance data points to a covariance map region of interest (ROI). A technique for using the covariance map ROI to measure and compensate for SHWFS misalignments is also introduced. We compare the accuracy of covariance matrix and map ROI optical turbulence profiling using both simulated and on-sky data from CANARY, an AO demonstrator on the 4.2 m William Herschel telescope, La Palma. On-sky CANARY results are compared to contemporaneous profiles from Stereo-SCIDAR - a dedicated high-resolution optical turbulence profiler. It is shown that the covariance map ROI optimizes the accuracy of AO telemetry optical turbulence profiling. In addition, we show that the covariance map ROI reduces the fitting time for an extremely large telescope-scale system by a factor of 72. The software package we developed to collect all of the presented results is now open source.

Knowledge of the atmospheric turbulence in the telescope line-of-sight is crucial for wide-field observations assisted by adaptive optics (AO), for which the Point Spread Function (PSF) becomes strongly elongated due to the anisoplanatism effect. This one must be modelled accurately to extrapolate the PSF anywhere across the Field of view (FOV) and improve the science exploitation. However, anisoplanatism is a function of the Cn2(h) profile, that is not directly accessible from single conjugate AO telemetry. One may rely on external profilers, but recent studies have highlighted more than 10% of discrepancies with AO internal measurements, while we aim better than 1% of accuracy for PSF modelling. To tackle this existing limitation, we present the Focal plane profiling (FPP), as a $C_n^2(h)$ profiling method that relies on post-AO focal plane images. We demonstrate such an approach complies with a 1%-level of accuracy on the $C_n^2(h)$ estimation and establish how this accuracy varies regarding the calibration stars magnitudes and positions in the field. We highlight that photometry and astrometry errors due to PSF mis-modelling reaches respectively 1% and 50{\mu}as using FPP on a Keck baseline, with a preliminary calibration using a star of magnitude H=14 at 20". We also validate this concept using Canadas NRC-Herzberg HENOS testbed images in comparing FPP retrieval with alternative $C_n^2(h)$ measurements on HeNOS. The FPP approach allows to profile the $C_n^2(h)$ using the SCAO systems and improve significantly the PSF characterisation. Such a methodology is also ELT-size compliant and will be extrapolated to tomographic systems in a near future.

Adaptive optics (AO) restore the angular resolution of ground-based telescopes, but at the cost of delivering a time- and space-varying point spread function (PSF) with a complex shape. PSF knowledge is crucial for breaking existing limits on the measured accuracy of photometry and astrometry in science observations. In this paper, we concentrate our analyses on anisoplanatism signature only onto PSF: for large-field observations (20") with single- conjugated AO, PSFs are strongly elongated due to anisoplanatism that manifests itself as three different terms for Laser-guide star (LGS) systems: angular, focal and tilt. We propose a generalized model that relies on a point-wise decomposition of the phase and encompasses the non-stationarity of LGS systems. We demonstrate it is more accurate and less computationally demanding than existing models: it agrees with end-to-end physical-optics simulations to within 0.1% of PSF measurables, such as Strehl-ratio, FWHM and fraction of variance unexplained.Secondly, we study off-axis PSF modelling is with respect to $C_n^2(h)$ profile (heights and fractional weights). For 10m class telescope, PSF morphology is estimated at 1% level as long as we model the atmosphere with at least 7 layers whose heights and weights are known respectively with 200m and 10% precision. As a verification test we used the Canada's NRC-Herzberg HeNOS testbed data, featuring four lasers. We highlight capability of retrieving off-axis PSF characteristics within 10% of fraction of variance unexplained, which complies with the expected range from the sensitivity analysis. Our new off-axis PSF modelling method lays the ground-work for testing on-sky in the near future.