OSA's Digital Library

Optics Express

Optics Express

  • Editor: J. H. Eberly
  • Vol. 7, Iss. 8 — Oct. 9, 2000
  • pp: 260–272
« Show journal navigation

Efficient and accurate numerical analysis of multilayer planar optical waveguides in lossy anisotropic media

Chengkun Chen, Pierre Berini, Dazeng Feng, Stoyan Tanev, and Velko P. Tzolov  »View Author Affiliations


Optics Express, Vol. 7, Issue 8, pp. 260-272 (2000)
http://dx.doi.org/10.1364/OE.7.000260


View Full Text Article

Acrobat PDF (164 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

This paper discusses a numerical method for computing the electromagnetic modes supported by multilayer planar optical waveguides constructed from lossy or active media, having in general a diagonal permittivity tensor. The method solves the dispersion equations in the complex plane via the Cauchy integration method. It is applicable to lossless, lossy and active waveguides, and to AntiResonant Reflecting Optical Waveguides (ARROW’s). Analytical derivatives for the dispersion equations are derived and presented for what is believed to be the first time, and a new algorithm that significantly reduces the time required to compute the derivatives is given. This has a double impact: improved accuracy and reduced computation time compared to the standard approach. A different integration contour, which is suitable for leaky modes is also presented. Comparisons are made with results found in the literature; excellent agreement is noted for all comparisons made.

© Optical Society of America

1. Introduction

The notions of guided and leaky modes are fundamental concepts in optical waveguiding theory. A knowledge of mode propagation characteristics is essential to the design of numerous guided-wave optoelectronic devices, passive and active components such as semiconductor lasers, electro-absorption and electro-optic modulators, switches, photodetectors, filters and couplers, to name but a few. Numerical methods that can efficiently and accurately model planar optical waveguides are thus of obvious importance since they are used as a basic tool in the design process.

The Transfer Matrix Method (TMM) [1

1. J. Chilwell and I. Hodgkinson, “Thin-films field-transfer matrix theory of planar multilayer waveguides and reflection from prism-loaded waveguide,” J. Opt. Soc. Am. A 1, 742–753, (1984) [CrossRef]

], [2

2. L. M. Walpita, “Solutions for planar optical waveguide equations by selecting zero elements in a characteristic matrix,” J. Opt. Soc. Am. A 2, 595–602, (1985) [CrossRef]

], as one of the primary tools for multilayer planar optical waveguide analysis, can generate the dispersion equation of the TE and TM modes supported by such structures in a straightforward manner. The waveguides can consist of any combination of lossless, lossy (dielectric, semiconductor, metallic) and active (including uniaxially anisotropic quantum wells) layers. Solving the dispersion equation yields the propagation constant and the electromagnetic field distribution for a mode. For the modes of lossy or active waveguides, and for leaky modes, the mode propagation constants, which are the roots of the dispersion equation, are complex numbers.

The novel features of the TMM and CIM that we present in this paper include: the derivation and use of analytical derivatives of the dispersion equations in anisotropic media, an algorithm for the rapid calculation of the derivatives, and a different integration contour for locating leaky modes. Using analytical derivatives in conjunction with the algorithm suggested improves the accuracy of the method and greatly reduces the CPU time required to find modes. The new integration contour improves the numerical efficiency of the integration calculation. The method is applicable to lossless, lossy and active waveguides, and to AntiResonant Reflecting Optical Waveguides (ARROW’s). The method can find leaky and guided modes. Anisotropic dielectrics described by a diagonal permittivity tensor are handled by the formulation.

2. Transfer Matrix Method

The transfer matrix method provides a straightforward formulation of the multilayer planar optical waveguide problem. A multilayer nonmagnetic anisotropic slab waveguide structure (µ=µ0), is shown in Figure 1. The refractive index tensor, n¯¯i, of the i-th layer can be in general complex, i.e., ñjji =njji-jkjji where njji and kjji are the refractive index and extinction coefficient along the jj direction of the i-th layer, and i=1,…, r is the layer number between the substrate and cover.

Figure 1. Structure of the multilayer planar optical waveguide

The TMM is a well-known approach for generating the dispersion equations governing the modes supported by a multilayer planar optical waveguide [1

1. J. Chilwell and I. Hodgkinson, “Thin-films field-transfer matrix theory of planar multilayer waveguides and reflection from prism-loaded waveguide,” J. Opt. Soc. Am. A 1, 742–753, (1984) [CrossRef]

]. The technique is only summarized here and the equations needed to obtain the derivatives analytically are given.

2.1 Maxwell’s Equations and TE Field Solutions

Maxwell’s curl equations [9

9. Hermann A. Haus, “Waves and Fields in Optoelectronics,” (New Jersey, Prentice-Hall Inc., 1984). Ch. 11

] for source-free, time-harmonic fields in anisotropic media are:

×E=jωμ0H
(1a)
×H=jωε0ε¯¯rE
(1b)

where ε0 is the free space permittivity, ω is the angular frequency, and ε̿r is a tensor of relative permittivity having the form:

ε¯¯r=[ε˜r,xxε˜r,yyε˜r,zz]=[n˜xx2n˜yy2n˜zz2]
(2)

For a TE mode (Ex, Ez, Hy =0) propagating in the + direction in the i-th layer, (xixx i+1), the non-zero electric and magnetic field components are:

Ei=ŷEyi(x)exp(jωtjγ˜z)
(3a)
Hi=[x̂Hxi(x)+ẑHzi(x)]exp(jωtjγ˜z)
(3b)

where x̂, ŷ, ẑ are the unit vectors in the x, y, z directions, respectively, γ˜=K 0(β-) is the complex propagation constant with β and α the normalized phase and attenuation constants, respectively, and k0 =ω/c=2π/λ0 , c is the speed of light in free space and λ0 is the free space wavelength. From equations (1) and (3), we can derive:

d2Eyi(x)dx+κ˜i2Eyi(x)=0
(4)

and

Hzi(x)=jωμ0dEyi(x)dx
(5a)
Hxi(x)=γ˜ωμ0Eyi(x)
(5b)

Eyi(x)=cos[κ˜i(xxi)]Eyi(xi)+1κ˜isin[κ˜i(xxi)]dEyi(xi)dx
(6a)
dEyi(xi)dx=κ˜isin[κ˜i(xxi)]Eyi(xi)+cos[κ˜i(xxi)]dEyi(xi)dx
(6b)

where xi defines the boundary between the i-th and (i+1)-th layer. Equation (6) implies that any sign for the square root of i is acceptable.

2.2 Transfer Matrix and the Dispersion Equation for the TE Modes

Using Equation (6), the tangential electric field and its derivative at the bottom of the i-th layer, (x=xi ), can be expressed as a function of the field and its derivative within that layer:

(Eyi(xi)dEyi(xi)dx)=(cos[κ˜i(xxi)]1κ˜isin[κ˜i(xxi)]κ˜isin[κ˜i(xxi)]cos[κ˜i(xxi)])(Eyi(x)dEyi(x)dx)
(7)

(EysdEysdx)=i=1rMi(EycdEycdx)
=(m11m12m21m22)(EycdEycdx)
(8)

where

Mi=(cos(κ˜idi)1κ˜isin(κ˜idi)κ˜isin(κ˜idi)cos(κ˜idi))fori=1,2,,r
(9)

The Mi are the transfer matrices for all of the r layers of thickness di . For guided modes in an open structure, the tangential fields in the substrate and cover must be exponentially decaying:

Eys(x)=Asexp(γ˜sx)dEys(x)dx=γ˜sAsexp(γ˜sx)}forx0
(10)
Eyc(x)=Bcexp[γ˜c(xxr+1)]dEyc(x)dx=γ˜cBcexp[γ˜c(xxr+1)]}forxxr+1
(11)

