AOtools is a Python package which is open-source and aimed at providing tools for adaptive optics users and researchers. We present version 1.0 which contains tools for adaptive optics processing, including analysing data in the pupil plane, images and point spread functions in the focal plane, wavefront sensors, modelling of atmospheric turbulence, physical optical propagation of wavefronts, and conversion between frequently used adaptive optics and astronomical units. The main drivers behind AOtools is that it should be easy to install and use. To achieve this the project features extensive documentation, automated unit testing and is registered on the Python Package Index. AOtools is under continuous active development to expand the features available and we encourage everyone involved in adaptive optics to become involved and contribute to the project.

Closed-loop adaptive optics systems which use minimum mean square error wavefront reconstruction require the computation of pseudo open loop wavefront slopes. These techniques incorporate a knowledge of atmospheric statistics which must therefore be represented within the wavefront slope measurements. These pseudo open loop slopes are computed from the sum of the measured residual slopes and the reconstructed slopes that would be given if the deformable mirror was flat, generally involving the multiplication of an interaction matrix with actuator demands from the previous time-step. When using dense algebra, this multiplication is computationally expensive for Extremely Large Telescopes, requiring a large memory bandwidth. Here we show that this requirement can be significantly reduced, maintaining mathematical correctness and significantly reducing system complexity. This therefore reduces the cost of these systems and increases robustness and reliability.

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 optical turbulence profile is important in adaptive optics (AO) systems, particularly tomographic AO systems such as those to be employed by the next generation of 40 m class extremely large telescopes (ELTs). Site characterisation and monitoring campaigns have produced large quantities of turbulence profiling data for sites around the world. However AO system design and performance characterisation is dependent on Monte-Carlo simulations that cannot make use of these large datasets due to long computation times. Here we address the question of how to reduce these large datasets into small sets of profiles that can feasibly be used in such Monte-Carlo simulations, whilst minimising the loss of information inherent in this effective compression of the data. We propose hierarchical clustering to partition the dataset according to the structure of the turbulence profiles and extract a single profile from each cluster. This method is applied to the Stereo-SCIDAR dataset from ESO Paranal containing over 10000 measurements of the turbulence profile from 83 nights. We present two methods of extracting turbulence profiles from the clusters, resulting in two sets of 18 profiles providing subtly different descriptions of the variability across the entire dataset. For generality we choose integrated parameters of the turbulence to measure the representativeness of our profiles and compare to others. Using these criterion we also show that such variability is difficult to capture with small sets of profiles associated with integrated turbulence parameters such as seeing.

Shack-Hartmann wavefront sensors for both solar and laser guide star adaptive optics (with elongated spots) need to observe extended objects. Correlation techniques have been successfully employed to measure the wavefront gradient in solar adaptive optics systems and have been proposed for laser guide star systems. In this paper we describe a method for synthesising reference images for correlation Shack-Hartmann wavefront sensors with a larger field of view than individual sub-apertures. We then show how these supersized reference images can increase the performance of correlation wavefront sensors in regimes where large relative shifts are induced between sub-apertures, such as those observed in open-loop wavefront sensors. The technique we describe requires no external knowledge outside of the wavefront-sensor images, making it available as an entirely "software" upgrade to an existing adaptive optics system. For solar adaptive optics we show the supersized reference images extend the magnitude of shifts which can be accurately measured from 12% to 50% of the field of view of a sub-aperture and in laser guide star wavefront sensors the magnitude of centroids that can be accurately measured is increased from 12% to 25% of the total field of view of the sub-aperture.

One of the most useful techniques in astronomical instrumentation is image slicing. It enables a spectrograph to have a more compact angular slit, whilst retaining throughput and increasing resolving power. Astrophotonic components like the photonic lanterns and photonic reformatters can be used to replace bulk optics used so far. This study investigates the performance of such devices using end-to-end simulations to approximate realistic on-sky conditions. It investigates existing components, tries to optimize their performance and aims to understand better how best to design instruments to maximize their performance. This work complements the recent work in the field and provides an estimation for the performance of the new components.

Knowledge of the Earth's atmospheric optical turbulence is critical for astronomical instrumentation. Not only does it enable performance verification and optimisation of existing systems but it is required for the design of future instruments. As a minimum this includes integrated astro-atmospheric parameters such as seeing, coherence time and isoplanatic angle, but for more sophisticated systems such as wide field adaptive optics enabled instrumentation the vertical structure of the turbulence is also required. Stereo-SCIDAR is a technique specifically designed to characterise the Earth's atmospheric turbulence with high altitude resolution and high sensitivity. Together with ESO, Durham University has commissioned a Stereo-SCIDAR instrument at Cerro Paranal, Chile, the site of the Very Large Telescope (VLT), and only 20~km from the site of the future Extremely Large Telescope (ELT). Here we provide results from the first 18 months of operation at ESO Paranal including instrument comparisons and atmospheric statistics. Based on a sample of 83 nights spread over 22 months covering all seasons, we find the median seeing to be 0.64" with 50% of the turbulence confined to an altitude below 2 km and 40% below 600 m. The median coherence time and isoplanatic angle are found as 4.18 ms and 1.75" respectively. A substantial campaign of inter-instrument comparison was also undertaken to assure the validity of the data. The Stereo-SCIDAR profiles (optical turbulence strength and velocity as a function of altitude) have been compared with the Surface-Layer SLODAR, MASS-DIMM and the ECMWF weather forecast model. The correlation coefficients are between 0.61 (isoplanatic angle) and 0.84 (seeing).

