Extremely Large Telescopes have overwhelmingly opted for the Pyramid wavefront sensor (PyWFS) over the more widely used Shack-Hartmann WaveFront Sensor (SHWFS) to perform their Single Conjugate Adaptive Optics (SCAO) mode. The PyWFS, a sensor based on Fourier filtering, has proven to be highly successful in many astronomy applications. However, it exhibits non-linearity behaviors that lead to a reduction of its sensitivity when working with non-zero residual wavefronts. This so-called Optical Gains (OG) effect, degrades the close loop performance of SCAO systems and prevents accurate correction of Non-Common Path Aberrations (NCPA). In this paper, we aim at computing the OG using a fast and agile strategy in order to control the PyWFS measurements in adaptive optics closed loop systems. Using a novel theoretical description of the PyFWS, which is based on a convolutional model, we are able to analytically predict the behavior of the PyWFS in closed-loop operation. This model enables us to explore the impact of residual wavefront error on particular aspects such as sensitivity and associated OG. The proposed method relies on the knowledge of the residual wavefront statistics and enables automatic estimation of the current OG. End-to-End numerical simulations are used to validate our predictions and test the relevance of our approach. We demonstrate, using on non-invasive strategy, that our method provides an accurate estimation of the OG. The model itself only requires AO telemetry data to derive statistical information on atmospheric turbulence. Furthermore, we show that by only using an estimation of the current Fried parameter r_0 and the basic system-level characteristics, OGs can be estimated with an accuracy of less than 10%. Finally, we highlight the importance of OG estimation in the case of NCPA compensation. The proposed method is applied to the PyWFS.

The performance of tomographic adaptive optics systems is intrinsically linked to the vertical profile of optical turbulence. Firstly, a sufficient number of discrete turbulent layers must be reconstructed to model the true continuous turbulence profile. Secondly over the course of an observation, the profile as seen by the telescope changes and the tomographic reconstructor must be updated. These changes can be due to the unpredictable evolution of turbulent layers on meteorological timescales as short as minutes. Here we investigate the effect of changing atmospheric conditions on the quality of tomographic reconstruction by coupling fast analytical adaptive optics simulation to a large database of 10 691 high resolution turbulence profiles measured over two years by the Stereo-SCIDAR instrument at ESO Paranal, Chile. This work represents the first investigation of these effects with a large, statistically significant sample of turbulence profiles. The statistical nature of the study allows us to assess not only the degradation and variability in tomographic error with a set of system parameters (e.g. number of layers, temporal update period) but also the required parameters to meet some error threshold. In the most challenging conditions where the profile is rapidly changing, these parameters must be far more tightly constrained in order to meet this threshold. By providing estimates of these constraints for a wide range of system geometries as well as the impact of different temporal optimisation strategies we may assist the designers of tomographic AO for the ELT to dimension their systems.

Current designs for all three extremely large telescopes show the overwhelming adoption of the pyramid wavefront sensor (P-WFS) as the WFS of choice for adaptive optics (AO) systems sensing on natural guide stars (NGS) or extended objects. The key advantages of the P-WFS over the Shack-Hartmann are known and are mainly provided by the improved sensitivity (fainter NGS) and reduced sensitivity to spatial aliasing. However, robustness and tolerances of the P-WFS for the ELTs are not currently well understood. In this paper, we present simulation results for the single-conjugate AO mode of HARMONI, a visible and near-infrared integral field spectrograph for the European Extremely Large Telescope. We first explore the wavefront sensing issues related to the telescope itself; namely the island effect (i.e. differential piston) and M1 segments phasing errors. We present mitigation strategies to the island effect and their performance. We then focus on some performance optimisation aspects of the AO design to explore the impact of the RTC latency and the optical gain issues, which will in particular affect the high-contrast mode of HARMONI. Finally, we investigate the influence of the quality of glass pyramid prism itself, and of optical aberrations on the final AO performance. By relaxing the tolerances on the fabrication of the prism, we are able to reduce hardware costs and simplify integration. We show the importance of calibration (i.e. updating the control matrix) to capture any displacement of the telescope pupil and rotation of the support structure for M4. We also show the importance of the number of pixels used for wavefront sensing to relax tolerances of the pyramid prism. Finally, we present a detailed optical design of the pyramid prism, central element of the P-WFS.

Context. The wind driven halo is a feature observed within the images delivered by the latest generation of ground-based instruments equipped with an extreme adaptive optics system and a coronagraphic device, such as SPHERE at the VLT. This signature appears when the atmospheric turbulence conditions are varying faster than the adaptive optics loop can correct. The wind driven halo shows as a radial extension of the point spread function along a distinct direction (sometimes referred to as the butterfly pattern). When present, it significantly limits the contrast capabilities of the instrument and prevents the extraction of signals at close separation or extended signals such as circumstellar disks. This limitation is consequential because it contaminates the data a substantial fraction of the time: about 30% of the data produced by the VLT/SPHERE instrument are affected by the wind driven halo.Aims. This paper reviews the causes of the wind driven halo and presents a method to analyze its contribution directly from the scientific images. Its effect on the raw contrast and on the final contrast after post-processing is demonstrated.Methods. We used simulations and on-sky SPHERE data to verify that the parameters extracted with our method are capable of describing the wind driven halo present in the images. We studied the temporal, spatial and spectral variation of these parameters to point out its deleterious effect on the final contrast.Results. The data driven analysis we propose does provide information to accurately describe the wind driven halo contribution in the images. This analysis justifies why this is a fundamental limitation to the final contrast performance reached.Conclusions. With the established procedure, we will analyze a large sample of data delivered by SPHERE in order to propose, in the future, post-processing techniques tailored to remove the wind driven halo.