where γ˜s=±γ˜2k02n˜yys2, γ˜c=±γ˜2k02n˜yyc2 and ñyys and ñyyc are the substrate and cover complex refractive indices along the y direction, respectively. Equations (10) and (11) in conjunction with (8) yield the dispersion equation for the TE modes:

F(γ˜)=γ˜sm11+γ˜cm22m21γ˜sγ˜cm12=0
(12)

The zeros of Equation (12) correspond to the complex propagation constants γ˜.

2.3 Transfer Matrix and the Dispersion Equation for the TM Modes

For the TM modes (Hx, Hz, Ey =0) the following equations hold:

d2Hyi(x)dx+κ˜i2Hyi(x)=0
(13)
Ezi(x)=jωε0n˜zzi2dHyi(x)dx
(14a)
Exi(x)=γ˜ωε0n˜xxi2Hyi(x)
(14b)

Following a derivation similar to the TE case the transfer matrices of each layer are determined:

Mi=(cos(κ˜idi)n˜zzi2κ˜isin(κ˜idi)κ˜in˜zzi2sin(κ˜idi)cos(κ˜idi))fori=1,2,,r
(15)

where κ˜i=±n˜zzin˜xxik02n˜xxi2γ˜2 and the dispersion equation is:

F(γ˜)=γ˜sn˜zzs2m11+γ˜cn˜zzc2m22m21γ˜sγ˜cn˜zzs2n˜zzc2m12=0
(16)

