Modeling second-order nonlinear effects in optical waveguides using a parallel-processing beam propagation method

In this work, we present a simple efficient numerical solution for the three-dimensional coupled wave equations containing a second-order nonlinearity, using an explicit finite difference beam propagation method (EFD BPM). The linear EFD-BPM is known to be very efficient and to gain large speed up when implemented on parallel computers. The new nonlinear version of the EFD-BPM has the same features of the linear counterpart in using two separate computational windows, one for the fundamental field and the other for the second-harmonic field. We demonstrate the implementation and discuss the application of this method to a nonlinear rib waveguide using the quasi-phase-matching technique