`ISRN OpticsVolume 2013 (2013), Article ID 108704, 8 pageshttp://dx.doi.org/10.1155/2013/108704`
Research Article

## Perturbation Analysis with Approximate Integration for Propagation Mode in Two-Dimensional Two-Slab Waveguides

Department of Integrated Information Technology, College of Science and Engineering, Aoyama Gauin University, Sagamihara-shi, Kanagawa 252-5258, Japan

Received 23 May 2013; Accepted 14 July 2013

Copyright © 2013 Naofumi Kitsunezaki. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

On the basis of perturbation expansion from a gapless system, we calculate the propagation constant and propagation mode wave function in two-dimensional two-slab waveguides with a core gap small enough that there is only one propagation mode. We also perform calculations without the approximation for comparison. Our result shows that first-order perturbation contains the first-order Taylor expansion of (core gap)/(core width), and when the integration of the perturbation is suitably approximated, the result of the first-order perturbation is the same as that of the first-order Taylor expansion of (core gap)/(core width).

#### 1. Introduction

When monochromatic light enters into one of the cores of a waveguide array in which each core has the same width and the same gap, the light propagates to adjacent cores in turn, and its trajectory becomes V-shaped. Such optical behavior has been theoretically and numerically analyzed. Theoretical analysis includes matrix method [1], the coupled mode equation, and the coupled power equation [29]. In the coupling mode equations, it is assumed that there is at least one propagation mode for each core and that all of the propagation modes are independent from each other.

This assumption is invalid in the case that the number of independent propagation modes is less than the number of cores. Such a situation can occur when the core gap becomes small enough that the index distribution is approximated by a gapless core. When such a situation occurs, analysis based on the coupled mode or the coupled power equation is invalid. Our interest is finding a simple method to analyze an optical behavior in such a situation.

An optical behavior for parallel slab waveguides with small core gap was theoretically analyzed using even and odd mode analysis, which was called supermode analysis later [35]. In this analysis, the exact solution for the Maxwell’s equations with the index distribution for two parallel slab waveguides, which becomes essentially one for Schrödinger equation for double finite wells in the field of quantum mechanics, is used to analyze the optical behavior. However, the method using supermodes is not useful when the number of cores increases, because the solution must be individually expressed in every core and clad.

In the field of quantum mechanics, perturbation theory is widely used as an approximation [10, 11]. In the field of electromagnetics and optics, it is also explained in [12], and applied to explain bend losses for a fiber [13, 14]. Coupled mode equation is also based on the perturbation theory [47]. The situation that we are interested in is suitable for using the perturbation theory.

In this paper, we calculate the supermode for a two-dimensional two-slab waveguide that has a small core gap so that the optical system has only one propagation mode, using first-order perturbation from gapless system. We also calculate the same without approximation for comparison. In this way, we show that first-order perturbation with suitable approximate integration and the first-order Taylor expansion of (core gap)/(core width) give the same result. Furthermore, we discuss higher-order relationships between perturbation with suitable approximate integration and Taylor expansion.

#### 2. Model

For our optical model of two cores with a gap, the spatial configuration of the cores and clads are shown in the left panel of Figure 1. There exist two cores with the same index and the same core width , and clads with index , where and are constants, and . The core width , core gap , and indices and are chosen so that there is only one propagation mode. The core is configured in and . In other words, the index distribution is set as which is constant except for and and is illustrated in the right panel of Figure 1. The index distribution is constant in the and directions. The direction is defined as the direction of light propagation, and is used to denote the angular frequency of the monochromatic light.

Figure 1: Optical system configuration (a) and its index distribution on . line (b).

