Hartmann-Shack Wavefront Sensor: Microlens Array Geometry, Spot Displacement & Phase Slope Reconstruction

A superluminescent diode ($\lambda \approx 840\ \text{nm}$) focuses a diffraction-limited spot onto the foveal retina. Light reflecting back out of the eye emerges as a planar wavefront in an ideal eye, or a warped, aberrated wavefront $W(x, y)$ in a real eye. The emerging beam passes through a microlens array (pitch $p \approx 100\text{ to } 200\ \mu\text{m}$, focal length $f_L \approx 4\text{ to } 8\ \text{mm}$):

$$\frac{\partial W(x, y)}{\partial x} = \frac{\Delta x}{f_L} \quad \text{and} \quad \frac{\partial W(x, y)}{\partial y} = \frac{\Delta y}{f_L}$$

Each lenslet focuses light onto a CCD/CMOS detector. The local gradient (slope) of the aberrated wavefront over that lenslet's sub-aperture causes the focused focal spot to displace by $(\Delta x, \Delta y)$ from its ideal optical axis coordinate.

With hundreds of discrete slope vectors measured simultaneously, mathematical integration reconstructs the continuous two-dimensional phase surface $W(x, y)$:

1. Zonal Reconstruction: Connects adjacent slope vectors numerically (spline interpolation).
2. Modal Reconstruction: Fits the measured slopes to the derivatives of orthogonal Zernike polynomials using least-squares linear regression:

$$\begin{bmatrix} \partial Z_i / \partial x \\ \partial Z_i / \partial y \end{bmatrix} \mathbf{c} = \begin{bmatrix} \Delta x / f_L \\ \Delta y / f_L \end{bmatrix} \implies \mathbf{c} = (A^T A)^{-1} A^T \mathbf{b}$$

Optical aberrometer design faces an intrinsic engineering limit: short microlens focal lengths ($f_L$) expand dynamic range to measure severe keratoconus without spot overlap, but reduce angular sensitivity. Modern aberrometers use adaptive spot tracking and variable software grids to measure from $-15.00\ \text{D}$ to $+10.00\ \text{D}$ and up to $6.00\ \text{D}$ of cylinder.