Image slicing is a powerful technique in astronomy. It allows the instrument designer to reduce the slit width of the spectrograph, increasing spectral resolving power whilst retaining throughput. Conventionally this is done using bulk optics, such as mirrors and prisms, however more recently astrophotonic components known as photonic lanterns (PLs) and photonic reformatters have also been used. These devices reformat the multi-mode (MM) input light from a telescope into single-mode (SM) outputs, which can then be re-arranged to suit the spectrograph. The photonic dicer (PD) is one such device, designed to reduce the dependence of spectrograph size on telescope aperture and eliminate modal noise. We simulate the PD, by optimising the throughput and geometrical design using Soapy and BeamProp. The simulated device shows a transmission between 8 and 20 %, depending upon the type of adaptive optics (AO) correction applied, matching the experimental results well. We also investigate our idealised model of the PD and show that the barycentre of the slit varies only slightly with time, meaning that the modal noise contribution is very low when compared to conventional fibre systems. We further optimise our model device for both higher throughput and reduced modal noise. This device improves throughput by 6.4 % and reduces the movement of the slit output by 50%, further improving stability. This shows the importance of properly simulating such devices, including atmospheric effects. Our work complements recent work in the field and is essential for optimising future photonic reformatters.

The performance of adaptive optics systems is partially dependant on the algorithms used within the real-time control system to compute wavefront slope measurements. We demonstrate use of a matched filter algorithm for the processing of elongated laser guide star (LGS) Shack-Hartmann images, using the CANARY adaptive optics instrument on the 4.2m William Herschel Telescope and the European Southern Observatory Wendelstein LGS Unit placed 40m away. This algorithm has been selected for use with the forthcoming Thirty Meter Telescope, but until now had not been demonstrated on-sky. From the results of a first observing run, we show that the use of matched filtering improves our adaptive optics system performance, with increases in on-sky H-band Strehl measured up to about a factor of 1.1 with respect to a conventional centre of gravity approach. We describe the algorithm used, and the methods that we implemented to enable on-sky demonstration.

As telescopes become larger, into the era of ~40 m Extremely Large Telescopes, the high- resolution vertical profile of the optical turbulence strength is critical for the validation, optimization and operation of optical systems. The velocity of atmospheric optical turbulence is an important parameter for several applications including astronomical adaptive optics systems. Here, we compare the vertical profile of the velocity of the atmospheric wind above La Palma by means of a comparison of Stereo-SCIntillation Detection And Ranging (Stereo- SCIDAR) with the Global Forecast System models and nearby balloon-borne radiosondes. We use these data to validate the automated optical turbulence velocity identification from the Stereo-SCIDAR instrument mounted on the 2.5 m Isaac Newton Telescope, La Palma. By comparing these data we infer that the turbulence velocity and the wind velocity are consistent and that the automated turbulence velocity identification of the Stereo-SCIDAR is precise. The turbulence velocities can be used to increase the sensitivity of the turbulence strength profiles, as weaker turbulence that may be misinterpreted as noise can be detected with a velocity vector. The turbulence velocities can also be used to increase the altitude resolution of a detected layer, as the altitude of the velocity vectors can be identified to a greater precision than the native resolution of the system. We also show examples of complex velocity structure within a turbulent layer caused by wind shear at the interface of atmospheric zones.

An important factor which affects performance of solar adaptive optics (AO) systems is the accuracy of tracking an extended object in the wavefront sensor. The accuracy of a centre-ofmass approach to image shift measurement depends on the parameters applied in thresholding the recorded image; however, there exists no analytical prediction for these parameters for extended objects. Motivated by this we present a new method for exploring the parameter space of image shift measurement algorithms, and apply this to optimize the parameters of the algorithm. Using a thresholded, windowed centre of mass, we are able to improve centroid accuracy compared to the typical parabolic fitting approach by a factor of 3 in a signal-to-noise regime typical for solar AO. Exploration of the parameters occurs after initial image crosscorrelation with a reference image, so does not require regeneration of correlation images. The results presented employ methods which can be used in real-time to estimate the error on centroids, allowing the system to use real data to optimize parameters, without needing to enter a separate calibration mode. This method can also be applied outside of solar AO to any field which requires the tracking of an extended object.