An efficient method for computing optical waveguides with discontinuous refractive index profiles using spectral collocation method with domain decomposition

An accurate and efficient solution method using spectral collocation method with domain decomposition is proposed for computing optical waveguides with discontinuous refractive index profiles. The use of domain decomposition divides the usual single domain into a few subdomains at the interfaces of discontinuous refractive index profiles. Each subdomain can be expanded by a suitable set of orthogonal basis functions and patched at these interfaces by matching the physical boundary conditions. In addition, a new technique incorporating the effective index method and the Wentzel-Kramers-Brillouin method for the a-priori determination of the scaling factor in Hermite-Gauss or Laguerre-Gauss basis functions is introduced to considerably save computational time without experimenting with it. This method shares the same desirable property of the spectral collocation method of providing a fast and accurate solution but avoids the oscillatory solutions and improves the poor convergence problem of the simple spectral collocation method with single domain where regions of discontinuous refractive index profiles exist. Applications to several two- and three-dimensional waveguide structures having exact or accurate approximate solutions are given to test the accuracy and efficiency; all the results are found to be in excellent agreement.