3. Cauchy Integration Method

3.1 Summary of the Method

The Cauchy integration method [5

5. L. M. Delves and J. N. Lyness, “A numerical method for locating the zeros of an analytic function,” Math. Comp. , 21, 543–560, (1967) [CrossRef]

] is based on the argument principle [10

10. J. W. Brown and R. V. Churchill, “Complex Variables and Applications,” (Sixth Edition, New York: McGraw-Hill, 1996)

] and the residue theorem of complex analysis. If a function f(z) is analytic and does not go to zero over a closed integral contour, then the argument principle is of the form:

S0=1j2πcf(z)f(z)dz=NzNp
(17)

where Nz is the number of zeros and Np is number of poles inside the region enclosed by the contour C. If there are no poles in the region enclosed by C, then S0 is the number of zeros, and from the residue theorem we have:

Sm=1j2πczmf(z)f(z)dz=i=1s0zimform=1,2,,S0
(18)

where zi, i=1, 2, …, S0 are the roots of f(z) inside C and Sm is the sum of zim with m=1, 2, …, S0 . Equation (18) leads to a system of equations that can be used to evaluate the coefficients of a polynomial p(z) of degree S 0, which has the same roots z 1, … zS0 as the function f(z) inside C. The approximation polynomial p(z) can be written as:

p(z)=i=1s0(zzi)=k=0s0Ckzk
(19)

with Cs0=1. The coefficients Ck are given via Newton’s recursive formula:

Ck=1(kS0)j=1s0kSjCk+jfork=S01,,0
(20)

The polynomial p(z) can be solved by standard techniques such as Laguerre’s method [11

11. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, “Numerical Recipes in C,” (Second Edition, Cambridge, 1994)

]. The problem of finding the zeros of an arbitrary function f(z) is thus transformed to the simpler problem of finding the zeros of the polynomial p(z), for which a variety of reliable and efficient numerical methods exist.

3.2 Numerical Implementation

The CIM can be used to solve the dispersion equation of the multilayer planar optical waveguide in a straightforward manner. The poles are first identified at β=ñc and β=ñs as shown in Figure 2, parts (a) and (b). The rectangular integral contour C1 is selected for the guided modes. The area between the minimum and the maximum real parts of all the refractive indices is enclosed by the contour C1 . The integral contour C2 is selected for leaky modes. These contours do not enclose any of the poles.

Figure 2. The integral contours in the complex plane. The ‘ 1 correspond to the poles, ñ max is the complex refractive index with the maximum real part and is enclosed by C1, for TE modes in anisotropic media: ñ max=ñ yymax, ñs =ñyys and ñc =ñyyc . For TM modes in anisotropic media: ñ max=ñ xxmax, ñs =ñxxs and ñc =ñxxc .

(a) ñc s , (b) ñc =ñs .

To solve Equations (17) and (18) the derivative of the dispersion equation is needed. If the derivative cannot be obtained analytically then it can be estimated numerically via Cauchy’s theorem [10

10. J. W. Brown and R. V. Churchill, “Complex Variables and Applications,” (Sixth Edition, New York: McGraw-Hill, 1996)

], which states that the first derivative of a function f(z) at z=z0 is given by:

f(z0)=1j2πDf(z)(zz0)2dz
(21)

where D is any closed path which encloses the point z0 and f(z) is analytic inside and on D.

Reference [7

7. E. Anemogiannis and E. N. Glytsis, “Multilayer waveguides: efficient numerical analysis of general structures”, J. Lightwave Tech. , 10, 1344–1351, (1992) [CrossRef]

] gives an algorithm for solving Equation (21) numerically. The algorithm is based on the parametric variable transformation z=z0 +Rej2πx where R is the radius of a circle centered at z0 defining the contour over which the function f(z) is evaluated. The m-point trapezoidal integration rule is then applied to Equation (21), with z transformed as described, yielding the following numerical approximation to the derivative:

Rf(z0)R[f(z0)]m=1ml=1mf(z0+Rej2πml)ej2πml
(22)

The computation of the numerical derivative according to Equation (22) is the most intensive computational task in the standard CIM [7

7. E. Anemogiannis and E. N. Glytsis, “Multilayer waveguides: efficient numerical analysis of general structures”, J. Lightwave Tech. , 10, 1344–1351, (1992) [CrossRef]

]. When z0 is very close to the zero of the function, the numerical derivative converges very slowly and inaccurately.

