OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 16, Iss. 21 — Oct. 13, 2008
  • pp: 17021–17030
« Show journal navigation

The complex Jacobi iterative method for three-dimensional wide-angle beam propagation

Khai Q. Le, R. Godoy-Rubio, Peter Bienstman, and G. Ronald Hadley  »View Author Affiliations


Optics Express, Vol. 16, Issue 21, pp. 17021-17030 (2008)
http://dx.doi.org/10.1364/OE.16.017021


View Full Text Article

Acrobat PDF (139 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

A new complex Jacobi iterative technique adapted for the solution of three-dimensional (3D) wide-angle (WA) beam propagation is presented. The beam propagation equation for analysis of optical propagation in waveguide structures is based on a novel modified Padé(1,1) approximant operator, which gives evanescent waves the desired damping. The resulting approach allows more accurate approximations to the true Helmholtz equation than the standard Padé approximant operators. Furthermore, a performance comparison of the traditional direct matrix inversion and this new iterative technique for WA-beam propagation method is reported. It is shown that complex Jacobi iteration is faster and better-suited for large problems or structures than direct matrix inversion.

© 2008 Optical Society of America

1. Introduction

The beam propagation method (BPM) has become one of the most widely used techniques for the study of optical waveguide devices [1

1. R. Scarmozzino, A. Gopinath, R. Pregla, and S. Helfert, “Numerical techniques for modeling guided-wave photonic devices,” IEEE J. Sel. Top. Quantum Electron. 6, 150–161 (2000). [CrossRef]

,2

2. C. Ma and E. V. Keuren, “A simple three dimensional wide-angle beam propagation method,” Opt. Express 14, 4668–4674 (2006). [CrossRef] [PubMed]

]. For the paraxial approximation or Fresnel equation used in BPM, this is usually accomplished by split-step methods, including the alternating-direction implicit (ADI) methods that are both fast and easy to implement [3

3. W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vettering, Numerical recipes: The art of scientific computing (Cambridge University Press, New York, 1986).

]. However, the split-step schemes are only first-order-accurate in the step size, and a direct matrix inversion is required for acceptable accuracy if WA propagation is needed. Besides the limitation to paraxial beams, these methods also restrict the simulations to a low refractive index contrast ratio between the core and cladding of the waveguide [4

4. M. D. Feit and J. A. Fleck Jr., “Analysis of rib waveguides and couplers by the propagation method,” J. Opt. Soc. Am. A 7, 73–79 (1990). [CrossRef]

]. Efforts have been taken to relax the limitations for WA simulations. Different treatments of WA-BPM based on the slowly varying envelope approximation have been developed, including the rational approximants of the square root operator [5

5. F. A. Milinazzo, C. A. Zala, and G. H. Brooke, “Rational square root approximations for parabolic equation algorithms,” J. Acoust. Soc. Am. 101, 760–766 (1997). [CrossRef]

], the one-way propagator [6

6. Y. Y. Lu and P. L. Ho, “Beam propagation method using a [(p-1)/p] Padé approximant of the propagator,” Opt. Lett. 27, 683–685 (2002). [CrossRef]

], the exponential of the square root operator [7

7. T. Anada, T. Hokazono, T. Hiraoka, J. P. Hsu, T. M. Benson, and P. Sewell, “Very-wide-angle beam propagation methods for integrated optical circuits,” IEICE Trans. Electron. E82-C, 1154–1158 (1999).

], and the Padé approximant operator [8

8. G. R. Hadley, “Wide-angle beam propagation using Padé approximant operators,” Opt. Lett. 17, 1426–1428 (1992). [CrossRef] [PubMed]

], for rectangular coordinates as well as an oblique coordinate system [9

9. S. Sujecki, “Wide-angle, finite-difference beam propagation in oblique coordinate system,” J. Opt. Soc. Am. 25, 138–145 (2007). [CrossRef]

]. The Padé-approximant-based WA-BPM is one of the most commonly used techniques for modeling optical waveguide structures. However, the method was originally limited to 2D structures due to the lack of efficient solvers. A 3D WA-BPM based on Hoekstra’s scheme was recently developed using the efficient Thomas algorithm and the splitting of the 3D Fresnel wave equation into three 2D wave equations [2

2. C. Ma and E. V. Keuren, “A simple three dimensional wide-angle beam propagation method,” Opt. Express 14, 4668–4674 (2006). [CrossRef] [PubMed]

]. However, this may cause splitting errors. Recently, C. Ma et al. [10

10. C. Ma and E. V. Keuren, “A three-dimensional wide-angle beam propagation method for optical waveguide structures,” Opt. Express 15, 402–407 (2007). [CrossRef] [PubMed]

] presented a new 3D WA-BPM also based on Hoekstra’s scheme that does not require the splitting of the Fresnel wave equation or the use of the ADI method. By using a technique for shifting the simulation window to reduce the dimension of the numerical equation and a threshold technique to further ensure its convergence, this approach shows accuracy and effectiveness. However, the resultant propagation scheme can be very slow if either the problem size is large or the structure or the boundary conditions are changing as the propagation proceeds, requiring frequent reinversions of the propagation matrix. Thus, it is imperative to find more efficient solution methods for 3D WA-BPM.

Recently, the complex Jacobi iterative method, a new iterative technique for solution of the indefinite Helmholtz equation, was introduced [11

11. G. R. Hadley, “A complex Jacobi iterative method for the indefinite Helmholtz equation,” J. Comp. Phys. 203, 358–370 (2005). [CrossRef]

]. It is based on a complex generalization of the point relaxation technique proposed by Jacobi in 1845, and has been shown to converge at a rate dependent only upon the grid size and effective absorption coefficient. For beam propagation of wave profiles with a 2D cross section, the wide angle beam propagation equation can be recast in terms of a Helmholtz equation with source term, and the effective absorption coefficient appearing in this equation is high, leading to rapid convergence of the method. In this work we apply the complex Jacobi iterative (CJI) method for beam propagation and show it to be highly efficient for the solution of large problems, compared with existing techniques.

The beam propagation equation for analysis of optical propagation in waveguide structures is based on a modified Padé(1,1) approximant operator we propose here, which gives evanescent waves the desired damping. Our modified Padé approximant propagation operators allow more accurate approximation to the true Helmholtz equation and faster convergence when the CJI method is used for the modified Padé approximant-based WA-BPM than those of the standard Padé approximant approach [8

8. G. R. Hadley, “Wide-angle beam propagation using Padé approximant operators,” Opt. Lett. 17, 1426–1428 (1992). [CrossRef] [PubMed]

]. Furthermore, since the utility of the CJI technique depends mostly upon its execution speed in comparison with the direct matrix inversion (DMI) method, we also present several speed comparisons. Numerical implementations are carried out for both 2D and 3D optical waveguide structures. The rest of this paper is structured as follows: In Section 2, we will describe the novel modified Padé approximant operators. Next, WA-BPM formulations using the CJI and DMI methods are presented in Section 3. In Section 4, the investigation of convergence rate of CJI is discussed. In Section 5, the benchmark results performed by WA-BPM using CJI and DMI are presented followed by conclusions in Section 6.

2. Modified Padé approximant operators

The scalar Helmholtz equation obtained by using the slowly varying envelope approximation is given by [8

8. G. R. Hadley, “Wide-angle beam propagation using Padé approximant operators,” Opt. Lett. 17, 1426–1428 (1992). [CrossRef] [PubMed]

],

Hzi2k2Hz2=iP2kH,
(1)

where P=2+k02(n2nref2)=2x2+2y2+k02(n2nref2) with k=k 0 nref, n the refractive index refractive index profile, nref the reference refractive index, k 0 the vacuum wavevector.

We may formally rewrite Eq. (1) in the form

Hz=iP2k1i2kzH.
(2)

Equation (2) suggests the recurrence relation

zn+1=iP2k1i2kzn.
(3)

Hadley [8

8. G. R. Hadley, “Wide-angle beam propagation using Padé approximant operators,” Opt. Lett. 17, 1426–1428 (1992). [CrossRef] [PubMed]

] proposed the rational approximation of WA beam propagation using Padé approximant operators with initial value of z0=0.Forz2, this gives us the well-known Padé (1,1) approximant-based WA beam propagation formula as follows:

HziP2k1+P4k2H.
(4)

If Eq. (4) is compared with a formal solution of Eq. (1) written in the well-known form.

Hz=i(P+k2k)H=ik(1+X1)H,
(5)

where X=Pk2, we obtain the approximation formula

1+X1P2k21+P4k2=X21+X4.
(6)

Since the operator X has a real spectrum, it is useful to consider the approximation of 1+X1 by the Padé approximant propagation operator. Figure 1 shows the absolute value of 1+X1 and its first-order Padé approximant called Hadley (1,1) as a function of X.

However, as the denominator of Hadley (1,1) gradually approaches zero, its absolute value approaches ∞ as can clearly be seen in Fig. 2. Physically this means that the standard Padé approximant incorrectly propagates the evanescent modes.

To circumvent this problem we propose a modified Padé approximant operator. First of all, following [12

12. Y. Y. Lu, “A complex coefficient rational approximation on of 1+x,” Appl. Numer. Math. 27, 141–156 (1998). [CrossRef]

], by multiplying both sides of Eq. (2) with ik, we obtain

ikz=Pk22ikz.
(7)

We may rewrite Eq. (7) as follows

f(X)=X2+f(X),
(8)

where f(X)=ikz.

Eq. (8) suggests the recurrence relation

fn+1(X)=X2+fn(X)forn=0,1,2,...
(9)

Y. Y. Lu [12

12. Y. Y. Lu, “A complex coefficient rational approximation on of 1+x,” Appl. Numer. Math. 27, 141–156 (1998). [CrossRef]

] has proved that Eq. (9) can provide a good approximation to 1+X1 with the initial value of

f0(X)=iβwhereβ>0.
(10)

Subsequently, we use this fact to go back to the original recurrence relation (3) and construct modified Padé approximant operators by using a different initial value of z0=kβ. For z2 the first-order modified Padé(1,1) approximant operator is given as follows:

HziP2k1+P4k2(1+iβ2)H,
(11)

The absolute value of the modified Padé(1,1) approximant of 1+X1 is also depicted in Fig. 1. It is seen that our modified Padé approximant operator (with β=2) allows more accurate approximations to the true Helmholtz equation than the standard Padé approximant operator. Furthermore, the standard rational Padé approximant incorrectly propagates the evanescent modes as their denominator gradually approaches zero while the modified Padé approximant gives the waves propagating in the evanescent region the desired damping as clearly seen in Fig. 2.

Fig. 1. The absolute values of (1+X)1/2-1 (solid line), its first-order standard Padé approximant (X/2)/(1+X/4) (solid line with circles) and modified Padé approximant (X/2)/(1+X/{4(1+ibeta/2)}) (dotted line).
Fig. 2. The absolute value of (1+X)1/2-1 (solid line), the first-order standard (X/2)/(1+X/4) (solid line with circles) and modified (X/2)/(1+X/{4(1+ibeta/2)}) (dotted line) Padé approximant of (1+X)1/2-1 with respect to X.

3. WA-beam propagation formulation

3.1 Basic equation

By using the modified Padé(1,1) approximant, the 3D semivectorial WA beam propagation equation can be written as follows [13

13. G. R. Hadley, “Low-truncation-error finite difference equations for photonics simulation I: beam propagation,” J. Lightwave Technol. 16, 134–141 (1998). [CrossRef]

]:

(1+ξP)ϕn+1=(1+ξ*P)ϕ,
(12)

where ξ=14k2(1+iβ/2)iΔz4k, ξ*the complex conjugate of ξ and Δz the propagation step.

3.2. WA-BPM using CJI

By dividing both sides of Eq. (12) by ξ, it may be written as an inhomogeneous Helmholtz equation

(2+k02(n2nref2)+1ξ)ϕn+1=(ξ*ξP+1ξ)ϕn,
(13)

or

(2+k02(n2nref2)+4k02nref21+(1+k2βΔz/2)β/2+ikΔz)(1+k2βΔz/2)2+k2Δz2)ϕn+1=sourceterm.
(14)

