Spherical Aberration (SA) is one key component missing from our MTF toolkit for modeling an ideal imaging system’s ‘sharpness’ in the center of the field of view in the frequency domain. In this article formulas will be presented to compute the two dimensional Point Spread and Modulation Transfer Functions of the combination of diffraction, defocus and third order Spherical Aberration for an otherwise perfect lens with a circular aperture.

Spherical Aberrations result because most photographic lenses are designed with quasi spherical surfaces that do not necessarily behave ideally in all situations. For instance, they may focus light on slightly different planes depending on whether the respective ray goes through the exit pupil closer or farther from the optical axis, as shown below:

If one imagines the sensing plane being inserted parallel to the lens at various distances from it, it can be seen that in the presence of SA light from a distant star projected by the lens onto the sensing plane would not converge to a single point but to a disk – and that the size and position along the optical axis of the smallest achievable disk (the circle of least confusion) would change as the lens is stopped down, resulting in focus shift. This effect interacts with axial chromatic aberrations (spherochromatism) and tends to be the key determinant of lens ‘sharpness’ loss in the center of most in-spec, well corrected lenses at wide apertures^{[1]}.

#### SA and Defocus Point Spread Function

As figure 1 shows, with Spherical Aberration there is no longer a single plane of focus so SA and defocus interact. Wyant and Creath^{[2]} suggest that the irradiance on the sensing plane (PSF) of the diffraction pattern of a circular aperture with a rotationally symmetric pupil function illuminated by a uniform light source can be written as

(1)

with representing the following function if defocus and third order SA are present

(2)

where:

- , the radial distance from the center of the circularly symmetric PSF on the sensing plane
- , the power of light contained within the pupil
- , the wavelength of light
- , the effective f-number
- , the radial distance on the exit pupil plane from its center, normalized to one
- and , the peak Optical Path Difference from the principal plane of third order SA and Defocus respectively, in wavelengths.

A slice through the center of the PSF looks like this when varying defocus coefficient only (no SA)

When the image is in focus the PSF has the shape of an Airy disk, as expected. In that case a slice of the PSF would look as follows when Spherical Aberration coefficient is varied:

These graphs thankfully look like the ones in Wyant so apparently my Octave/Matlab implementation of the formula works^{[3]}. They represent a slice through the center of the 2D PSFs, so to obtain the full 2D PSF all one needs to do is rotate the relative radial slice around zero through 180 degrees.

#### SA and Defocus Modulation Transfer Functions

We are after a radial slice through the center of the 2D MTF, which yields the Transfer Function of the system in the direction of the slice. Taking the discrete Fourier Transform of a one dimensional slice as depicted above will not produce the desired result, even though the PSF is circularly symmetric. We have two options: either expand the radial PSF above to two dimensions by rotating it and taking the normalized modulus of the 2D DFT of that; or obtain the MTF directly as the normalized modulus of the Fourier-Bessel Transform of equation (1) as is. I opted for the latter in order to stay in FT, instead of DFT, math.

The Fourier Transform of a circularly symmetric function is itself circularly symmetric and can be obtained by applying the Fourier-Bessel Transformation (also known as the Hankel transform of zero order) to equation (1)^{[4]}. The result is a radial slice through the 2D MTF of the 2D PSF of the combination of diffraction, third order Spherical Aberrations and defocus:

(3)

with ^{[5]} the spatial frequencies of interest, indicating normalization to one at the origin, and as defined earlier. I don’t know how to solve equation (3) analytically but it’s easy enough to do it numerically in Octave/Matlab^{[3]}. Doing so produces the following graph when varying SA with zero defocus:

The first curve with zero SA and defocus is just the MTF of the diffraction pattern caused by a circular aperture. When varying defocus without SA these are the resulting curves instead:

What do you know, this last graph matches MTFs generated by Hopkins’ dedicated defocus equation explored in the last post:

#### SA and Defocus Interact to Produce the CLC

The signs of the and coefficients refer to whether the wavefront is ahead or behind the reference sphere. They interact, so when SA is present the plane of best focus is not necessarily with at zero, as can be gleaned from figure 1. The plane of best focus is the one where the diameter of the PSF on the sensing plane is smallest, the resulting ‘disk’ denoted the circle of least confusion (CLC). Wyant suggests that, geometrically, the CLC is smallest when is equal to -1.5 * . Indeed a curve with the coefficients in those proportions is the best one below:

Taking diffraction into account instead I am told the least square wavefront error should be produced when the combined optical path differences sum to zero. In terms of peak MTF50, performing an empirical simulation using the code described in the post scriptum article yields about a factor of -0.9. That is peak MTF50 is achieved when W020 = -0.9*W040.

#### Bringing It All Together

So equation (3) allows us to calculate the combined two dimensional Modulation Transfer Function of diffraction, defocus and third order Spherical Aberration for an otherwise perfect lens with a circular aperture. Together with the Transfer Functions of pixel aperture, sampling pitch and AA discussed in previous articles we should now have all the principal ingredients necessary to roughly model the linear spatial resolution (‘sharpness’) of a photographic digital camera matched to an excellent lens in the center of the field of view:

(4)

with indicating element-wise multiplication and convolution. We could for instance simulate the effect of varying f-number, pixel size and shape while taking into account these basic lens imperfections. And by analyzing the raw data of an image captured with good technique we should be able to learn a fair amount about the imaging system that took it.

PS. This article allows the simulation of complex lens imperfections and aberrations from a numerical perspective.

#### Notes and References

_{1. See this excellent paper for for a more in-depth treatment of SA.
2. Basic Wavefront Aberration Theory for Optical Metrology, James C. Wyant and Katherine Creath.
3. The Matlab/Octave scripts used to generate these plots can be downloaded from here.
4. Introduction to Fourier Optics, 3rd Edition. Joseph W Goodman, equation 2-31.
5. See here for a description of linear spatial frequency .
}