Norm-conserving finite-difference beam-propagation method for TM wave analysis in step-index optical waveguides

Nonconservation of power is a perplexing problem in the propagating beam analysis of transverse magnetic (TM) waves in a z-variant step-index optical waveguide. To conserve the power in terms of a squared norm, a modified finite-difference (FD) formula is introduced that allows a general position of a core-cladding interface. The use of the modified formula contributes to a reduction in a field profile error caused by a staircase approximation with subsequent conservation of power, particularly for a symmetrical waveguide. To obtain the power conservation even in the analysis of an asymmetrical waveguide, a z-derivative of the refractive index is taken into account. An asymmetrical taper and tilted waveguides placed in parallel are investigated to validate the present technique.

[1]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[2]  Giuseppe Pelosi,et al.  Finite Element Software for Microwave Engineering , 1996 .

[3]  J. Yamauchi,et al.  Modified finite-difference formula for semivectorial H-field solutions of optical waveguides , 1998, IEEE Photonics Technology Letters.

[4]  Gijsbertus J.M. Krijnen,et al.  A comparison between different propagative schemes for the simulation of tapered step index slab waveguides , 1996 .

[5]  Thomas Weiland,et al.  Using the Maxwell grid equations to solve large problems , 1993 .

[6]  Hugo Hoekstra,et al.  On beam propagation methods for modelling in integrated optics , 1997 .

[7]  L. Poladian,et al.  Unification of TE and TM beam propagation algorithms , 1998, IEEE Photonics Technology Letters.

[8]  J. Yamauchi,et al.  Modified finite-difference formula for the analysis of semivectorial modes in step-index optical waveguides , 1997, IEEE Photonics Technology Letters.

[9]  T. Weiland,et al.  Maxwell's Grid Equations , 1990 .

[10]  C. Vassallo,et al.  Difficulty with vectorial BPM , 1997 .

[11]  Wei-Ping Huang,et al.  Efficient and accurate vector mode calculations by beam propagation method , 1993 .

[12]  Frank Schmidt,et al.  An adaptive approach to the numerical solution of Fresnel's wave equation , 1993 .

[13]  Charles Vassallo.  Interest of improved three-point formulas for finite-difference modeling of optical devices , 1997 .

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

[15]  M. S. Stern.  Semivectorial polarised H˜ field solutions for dielectric waveguides with arbitrary index profiles , 1988 .

[16]  Junji Yamauchi,et al.  Full-vectorial beam-propagation method based on the McKee-Mitchell scheme with improved finite-difference formulas , 1998 .

[17]  D Yevick,et al.  Correspondence of variational finite-difference (relaxation) and imaginary-distance propagation methods for modal analysis. , 1992, Optics letters.

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