To solve Equations (17) and (18), integration is performed numerically by applying an adaptive integration method [12

12. J. R. Rice, “Numerical Methods, Software, and Analysis,” (IMSL Reference Edition. New York: McGraw-Hill, 1983)

] in a manner similar to reference [7

7. E. Anemogiannis and E. N. Glytsis, “Multilayer waveguides: efficient numerical analysis of general structures”, J. Lightwave Tech. , 10, 1344–1351, (1992) [CrossRef]

]. This approach is summarized here for completeness. The quadrature integration of a function F(z) is represented by the summation:

abF(z)=i=1KWiF(zi)
(23)

where Wi are the weights and azib are the integration nodes. A 7-point Gaussian rule G7f and 15-point Kronrod [13

13. A. S. Kronrod, “Nodes and Weights of Quadrature Formulas,” (NewYork: Consultants Bureau, 1965)

] rule K15f are used to estimate integration along the straight line connecting the end-points of the interval [a, b] in the complex plane. The local estimate is taken as K15f since the Kronrod rule is more accurate than the Gaussian rule. If the local error estimate |K15f-G7f|/|K15f| is less than a specified tolerance, the local estimate is accepted as the integration result for the specific interval. Otherwise the interval is bisected and the same procedure is applied to the two new subintervals. In the same integration routine the calculation of the products zimf,(zi )/f(zi ) for m=1, …, S 0 can be incorporated such that f(zi) and f’(zi ) are evaluated only once for the specific node zi for all the summations S 1, …, SS0.

The polynomial p(z) can be solved by Laguerre’s method. The roots of polynomial p(z) and f(z) do not coincide exactly due to the integration errors introduced in Equation (23). After the roots of the polynomial p(z) are obtained, a further refinement must be performed to find the roots of f(z) by applying Muller’s Method [11

11. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, “Numerical Recipes in C,” (Second Edition, Cambridge, 1994)

] with the initial guess being the roots of p(z).

As discussed in reference [7

7. E. Anemogiannis and E. N. Glytsis, “Multilayer waveguides: efficient numerical analysis of general structures”, J. Lightwave Tech. , 10, 1344–1351, (1992) [CrossRef]

], the region enclosed by the contour of interest C ought to be subdivided into smaller regions such that S0 ≤4. This ensures that the degree and the coefficients of the polynomial p(z) remain small, thus helping to reduce numerical errors when locating the roots. Considering all sub-regions enclosed by the original contour leads to all roots of interest.

4. Analytical Derivatives of the Transfer Matrices and Dispersion Equations

4.1. TE Modes

For the TE modes, let γ˜=k0ũ, κ˜i=±k0n˜yyi2u˜2, dκ˜idu˜=k02κ˜iu˜ and ddu˜(1κ˜i)=k02κ˜i3u˜. Then from Equation (9) we derive for each layer i:

dMidu˜=u˜(k02diκ˜isin(κ˜idi)k02diκ˜i2cos(κ˜idi)k0κ˜i3sin(κ˜idi)k02κ˜isin(κ˜idi)k02dicos(κ˜idi)k02diκ˜isin(κ˜idi))fori=1,,r
(24)

The derivative of the transfer matrix describing the multilayer structure for i=1, …, r is then:

dMdu˜=(dm11du˜dm12du˜dm21du˜dm22du˜)=j=1r(dMjdu˜·i=1ijrMi)
(25)

The derivatives of γ˜s and γ˜c with respect to ũ are: dγ˜sdu˜=k02u˜γ˜s and dγ˜cdu˜=k02u˜γ˜c, thus from Equation (12) we obtain the derivative of the dispersion equation as:

dF(u˜)du˜=k02u˜γ˜sm11+k02u˜γ˜cm22+γ˜sdm11du˜+γ˜cdm22du˜dm21du˜k02u˜γ˜cγ˜sm12k02u˜γ˜sγ˜cm12γ˜sγ˜cdm12du˜
(26)

4.2 TM Modes

Similarly, for the TM modes κ˜i=±n˜zzin˜xxik0n˜xxi2u˜2=±kn˜xxi2u˜2 where k=n˜zzin˜xxik0. Then dκ˜idu˜=k2κ˜iu˜ and ddu˜(1κ˜i)=k2κ˜i3u˜, and from Equation (15) we derive for each layer i:

dMidu˜=u˜(k2diκ˜isin(κ˜idi)k2n˜zzi2diκ˜i2cos(κ˜idi)k2n˜zzi2κ˜i3sin(κ˜idi)k2κ˜in˜zzi2sin(κ˜idi)k2din˜zzi2cos(κ˜idi)k2diκ˜isin(κ˜idi))fori=1,,r
(27)