Thus the beam propagation can be cast as a 2D Helmholtz equation with source term in an effective medium with loss of 4k02nref2kΔz(1+k02nref2βΔz/2)2+k2Δz2. This loss is high for a typical choice of k Δz. This is a condition that favors rapid convergence for the CJI method.

3.3 WA-BPM using DMI

By discretizing Eq. (13), we find that

An+1ϕi+1,jn+1+Bn+1ϕi1,jn+1+Cn+1ϕi,jn+1+Dn+1ϕi,j1n+1+En+1ϕi,j+1n+1
=Anϕi+1,jn+Bnϕi1,jn+Cnϕi,jn+Dnϕi,j1n+Enϕi,j+1n.
(15)

where

An+1=Bn+1=1Δx2,Dn+1=En+1=1Δy2,
Cn+1=k02(n2nref2)+1ξ2(1Δx2+1Δy2),
An=Bn=ξ*ξ1Δx2,Dn=En=ξ*ξ1Δy2,
Cn=ξ*ξ(k02(n2nref2)+1ξ*2(1Δx2+1Δy2).
(16)

Actually, the coefficients in Eq. (16) are the same as the result when applying Hoekstra’s scheme to the 3D wave equation [10

10. C. Ma and E. V. Keuren, “A three-dimensional wide-angle beam propagation method for optical waveguide structures,” Opt. Express 15, 402–407 (2007). [CrossRef] [PubMed]

]. Eq. (15) is an M 2 by M 2matrix equation for an M by M mesh grid. However, each row of the coefficient matrix has no more than five non-zero values. As a result, this sparse matrix equation can be efficiently solved using various methods [3

3. W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vettering, Numerical recipes: The art of scientific computing (Cambridge University Press, New York, 1986).

]. In our calculations, the sparse matrix solver-UMFPACK package has been used [14

14. T.A. Davis, Direct Methods for Sparse Linear Systems (SIAM, 2006). [CrossRef]

].

4. Convergence studies of CJI

In this section, we investigate the convergence rate of the CJI method for 3D WA-BPM based on the modified Padé(1,1) approximant and Hadley(1,1). We simulate the propagation of a 3D Gaussian beam through a symmetric Y-branch waveguide. The structure parameters of the Y-branch waveguide are the same as in [17

17. Z. Ju, J. Fu, and E. Feng, “A simple wide-angle beam-propagation method for integrated optics,” Microwave Opt. Technol. Lett. 14, 345–347 (1997). [CrossRef]

] with d 1=2µm and d 2=1µm. In Fig. 3 we show that the CJI method converges faster with WA-BPM based on the modified Padé(1,1) approximant than those of Hadley(1,1). However, it suffers from the fact that the iteration count between two successive 2D cross sections increases throughout the propagation direction. To overcome this problem we propose the use of a perfectly matched layer (PML), which can absorb incident radiation without any additional parasitic reflections, regardless of wavelength, incident angle or polarization as boundary conditions [15

15. J. P. Bérenger, “A perfectly matched layer for the absorption of electromagnetic waves,” J. Comput. Phys. 114, 185–200 (1994). [CrossRef]

16

16. P. Vandersteengen, B. Maes, P. Bienstman, and R. Baets, “Using the complex Jacobi method to simulate Kerr non-linear photonic components,” Opt. Quantum Electron. 38, 35–44 (2006). [CrossRef]

]. It is seen that the CJI technique becomes more stable with the use of a PML as shown in Fig. 3.

Fig. 3. The iteration count per propagation step for propagation through a symmetric Y-branch waveguide obtained using the CJI method for standard Padé approximant-based WA BPM without (1) and with (2) PML and the CJI method for the modified Padé without (3) and with (4) PML. (um is defined as µm)

5. Benchmark results

We now employ the WA-BPM using the new CJI and the traditional DMI methods to perform benchmark tests on both 2D and 3D optical waveguide structures. The simulations were all run on a notebook PC using Matlab. For the 2D case, we consider a 10-degree tilted waveguide [7

7. T. Anada, T. Hokazono, T. Hiraoka, J. P. Hsu, T. M. Benson, and P. Sewell, “Very-wide-angle beam propagation methods for integrated optical circuits,” IEICE Trans. Electron. E82-C, 1154–1158 (1999).

] and Y-branch waveguide [17

17. Z. Ju, J. Fu, and E. Feng, “A simple wide-angle beam-propagation method for integrated optics,” Microwave Opt. Technol. Lett. 14, 345–347 (1997). [CrossRef]

]. In the tilted waveguide the fundamental mode for the slab of width w=0.2µm is propagated through 30µm at wavelength λ=1.55µm in a medium of refractive index n=3.4 with the propagation step size of Δz=0.1µm. With a very strict propagation error tolerance of 10-9 the CJI method only took 84.8 seconds, whereas the DMI method took 296.5 seconds.

In the Y-branch waveguide the parameters needed for calculation are the same as in [17

17. Z. Ju, J. Fu, and E. Feng, “A simple wide-angle beam-propagation method for integrated optics,” Microwave Opt. Technol. Lett. 14, 345–347 (1997). [CrossRef]

]. With a small propagation step size Δz=0.01µm (requiring frequent matrix inversion) the DMI method performed the propagation in 1414.8 seconds while the CJI method took only 584.1 seconds. It is obvious that for these 2D waveguide structures the CJI method is faster than DMI.

For the 3D case, we consider Gaussian beam propagation in a straight rib waveguide [18

18. P. C. Lee and E. Voges, “Three-dimensional semi-vectorial wide-angle beam propagation method,” J. Lightwave Technol. 12, 215–225 (1994) [CrossRef]

- 19

19. Y. Tsuji, M. Koshiba, and T. Shiraishi, “Finite element beam propagation method for three-dimensional optical waveguide structures,” J. Lightwave Technol. 15, 1728–1734 (1997). [CrossRef]

] and guided-mode propagation in a Y-branch rib waveguide. The width and height of the straight rib waveguide are w=2µm and h=1.1µm, as seen in Fig. 4 of [19

19. Y. Tsuji, M. Koshiba, and T. Shiraishi, “Finite element beam propagation method for three-dimensional optical waveguide structures,” J. Lightwave Technol. 15, 1728–1734 (1997). [CrossRef]

]. The guiding core has an index nf=3.44 and a thickness t=0.2µm while the refractive index of substrate and cover is ns=3.4 and nc=1, respectively. The Gaussian beam with a waist radius w 0=0.3µm has been injected into the rib waveguide at wavelength λ=1.55µm. Due to the large memory required for DMI, the small computational window of 2x2µm is discretized with a grid size of Δx= Δy=0.1µm, and the short path length of 2µm is discretized with a propagation step size Δz=0.1µm. The resulting runtime of DMI is 177.9 seconds while runtime for CJI is only 4.7 seconds.

For a Y-branch, the initial rib waveguide is split into two 5-degree tilted waveguides as shown in Fig. 4. The longitudinal dimension is h1=1µm. The other structure parameters are the same as the above straight rib waveguide. The fundamental TE mode of the ridge waveguide of width w=2µm at 1.55µm wavelength is used as the excited field at z=0. The propagation step size is Δz=0.1µm. The field pattern at z=3µm calculated by DMI and CJI are depicted in Fig. 5(a) and 5(b), respectively. Due to the high effective loss in the propagation medium the complex Jacobi method performed the propagation in only 5.9 seconds while DMI required 268.9 seconds.

Table 1. Quantitative comparison of runtimes of the direct matrix inversion and the complex Jacobi iteration for WA beam propagation in waveguide (WG) structures

table-icon
View This Table

Table 1 shows the performance of the two methods for the optical waveguide structures chosen here. It is clearly seen that the runtime of the iterative method is substantially lower than that of the DMI method. For large problems requiring very large storage space and also for structures with a long path length with small propagation step size that require frequent matrix inversions, the DMI technique is numerically very intensive. In contrast, for typical choices of k Δz the CJI technique offers rapid convergence and shorter runtimes.

Fig. 4. Y-branch optical rib waveguide
Fig. 5. Intensity contours for TE mode propagating in a 3D Y-branch rib waveguide at z=3µm calculated by (a) DMI and (b) the CJI method.

4. Conclusions

A new complex Jacobi iterative method adapted for the solution of 3D WA beam propagation has been presented. We proposed the modified Padé approximant of the WA beam propagation operator that gives the waves propagating in the evanescent region the desired damping. The resulting approach allows accurate approximations to the true Helmholtz equation. Furthermore, a quantitative comparison of runtimes between the traditional direct matrix inversion and the new complex Jacobi iterative method for both 2D and 3D WA beam propagation demonstrates convincingly that the complex Jacobi iterative method is very competitive for demanding problems. This solution technique will also enable a development of higher order 3D Padé approximant-based WA-beam propagation algorithms using the multistep method [20

20. G. R. Hadley, “Multistep method for wide-angle beam propagation,” Opt. Lett. 17, 1743–1745 (1992). [CrossRef] [PubMed]

], which will be presented in a future publication.

Acknowledgments

Parts of this work were performed within the context of the Belgian IAP project Photonics@Be. R. Godoy-Rubio would like to thank the Spanish CICYT Project TEC2006-02868, the Andalusian CICYE Project TIC-02946 and the “Juan de la Cierva” National Fellowship program.

References and Links

1.

R. Scarmozzino, A. Gopinath, R. Pregla, and S. Helfert, “Numerical techniques for modeling guided-wave photonic devices,” IEEE J. Sel. Top. Quantum Electron. 6, 150–161 (2000). [CrossRef]

2.

C. Ma and E. V. Keuren, “A simple three dimensional wide-angle beam propagation method,” Opt. Express 14, 4668–4674 (2006). [CrossRef] [PubMed]

3.

W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vettering, Numerical recipes: The art of scientific computing (Cambridge University Press, New York, 1986).

4.

M. D. Feit and J. A. Fleck Jr., “Analysis of rib waveguides and couplers by the propagation method,” J. Opt. Soc. Am. A 7, 73–79 (1990). [CrossRef]

5.

F. A. Milinazzo, C. A. Zala, and G. H. Brooke, “Rational square root approximations for parabolic equation algorithms,” J. Acoust. Soc. Am. 101, 760–766 (1997). [CrossRef]

6.

Y. Y. Lu and P. L. Ho, “Beam propagation method using a [(p-1)/p] Padé approximant of the propagator,” Opt. Lett. 27, 683–685 (2002). [CrossRef]

7.

T. Anada, T. Hokazono, T. Hiraoka, J. P. Hsu, T. M. Benson, and P. Sewell, “Very-wide-angle beam propagation methods for integrated optical circuits,” IEICE Trans. Electron. E82-C, 1154–1158 (1999).

8.

G. R. Hadley, “Wide-angle beam propagation using Padé approximant operators,” Opt. Lett. 17, 1426–1428 (1992). [CrossRef] [PubMed]

9.

S. Sujecki, “Wide-angle, finite-difference beam propagation in oblique coordinate system,” J. Opt. Soc. Am. 25, 138–145 (2007). [CrossRef]

10.

C. Ma and E. V. Keuren, “A three-dimensional wide-angle beam propagation method for optical waveguide structures,” Opt. Express 15, 402–407 (2007). [CrossRef] [PubMed]

11.

G. R. Hadley, “A complex Jacobi iterative method for the indefinite Helmholtz equation,” J. Comp. Phys. 203, 358–370 (2005). [CrossRef]

12.

Y. Y. Lu, “A complex coefficient rational approximation on of 1+x,” Appl. Numer. Math. 27, 141–156 (1998). [CrossRef]

13.

G. R. Hadley, “Low-truncation-error finite difference equations for photonics simulation I: beam propagation,” J. Lightwave Technol. 16, 134–141 (1998). [CrossRef]

14.

T.A. Davis, Direct Methods for Sparse Linear Systems (SIAM, 2006). [CrossRef]

15.

J. P. Bérenger, “A perfectly matched layer for the absorption of electromagnetic waves,” J. Comput. Phys. 114, 185–200 (1994). [CrossRef]

16.

P. Vandersteengen, B. Maes, P. Bienstman, and R. Baets, “Using the complex Jacobi method to simulate Kerr non-linear photonic components,” Opt. Quantum Electron. 38, 35–44 (2006). [CrossRef]

17.

Z. Ju, J. Fu, and E. Feng, “A simple wide-angle beam-propagation method for integrated optics,” Microwave Opt. Technol. Lett. 14, 345–347 (1997). [CrossRef]

18.

P. C. Lee and E. Voges, “Three-dimensional semi-vectorial wide-angle beam propagation method,” J. Lightwave Technol. 12, 215–225 (1994) [CrossRef]

19.

Y. Tsuji, M. Koshiba, and T. Shiraishi, “Finite element beam propagation method for three-dimensional optical waveguide structures,” J. Lightwave Technol. 15, 1728–1734 (1997). [CrossRef]

20.

G. R. Hadley, “Multistep method for wide-angle beam propagation,” Opt. Lett. 17, 1743–1745 (1992). [CrossRef] [PubMed]

OCIS Codes
(000.4430) General : Numerical approximation and analysis
(220.2560) Optical design and fabrication : Propagating methods
(350.5500) Other areas of optics : Propagation

ToC Category:
Physical Optics

History
Original Manuscript: August 1, 2008
Revised Manuscript: September 17, 2008
Manuscript Accepted: October 1, 2008
Published: October 10, 2008

Citation
Khai Q. Le, R. Godoy-Rubio, Peter Bienstman, and G. Ronald Hadley, "The complex Jacobi iterative method for three-dimensional wide-angle beam propagation," Opt. Express 16, 17021-17030 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-21-17021


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. R. Scarmozzino, A. Gopinath, R. Pregla, and S. Helfert, "Numerical techniques for modeling guided-wave photonic devices," IEEE J. Sel. Top. Quantum Electron. 6, 150-161 (2000). [CrossRef]
  2. C. Ma and E. V. Keuren, "A simple three dimensional wide-angle beam propagation method," Opt. Express 14, 4668-4674 (2006). [CrossRef] [PubMed]
  3. W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vettering, Numerical recipes: The art of scientific computing (Cambridge University Press, New York, 1986).
  4. M. D. Feit and J. A. Fleck Jr., "Analysis of rib waveguides and couplers by the propagation method," J. Opt. Soc. Am. A 7, 73-79 (1990). [CrossRef]
  5. F. A. Milinazzo, C. A. Zala, and G. H. Brooke, "Rational square root approximations for parabolic equation algorithms," J. Acoust. Soc. Am. 101, 760-766 (1997). [CrossRef]
  6. Y. Y. Lu and P. L. Ho, "Beam propagation method using a [(p-1)/p] Padé approximant of the propagator," Opt. Lett. 27, 683-685 (2002). [CrossRef]
  7. T. Anada, T. Hokazono, T. Hiraoka, J. P. Hsu, T. M. Benson, and P. Sewell, "Very-wide-angle beam propagation methods for integrated optical circuits," IEICE Trans. Electron.E 82-C, 1154-1158 (1999).
  8. G. R. Hadley, "Wide-angle beam propagation using Padé approximant operators," Opt. Lett. 17, 1426-1428 (1992). [CrossRef] [PubMed]
  9. S. Sujecki, "Wide-angle, finite-difference beam propagation in oblique coordinate system," J. Opt. Soc. Am. 25, 138-145 (2007). [CrossRef]
  10. C. Ma and E. V. Keuren, "A three-dimensional wide-angle beam propagation method for optical waveguide structures," Opt. Express 15, 402-407 (2007). [CrossRef] [PubMed]
  11. G. R. Hadley, "A complex Jacobi iterative method for the indefinite Helmholtz equation," J. Comp. Phys. 203, 358-370 (2005). [CrossRef]
  12. Y. Y. Lu, "A complex coefficient rational approximation of1+x ," Appl. Numer. Math. 27, 141-156 (1998). [CrossRef]
  13. G. R. Hadley, "Low-truncation-error finite difference equations for photonics simulation I: beam propagation," J. Lightwave Technol. 16, 134-141 (1998). [CrossRef]
  14. T. A. Davis, Direct Methods for Sparse Linear Systems (SIAM, 2006). [CrossRef]
  15. J. P. Bérenger, "A perfectly matched layer for the absorption of electromagnetic waves," J. Comput. Phys. 114, 185-200 (1994). [CrossRef]
  16. P. Vandersteengen, B. Maes, P. Bienstman, and R. Baets, "Using the complex Jacobi method to simulate Kerr non-linear photonic components," Opt. Quantum Electron. 38, 35-44 (2006). [CrossRef]
  17. Z. Ju, J. Fu, and E. Feng, "A simple wide-angle beam-propagation method for integrated optics," Microwave Opt. Technol. Lett. 14, 345-347 (1997). [CrossRef]
  18. P. C. Lee and E. Voges, "Three-dimensional semi-vectorial wide-angle beam propagation method," J. Lightwave Technol. 12, 215-225 (1994) [CrossRef]
  19. Y. Tsuji, M. Koshiba, and T. Shiraishi, "Finite element beam propagation method for three-dimensional optical waveguide structures," J. Lightwave Technol. 15, 1728-1734 (1997). [CrossRef]
  20. G. R. Hadley, "Multistep method for wide-angle beam propagation," Opt. Lett. 17, 1743-1745 (1992). [CrossRef] [PubMed]

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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited