Efficient computation of 2D point-spread functions for diffractive lenses

Diffractive lenses, such as Fresnel zone plates, photon sieves, and their modified versions, have been of significant recent interest in high-resolution imaging applications. As the advent of diffractive lens systems with different configurations expands, the fast and accurate simulation of these systems becomes crucial for both the design and image reconstruction tasks. Here we present a fast and accurate method for computing the 2D point-spread function (PSF) of an arbitrary diffractive lens. The method is based on the recently derived closed-form mathematical formula for the PSF and the transfer function of a diffractive lens. In the method, first, the samples of the transfer function are computed using the transmittance function of the diffractive lens, and then the inverse Fourier transform of this transfer function is computed to obtain the PSF. For accurate computation, the selection of the sampling parameters is handled with care, and simple selection rules are provided for this purpose. The developed method requires a single fast Fourier transform, and, therefore, has little computational complexity. Moreover, it is also applicable to any diffractive lens configuration with arbitrary-shaped structures and modulation. As a result, this fast and accurate PSF computation method enables efficient simulation, analysis, and development of diffractive lens systems under both focused and defocused settings.