From Equation (16), the derivative of the dispersion equation is:

dF(u˜)du˜=k02u˜γ˜sn˜xxs2m11+k02u˜γ˜cn˜xxc2m22+γ˜sn˜zzs2dm11du˜+γ˜cn˜zzc2dm22du˜dm21du˜γ˜sγ˜cn˜zzs2n˜zzc2dm12du˜
k02u˜γ˜cn˜xxs2n˜zzc2γ˜sm12k02u˜γ˜sn˜xxc2n˜zzs2γ˜cm12
(28)

4.3 Algorithm for Computing the Derivative of the Transfer Matrices

The transfer matrices, and their derivatives as expressed by Equation (25), can be computed in an efficient manner. The idea is to accumulate the derivative of the transfer matrix as the latter is being constructed. Consider accumulation over the first 3 layers of a waveguide structure. For the first layer (i=1), M 1 and dM1du˜ are computed. For the second layer (i=2), M 2 and dM2du˜ are computed, then dMdu˜=d(M1M2)du˜=M1dM2du˜+dM1du˜M2 and M=M 1 M 2 are accumulated. For the third layer (i=3), M 3 and dM3du˜ are computed, then dMdu˜=d(M1M2M3)du˜=d(M1M2)du˜M3+(M1M2)dM3du˜ and M=M 1 M 2 M 3 are accumulated. The accumulation process continues until all layers in the structure have been handled. Clearly, one matrix product is required per layer to accumulate the transfer matrix and two products are required to accumulate its derivative. This algorithm thus accumulates the transfer matrix for r layers in r matrix products and the derivative of the transfer matrix in 2r matrix products. The CIM based on our analytical derivatives and algorithm is thus an order r method, which is a substantial improvement over a similar approach of order r! discussed in [14].

5. Numerical Results and Discussion

5.1 Guided Modes in a Lossless Waveguide

An inhomogeneous waveguide having an index of refraction that follows an exponential distribution was analyzed in order to compare the computational efficiency of the CIM using numerical and analytical derivatives. The ADR method [14

14. E. Anemogiannis, E. N. Glytsis, and T. K. Gaylord, “Efficient solution of eigenvalue equations of optical waveguiding structures”, J. Lightwave Tech. , 12, 2080–2084, (1994) [CrossRef]

] was also implemented and its computation times are included for comparison. The parameters of the structure analyzed were taken from reference [14

14. E. Anemogiannis, E. N. Glytsis, and T. K. Gaylord, “Efficient solution of eigenvalue equations of optical waveguiding structures”, J. Lightwave Tech. , 12, 2080–2084, (1994) [CrossRef]

]:

n(x)=ns2+2Δnse(xα)forx<0
(29a)
n(x)=ncforx>0
(29b)

where ns =2.177, nc =1.0, Δ=0.043, α=0.931 µm and the thickness of the inhomogeneous layer is set to d=4 µm. The normalized phase constant (or effective index) of the two TE modes supported by this structure can be obtained analytically. They are β=2.19075 for the TE0 mode and β=2.17930 for the TE1 mode [14

14. E. Anemogiannis, E. N. Glytsis, and T. K. Gaylord, “Efficient solution of eigenvalue equations of optical waveguiding structures”, J. Lightwave Tech. , 12, 2080–2084, (1994) [CrossRef]

]. The effective indices of both modes, computed using our approach, are shown in Table I. The computed values are seen to converge to the analytical values as the number of layers used to approximate the continuous profile is increased. Agreement with the analytical values is excellent for a large number of layers.

TABLE I. Effective Index of Guided Modes Supported by an Inhomogeneous Waveguide.

table-icon
View This Table
| View All Tables

Figure 3 shows the execution time required to compute the effective index of the TE0 mode supported by the structure as a function of the number of layers used to approximate Equation (29a). The numerical derivative was computed to a relative accuracy of 10-3 and the tolerance in assessing the local error estimate when evaluating Equation (23) was set to 10-5 for all three methods. From Figure 3, it is apparent that the CIM based on the analytical derivative is order r in matrix operations, like the CIM based on the numerical derivative and the ADR method. Using an analytical derivative is more accurate and clearly much more computationally efficient than using a numerical one. The computational cost of the CIM with the analytical derivative is only slightly higher than the ADR method, which is less robust and more complex.

Figure 3. Comparison of the normalized execution times required for the ADR and the CIM based on a numerical derivative (ND) and an analytical derivative (AD).

