# The threshold mode structure analysis of the two-dimensional photonic crystal lasers.

1. INTRODUCTIONRecently, there has been increasing interest in two-dimensional (2D) photonic crystal (PC) lasers utilizing a 2D distributed feedback (DFB) mechanism. They have been attracting much attention because of their important characteristics, such as single mode operation in a large cavity area and a diffraction-limited circular-shape output beam [1,2]. The PC laser demonstrations were obtained under pulse excitation, weather optical or electrical, as well as continues wave regime [3-11]. Simultaneously, theoretical models of PC laser structures have been developed [12-19]. Particularly, coupled-wave models for square and triangular lattices 2D PC lasers with transverse electric and transverse magnetic like modes have been presented. In [16, 17], the coupled-wave equations and field distributions for square lattice PC laser were shown for TE and TM modes, respectively.

In [15], the expressions for the resonant frequencies of mode derived from the coupled-mode equations describing the characteristics of experimental results for the band-edge frequencies of the 2D PC laser are shown. In the case of TM-like modes [17], the threshold characteristics for this kind of laser structures have been obtained. It is also shown that this lasing mode can be selected by manipulating the filling factor or the boundary reflection. In order to further develop the PC laser, it seems to be important to investigate different crystal geometries such as triangular lattice. A PC laser with triangular lattice geometry has been studied [3, 4] and was found to have six resonant modes at the edge of photonic band structure. Some of the resonant properties have been studied in further works [7, 12]. Recently, coupled-wave model for triangular-lattice photonic crystal with transverse electric polarization has been presented [18]. The analytic expressions that describe the relations between the resonant mode frequencies and coupling constant have been derived. However, it is also important to understand the resonant modes with TM polarization, since for example the promising light sources for TM region (e.g., quantum cascade lasers) have their gain in TM polarization [20]. In this paper, two cases of photonic crystal symmetries (square and triangular) and TE and TM polarizations, for each structure, are analyzed. In comparison with earlier works, this one gathers all four possible scenarios and analyzes them in the wide range of the coupling coefficient values, extending the two-dimensional coupled-wave model for triangular lattice photonic crystal laser to describe threshold behavior of the TM-like modes.

2. COUPLED-WAVE MODELS

Throughout this work a PC structure in a fixed square area with circular holes arranged in square or triangular lattice is considered. Following works of, e.g., Plihal and Maradudin [14] and Sakai et al. [2,16-18], four cases are distinguished: PC cavity with square and triangular lattices with, in each case, TE and TM modes. The discussed structures are schematically shown in Figure 1. These structures are assumed to be uniform and not confined in the z direction.

[FIGURE 1 OMITTED]

The scalar wave equations for the electric and magnetic fields, [E.sub.z] and [H.sub.z], respectively, are written in the following form [14]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

In Equations (1) and (2), the constant fc is given by:

[k.sup.2] = [[beta].sup.2] + 2i ([alpha] - [[alpha].sub.L])[beta] + 2[beta] [summation over (G [not equal to] 0)] k(G) exp (i(G x r)) (3)

in case of TM mode [17] and by

1/[k.sup.2] = 1/[[beta].sup.4] ([[beta].sup.2] - i2([alpha] - [[alpha].sub.L])[beta] + 2[beta] [summation over (G [not equal to] 0)] K(G) exp (i(G x r)))(4)