Since the index distribution has no variation in the direction, electromagnetic fields do not depend on . For the electric field and the magnetic field , the time variable can be separated from the spatial variables as and , because we are interested in monochromatic light with angular frequency . Then, from Faraday’s and Ampére-Maxwell’s laws, components of the electric and magnetic fields have the following relationships: It is straightforward to see that (2) to (7) are also solutions to the chargeless Gauss’s law for electric and magnetic fields, which are expressed by the rest of Maxwell’s equations. Equations (3) and (6) can be written as the second-order differential equations by using (5) and (7) and (2) and (4), respectively, where the derivative of is neglected in (9) because is constant except at some discrete points. Since (8) and (9) are basically the same differential equation, we introduce the field instead of and . Equations (8) and (9) are usually solved by separation of variables: So where is a constant. In (11), obviously cannot be positive if the region of is , because positive makes diverge in or . For any negative , there exists a positive such that where is the propagation constant, and the solution of (11) represents propagation in the positive or negative direction with .

Equation (12) is essentially the one-dimensional Schrödinger equation [10, 11], which is an eigenfunction problem, where is the eigenvalue and is the eigenfunction of the operator: Equation (12) has nontrivial solutions for special values in (propagation mode solution), and two nontrivial solutions exist for each in (radiation mode solution).

In the following section, we solve (12) by using first-order perturbation with a suitable approximation integration in order to calculate the difference of its propagation constant and propagation mode solution from that of a gapless system. In Section 4, we solve (12) without approximation in order to obtain the first-order Taylor expansion of , and then we compared the result with those obtained in Section 3. The final section is devoted to summary and discussion.

#### 3. Solution with First-Order Perturbation

To solve (12), we decompose the square of the index distribution into the square of the index distribution for () and the rest () as where Figure 2 shows the distribution of (solid line) and (dotted line) on the line. The index distribution of represents the core that is laid on . Instead of directly solving (12), we use the complete set of solutions for in order to solve (12) perturbatively. All independent solutions for (17) are where and (or the pair of and ) in (21) satisfy where In (22), is in the region

Figure 2: Decomposed index distributions. Solid line represents , and dotted line represents .

represents the propagation mode, and and represent symmetric and antisymmetric radiation modes according to , respectively. For the reasons stated in the introduction, we are interested in the case that the propagation mode of the optical system has only one independent mode. Therefore, we examine the case where the normalized propagation mode solution of (18) is only.

Here, let us define the inner product of fields as Then, , , and are defined so that , , and are normalized to where (note that and can be regarded as eigenfunctions of the eigenequation (17), and their eigenvalues are , , and , respectively. In addition to the relationships in (28), and [10, 11]). For instance, is In perturbation theory, the propagation mode solution for (12) is determined by using (18) to (20), which results in The results of the first-order perturbation for and are as follows [10, 11]: where In (33), and mean the and components of , respectively. Note that for any field , the inner product with is which means that it contains the first-order Taylor expansion of . In the deformation from (35) to (36), we use the Euler method for the integration. Using (18) through (20), we obtain the first-order Taylor expansion of for , , and from (31) and (33), respectively, as follows:

#### 4. Solution without Approximation

When the exact eigenvalue of the propagation mode of (12) is , its eigenfunction is where and satisfy the relationships such that and its derivative are continuous for any . Note that agrees with in the limit of , and becomes when . When , the relationships between and and between and become the following: Substituting (44) into (42) by using (41), and applying (23), we obtain Comparing the lowest order of and in (45) using (23), we obtain (37) from (45). Because of (13), the first-order perturbation and the first-order Taylor expansion of give the same propagation constant for a two-core optical system with gap .

and are directly derived, which results in

From (46), we obtain the first-order Taylor expansion of of and after a straightforward but lengthy calculation: and are equal to and , respectively, because of (32) and (43). As shown above, (38) and (39) completely coincide with (47) and (48), respectively.

#### 5. Conclusion and Discussion

We have shown that in a two-dimensional two-slab waveguide system that has a small core gap so that the optical system has only one propagation mode, the supermode and its propagation constant calculated by first-order perturbation by the Euler method and the first-order Taylor expansion of for the exact solution give the same result up to order.