5.2 Guided Modes in Quantum Well Active Waveguides

Quantum well (QW) layers used in active photonic devices are in fact anisotropic. Their optical response is different for polarizations parallel and perpendicular to the direction of growth. The anisotropy is uniaxial and appears even though quantum wells are made up from cubic semiconductor materials. Modeling uniaxial media is essential if accurate waveguide propagation constants are desired. We applied our method to a novel InP-based QW laser structure as defined in Table II in order to illustrate its ability to handle anisotropy and gain. The second layer is the active anisotropic quantum well layer. The results are shown in Table II and gain is observed for the first TE and TM modes.

TABLE II. Effective Index of Guided Modes Supported by a 6-layer Quantum Well Active Waveguide.

table-icon
View This Table
| View All Tables

5.3 Leaky Modes in Lossless Waveguides

The CIM can be used to find leaky modes. The leaky mode propagation constants satisfy the same dispersion equation as the guided modes, but the appropriate sign of γ˜s and γ˜c must be chosen. If we assume γ˜s=γsr +jγsi and γ˜c=γcr +jγci , for β<ns (β<nc ), we should choose γsr <0 and γsi >0 (γcr <0 and γci >0). The integral contour is selected as C2 , shown in Fig. 2(a) for ns > nc . The reason for selecting the top part of the integral contour along the β axis is that the imaginary part of the leaky mode propagation constants is negative. The integral contour above the β axis will lead to large integration errors and generate roots for the polynomial p(z) far away from the roots of the dispersion equation. For the TE modes supported by the waveguide structure described in Table III, our results are complete agreement with those reported in [1

1. J. Chilwell and I. Hodgkinson, “Thin-films field-transfer matrix theory of planar multilayer waveguides and reflection from prism-loaded waveguide,” J. Opt. Soc. Am. A 1, 742–753, (1984) [CrossRef]

], as shown.

TABLE III. Effective Index of Leaky Modes Supported by a 4-layer Lossless Waveguide.

table-icon
View This Table
| View All Tables

5.4 ARROW Waveguides

ARROW (AntiResonant Reflecting Optical Waveguide) waveguides [15

15. T. Baba and Y. Kokubun, “Dispersion and radiation loss characteristics of antiresonsnt reflecting optical waveguides-Numerical Results and Analytical Expressions”, IEEE J. Quantum Electron. , 28, 1689–1700, (1992) [CrossRef]

], [16

16. W. Huang, R. M. Shubiar, A. Nathan, and Y. L. Chow, “The modal characteristics of ARROW structures”, J. Lightwave Tech. , 10, 1015–1022. (1992) [CrossRef]

] are based on Fabry-Perot reflection instead of total internal reflection. The modes supported by such structures can be particularly difficult to find due to the proximity of the corresponding zeros in the complex plane. An example structure is described in Table IV [17

17. J. Deng and Y. Huang,“A novel hybrid coupler based on antiresonant reflecting optical waveguides,” J. Lightwave Tech. , 16, 1062–1069, (1998) [CrossRef]

], where the effective index of the first 6 TE and TM ARROW modes found using our CIM are also given.

TABLE IV. Effective Index of ARROW Modes Supported by a 9-layer ARROW Waveguide

table-icon
View This Table
| View All Tables

Figure 4. (a) Field distribution for the symmetric TM mode: TM1. (b) Field distribution for the first anti-symmetric TM mode: TM2.

5.5 Anisotropic ARROW Waveguides

Many anisotropic materials have desirable properties at optical wavelengths, including low losses, and large electro-optic or photo-elastic effects. In addition, anisotropy may be introduced during material growth or device processing. It is thus important for a mutilayer waveguide modeling tool to be able to handle effectively such materials.

As the last example, we have applied our method to obtain the propagation constant of modes supported by an anisotropic ARROW waveguide [18

18. B. Ray and G W. Hanson, “Some effects of anisotropy on planar antiresonant reflecting optical waveguides”, J. Lightwave Tech. , 14, 202–208, (1996) [CrossRef]

]. The structure of interest is defined in Table V and the first 4 ARROW modes for the isotropic (AR=1) and anisotropic cases (AR=nxxi/nzzi=1.03) are given. Our results for the TE modes are in complete agreement with those reported in [19

19. E. Anemogiannis, E. N. Glytsis, and T. K. Gaylord, “Determination of guided and leaky modes in lossless and lossy planar multiplayer optical waveguides: reflection pole method and wavevector density method”, J. Lightwave Tech. , 17, 929–941, (1999) [CrossRef]

] for the case AR=1.