in case of TE mode [18]. In Equations (3) and (4) [beta] = 2[pi][[epsilon].sup.1/2.sub.0]/[lambda] where [[epsilon].sub.0] = [epsilon](G = 0) is the averaged dielectric permittivity ([[epsilon].sup.1/2.sub.0] corresponds to averaged refractive index n), [alpha] an averaged gain in the medium, K(G) the coupling constant, [lambda] the Bragg wavelength, G = (m[b.sub.1], n[b.sub.2]) the reciprocal lattice vector, and m and n the arbitrary integers. [b.sub.1] and [b.sub.2] vary depending on the structure symmetry. Therefore, these vectors are expressed in the following forms [b.sub.1] = ([[beta].sup.s.sub.0], 0) and [b.sub.2] = (0, [[beta].sup.s.sub.0]) for square lattice, and [b.sub.1] = ([[beta].sup.t.sub.0], 0) and [b.sub.2] = (- [[beta].sup.t.sub.0]/2, [square root of 3[[beta].sup.t.sub.0]/2) for triangular lattice structure, where [[beta].sup.s.sub.0] = 2[pi]/a and [[beta].sup.t.sub.0] = 4[pi]/[square root of 3a]. In the derivation of Equations (3) and (4) following e.g., [17, 21], we set [alpha] [much less than] [beta] = 2[pi][[epsilon].sup.1/2.sub.0]/[lambda], [[epsilon].sub.G [not equal to] 0] [much less than] [[epsilon].sub.0], and [[alpha].sub.G] [much less than] [beta]. In these equations, the periodic variation in the refractive index is included as a small perturbation and appears in the third term through the coupling constant K(G) of the form:

K(G) = -[pi]/[lambda][[epsilon].sup.1/2.sub.0] [epsilon](G) [+ or -] i [alpha](G)/2 (5)

In (5), plus sign refers to TM polarization (3), while minus sign refers to TE polarization (4). Further, we set [[alpha](G).sub.G[not equal to] 0] = 0 neglecting spatial periodicity of gain. In the vicinity of the Bragg wavelength, some of the diffraction orders contribute to the coupling of waves in more significant way than the others. In general, a periodic perturbation produces an infinite set of diffraction orders. The Bragg frequency corresponding to the [GAMMA] point in the photonic band structure [12] is chosen for the purpose of this paper, and the most significantly contributing coupling constants are expressed as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

In Equations (1) and (2), electric and magnetic fields for the infinite periodic structure are given by the Bloch mode [14]:

[E.sub.z] (r) = [summation over (G)] e(G) exp (i(k + G) x r) (7)

[H.sub.z] (r) = [summation over (G)] h(G) exp (i(k + G) x r) (8)

where functions e(G) and h(G) correspond to the plane wave amplitudes, whereas wave vector is denoted by k. It is worth noting here that for the more general three-dimensional case with vertical confinement (i.e., in z direction) only the waves diffracted in the plane should be included, e.g., [24]. In the first Brillouin zone k = 0 at the [GAMMA] point [16-18]. For a finite structure, the amplitude of each plane wave is not a constant, so e(G) and h(G) become functions of space. At the r point, we consider only the amplitudes (e(G), h(G)) which are the most significant and let include one- and two-dimensional coupling effects, i.e., in most cases with [absolute value of G] = [[beta].sup.s,t.sub.0], except for square lattice with TE polarization where additional h(G) amplitudes with [absolute value of G] = [square root of 2[[beta].sup.s.sub.0] have to be included [15]. For the purpose of this article, the contributions of higher order waves in the Bloch mode are considered negligible, which in general is not true, especially for non circular holes, e.g., [24,25].

2.1. Square Lattice

2.1.1. TM Polarization

Considering square lattice photonic crystal with TM polarization, it is assumed that at the r point, the most significant contribution to coupling is given by the electric waves, which fulfill the condition ([absolute value of G] = [[beta].sup.s.sub.0]). Thus, all higher order electric wave expansion coefficients ([absolute value of G] [greater than or equal to] [square root of 2[[beta].sub.0])] are neglected. Four basic waves most significantly contributing to coupling are depicted in Figure 2.

Equation (7) describes infinite structures. It is possible to take into account the fact that the structure is finite by using the space dependent amplitudes, e.g., [15]. Thus, the electric field (7) in the finite periodic structure can be expressed in the following way:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

In (9) [E.sup.s.sub.1], i = 1... 4 are the four basic electric field amplitudes of waves propagating in four directions +x, - x, +y, y. These amplitudes correspond to e(G) in Equation (7). In further analysis, the space dependence notation is omitted.

Putting Equations (3) and (9) into Equation (1), and assuming the slow varying electromagnetic field, one can get the set of coupled mode equations [17]:

-[partial derivative]/[partial derivative]x [E.sup.s.sub.1] + ([alpha] - [[alpha].sub.L] - [k.sub.0] - i[xi]) [E.sup.s.sub.1] = (i[k.sub.3] - [k.sub.0]) [E.sup.s.sub.3] + i[k.sub.2]([E.sup.s.sub.2] + [E.sup.s.sub.4] (10)

[partial derivative]/[partial derivative]x [E.sup.s.sub.3] + ([alpha] - [[alpha].sub.L] - [k.sub.0] - i[xi]) [E.sup.s.sub.3] = (i[k.sub.2] - [k.sub.0]) [E.sup.s.sub.1] + i[k.sub.2]([E.sup.s.sub.2] + [E.sup.s.sub.4]) (11)

-[partial derivative]/[partial derivative]y [E.sup.s.sub.2] + ([alpha] - [[alpha].sub.L] - [k.sub.0] - i[xi]) [E.sup.s.sub.2] = (i[k.sub.3] - [k.sub.0]) [E.sup.s.sub.4] + i[k.sub.2]([E.sup.s.sub.1] + [E.sup.s.sub.3]) (12)

[partial derivative]/[partial derivative]y [E.sup.s.sub.4] + ([alpha] - [[alpha].sub.L] - [k.sub.0] - i[xi]) [E.sup.s.sub.4] = (i[k.sub.3] - [k.sub.0]) [E.sup.s.sub.2] + i[k.sub.2]([E.sup.s.sub.1] + [E.sup.s.sub.3]) (13)

where

[xi] = ([[beta].sup.2] - [[beta].sup.s2.sub.0])/2[beta] [approximately equal to] [beta] - [[beta].sup.s.sub.0] (14)

is the Bragg frequency deviation, [k.sub.2] and [k.sub.3] are coupling coefficients [18]. The [k.sub.2] coefficient is responsible for orthogonal coupling (e.g., the coupling of [E.sup.s.sub.1] to [E.sup.s.sub.2] and [E.sup.s.sub.4]), and [k.sub.2] corresponds to backward coupling (e.g., the coupling of [E.sup.s.sub.1] to [E.sup.s.sub.3]). The additional coefficient [k.sub.0] denotes surface emission losses [16,17]. Solution of Equations (10)-(13) for the boundary conditions:

[E.sup.s.sub.1] (-L/2, y) = [E.sup.s.sub.3] (L/2, y) = 0, [E.sup.s.sub.2] (x, -L/2) = [E.sup.s.sub.4] (x, L/2) = 0 (15)

defines eigenmodes of the photonic structure. The analysis of this solution will be shown in Section 3.

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

2.1.2. TE Polarization

In the square lattice photonic crystal cavity with TE polarization, as mentioned before, the coupling process involves magnetic waves satisfying following conditions ([absolute value of G] = [[beta].sub.0]) and ([absolute value of G] = [square root of 2[[beta].sub.0]) [16], neglecting higher order Bloch modes. Eight basic waves most significantly contributing to coupling are depicted in Figure 3.

Similar as in the case of TM polarization, the equation for magnetic field (8) describes modes for infinite structure. Thus, the finite dimensions of the structure are described by spatial dependence of magnetic field amplitudes [16], and the magnetic field (8) is written in the following form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

In Equation (16), [H.sup.s.sub.i], i = 1... 8 are the eight basic magnetic field amplitudes of waves propagating in directions schematically shown in Figure 3. These amplitudes correspond to h(G) in Equation (8). Joining Equations (4), (16), and (2), and assuming slowly varying amplitudes, the coupled-wave equations for TE modes in square lattice PC are obtained:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

In the derivation of the above equations, the derivatives of slowly varying amplitudes [H.sup.s.sub.i], i = 2, 4, 6, 8 and terms much smaller than [xi] were neglected. In Equations (17)-(20), S is the Bragg frequency deviation, given by (14). The expressions for the coupling coefficients [k.sub.1], [k.sub.2], and [k.sub.3] can be found in [16,18].

In contrast to TM polarization, in Equations (17)-(20), the coupling coefficient [k.sub.2] responsible for coupling in perpendicular direction vanishes. The coupling coefficient [k.sub.3] has the same meaning as described in the previous case, whereas the coupling coefficient [k.sub.1] describes, for example, the coupling of waves [H.sup.s.sub.1], [H.sup.s.sub.2], and [H.sup.s.sub.8].

Solution of Equations (17)-(20) while taking into account boundary conditions:

[H.sup.s.sub.7] (- L/2, y) = [H.sup.s.sub.5] (L/2, y) = 0, [H.sup.s.sub.3] (x, -L/2) = [H.sup.s.sub.7] (x, L/2) = 0 (21)

defines structure eigenmodes at lasing threshold, i.e., in the linear case.

2.2. Triangular Lattice

2.2.1. TM Polarization

In the triangular lattice photonic crystal cavity with TM polarization, the coupling process involves waves satisfying following conditions ([absolute value of G] = [[beta].sub.0]) [18,19], neglecting higher order Bloch modes. Six basic waves most significantly contributing to coupling are depicted in Figure 4.

The space dependent amplitudes for electric field e(G) (Equation (7)) in triangular lattice photonic crystal cavity are written in the following form [19]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

In Equation (22), [E.sup.t.sub.i], i = 1... 6, are the six electric field amplitudes propagating in the symmetry directions (Figure 4). Combining Equations (3), (22) and (1), and assuming slowly varying amplitudes, the coupled-wave equations for TM modes in triangular lattice PC are obtained:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (27)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (28)

In Equations (23)-(28), as in the case of square lattice, [xi] is the Bragg frequency deviation, given by (17), while [k.sub.1], [k.sub.2], and [k.sub.3] are the coupling coefficients, defined by (6), and expressed in, e.g., [19]. These coefficients describe strength and direction of the coupling of the waves, e.g., the coupling of [E.sup.t.sub.1] and [E.sup.t.sub.4] is described by [k.sub.3], coupling of [E.sup.t.sub.1], [E.sup.t.sub.2], and [E.sup.t.sub.6] by [k.sub.1], and coupling of [E.sup.t.sub.1], [E.sup.t.sub.3], and [E.sup.t.sub.5] by [k.sub.2]. In Equations (23)-(28), there is an additional coefficient [k.sub.0] which, as in the square lattice case, is responsible for surface emission losses [7, 22].

Solution of Equations (23)-(28) for the boundary conditions

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (29)

defines structure eigenmodes at lasing threshold. 2.2.2. TE Polarization

In the triangular lattice photonic crystal cavity with TE polarization, the coupling process involves waves satisfying the same condition as it was stated in TM polarization case, i.e., ([absolute value of G] = [[beta].sub.0]) [18], neglecting higher order Bloch modes. Six basic waves most significantly contributing to coupling are depicted in Figure 5.

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

The magnetic field amplitudes h(G) (Equation (8)) in the triangular lattice photonic crystal cavity are written as follows [18]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (30)

In Equation (43), [H.sup.t.sub.i], i = 1... 6, are the six magnetic field amplitudes propagating in the symmetry directions (Figure 5). Combining Equations (4), (30) and (2), and assuming slowly varying magnetic field amplitudes, the coupled-wave equations for TE modes in triangular lattice PC are obtained:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (32)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (33)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (34)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (35)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (36)

where the coupling coefficients [k.sub.1], [k.sub.2], and [k.sub.3] have the same physical meaning as described in the TM polarization case. The boundary conditions for the square region of PC with triangular symmetry are written as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (37)

3. NUMERICAL ANALYSIS OF THE PC LASER THRESHOLD OPERATION

3.1. Square Lattice

3.1.1. TM Polarization

In Figure 6, an enlarged area of a square lattice photonic crystal dispersion curves for the first four modes (A, B, C, D) in the vicinity of [GAMMA] point is shown. The plane wave method [23] was used to plot the dispersion characteristic for the infinite two-dimensional PC structure with circular holes [[epsilon].sub.b] = 9.8 arranged in square lattice with background material [[epsilon].sub.a] = 12.0. The ratio of rods radius to lattice constant was set to 0.24, and the number of plane waves was set at about 1024.

In Figure 6, each of the selected (A, B, C, D) points represents a mode, characterized by Bragg frequency deviation [xi], threshold gain [alpha], and threshold field distribution. These characteristic values were calculated by the numerical solution of Equations (10)-(13). In order to assign appropriate points A, B, C, D to the obtained numerical values, it was necessary to use the analytic expressions for the Bragg frequency deviation [15]:

[[xi].sub.A] = -2[k.sub.2] - [k.sub.3], [[xi].sub.B,C] = [K.sub.3], [[xi].sub.D] = 2[k.sub.2] - [k.sub.3]. (38)

[FIGURE 6 OMITTED]

These expressions were obtained from Equations (10)-(13) where no gain ([alpha] = 0), no loss ([k.sub.0] = 0, [[alpha].sub.L] = 0), and no spatial dependence of electric field amplitude were assumed.

Equations (10)-(13) were solved numerically for the wide range of coupling coefficients ([k.sub.2], [k.sub.3]), obtaining solutions [([([xi], [alpha], [E.sup.s.sub.l]).sup.j]).sub.K3i], where l = 1, 2, 3, 4, [k.sub.3i] corresponds to subsequent values of coupling coefficient for different modes j = A, B, C, D. Assigning numerical values of [[xi].sub.j] to analytical solutions (38), the mode structure of 2D square lattice PC laser with TM polarization has been obtained.

Figure 7 shows the field distributions [[absolute value of [E.sup.s.sub.1]].sup.2] + [[absolute value of [E.sup.s.sub.2]].sup.2] + [[absolute value of [E.sup.s.sub.3]].sup.2] + [[absolute value of [E.sup.s.sub.4]].sup.2] corresponding to modes: A--Figure 7(a), D--Figure 7(b), B, C --Figures 7(b), (c). They were made for the normalized coupling coefficients [absolute value of [k.sub.2]L] = 8, [absolute value of [k.sub.3]L] = 4 and filling factor f = 0.16.

[FIGURE 7 OMITTED]

In Figure 8, the normalized threshold gain al was plotted as a function of Bragg frequency deviation [xi], for various values of the normalized coupling coefficient [k.sub.3]L (which takes values from 0.01 to 50).

Figure 8 shows that increasing values of coupling coefficient refer to the Bragg frequency deviation increment and threshold gain decrement. Simultaneously for larger values of coupling coefficient, the threshold gain tends to be similar values. It is also worth noting that the threshold gain values for mode A are the lowest in vide range of coupling coefficient.

[FIGURE 8 OMITTED]

[FIGURE 9 OMITTED]

[FIGURE 10 OMITTED]

Figure 9 illustrates mode A (shown in Figure 6) development on the threshold gain a vs. coupling coefficient [k.sub.3] curve.

3.1.2. TE Polarization

In Figure 10, similar as in TM polarization case, an enlarged area of a square lattice photonic crystal dispersion curves for the first four modes (A, B, C, D) in the vicinity of r point is shown. The plane wave method [23] was used with the same parameters as done for TM polarization.

Each of the selected (A, B, C, D) points represents a mode, characterized by Bragg frequency deviation S, threshold gain a, and threshold field distribution. These characteristic values were calculated by the numerical solution of Equations (17)-(20). In order to assign appropriate points A, B, C, D to the obtained numerical values, as for TM polarization, it was necessary to use the analytic expressions for the Bragg frequency deviation [17]:

[[xi].sub.A] = -8[K.sup.2.sub.1]/[[beta].sub.0] - [K.sub.3], [[xi].sub.B] = -[K.sub.3], [[xi].sub.C,D] = - 4[K.sup.2.sub.1]/[[beta].sub.0] + [K.sub.3]. (39)

These expressions were obtained from Equations (17)-(20) where no gain ([alpha] = 0), no loss ([k.sub.0] = 0, [[alpha].sub.L] = 0), and no spatial dependence of magnetic field amplitude were assumed.

Equations (17)-(20) were solved numerically for the wide range of coupling coefficients ([k.sub.1], [k.sub.3]), obtaining solutions [([([xi], [alpha], [H.sup.s.sub.l]).sup.j]).sub.K3i], where l = 1, 3, 5, 7, [K.sub.3i] corresponds to subsequent values of coupling coefficient for different modes j = A, B, C, D. Assigning numerical values of [[xi].sub.j] to analytical solutions (39), the mode structure of 2D square lattice PC laser with TE polarization has been obtained. Figure 11 shows the field distributions [[absolute value of [H.sup.s.sub.1]].sup.2] + [[absolute value of [H.sup.s.sub.3]].sup.2] + [[absolute value of [H.sup.s.sub.5]].sup.2] + [[absolute value of [H.sup.s.sub.7]].sup.2] corresponding to modes: A--Figure 11(a), B--Figure 11(b), C, D --Figures 11(c), (d). They were made for the normalized coupling coefficients [absolute value of [k.sub.1]L] = 10.96, [absolute value of [k.sub.3]L] = 4 and filling factor f = 0.16.

In Figure 12, the normalized threshold gain aL was plotted as a function of Bragg frequency deviation S, for various values of the normalized coupling coefficient [k.sub.3]L (which takes values from 0.01 to 50).

Figure 12 shows that increasing values of coupling coefficient refer to the Bragg frequency deviation increment and the threshold gain decrement. Simultaneously, for larger values of coupling coefficient the threshold gain tends to be similar values.

Figure 13 depicts mode A (shown in Figure 10) development with increasing values of coupling coefficient [k.sub.3] on the threshold gain [alpha] curve. Compared to square lattice transverse magnetic polarization case, it is noticed that TE modes develop higher field confinement in the structure for the greater difference of refractive indices (i.e., higher values of coupling coefficient).

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

[FIGURE 13 OMITTED]

3.2. Triangular Lattice

3.2.1. TM Polarization

Similar as in case of square lattice, in Figure 14 an enlarged area of a triangular lattice photonic crystal dispersion curves for the first six modes (A, B, C, D, E, F) in the vicinity of r point is shown. The plane wave method [23] was used to plot the dispersion characteristic for the infinite two-dimensional PC structure with circular holes [[epsilon].sub.b] = 9.8 arranged in triangular lattice with background material [[epsilon].sub.a] = 12.0. The ratio of rods radius to lattice constant was set to 0.24, and the number of plane waves to 1024.

In Figure 14, each of the selected (A, B, C, D, E, F) points represents a mode, characterized by Bragg frequency deviation [xi], threshold gain a, and threshold field distribution. These characteristic values were calculated by the numerical solution of Equations (23)-(28). In order to assign appropriate points A, B, C, D, E, F to the obtained numerical values, it was necessary to express the Bragg frequency deviation through the coupling coefficients:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (40)

Similar to previous cases, expressions (40) were obtained from coupled-mode Equations (23)-(28) where no gain ([alpha] = 0), no loss ([k.sub.0] = 0, [[alpha].sub.L] = 0), and no spatial dependence of electric field amplitude were assumed. Equations (23)-(28) were solved numerically for the wide range of coupling coefficients ([k.sub.1], [k.sub.2], [k.sub.3]), obtaining solutions [([([xi], [alpha], [E.sup.t.sub.k]).sup.j]).sub.K3i], where [K.sub.3i] corresponds to subsequent values of coupling coefficient for different modes j = A, B, C, D, E, F; k = 1... 6. Assigning numerical values of [[xi].sub.j] to the analytical solutions (40), the mode structure of 2D triangular lattice PC laser with TM polarization has been obtained. Figure 15 shows the field distributions [6.summation over (i = 1)] [[absolute value of [E.sup.t.sub.i]].sup.2] corresponding to the modes: A--Figure 15(a), F--Figure 15(b), B, C--Figures 15(c), (d), D, E--Figures 15(e), (f). They were made for the normalized coupling coefficients [absolute value of [k.sub.1]L] = 13.96, [absolute value of [k.sub.2]L] = 6.6, [absolute value of [k.sub.3]L] = 4, and filling factor f = 0.16.

[FIGURE 14 OMITTED]

In Figure 16, the normalized threshold gain [[alpha].sub.L] was plotted as a function of Bragg frequency deviation [xi], for various values of the normalized coupling coefficient [k.sub.3]L (which takes values from 0.01 to 50).

Figure 16 shows that increasing values of coupling coefficient refer to the Bragg frequency deviation increment and the threshold gain decrement. Simultaneously, for larger values of coupling coefficient the threshold gain tends to be similar values. The difference in the values of degenerate modes' threshold gains stem from numerical inaccuracy.

[FIGURE 15 OMITTED]

[FIGURE 16 OMITTED]

Similar as shown in square lattice case, here in Figure 17 mode A (showed in Figure 14) development is shown on the threshold gain a vs. coupling coefficient [k.sub.3] curve.

3.2.2. TE Polarization

In Figure 18, an enlarged area of a triangular lattice photonic crystal dispersion curves for the first six modes (A, B, C, D, E, F) in the vicinity of [GAMMA] point is shown. The plane wave method was used to plot the dispersion characteristic with the same parameters as done for TM polarization.

In Figure 18, each of the selected points (A, B, C, D, E, F) represents a mode, characterized by Bragg frequency deviation [xi], threshold gain [alpha], and threshold field distribution. These characteristic values were calculated by the numerical solution of Equations (31)(36). In order to assign appropriate points A, B, C, D, E, F to the obtained numerical values, it was necessary to express Bragg frequency deviation through the coupling coefficients:

[[xi].sub.A] = -[K.sub.1] - [K.sub.2] + [K.sub.3], [[xi].sub.B,C] = -[K.sub.1]/2 + [K.sub.2]/2 - [K.sub.3], [[xi].sub.D] = [K.sub.1] - [K.sub.2] - [K.sub.3], [[xi].sub.E,F] = [K.sub.1]/2 + [K.sub.2]/2 + [K.sub.3]. (41)

[FIGURE 17 OMITTED]

[FIGURE 18 OMITTED]

Expressions (41) were obtained from coupled mode Equations (31)-(36) where no gain ([alpha] = 0), no loss ([k.sub.0] = 0, [[alpha].sub.L] = 0), and no spatial dependence of magnetic field amplitude were assumed. Equations (31)-(36) were solved numerically for the wide range of coupling coefficients ([k.sub.1], [k.sub.2], [k.sub.3]), obtaining solutions [([([xi], [alpha], [H.sup.t.sub.k]).sup.j]).sub.K3i], where [k.sub.3i] corresponds to subsequent values of coupling coefficient for different modes j = A, B, C, D, E, F; k = 1 ... 6. Assigning numerical values of [[xi].sub.j] to the analytical solutions (41), the mode structure of 2D triangular lattice PC laser with TE polarization has been obtained. Figure 19 shows the field distributions [6.summation over i = 1] [[absolute value of [H.sup.t.sub.i]].sup.2] corresponding to the modes: A --Figure 19(a), D--Figure 19(b), B, C--Figures 19(c), (d), E, F --Figures 19(e), (f).

Plots in Figure 19 were made for the same normalized coupling coefficients as in the case of TM polarization, i.e., [absolute value of [K.sub.1]L] = 13.96, [absolute value of [k.sub.2]L] = 6.6, [absolute value of [k.sub.3]L] = 4, and filling factor f = 0.16.

In Figure 20, the normalized threshold gain [[alpha].sub.L] was plotted as a function of Bragg frequency deviation [xi], for various values of the normalized coupling coefficient [k.sub.3]L (which takes values from 0.01 to 50).

[FIGURE 19 OMITTED]

[FIGURE 20 OMITTED]

[FIGURE 21 OMITTED]

Figure 20 shows that increasing values of coupling coefficient refer to the Bragg frequency deviation increment and threshold gain decrement. Simultaneously, for larger values of coupling coefficient the threshold gain tends to be similar values. The difference in the values of degenerate modes' threshold gains stem from numerical inaccuracy.

Figure 21 illustrates mode A (showed in Figure 18) development with increasing values of coupling coefficient [k.sub.3] on the threshold gain [alpha] curve. Compared to TM polarization case, as in square symmetry structures, it can be noticed that TE modes develop higher field confinement in the structure for the greater difference of refractive indices (i.e., higher values of coupling coefficient). Thus in general it can be stated that TM modes are better confined in the structures with small refractive indices difference.

4. CONCLUSIONS

We have presented a coupled-wave analysis for square and triangular lattices photonic crystal laser with TE and TM polarization. We used coupled-mode equations to analyze PC structures in the wide range of coupling coefficient values distinguishing basic modes. We have plotted threshold characteristics showing lowest threshold modes and their frequency deviations. We also illustrate the lowest frequency mode development with increasing values of coupling coefficient and point out that TM modes are most likely to occur in structures with low dielectric constants contrast.

Received 12 January 2012, Accepted 17 February 2012, Scheduled 5 March 2012

REFERENCES

[1.] Miyai, E., K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, "Lasers producing tailored beams," Nature, Vol. 441, 946, 2006.

[2.] Sakai, K., E. Miyai, T. Sakaguchi, D. Ohnishi, T. Okano, and S. Noda, "Lasing band-edge identification for a surface-emitting photonic crystal laser," IEEE J. Sel. Areas Commun., Vol. 23, No. 7, 1335-1340, 2005.

[3.] Imada, M., S. Noda, A. Chutinan, T. Tokuda, M. Murata, and G. Sasaki, "Coherent two-dimensional lasing action in surface-emitting laser with triangular-lattice photonic crystal structure," Appl. Phys. Lett., Vol. 75, No. 3, 316-318, 1999.

[4.] Meier, M., A. Mekis, A. Dodabalapur, A. Timko, R. E. Slusher, J. D. Joannopoulos, and O. Nalamasu, "Laser action from two-dimensional distributed feedback in photonic crystals," Appl. Phys. Lett., Vol. 74, No. 1, 7-9, 1999.

[5.] Noda, S., M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, "Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design," Science, Vol. 293, No. 5532, 1123-1125, 2001.

[6.] Turnbull, G. A., P. Andrew, W. L. Barnes, and I. D. W. Samuel, "Operating characteristics of a semiconducting polymer laser pumped by a microchip laser," Appl. Phys. Lett., Vol. 82, No. 3, 313-315, 2003.

[7.] Vurgaftman, I. and J. R. Meyer, "Design optimization for high-brightness surface-emitting photonic-crystal distributed-feedback lasers," IEEE J. Quantum Electron., Vol. 39, No. 6, 689-700, 2003.

[8.] Ohnishi, D., T. Okano, M. Imada, and S. Noda, "Room temperature continuous wave operation of a surface-emitting two-dimensional photonic crystal diode laser," Opt. Express, Vol. 12, No. 8, 1562-1568, 2004.

[9.] Matsubara, H., S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, "GaN photonic-crystal surface-emitting laser at blue-violet wavelengths," Science, Vol. 319, No. 5862, 445-447, 2008.

[10.] Lu, T. C., S. W. Chen, L. F. Lin, T. T. Kao, C. C. Kao, P. Yu, H. C. Kuo, and S. C. Wang, "GaN-based two-dimensional surface-emitting photonic crystal lasers with AlN/GaN distributed Bragg reflector," Appl. Phys. Lett., Vol. 92, No. 1, 011129 1-3, 2008.

[11.] Kim, M., C. S. Kim, W. W. Bewley, J. R. Lindle, C. L. Canedy, I. Vurgaftman, and J. R. Meyer, "Surface emitting photonic-crystal distributed-feedback laser for the midinfrared," Appl. Phys. Lett., Vol. 88, No. 19, 191105 1-3, 2006.

[12.] Imada, M., A. Chutinan, S. Noda, and M. Mochizuki, "Multidirectionally distributed feedback photonic crystal lasers," Phys. Rev. B, Vol. 65, No. 19, 195306 1-8, 2002.

[13.] Yokoyama, M. and S. Noda, "Finite-difference time-domain simulation of two-dimensional photonic crystal surface-emitting laser," Opt. Express, Vol. 13, No. 8, 2869-2880, 2005.

[14.] Plihal, M. and A. A. Maradudin, "Photonic band structure of two-dimensional systems: The triangular lattice," Phys. Rev. B, Vol. 44, No. 16, 8565-8571, 1991.

[15.] Sakai, K., E. Miyai, and S. Noda, "Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse electric-like mode," Appl. Phys. Lett., Vol. 89, No. 2, 021101 1-3, 2006.

[16.] Sakai, K., E. Miyai, and S. Noda, "Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization," IEEE J. Quantum Electron., Vol. 46, No. 5, 788-795, 2010.

[17.] Sakai, K., E. Miyai, and S. Noda, "Two-dimensional coupled wave theory for square-lattice photonic-crystal lasers with TM polarization," Opt. Express, Vol. 15, 3981-3990, 2007.

[18.] Sakai, K., J. Yue, and S. Noda, "Coupled-wave model for triangular-lattice photonic crystal with transverse electric polarization," Opt. Express, Vol. 16, No. 9, 6033-6040, 2008.

[19.] Koba, M., P. Szczepanski, and T. Kossek, "Nonlinear operation of a 2D triangular lattice photonic crystal laser," IEEE J. Quantum Electron., Vol. 47, No. 1, 13-19, 2011.

[20.] Scamarcio, G., F. Capasso, C. Sirtori, J. Faist, A. L. Hutchinson, D. L. Sivco, and A. Y. Cho, "High-power infrared (8-micrometer wavelength) superlattice lasers," Science, Vol. 276, No. 5313, 773776, 1997.

[21.] Kogelnik, H., "Coupled wave theory for thick hologram gratings," Bell Syst. Tech. J., Vol. 48, 2909-2947, 1969.

[22.] Kazarinov, R. and C. Henry, "Second-order distributed feedback lasers with mode selection provided by first-order radiation losses," IEEE J. Quantum Electron., Vol. 21, No. 2, 144-150, 1985.

[23.] Johnson, S. and J. Joannopoulos, "Block-iterative frequency-domain methods for Maxwell's equations in a planewave basis," Opt. Express, Vol. 8, No. 10, 173-190, 2001.

[24.] Liang, Y., C. Peng, K. Sakai, S. Iwahashi, and S. Noda, "Three-dimensional coupled-wave model for square-lattice photonic crystal lasers with transverse electric polarization: A general approach," Phys. Rev. B, Vol. 84, No. 19, 195119 1-11, 2011.

[25.] Peng, C., Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, "Coupled-wave analysis for photonic-crystal surface-emitting lasers on air holes with arbitrary sidewalls," Opt. Express, Vol. 19, No. 24, 24672-24686, 2011.

M. Koba (1,2,3), *, P. Szczepanski (2,3)

(1) Institute of Experimental Physics, University of Warsaw, 69 Hoza St., Warsaw 00-681, Poland

(2) Institute of Microelectronics and Optoelectronics, Warsaw University of Technology, 75 Koszykowa St., Warsaw 00-662, Poland

(3) National Institute of Telecommunications, 1 Szachowa St., Warsaw 04-894, Poland

* Corresponding author: Marcin Koba (mkoba@elka.pw.edu.pl).

Printer friendly Cite/link Email Feedback | |

Author: | Koba, M.; Szczepanski, P. |
---|---|

Publication: | Progress In Electromagnetics Research |

Article Type: | Report |

Geographic Code: | 4EXPO |

Date: | Mar 1, 2012 |

Words: | 5875 |

Previous Article: | Simple traffic surveillance system based on range-Doppler radar images. |

Next Article: | Design and full-wave analysis of piezoelectric micro-needle antenna sensors for enhanced near-field detection of skin cancer. |

Topics: |