As for higher-order perturbation, the correction of by second-order perturbation is given by [10] where represents mode density function, and the correction of is The lowest order of for (49) and (50) is , because of the argument given in Section 3. Therefore, second-order perturbation with the Euler method is a part of the second-order Taylor expansion of . As shown in (36), the first-order perturbation can include the second order of . To obtain the second order of from the first-order perturbation, integration over and in the perturbation must be approximated by the trapezoidal rule instead of the Euler method. Thus, first-order perturbation by the trapezoidal rule in combination with a second-order perturbation by the Euler method gives the same result as the exact calculation given by the second-order Taylor expansion of . In a similar way, the th-order coefficient of the Taylor expansion of is obtained from all orders of perturbation less than the th order. In each order of perturbation, exact integration over and is not necessary; instead, a suitable approximation of integration including the order term is necessary. Thus, when an optical system has two cores with a small gap, it can be analyzed by perturbation using modes of an optical system with a gapless core with twice the width.

Obtaining the exact solution becomes time consuming as the number of cores increases. Even in such situations, perturbation analysis can be executed easily. Hence, perturbation analysis using suitable approximate integration is a promising candidate for approximation in the region when the coupled mode or power mode equation is invalid because is very small.

#### Acknowledgment

The authors thank all members of the Tobe Laboratory in the Department of Integrated Information Technology, College of Science and Engineering, Aoyama Gakuin University, for comments and discussion.

#### References

1. A. Ghatak and V. Lakshminarayanan, “Propagation characteristics of planar waveguides,” in Optical Waveguides from Theory to Applied Technologies, M. L. Calvo and V. Lakshminarayanan, Eds., p. CRC press, New York, NY, USA, 2007.
2. A. Yariv, “Coupling-mode theory for guided-wave optics,” IEEE Journal of Quantum Electronics, vol. QE-9, no. 9, pp. 919–933, 1973.
3. Y. Suematsu and K. Kishino, “Coupling coefficient in strongly coupled dielectric waveguides,” Radio Science, vol. 12, no. 4, pp. 587–592, 1977.
4. A. Hardy and W. Streifer, “Coupled mode theory of parallel waveguides,” Journal of Lightwave Technology, vol. 3, no. 5, pp. 1135–1146, 1985.
5. A. Hardy and W. Streifer, “Coupled modes of multiwaveguide system and phased arrays,” Journal of Lightwave Technology, vol. 4, no. 1, pp. 90–99, 1986.
6. H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” Journal of Lightwave Technology, vol. 5, no. 1, pp. 16–23, 1987.
7. K. Yasumoto, “Coupled-mode formulation of parallel dielectric waveguides using singular perturbation technique,” Microwave and Optical Technology Letters, vol. 4, no. 11, pp. 486–491, 1991.
8. H. Yoshida, K. Hidaka, and N. Tsukada, “Proposal of optical waveguide array beam splitter and Mach-Zehnder optical interferometer,” IEICE Technical Report OCS2010-36, 2010.
9. A. Komiyama, “Energy conservation in a waveguide system with an imperfection core,” The Papers of Technical Meeting on Electromagnetic Theory, IEEJ, EMT-11-131, 2011.
10. L. I. Shiff, Quantum Mechanics, Mc-Graw Hill, Singapore, 1969.
11. W. Greiner, Quantum Mechanics—An Introduction, Springer, Berlin, Germany, 1989.
12. R. F. Harrington, Time-Harmonic Electromagnetic Fields, Wiley-Interscience, New York, NY, USA, 2001.
13. C. Vassallo, “Perturbation of a LP mode of an optical fibre by a quasi-degenerate field: a simple formula,” Optical and Quantum Electronics, vol. 17, no. 3, pp. 201–205, 1985.
14. Q. Wang, G. Farrell, and T. Freir, “Theoretical and experimental investigations of macro-bend Losses for standard single mode fibers,” Optics Express, vol. 13, no. 12, pp. 4476–4484, 2005.