Advanced AO systems will likely utilise Pyramid wave-front sensors (PWFS) over the traditional Shack-Hartmann sensor in the quest for increased sensitivity, peak performance and ultimate contrast. Here, we wish to bring knowledge and quantify the PWFS theoretical limits as a means to highlight its properties and use cases. We explore forward models for the PWFS in the spatial-frequency domain for they prove quite useful since a) they emanate directly from physical-optics (Fourier) diffraction theory; b) provide a straightforward path to meaningful error breakdowns, c) allow for reconstruction algorithms with $O (n\,log(n))$ complexity for large-scale systems and d) tie in seamlessly with decoupled (distributed) optimal predictive dynamic control for performance and contrast optimisation. All these aspects are dealt with here. We focus on recent analytical PWFS developments and demonstrate the performance using both analytic and end-to-end simulations. We anchor our estimates with observed on-sky contrast on existing systems and then show very good agreement between analytical and Monte-Carlo estimates for the PWFS. For a potential upgrade of existing high-contrast imagers on 10\,m-class telescopes with visible or near-infrared PWFS, we show under median conditions at Paranal a contrast improvement (limited by chromatic and scintillation effects) of 2x-5x by replacing the wave-front sensor alone at large separations close to the AO control radius where aliasing dominates, and factors in excess of 10x by coupling distributed control with the PWFS over most of the AO control region, from small separations starting with the Inner Working Angle of typically 1-2 $\lambda/D$ to the AO correction edge (here 20 $\lambda/D$).

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.

As part of the National Science Foundation funded "Gemini in the Era of MultiMessenger Astronomy" (GEMMA) program, Gemini Observatory is developing GNAO, a widefield adaptive optics (AO) facility for Gemini-North on Maunakea, the only 8m-class open-access telescope available to the US astronomers in the northern hemisphere. GNAO will provide the user community with a queue-operated Multi-Conjugate AO (MCAO) system, enabling a wide range of innovative solar system, Galactic, and extragalactic science with a particular focus on synergies with JWST in the area of time-domain astronomy. The GNAO effort builds on institutional investment and experience with the more limited block-scheduled Gemini Multi-Conjugate System (GeMS), commissioned at Gemini South in 2013. The project involves close partnerships with the community through the recently established Gemini AO Working Group and the GNAO Science Team, as well as external instrument teams. The modular design of GNAO will enable a planned upgrade to a Ground Layer AO (GLAO) mode when combined with an Adaptive Secondary Mirror (ASM). By enhancing the natural seeing by an expected factor of two, GLAO will vastly improve Gemini North's observing efficiency for seeing-limited instruments and strengthen its survey capabilities for multi-messenger astronomy.

In this paper, we describe Fourier-based Wave Front Sensors (WFS) as linear integral operators, characterized by their Kernel. In a first part, we derive the dependency of this quantity with respect to the WFS's optical parameters: pupil geometry, filtering mask, tip/tilt modulation. In a second part we focus the study on the special case of convolutional Kernels. The assumptions required to be in such a regime are described. We then show that these convolutional kernels allow to drastically simplify the WFS's model by summarizing its behavior in a concise and comprehensive quantity called the WFS's Impulse Response. We explain in particular how it allows to compute the sensor's sensitivity with respect to the spatial frequencies. Such an approach therefore provides a fast diagnostic tool to compare and optimize Fourier-based WFSs. In a third part, we develop the impact of the residual phases on the sensor's impulse response, and show that the convolutional model remains valid. Finally, a section dedicated to the Pyramid WFS concludes this work, and illustrates how the slopes maps are easily handled by the convolutional model.

For extremely large telescopes, adaptive optics will be required to correct the Earth's turbulent atmosphere. The performance of tomographic adaptive optics is strongly dependent on the vertical distribution (profile) of this turbulence. An important way in which this manifests is the tomographic error, arising from imperfect measurement and reconstruction of the turbulent phase at altitude. Conventionally, a small number of reference profiles are used to obtain this error in simulation however these profiles are not constructed to be representative in terms of tomographic error. It is therefore unknown whether these simulations are providing realistic performance estimates. Here, we employ analytical adaptive optics simulation that drastically reduces computation times to compute tomographic error for 10 691 measurements of the turbulence profile gathered by the Stereo-SCIDAR instrument at ESO Paranal. We assess for the first time the impact of the profile on tomographic error in a statistical manner. We find, in agreement with previous work, that the tomographic error is most directly linked with the distribution of turbulence into discrete, stratified layers. Reference profiles are found to provide mostly higher tomographic error than expected, which we attribute to the fact that these profiles are primarily composed of averages of many measurements resulting in unrealistic, continuous distributions of turbulence. We propose that a representative profile should be defined with respect to a particular system, and that as such simulations with a large statistical sample of profiles must be an important step in the design process.

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.

