Abstract

A new beam-propagation method is presented whereby the exact scalar Helmholtz propagation operator is replaced by any one of a sequence of higher-order (n, n) Padé approximant operators. The resulting differential equation may then be discretized to obtain (in two dimensions) a matrix equation of bandwidth 2n + 1 that is solvable by using Standard implicit solution techniques. The final algorithm allows (for n = 2) accurate propagation at angles of greater than 55 deg from the propagation axis as well as propagation through materials with widely differing indices of refraction.

© 1992 Optical Society of America

Multistep method for wide-angle beam propagation

G. Ronald Hadley
Opt. Lett. 17(24) 1743-1745 (1992)

Wide-angle beam propagation method without using slowly varying envelope approximation

Khai Q. Le and Peter Bienstman
J. Opt. Soc. Am. B 26(2) 353-356 (2009)

Beam propagation method using a [(p–1)/p] Padé approximant of the propagator

Ya Yan Lu and Pui Lin Ho
Opt. Lett. 27(9) 683-685 (2002)

Numerical studies of split-operator finite-difference alternating-direction implicit propagation techniques based on generalized Padé approximants

David Yevick, Moses Glasner, and Björn Hermansson
Opt. Lett. 17(10) 725-727 (1992)

Application of modified Padé approximant operators to time-domain beam propagation methods

Khai Q. Le, Trevor Benson, and Peter Bienstman
J. Opt. Soc. Am. B 26(12) 2285-2289 (2009)