Computation of wave propagation in integrated optical devices using z-transient variational principles

As an alternative to the classical beam propagation method (BPM), a variational method is presented to solve the TE and TM Helmholtz equations in the paraxial approximation for the propagation of polarized beams through optical waveguides. Using the method of local potentials, the paraxial wave equations are first converted into equivalent z-transient variational principles. These functionals are minimized using a combination of the Rayleigh-Ritz finite-element procedure and a Crank-Nicholson-like finite-difference scheme. Solutions for anisotropic materials are obtained by applying standard Galerkin finite-element and finite-difference methods to a variational formulation derived from the coupled TE/TM paraxial Helmholtz equations.>