Spatial-temporal-covariance-based modeling, analysis, and simulation of aero-optics wavefront aberrations.

We introduce a framework for modeling, analysis, and simulation of aero-optics wavefront aberrations that is based on spatial-temporal covariance matrices extracted from wavefront sensor measurements. Within this framework, we present a quasi-homogeneous structure function to analyze nonhomogeneous, mildly anisotropic spatial random processes, and we use this structure function to show that phase aberrations arising in aero-optics are, for an important range of operating parameters, locally Kolmogorov. This strongly suggests that the d5/3 power law for adaptive optics (AO) deformable mirror fitting error, where d denotes actuator separation, holds for certain important aero-optics scenarios. This framework also allows us to compute bounds on AO servo lag error and predictive control error. In addition, it provides us with the means to accurately simulate AO systems for the mitigation of aero-effects, and it may provide insight into underlying physical processes associated with turbulent flow. The techniques introduced here are demonstrated using data obtained from the Airborne Aero-Optics Laboratory.

Design of an integrated hardware interface for AOSLO image capture and cone-targeted stimulus delivery.

We demonstrate an integrated FPGA solution to project highly stabilized, aberration-corrected stimuli directly onto the retina by means of real-time retinal image motion signals in combination with high speed modulation of a scanning laser. By reducing the latency between target location prediction and stimulus delivery, the stimulus location accuracy, in a subject with good fixation, is improved to 0.15 arcminutes from 0.26 arcminutes in our earlier solution. We also demonstrate the new FPGA solution is capable of delivering stabilized large stimulus pattern (up to 256 x 256 pixels) to the retina.

Retinally stabilized cone-targeted stimulus delivery.

We demonstrate projection of highly stabilized, aberration-corrected stimuli directly onto the retina by means of real-time retinal image motion signals in combination with high speed modulation of a scanning laser. In three subjects with good fixation stability, stimulus location accuracy averaged 0.26 arcminutes or approximately 1.3 microns, which is smaller than the cone-to-cone spacing at the fovea. We also demonstrate real-time correction for image distortions in adaptive optics scanning laser ophthalmoscope (AOSLO) with an intraframe accuracy of about 7 arcseconds.

Fourier domain preconditioned conjugate gradient algorithm for atmospheric tomography.

By 'atmospheric tomography' we mean the estimation of a layered atmospheric turbulence profile from measurements of the pupil-plane phase (or phase gradients) corresponding to several different guide star directions. We introduce what we believe to be a new Fourier domain preconditioned conjugate gradient (FD-PCG) algorithm for atmospheric tomography, and we compare its performance against an existing multigrid preconditioned conjugate gradient (MG-PCG) approach. Numerical results indicate that on conventional serial computers, FD-PCG is as accurate and robust as MG-PCG, but it is from one to two orders of magnitude faster for atmospheric tomography on 30 m class telescopes. Simulations are carried out for both natural guide stars and for a combination of finite-altitude laser guide stars and natural guide stars to resolve tip-tilt uncertainty.

Modeling, simulation, and open-loop control of a continuous facesheet MEMS deformable mirror.

We introduce a coupled system of nonlinear partial differential equations (PDEs) to model a particular microelectromechanical systems (MEMS) deformable mirror consisting of a continuous facesheet coupled with electrostatic plate-driven actuators. To make the problem computationally tractable, we reduced the system to a single linear PDE for the facesheet coupled with nonlinear algebraic constraints for each of the actuators. We also introduce a nonlinearly constrained quadratic optimization problem for open-loop control of the MEMS mirror. Numerical simulation and control results are presented, and shortcomings of the model are discussed.

Retinal motion estimation in adaptive optics scanning laser ophthalmoscopy.

We apply a novel computational technique known as the map-seeking circuit algorithm to estimate the motion of the retina of eye from a sequence of frames of data from a scanning laser ophthalmoscope. We also present a scheme to dewarp and co-add frames of retinal image data, given the estimated motion. The motion estimation and dewarping techniques are applied to data collected from an adaptive optics scanning laser ophthalmoscopy.

Fast optimal wavefront reconstruction for multi-conjugate adaptive optics using the Fourier domain preconditioned conjugate gradient algorithm.

We present two different implementations of the Fourier domain preconditioned conjugate gradient algorithm (FD-PCG) to efficiently solve the large structured linear systems that arise in optimal volume turbulence estimation, or tomography, for multi-conjugate adaptive optics (MCAO). We describe how to deal with several critical technical issues, including the cone coordinate transformation problem and sensor subaperture grid spacing. We also extend the FD-PCG approach to handle the deformable mirror fitting problem for MCAO.

Preconditioned conjugate gradient wave-front reconstructors for multiconjugate adaptive optics.

Multiconjugate adaptive optics (MCAO) systems with 10(4)-10(5) degrees of freedom have been proposed for future giant telescopes. Using standard matrix methods to compute, optimize, and implement wavefront control algorithms for these systems is impractical, since the number of calculations required to compute and apply the reconstruction matrix scales respectively with the cube and the square of the number of adaptive optics degrees of freedom. We develop scalable open-loop iterative sparse matrix implementations of minimum variance wave-front reconstruction for telescope diameters up to 32 m with more than 10(4) actuators. The basic approach is the preconditioned conjugate gradient method with an efficient preconditioner, whose block structure is defined by the atmospheric turbulent layers very much like the layer-oriented MCAO algorithms of current interest. Two cost-effective preconditioners are investigated: a multigrid solver and a simpler block symmetric Gauss-Seidel (BSGS) sweep. Both options require off-line sparse Cholesky factorizations of the diagonal blocks of the matrix system. The cost to precompute these factors scales approximately as the three-halves power of the number of estimated phase grid points per atmospheric layer, and their average update rate is typically of the order of 10(-2) Hz, i.e., 4-5 orders of magnitude lower than the typical 10(3) Hz temporal sampling rate. All other computations scale almost linearly with the total number of estimated phase grid points. We present numerical simulation results to illustrate algorithm convergence. Convergence rates of both preconditioners are similar, regardless of measurement noise level, indicating that the layer-oriented BSGS sweep is as effective as the more elaborated multiresolution preconditioner.

Numerical simulations of multiconjugate adaptive optics wave-front reconstruction on giant telescopes.

We present sample Monte Carlo simulation results to illustrate the trends in multiconjugate adaptive optics (MCAO) performance as the telescope aperture diameter increases from 8 to 32 m with all other first-order system parameters held constant. The MCAO system considered includes three deformable mirrors, a 1-arc min square field of view, and five wave-front-sensing references consisting of either natural guide stars or laser guide stars at a range of either 30 or 90 km. The rms residual wave-front error decreases slowly with increasing aperture diameter with natural guide stars, whereas performance degrades significantly with increasing aperture diameter for laser guide stars at 30 km if the number of guide stars is held fixed. Performance with laser guide stars at 90 km is a weak function of telescope aperture diameter in the range from 8 to 32 m, with rms wave-front errors no more than 20% greater than the corresponding natural guide-star case for the same level of wave-front sensor's measurement noise.

Multigrid preconditioned conjugate-gradient method for large-scale wave-front reconstruction.

We introduce a multigrid preconditioned conjugate-gradient (MGCG) iterative scheme for computing open-loop wave-front reconstructors for extreme adaptive optics systems. We present numerical simulations for a 17-m class telescope with n = 48756 sensor measurement grid points within the aperture, which indicate that our MGCG method has a rapid convergence rate for a wide range of subaperture average slope measurement signal-to-noise ratios. The total computational cost is of order n log n. Hence our scheme provides for fast wave-front simulation and control in large-scale adaptive optics systems.