Wavefront sensors encode phase information of an incoming wavefront into an intensity pattern that can be measured on a camera. Several kinds of wavefront sensors (WFS) are used in astronomical adaptive optics. Amongst them, Fourier-based wavefront sensors perform a filtering operation on the wavefront in the focal plane. The most well known example of a WFS of this kind is the Zernike wavefront sensor, and the pyramid wavefront sensor (PWFS) also belongs to this class. Based on this same principle, new WFSs can be proposed such as the n-faced pyramid (which ultimately becomes an axicone) or the flattened pyramid, depending on whether the image formation is incoherent or coherent. In order to test such novel concepts, the LOOPS adaptive optics testbed hosted at the Laboratoire d'Astrophysique de Marseille has been upgraded by adding a Spatial Light Modulator (SLM). This device, placed in a focal plane produces high-definition phase masks that mimic otherwise bulk optic devices. In this paper, we first present the optical design and upgrades made to the experimental setup of the LOOPS bench. Then, we focus on the generation of the phase masks with the SLM and the implications of having such a device in a focal plane. Finally, we present the first closed-loop results in either static or dynamic mode with different WFS applied on the SLM.

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.

The European Extremely Large Telescope will see first lights by the end of 2024. With a diameter of almost 40 meters, it will be the biggest optical telescope ever built from the ground. The ELT will open a brand new window in a sensitivity / spatial angular resolution parameter space. To take the full benefit of the scientific potential of this giant, all the instruments will be equipped with Adaptive Optics (AO), providing the sharpest images. This paper provides a quick overview of the AO capabilities of the future instruments to be deployed at the ELT, highlight some of the expected performance and describe a couple of technical challenges that are still to tackle for an optimal scientific use. This paper has been presented at the "Societe Francaise d'Astronomie et Astrophysique" symposium in Bordeaux 2018, it is then naturally biased toward the French contribution for the ELT.

Geared by the increasing need for enhanced performance, both optical and computational, new dynamic control laws have been researched in recent years for next generation adaptive optics systems on current 10 m-class and extremely large telescopes up to 40 m. We provide an overview of these developments and point out prospects to making such controllers drive actual systems on-sky.

The estimation of atmospheric turbulence parameters is of relevance for: a) site evaluation & characterisation; b) prediction of the point spread function; c) live assessment of error budgets and optimisation of adaptive optics performance; d) optimisation of fringe trackers for long baseline optical interferometry. The ubiquitous deployment of Shack-Hartmann wavefront sensors in large telescopes makes them central for atmospheric turbulence parameter estimation via adaptive optics telemetry. Several methods for the estimation of the Fried parameter and outer scale have been developed, most of which are based on the fitting of Zernike polynomial coefficients variances reconstructed from the telemetry. The non-orthogonality of Zernike polynomial derivatives introduces modal cross coupling, which affects the variances. Furthermore, the finite resolution of the sensor introduces aliasing. In this article the impact of these effects on atmospheric turbulence parameter estimation is addressed with simulations. It is found that cross coupling is the dominant bias. An iterative algorithm to overcome it is presented. Simulations are conducted for typical ranges of the outer scale (4 to 32m), Fried parameter (10 cm) and noise in the variances (signal-to-noise ratio of 10 and above). It is found that, using the algorithm, both parameters are recovered with sub-percent accuracy.

The latest generation of high contrast instruments dedicated to exoplanets and circumstellar disks imaging are equipped with extreme adaptive optics and coronagraphs to reach contrasts of up to 10^--4 at a few tenths of arc-seconds in the near infrared. The resulting image shows faint features, only revealed with this combination, such as the wind driven halo. The wind driven halo is due to the lag between the adaptive optics correction and the turbulence speed over the telescope pupil. However we observe an asymmetry of this wind driven halo that was not expected when the instrument was designed. In this letter, we describe and demonstrate the physical origin of this asymmetry and support our explanation by simulating the asymmetry with an end-to-end approach. From this work, we found out that the observed asymmetry is explained by the interference between the AO-lag error and scintillation effects, mainly originating from the fast jet stream layer located at about 12 km in altitude. Now identified and interpreted, this effect can be taken into account for further design of high contrast imaging simulators, next generation or upgrade of high contrast instruments, predictive control algorithms for adaptive optics or image post-processing techniques.

The segmentation of the telescope pupil (by spiders & the segmented M4) create areas of phase isolated by the width of the spiders on the wavefront sensor (WFS), breaking the spatial continuity of the wavefront. The poor sensitivity of the Pyramid WFS (PWFS) to differential piston leads to badly seen and therefore uncontrollable differential pistons. In close loop operation, differential pistons between segments will settle around integer values of the average sensing wavelength. The differential pistons typically range from one to ten times the sensing wavelength and vary rapidly over time, leading to extremely poor performance. In addition, aberrations created by atmospheric turbulence will contain large amounts of differential piston between the segments. Removing piston contribution over each of the DM segments leads to poor performance. In an attempt to reduce the impact of unwanted differential pistons that are injected by the AO correction, we compare three different approaches. We first limit ourselves to only use the information measured by the PWFS, in particular by reducing the modulation. We show that using this information sensibly is important but will not be sufficient. We discuss possible ways of improvement by using prior information. A second approach is based on phase closure of the DM commands and assumes the continuity of the correction wavefront over the entire unsegmented pupil. The last approach is based on the pair-wise slaving of edge actuators and shows the best results. We compare the performance of these methods using realistic end-to-end simulations. We find that pair-wise slaving leads to a small increase of the total wavefront error, only adding between 20-45 nm RMS in quadrature for seeing conditions between 0.45-0.85 arcsec. Finally, we discuss the possibility of combining the different proposed solutions to increase robustness.

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.

