A novel Fourier based 3D full-vectorial beam propagation method

In this paper a Fourier based full-vectorial beam propagation method (BPM) with perfectly matched layer (PML) absorbing boundary conditions is presented. The method starts with the slowly varying envelope coupled wave equations, including coordinate-stretching type PMLs, which are discretised using the Galerkin method with Fourier basis functions in the transversal direction. The proposed method is then assessed in a wide variety of situations, including 2D propagation of TE and TM waves and 3D full-vectorial situations. The obtained results show that good accuracy can be obtained even for complex strongly-guiding waveguides with a small number of basis functions.