Efficient interface conditions for the semi-vectorial finite-difference beam propagation method

Efficient interface conditions (EICs) are derived for the propagation equation using the slowly varying envelope approximation for the dominant electric field component. At the interface between two different media, the two lateral second derivatives in the discretized propagation equation are adapted such that the discretized modal field equation is correct up to second order in the lateral grid spacing. Since the error term is then of the order of the lateral grid spacing, our EICs are first-order EICs. These interface conditions are compared with well-known zero-order EICs derived by Stern and Kim and Ramaswamy. It is shown that the first-order EICs yield faster convergence to the exact effective index value as the lateral grid spacing is decreased than do the zero-order EICs. It turns out that our EICs are very much like those derived by Vassallo. Using essentially the same method, he derived EICs of second and first order for the field component respectively parallel and perpendicular, to the interface. Hence the accuracy of his EICs is one order higher for the field component parallel to the interface, although it introduces an extra asymmetry in the propagation matrix.

[1]  Piero Ravetto,et al.  Finite difference methods for the analysis of integrated optical waveguides , 1990 .

[2]  A. R. Mitchell,et al.  The Finite Difference Method in Partial Differential Equations , 1980 .

[3]  Ramu V. Ramaswamy,et al.  Modeling of graded-index channel waveguides using nonuniform finite difference method , 1989 .

[4]  R. Ridder,et al.  Modal fields calculation using the finite difference beam propagation method , 1994 .

[5]  J. C. Chen,et al.  Computation of higher-order waveguide modes by imaginary-distance beam propagation method , 1994 .

[6]  K. L. Andrew,et al.  Polarization, penetration, and exchange effects in the hydrogenlike nf and ng terms of cesium , 1981 .

[7]  A. Ralston.  A first course in numerical analysis , 1965 .

[8]  Lars Thylén,et al.  The beam propagation method: an analysis of its applicability , 1983 .

[9]  G. R. Hadley,et al.  Transparent boundary condition for beam propagation. , 1991, Optics letters.

[10]  B. Hermansson,et al.  Efficient beam propagation techniques , 1990 .

[11]  Gijsbertus J.M. Krijnen,et al.  Numerical solution method of nonlinear guided modes with a finite difference complex axis beam propagation method , 1995 .

[12]  Gijsbertus J.M. Krijnen,et al.  Efficient interface conditions for the finite difference beam propagation method , 1992 .

[13]  S. Jungling,et al.  A study and optimization of eigenmode calculations using the imaginary-distance beam-propagation method , 1994 .

[14]  W. Burns,et al.  Mode dispersion in diffused channel waveguides by the effective index method. , 1977, Applied optics.

[15]  Pao-Lo Liu,et al.  The semivectorial beam propagation method , 1993 .

[16]  Paul Lagasse,et al.  Beam-propagation method: analysis and assessment , 1981 .

[17]  William H. Press,et al.  Numerical recipes , 1990 .

[18]  Paul Lagasse,et al.  Application of propagating beam methods to electromagnetic and acoustic wave propagation problems: A review , 1987 .

[19]  Youngchul Chung,et al.  An assessment of finite difference beam propagation method , 1990 .

[20]  David Yevick,et al.  New formulations of the matrix beam propagation method: application to rib waveguides , 1989 .

[21]  A. Bjarklev,et al.  Accurate finite difference beam propagation method for complex integrated optical structures , 1993, IEEE Photonics Technology Letters.