Magnetohydrodynamic(MHD) turbulence displays velocity anisotropies which reflect the direction of the magnetic field. This anisotropy has led to the development of a number of statistical techniques for studying magnetic fields in the interstellar medium. In this paper, we review and compare three techniques that use radio position-position-velocity data for determining magnetic field strength and morphology : the correlation function anisotropy (CFA), Principal Component Analysis of Anisotropies (PCAA), and the more recent Velocity Gradient Technique (VGT). We compare these three techniques and suggest improvements to the CFA and PCAA techniques to increase their accuracy and versatility. In particular, we suggest and successfully implement a much faster way of calculating non-periodic correlation functions for the CFA. We discuss possible improvements to the current implementation of the PCAA. We show the advantages of the VGT in terms of magnetic field tracing and stress the complementary nature with the other two techniques.

The Gemini Infrared Multi-Object Spectrograph (GIRMOS) is a powerful new instrument being built to facility-class standards for the Gemini telescope. It takes advantage of the latest developments in adaptive optics and integral field spectrographs. GIRMOS will carry out simultaneous high-angular-resolution, spatially-resolved infrared ($1-2.4$ $\mu$m) spectroscopy of four objects within a two-arcminute field-of-regard by taking advantage of multi-object adaptive optics. This capability does not currently exist anywhere in the world and therefore offers significant scientific gains over a very broad range of topics in astronomical research. For example, current programs for high redshift galaxies are pushing the limits of what is possible with infrared spectroscopy at $8-10$-meter class facilities by requiring up to several nights of observing time per target. Therefore, the observation of multiple objects simultaneously with adaptive optics is absolutely necessary to make effective use of telescope time and obtain statistically significant samples for high redshift science. With an expected commissioning date of 2023, GIRMOS's capabilities will also make it a key followup instrument for the James Webb Space Telescope when it is launched in 2021, as well as a true scientific and technical pathfinder for future Thirty Meter Telescope (TMT) multi-object spectroscopic instrumentation. In this paper, we will present an overview of this instrument's capabilities and overall architecture. We also highlight how this instrument lays the ground work for a future TMT early-light instrument.

Tomographic wave-front reconstruction is the main computational bottleneck to realize real-time correction for turbulence-induced wave-front aberrations in future laser-assisted tomographic adaptive-optics (AO) systems for ground-based Giant Segmented Mirror Telescopes (GSMT), because of its unprecedented number of degrees of freedom, $N$, i.e. the number of measurements from wave-front sensors (WFS). In this paper, we provide an efficient implementation of the minimum-mean-square error (MMSE) tomographic wave-front reconstruction mainly useful for some classes of AO systems not requiring a multi-conjugation, such as laser-tomographic AO (LTAO), multi-object AO (MOAO) and ground-layer AO (GLAO) systems, but also applicable to multi-conjugate AO (MCAO) systems. This work expands that by R. Conan [ProcSPIE, 9148, 91480R (2014)] to the multi-wave-front, tomographic case using natural and laser guide stars. The new implementation exploits the Toeplitz structure of covariance matrices used in a MMSE reconstructor, which leads to an overall $O(N\log N)$ real-time complexity compared to $O(N^2)$ of the original implementation using straight vector-matrix multiplication. We show that the Toeplitz-based algorithm leads to 60\,nm rms wave-front error improvement for the European Extremely Large Telescope Laser-Tomography AO system over a well-known sparse-based tomographic reconstruction, but the number of iterations required for suitable performance is still beyond what a real-time system can accommodate to keep up with the time-varying turbulence

HARMONI is a visible and NIR integral field spectrograph, providing the E-ELT's core spectroscopic capability at first light. HARMONI will work at the diffraction limit of the E-ELT, thanks to a Classical and a Laser Tomographic AO system. In this paper, we present the system choices that have been made for these SCAO and LTAO modules. In particular, we describe the strategy developed for the different Wave-Front Sensors: pyramid for SCAO, the LGSWFS concept, the NGSWFS path, and the truth sensor capabilities. We present first potential implementations. And we asses the first system performance.

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.

Shack-Hartmann wavefront sensing relies on accurate spot centre measurement. Several algorithms were developed with this aim, mostly focused on precision, i.e. minimizing random errors. In the solar and extended scene community, the importance of the accuracy (bias error due to peak-locking, quantisation or sampling) of the centroid determination was identified and solutions proposed. But these solutions only allow partial bias corrections. To date, no systematic study of the bias error was conducted. This article bridges the gap by quantifying the bias error for different correlation peak-finding algorithms and types of sub-aperture images and by proposing a practical solution to minimize its effects. Four classes of sub-aperture images (point source, elongated laser guide star, crowded field and solar extended scene) together with five types of peak-finding algorithms (1D parabola, the centre of gravity, Gaussian, 2D quadratic polynomial and pyramid) are considered, in a variety of signal-to-noise conditions. The best performing peak-finding algorithm depends on the sub-aperture image type, but none is satisfactory to both bias and random errors. A practical solution is proposed that relies on the anti-symmetric response of the bias to the sub-pixel position of the true centre. The solution decreases the bias by a factor of ~7 to values of < 0.02 pix. The computational cost is typically twice of current cross-correlation algorithms.

We build on a long-standing tradition in astronomical adaptive optics (AO) of specifying performance metrics and error budgets using linear systems modeling in the spatial-frequency domain. Our goal is to provide a comprehensive tool for the calculation of error budgets in terms of residual temporally filtered phase power spectral densities and variances. In addition, the fast simulation of AO-corrected point spread functions (PSFs) provided by this method can be used as inputs for simulations of science observations with next-generation instruments and telescopes, in particular to predict post-coronagraphic contrast improvements for planet finder systems. We extend the previous results and propose the synthesis of a distributed Kalman filter to mitigate both aniso-servo-lag and aliasing errors whilst minimizing the overall residual variance. We discuss applications to (i) analytic AO-corrected PSF modeling in the spatial-frequency domain, (ii) post-coronagraphic contrast enhancement, (iii) filter optimization for real-time wavefront reconstruction, and (iv) PSF reconstruction from system telemetry. Under perfect knowledge of wind velocities, we show that $\sim$60 nm rms error reduction can be achieved with the distributed Kalman filter embodying anti- aliasing reconstructors on 10 m class high-order AO systems, leading to contrast improvement factors of up to three orders of magnitude at few ${\lambda}/D$ separations ($\sim1-5{\lambda}/D$) for a 0 magnitude star and reaching close to one order of magnitude for a 12 magnitude star.

