Eigenfunction expansion of the electric fields in the focal region of a high numerical aperture focusing system

The Debye-Wolf electromagnetic diffraction integral is now routinely used to describe focusing by high numerical (NA) lenses. We obtain an eigenfunction expansion of the electric vector field in the focal region in terms of Bessel and generalized prolate spheroidal functions. Our representation has many optimal and desirable properties which offer considerable simplification to the evaluation and analysis of the DebyeWolf integral. It is potentially also useful in implementing two-dimensional apodization techniques to synthesize electromagnetic field distributions in the focal region of a high NA lenses. Our work is applicable to many areas, such as optical microscopy, optical data storage and lithography.