TABLE V. Effective Index of ARROW Modes Supported by a 3-layer Anisotropic ARROW Waveguide

table-icon
View This Table
| View All Tables

6. Conclusion

A numerical method based on integration in the complex plane has been developed to characterize the modes supported by multilayer planar optical waveguides constructed from lossy or active anisotropic media. The analytical derivative of the transfer matrices and dispersion equations were obtained and presented for the first time. An algorithm for efficiently computing the derivative was given and a different integral contour to be used for finding leaky modes was presented. The method reported has several important advantages over other methods for finding modes of planar waveguides. An advantage is the ability to find the total number of modes in a region of interest in the complex plane. Another advantage is that the method can be used to characterize lossless, lossy, active and ARROW waveguides in anisotropic media. The method can generate guided modes as well as leaky modes. The combination of high accuracy, low computation time and broad range of applicability makes the method very attractive for integration into commercial computer-aided design and modeling tools for integrated optics.

7. Acknowledgements

The authors gratefully acknowledge Valery I. Tolstikhin, for providing the second example structure (an active anisotropic waveguide) investigated in this work.

References and links

1.

J. Chilwell and I. Hodgkinson, “Thin-films field-transfer matrix theory of planar multilayer waveguides and reflection from prism-loaded waveguide,” J. Opt. Soc. Am. A 1, 742–753, (1984) [CrossRef]

2.

L. M. Walpita, “Solutions for planar optical waveguide equations by selecting zero elements in a characteristic matrix,” J. Opt. Soc. Am. A 2, 595–602, (1985) [CrossRef]

3.

K. H. Schlereth and M. Tacke, “The complex propagation constant of multilayer waveguides: An algorithm for a personal computer,” IEEE J. Quantum Electron. , 26, 627–630, (1990) [CrossRef]

4.

L. Sun and E. Marhic, “Numerical study of attenuation in multilayer infrared waveguides by the circle-chain convergence method”, J. Opt. Soc. Am. B 8, 478–483, (1991) [CrossRef]

5.

L. M. Delves and J. N. Lyness, “A numerical method for locating the zeros of an analytic function,” Math. Comp. , 21, 543–560, (1967) [CrossRef]

6.

L. C. Botten and M. S. Craig, “Complex zeros of analytic functions”, Comput. Phys. Commun. , 29, 245–259, (1983) [CrossRef]

7.

E. Anemogiannis and E. N. Glytsis, “Multilayer waveguides: efficient numerical analysis of general structures”, J. Lightwave Tech. , 10, 1344–1351, (1992) [CrossRef]

8.

R. E. Smith, S. N. Houde-Walter, and G. W. Forbes, “Mode determination for planar waveguides using the four-sheeted dispersion relation,” IEEE J. Quantum Electron. , 28, 1520–1526, (1992) [CrossRef]

9.

Hermann A. Haus, “Waves and Fields in Optoelectronics,” (New Jersey, Prentice-Hall Inc., 1984). Ch. 11

10.

J. W. Brown and R. V. Churchill, “Complex Variables and Applications,” (Sixth Edition, New York: McGraw-Hill, 1996)

11.

W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, “Numerical Recipes in C,” (Second Edition, Cambridge, 1994)

12.

J. R. Rice, “Numerical Methods, Software, and Analysis,” (IMSL Reference Edition. New York: McGraw-Hill, 1983)

13.

A. S. Kronrod, “Nodes and Weights of Quadrature Formulas,” (NewYork: Consultants Bureau, 1965)

14.

E. Anemogiannis, E. N. Glytsis, and T. K. Gaylord, “Efficient solution of eigenvalue equations of optical waveguiding structures”, J. Lightwave Tech. , 12, 2080–2084, (1994) [CrossRef]

15.

T. Baba and Y. Kokubun, “Dispersion and radiation loss characteristics of antiresonsnt reflecting optical waveguides-Numerical Results and Analytical Expressions”, IEEE J. Quantum Electron. , 28, 1689–1700, (1992) [CrossRef]

16.

W. Huang, R. M. Shubiar, A. Nathan, and Y. L. Chow, “The modal characteristics of ARROW structures”, J. Lightwave Tech. , 10, 1015–1022. (1992) [CrossRef]

17.

J. Deng and Y. Huang,“A novel hybrid coupler based on antiresonant reflecting optical waveguides,” J. Lightwave Tech. , 16, 1062–1069, (1998) [CrossRef]

18.

B. Ray and G W. Hanson, “Some effects of anisotropy on planar antiresonant reflecting optical waveguides”, J. Lightwave Tech. , 14, 202–208, (1996) [CrossRef]