We propose and apply two methods to estimate pupil plane phase discontinuities for two realistic scenarios on VLT and Keck. The methods use both Phase Diversity and a form of image sharpening. For the case of VLT, we simulate the `low wind effect' (LWE) which is responsible for focal plane errors in the SPHERE system in low wind and good seeing conditions. We successfully estimate the simulated LWE using both methods, and show that they are complimentary to one another. We also demonstrate that single image Phase Diversity (also known as Phase Retrieval with diversity) is also capable of estimating the simulated LWE when using the natural de-focus on the SPHERE/DTTS imager. We demonstrate that Phase Diversity can estimate the LWE to within 30 nm RMS WFE, which is within the allowable tolerances to achieve a target SPHERE contrast of 10$^{-6}$. Finally, we simulate 153 nm RMS of piston errors on the mirror segments of Keck and produce NIRC2 images subject to these effects. We show that a single, diverse image with 1.5 waves (PV) of focus can be used to estimate this error to within 29 nm RMS WFE, and a perfect correction of our estimation would increase the Strehl ratio of a NIRC2 image by 12\%

The use of Fourier methods in wave-front reconstruction can significantly reduce the computation time for large telescopes with a high number of degrees of freedom. However, Fourier algorithms for discrete data require a rectangular data set which conform to specific boundary requirements, whereas wave-front sensor data is typically defined over a circular domain (the telescope pupil). Here we present an iterative Gerchberg routine modified for the purposes of discrete wave-front reconstruction which adapts the measurement data (wave-front sensor slopes) for Fourier analysis, fulfilling the requirements of the Fast Fourier Transform (FFT) and providing accurate reconstruction. The routine is used in the adaptation step only and can be coupled to any other Wiener-like or least-squares method. We compare simulations using this method with previous Fourier methods and show an increase in performance in terms of Strehl ratio and a reduction in noise propagation for a 40x40 SPHERE-like adaptive optics system. For closed loop operation with minimal iterations the Gerchberg method provides an improvement in Strehl, from 95.4% to 96.9% in K-band. This corresponds to ~40 nm improvement in rms, and avoids the high spatial frequency errors present in other methods, providing an increase in contrast towards the edge of the correctable band.

The chemical abundances for five metal-poor stars in and towards the Galactic bulge have been determined from H-band infrared spectroscopy taken with the RAVEN multi-object adaptive optics science demonstrator and the IRCS spectrograph at the Subaru 8.2-m telescope. Three of these stars are in the Galactic bulge and have metallicities between -2.1 < [Fe/H] < -1.5, and high [alpha/Fe] ~+0.3, typical of Galactic disk and bulge stars in this metallicity range; [Al/Fe] and [N/Fe] are also high, whereas [C/Fe] < +0.3. An examination of their orbits suggests that two of these stars may be confined to the Galactic bulge and one is a halo trespasser, though proper motion values used to calculate orbits are quite uncertain. An additional two stars in the globular cluster M22 show [Fe/H] values consistent to within 1 sigma, although one of these two stars has [Fe/H] = -2.01 +/- 0.09, which is on the low end for this cluster. The [alpha/Fe] and [Ni/Fe] values differ by 2 sigma, with the most metal-poor star showing significantly higher values for these elements. M22 is known to show element abundance variations, consistent with a multi-population scenario (i.e. Marino et al. 2009, 2011; Alves-Brito et al. 2012) though our results cannot discriminate this clearly given our abundance uncertainties. This is the first science demonstration of multi-object adaptive optics with high resolution infrared spectroscopy, and we also discuss the feasibility of this technique for use in the upcoming era of 30-m class telescope facilities.

In preparation of future Multi-Object Spectrographs (MOS) whose one of the major role is to provide an extensive statistical studies of high redshifted galaxies surveyed, the demonstrator Canary has been designed to tackle technical challenges related to open-loop Adaptive-Optics (AO) control with jointed Natural Guide Star (NGS) and Laser Guide Star (LGS) tomography. We have developed a Point Spread Function (PSF)-Reconstruction algorithm dedicated to MOAO systems using system telemetry to estimate the PSF potentially anywhere in the observed field, a prerequisite to post- process AO-corrected observations in Integral Field Spectroscopy (IFS). In this paper we show how to handle off-axis data to estimate the PSF using atmospheric tomography and compare it to a classical approach that uses on-axis residual phase from a truth sensor observing a natural bright source. We have reconstructed over 450 on-sky Canary PSFs and we get bias/1-${\sigma}$ standard-deviation (std) of 1.3/4.8 on the H-band Strehl ratio (SR) with 92.3% of correlation between reconstructed and sky SR. On the Full Width at Half Maximum (FWHM), we get respectively 2.94 mas, 19.9 mas and 88.3% for the bias, std and correlation. The reference method achieves 0.4/3.5/95% on the SR and 2.71 mas/14.9 mas/92.5% on the FWHM for the bias/std/correlation.

Canary is the Multi-Object Adaptive Optics (MOAO) pathfinder for the future MOAO-assisted Integral-Field Units (IFU) proposed for Extremely Large Telescopes (ELT). The MOAO concept relies on tomographically reconstructing the turbulence using multiple measurements along different lines of sight. Tomography requires the knowledge of the statistical turbulence parameters, commonly recovered from the system telemetry using a dedicated profiling technique. For demonstration purposes with the MOAO pathfinder Canary , this identification is performed thanks to the Learn & Apply (L&A) algorithm, that consists in model- fitting the covariance matrix of WFS measurements dependent on relevant parameters: $C_n^2(h)$ profile, outer scale profile and system mis-registration. We explore an upgrade of this algorithm, the Learn 3 Steps (L3S) approach, that allows one to dissociate the identification of the altitude layers from the ground in order to mitigate the lack of convergence of the required empirical covariance matrices therefore reducing the required length of data time-series for reaching a given accuracy. For nominal observation conditions, the L3S can reach the same level of tomographic error in using five times less data frames than the L&A approach. The L3S technique has been applied over a large amount of Canary data to characterize the turbu- lence above the William Herschel Telescope (WHT). These data have been acquired the 13th, 15th, 16th, 17th and 18th September 2013 and we find 0.67"/8.9m/3.07m/s of total seeing/outer scale/wind-speed, with 0.552"/9.2m/2.89m/s below 1.5 km and 0.263"/10.3m/5.22m/s between 1.5 and 20 km. We have also de- termined the high altitude layers above 20 km, missed by the tomographic reconstruction on Canary , have a median seeing of 0.187" and have occurred 16% of observation time.

CANARY is the multi-object adaptive optics (MOAO) on-sky pathfinder developed in the perspective of Multi-Object Spectrograph on Extremely Large Telescopes~(ELTs). In 2013, CANARY was operated on-sky at the William Herschel telescope~(WHT), using three off-axis natural guide stars~(NGS) and four off-axis Rayleigh laser guide stars~(LGS), in open-loop, with the on-axis compensated turbulence observed with a H-band imaging camera and a Truth wave-front sensor~(TS) for diagnostic purposes. Our purpose is to establish a reliable and accurate wave-front error breakdown for LGS MOAO. This will enable a comprehensive analysis of \cana on-sky results and provide tools for validating simulations of MOAO systems for ELTs. To evaluate the MOAO performance, we compared the CANARY on-sky results running in MOAO, in Single Conjugated Adaptive Optics~(SCAO) and in Ground Layer Adaptive Optics~(GLAO) modes, over a large set of data acquired in 2013. We provide a statistical study of the seeing. We also evaluated the wave-front error breakdown from both analytic computations, one based on a MOAO system modelling and the other on the measurements from the CANARY TS. We have focussed especially on the tomographic error and we detail its vertical error decomposition~(VED). We show that CANARY obtained 30.1\%, 21.4\% and 17.1\% H-band Strehl ratios in SCAO, MOAO and GLAO respectively, for median seeing conditions with 0.66" of total seeing including 0.59" at the ground. Moreover, we get 99\% of correlation over 4,500 samples, for any AO modes, between two analytic computations of residual phase variance. Based on these variances, we obtain a reasonable Strehl-ratio~(SR) estimation when compared to the measured IR image SR. We evaluate the gain in compensation for the altitude turbulence brought by MOAO when compared to GLAO.

Prior statistical knowledge of the atmospheric turbulence is essential for designing, optimizing and evaluating tomographic adaptive optics systems. We present the statistics of the vertical profiles of $C_N^2$ and the outer scale at Maunakea estimated using a Slope Detection And Ranging (SLODAR) method from on-sky telemetry taken by RAVEN, which is a MOAO demonstrator in the Subaru telescope. In our SLODAR method, the profiles are estimated by a fit of the theoretical auto- and cross-correlation of measurements from multiple Shack-Haltmann wavefront sensors to the observed correlations via the non-linear Levenberg-Marquardt Algorithm (LMA), and the analytic derivatives of the spatial phase structure function with respect to its parameters for the LMA are also developed. The estimated profile has the median total seeing of 0.460$^{\prime\prime}$ and large $C_N^2$ fraction of the ground layer of 54.3%. The $C_N^2$ profile has a good agreement with the result from literatures, except for the ground layer. The median value of the outer scale is 25.5m and the outer scale is larger at higher altitudes, and these trends of the outer scale are consistent with findings in literatures.

Sky-coverage in laser-assisted AO observations largely depends on the system's capability to guide on the faintest natural guide-stars possible. Here we give an up-to-date status of our natural guide-star processing tailored to the European-ELT's visible and near-infrared (0.47 to 2.45 {\mu}m) integral field spectrograph - Harmoni. We tour the processing of both the isoplanatic and anisoplanatic tilt modes using the spatio-angular approach whereby the wave-front is estimated directly in the pupil plane avoiding a cumbersome explicit layered estimation on the 35-layer profiles we're currently using. Taking the case of Harmoni, we cover the choice of wave-front sensors, the number and field location of guide-stars, the optimised algorithms to beat down angular anisoplanatism and the performance obtained with different temporal controllers under split high-order/low-order tomography or joint tomography. We consider both atmospheric and far greater telescope wind buffeting disturbances. In addition we provide the sky-coverage estimates thus obtained.

In tomographic adaptive-optics (AO) systems, errors due to tomographic wave-front reconstruction limit the performance and angular size of the scientific field of view (FoV), where AO correction is effective. We propose a multi time-step tomographic wave-front reconstruction method to reduce the tomographic error by using the measurements from both the current and the previous time-steps simultaneously. We further outline the method to feed the reconstructor with both wind speed and direction of each turbulence layer. An end-to-end numerical simulation, assuming a multi-object AO (MOAO) system on a 30 m aperture telescope, shows that the multi time-step reconstruction increases the Strehl ratio (SR) over a scientific FoV of 10 arcminutes in diameter by a factor of 1.5--1.8 when compared to the classical tomographic reconstructor, depending on the guide star asterism and with perfect knowledge of wind speeds and directions. We also evaluate the multi time-step reconstruction method and the wind estimation method on the RAVEN demonstrator under laboratory setting conditions. The wind speeds and directions at multiple atmospheric layers are measured successfully in the laboratory experiment by our wind estimation method with errors below 2 \ms. With these wind estimates, the multi time-step reconstructor increases the SR value by a factor of 1.2--1.5, which is consistent with a prediction from end-to-end numerical simulation.

In this work we investigate the Principal Component Analysis (PCA) sensitivity to the velocity power spectrum in high opacity regimes of the interstellar medium (ISM). For our analysis we use synthetic Position-Position-Velocity (PPV) cubes of fractional Brownian motion (fBm) and magnetohydrodynamics (MHD) simulations, post processed to include radiative transfer effects from CO. We find that PCA analysis is very different from the tools based on the traditional power spectrum of PPV data cubes. Our major finding is that PCA is also sensitive to the phase information of PPV cubes and this allows PCA to detect the changes of the underlying velocity and density spectra at high opacities, where the spectral analysis of the maps provides the universal -3 spectrum in accordance with the predictions of Lazarian \& Pogosyan (2004) theory. This makes PCA potentially a valuable tool for studies of turbulence at high opacities provided that the proper gauging of the PCA index is made. The later, however, we found to be not easy, as the PCA results change in an irregular way for data with high sonic Mach numbers. This is in contrast to synthetic Brownian noise data used for velocity and density fields that show monotonic PCA behavior. We attribute this difference to the PCA's sensitivity to Fourier phase information.

Aims:To support the computation and evolutionary interpretation of periods associated with the rotational modulation, oscillations, and variability of stars located in the CoRoT fields, we are conducting a spectroscopic survey for stars located in the fields already observed by the satellite. These observations allow us to compute physical and chemical parameters for our stellar sample. Method: Using spectroscopic observations obtained with UVES/VLT and Hydra/Blanco, and based on standard analysis techniques, we computed physical and chemical parameters ($T_{\rm{eff}}$, $\log \,(g)$, $\rm{[Fe/H]}$, $v_{\rm{mic}}$, $v_{\rm{rad}}$, $v \sin \,(i)$, and $A(\rm{Li})$) for a large sample of CoRoT targets. Results: We provide physical and chemical parameters for a sample comprised of 138 CoRoT targets. Our analysis shows the stars in our sample are located in different evolutionary stages, ranging from the main sequence to the red giant branch, and range in spectral type from F to K. The physical and chemical properties for the stellar sample are in agreement with typical values reported for FGK stars. However, we report three stars presenting abnormal lithium behavior in the CoRoT fields. These parameters allow us to properly characterize the intrinsic properties of the stars in these fields. Our results reveal important differences in the distributions of metallicity, $T_{\rm eff}$, and evolutionary status for stars belonging to different CoRoT fields, in agreement with results obtained independently from ground-based photometric surveys. Conclusions: Our spectroscopic catalog, by providing much-needed spectroscopic information for a large sample of CoRoT targets, will be of key importance for the successful accomplishment of several different programs related to the CoRoT mission, thus it will help further boost the scientific return associated with this space mission.

Multi-object astronomical adaptive-optics (MOAO) is now a mature wide-field observation mode to enlarge the adaptive-optics-corrected field in a few specific locations over tens of arc-minutes. The work-scope provided by open-loop tomography and pupil conjugation is amenable to a spatio-angular Linear-Quadratic Gaussian (SA-LQG) formulation aiming to provide enhanced correction across the field with improved performance over static reconstruction methods and less stringent computational complexity scaling laws. Starting from our previous work [1], we use stochastic time-progression models coupled to approximate sparse measurement operators to outline a suitable SA-LQG formulation capable of delivering near optimal correction. Under the spatio-angular framework the wave-fronts are never explicitly estimated in the volume,providing considerable computational savings on 10m-class telescopes and beyond. We find that for Raven, a 10m-class MOAO system with two science channels, the SA-LQG improves the limiting magnitude by two stellar magnitudes when both Strehl-ratio and Ensquared-energy are used as figures of merit. The sky-coverage is therefore improved by a factor of 5.

We use a theoretical frame-work to analytically assess temporal prediction error functions on von-Karman turbulence when a zonal representation of wave-fronts is assumed. Linear prediction models analysed include auto-regressive of order up to three, bilinear interpolation functions and a minimum mean square error predictor. This is an extension of the authors' previously published work (see ref. 2) in which the efficacy of various temporal prediction models was established. Here we examine the tolerance of these algorithms to specific forms of model errors, thus defining the expected change in behaviour of the previous results under less ideal conditions. Results show that +/- 100pc wind-speed error and +/- 50 deg are tolerable before the best linear predictor delivers poorer performance than the no-prediction case.

Context. Frequency analyses are very important in astronomy today, not least in the ever-growing field of exoplanets, where short-period signals in stellar radial velocity data are investigated. Periodograms are the main (and powerful) tools for this purpose. However, recovering the correct frequencies and assessing the probability of each frequency is not straightforward. Aims. We provide a formalism that is easy to implement in a code, to describe a Bayesian periodogram that includes weights and a constant offset in the data. The relative probability between peaks can be easily calculated with this formalism. We discuss the differences and agreements between the various periodogram formalisms with simulated examples. Methods. We used the Bayesian probability theory to describe the probability that a full sine function (including weights derived from the errors on the data values and a constant offset) with a specific frequency is present in the data. Results. From the expression for our Baysian generalised Lomb-Scargle periodogram (BGLS), we can easily recover the expression for the non-Bayesian version. In the simulated examples we show that this new formalism recovers the underlying periods better than previous versions. A Python-based code is available for the community.

Computationally-efficient wave-front reconstruction techniques for astronomical adaptive optics systems have seen a great development in the past decade. Algorithms developed in the spatial-frequency (Fourier) domain have gathered large attention specially for high-contrast imaging systems. In this paper we present the Wiener filter (resulting in the maximization of the Strehl-ratio) and further develop formulae for the anti-aliasing Wiener filter that optimally takes into account high-order wave-front terms folded in-band during the sensing (i.e. discrete sampling) process. We employ a continuous spatial-frequency representation for the forward measurement operators and derive the Wiener filter when aliasing is explicitly taken into account. We further investigate and compare to classical estimates using least-squares filters the reconstructed wave-front, measurement noise and aliasing propagation coefficients as a function of the system order. Regarding high-contrast systems, we provide achievable performance results as a function of an ensemble of for ward models for the Shack-Hartmann wave-front sensor (using sparse and non-sparse representations) and compute point-spread function raw intensities. We find that for a 32x32 single-conjugated adaptive optics system the aliasing propagation coefficient is roughly 60% of the least-squares filters whereas the noise propagation is around 80%. Contrast improvements of factors of up to 2 are achievable across the field in H-band. For current and next generation high-contrast imagers, despite better aliasing mitigation, anti-aliasing Wiener filtering cannot be used as a stand-alone method and must therefore be used in combination with optical spatial filters deployed before image formation takes actual place.

To directly image exoplanets and faint circumstellar disks, the noisy stellar halo must be suppressed to a high level. To achieve this feat, the angular differential imaging observing technique and the least-squares Locally Optimized Combination of Images (LOCI) algorithm have now become the standard in single band direct imaging observations and data reduction. With the development and commissioning of new high-order high-contrast adaptive optics equipped with integral field units, the image subtraction algorithm needs to be modified to allow the optimal use of polychromatic images, field-rotated images and archival data. A new algorithm, TLOCI (for Template LOCI), is designed to achieve this task by maximizing a companion signal-to-noise ratio instead of simply minimizing the noise as in the original LOCI algorithm. The TLOCI technique uses an input spectrum and template Point Spread Functions (PSFs, generated from unocculted and unsaturated stellar images) to optimize the reference image least-squares coefficients to minimize the planet self-subtraction, thus maximizing its throughput per wavelength, while simultaneously providing a maximum suppression of the speckle noise. The new algorithm has been developed using on-sky GPI data and has achieved impressive contrast. This paper presents the TLOCI algorithm, on-sky performance, and will discuss the challenges in recovering the planet spectrum with high fidelity.

We study how the estimation of the sonic Mach number ($M_s$) from $^{13}$CO linewidths relates to the actual 3D sonic Mach number. For this purpose we analyze MHD simulations which include post-processing to take radiative transfer effects into account. As expected, we find very good agreement between the linewidth estimated sonic Mach number and the actual sonic Mach number of the simulations for optically thin tracers. However, we find that opacity broadening causes $M_s$ to be overestimated by a factor of ~ 1.16 - 1.3 when calculated from optically thick $^{13}$CO lines. We also find that there is a dependency on the magnetic field: super-Alfv\'enic turbulence shows increased line broadening as compared with sub-Alfv\'enic turbulence for all values of optical depth for supersonic turbulence. Our results have implications for the observationally derived sonic Mach number--density standard deviation ($\sigma_{\rho/<\rho>}$) relationship, $\sigma^2_{\rho/<\rho>}=b^2M_s^2$, and the related column density standard deviation ($\sigma_{N/<N>}$) sonic Mach number relationship. In particular, we find that the parameter b, as an indicator of solenoidal vs. compressive driving, will be underestimated as a result of opacity broadening. We compare the $\sigma_{N/<N>}-M_s$ relation derived from synthetic dust extinction maps and $^{13}$CO linewidths with recent observational studies and find that solenoidally driven MHD turbulence simulations have values of $\sigma_{N/<N>}$ which are lower than real molecular clouds. This may be due to the influence of self-gravity which should be included in simulations of molecular cloud dynamics.