The fast and accurate propagation of general optical fields in free space is still a challenging task. Most of the standard algorithms are either fast or accurate. In the paper we introduce without further physical approximations three new algorithms for the fast propagation of non-paraxial vectorial optical fields containing smooth but strong phase terms. Dependent on the shape of the smooth phase term different propagation operators are applied. The first method for the efficient propagation of fields, which are containing smooth spherical phase terms, is based on Mansuripur's extended Fresnel diffraction integral1 using fast Fourier Transformations. This concept is improved by Avoort's parabolic fitting technique2 and the parameter space, for which the extended Fresnel operator is numerically feasible, is discussed in detail. Furthermore we introduce the inversion of the extended Fresnel operator for the fast propagation of non-paraxial fields into the focal region. Secondly we discuss a new semi-analytical spectrum of plane waves (SPW) operator for the quick propagation of fields with smooth linear phase terms. The method is based on the analytical handling of the linear phase term and the lateral offset, which reduces the required computational window sizes in the target plane. Finally we generalize the semi-analytical SPW operator concept to universal shapes of smooth phases by de- composing non-paraxial fields into subfields with smooth linear phase terms. In the target plane, all propagated subfields are added coherently where the analytical known smooth linear phase terms are handled numerical efficient by a new inverse parabasal decomposition technique (PDT). Numerical results are presented for examples, demonstrating the efficiency and the accuracy of the three new propagation methods. All simulations were done with the optics software VirtualLabTM.3
[1]
Frank Wyrowski,et al.
Introduction to field tracing
,
2011
.
[2]
Frank Wyrowski,et al.
Parabasal field decomposition and its application to non-paraxial propagation.
,
2012,
Optics express.
[3]
E. Brigham,et al.
The fast Fourier transform and its applications
,
1988
.
[4]
P. Roman.
A Scalar Representation of Electromagnetic Fields: III
,
1959
.
[5]
M..
Certain computational aspects of vector diffraction problems
,
2002
.
[6]
J. Tukey,et al.
An algorithm for the machine calculation of complex Fourier series
,
1965
.
[7]
L. Mandel,et al.
Optical Coherence and Quantum Optics
,
1995
.
[8]
J. J. M. Braat,et al.
Aberration retrieval from the intensity point-spread function in the focal region using the extended Nijboer–Zernike approach
,
2005
.
[9]
K. Matsushima.
Shifted angular spectrum method for off-axis numerical propagation.
,
2010,
Optics express.
[10]
Shlomo Ruschin,et al.
Fast method for physical optics propagation of high-numerical-aperture beams.
,
2004,
Journal of the Optical Society of America. A, Optics, image science, and vision.
[11]
Colin J. R. Sheppard,et al.
Diffraction by a circular aperture: a generalization of Fresnel diffraction theory
,
1992
.
[12]
Jari Turunen,et al.
Elementary-field representations in partially coherent optics
,
2011
.
[13]
Charles L. Lawson,et al.
Solving least squares problems
,
1976,
Classics in applied mathematics.