19.

E. Anemogiannis, E. N. Glytsis, and T. K. Gaylord, “Determination of guided and leaky modes in lossless and lossy planar multiplayer optical waveguides: reflection pole method and wavevector density method”, J. Lightwave Tech. , 17, 929–941, (1999) [CrossRef]

OCIS Codes
(130.2790) Integrated optics : Guided waves
(230.7390) Optical devices : Waveguides, planar

ToC Category:
Research Papers

History
Original Manuscript: September 5, 2000
Published: October 9, 2000

Citation
Chengkun Chen, Pierre Berini, Dazeng Feng, Stoyan Tanev, and Velko Tzolov, "Efficient and accurate numerical analysis of multilayer planar optical waveguides in lossy anisotropic media," Opt. Express 7, 260-272 (2000)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-7-8-260


Sort:  Journal  |  Reset  

References

  1. J. Chilwell and I. Hodgkinson, " Thin-films field-transfer matrix theory of planar multilayer waveguides and reflection from prism-loaded waveguide," J. Opt. Soc. Am. A 1, 742-753, (1984) [CrossRef]
  2. L. M. Walpita, " Solutions for planar optical waveguide equations by selecting zero elements in a characteristic matrix," J. Opt. Soc. Am. A 2, 595-602, (1985) [CrossRef]
  3. K. H. Schlereth and M. Tacke, "The complex propagation constant of multilayer waveguides: An algorithm for a personal computer," IEEE J. Quantum Electron., 26, 627-630, (1990) [CrossRef]
  4. L. Sun and E. Marhic, "Numerical study of attenuation in multilayer infrared waveguides by the circle-chain convergence method," J. Opt. Soc. Am. B 8, 478-483, (1991) [CrossRef]
  5. L. M. Delves and J. N. Lyness, "A numerical method for locating the zeros of an analytic function," Math. Comp., 21, 543-560, (1967) [CrossRef]
  6. L. C. Botten and M. S. Craig, "Complex zeros of analytic functions," Comput. Phys. Commun. 29, 245- 259, (1983) [CrossRef]
  7. E. Anemogiannis, and E. N. Glytsis, "Multilayer waveguides: efficient numerical analysis of general structures," J. Lightwave Tech. 10, 1344-1351, (1992) [CrossRef]
  8. R. E. Smith, S. N. Houde-Walter, and G. W. Forbes, " Mode determination for planar waveguides using the four-sheeted dispersion relation," IEEE J. Quantum Electron. 28, 1520-1526, (1992) [CrossRef]
  9. Hermann A. Haus, Waves and Fields in Optoelectronics, (New Jersey, Prentice-Hall Inc., 1984). Ch. 11
  10. J. W. Brown and R. V. Churchill, Complex Variables and Applications, (Sixth Edition, New York: McGraw-Hill, 1996)
  11. W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in C, (Second Edition, Cambridge, 1994).
  12. J. R. Rice, Numerical Methods, Software, and Analysis, (IMSL Reference Edition. New York: McGraw-Hill, 1983)
  13. A. S. Kronrod, Nodes and Weights of Quadrature Formulas, (NewYork: Consultants Bureau, 1965)
  14. E. Anemogiannis, E. N. Glytsis and T. K. Gaylord, "Efficient solution of eigenvalue equations of optical waveguiding structures," J. Lightwave Tech. 12, 2080-2084, (1994) [CrossRef]
  15. T. Baba and Y. Kokubun, "Dispersion and radiation loss characteristics of antiresonsnt reflecting optical waveguides-Numerical Results and Analytical Expressions," IEEE J. Quantum Electron., 28, 1689-1700, (1992) [CrossRef]
  16. W. Huang, R. M. Shubiar, A. Nathan and Y. L. Chow, "The modal characteristics of ARROW structures," J. Lightwave Tech. 10, 1015-1022. (1992) [CrossRef]
  17. J. Deng and Y. Huang,"A novel hybrid coupler based on antiresonant reflecting optical waveguides," J. Lightwave Tech. 16, 1062-1069, (1998) [CrossRef]
  18. B. Ray and G W. Hanson, "Some effects of anisotropy on planar antiresonant reflecting optical waveguides," J. Lightwave Tech. 14, 202-208, (1996) [CrossRef]
  19. E. Anemogiannis, E. N. Glytsis and T. K. Gaylord, "Determination of guided and leaky modes in lossless and lossy planar multiplayer optical waveguides: reflection pole method and wavevector density method," J. Lightwave Tech. 17, 929-941, (1999) [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited