Abstract
An efficient algorithm for calculating nonparaxial scalar field distributions in the focal region of a lens is discussed. The algorithm is based on fast Fourier transform implementations of the first Rayleigh–Sommerfeld diffraction integral and assumes that the input field at the pupil plane has a larger extent than the field in the focal region. A sampling grid is defined over a finite region in the output plane and referred to as a tile. The input field is divided into multiple separate spatial regions of the size of the output tile. Finally, the input tiles are added coherently to form a summed tile, which is propagated to the output plane. Since only a single tile is propagated, there are significant reductions of computational load and memory requirements. This method is combined either with a subpixel sampling technique or with a chirp z-transform to realize smaller sampling intervals in the output plane than in the input plane. For a given example the resulting methods enable a speedup of approximately in comparison to the normal angular spectrum method, while the memory requirements are reduced by more than 99%.
© 2014 Optical Society of America
Full Article | PDF ArticleMore Like This
Matthias Hillenbrand, Damien P. Kelly, and Stefan Sinzinger
J. Opt. Soc. Am. A 31(8) 1832-1841 (2014)
Giang-Nam Nguyen, Kevin Heggarty, Philippe Gérard, Bruno Serio, and Patrick Meyrueis
Appl. Opt. 53(10) 2196-2205 (2014)
Nicole J. Moore, Miguel A. Alonso, and Colin J. R. Sheppard
J. Opt. Soc. Am. A 24(7) 2057-2064 (2007)