Fast method for physical optics propagation of high-numerical-aperture beams.

A method is presented that expands the scheme of physical optics propagation beyond the Fresnel approximation to include beams that are nonparaxial. The formalism retains most of the calculation advantages of the Fresnel approach; i.e., it is based on a single Fourier transform step. The kernel of the new transformation is no longer separable in Cartesian coordinates; thus the formalism can account for astigmatic coupling effects originating purely from diffraction. The validity limits of the proposed algorithm are explored. Analytical expressions, numerical simulation results, and experimental data are compared.