Abstract

Three-dimensional coupled-wave theory is extended to model triangular-lattice photonic-crystal surface-emitting lasers with transverse-electric polarization. A generalized coupled-wave equation is derived to describe the sixfold symmetry of the eigenmodes in a triangular lattice. The extended theory includes the effects of both surface radiation and in-plane losses in a finite-size laser structure. Modal properties of interest including the band structure, radiation constant, threshold gain, field intensity profile, and far-field pattern (FFP) are calculated. The calculated band structure and FFP, as well as the predicted lasing mode, agree well with experimental observations. The effect of air-hole size on mode selection is also studied and confirmed by experiment.

© 2013 OSA

1. Introduction

Two-dimensional (2D) photonic-crystal surface-emitting lasers (PC-SELs) are becoming increasingly important owing to their promising functionality and improved performance compared to conventional semiconductor lasers. A number of successful demonstrations underpinning this promise have already been conducted [113]. By utilizing the band edge of the photonic band structure, single longitudinal and transverse mode oscillation in two dimensions has been achieved with a large lasing area, which enables high-power, single-mode operation [5,6]. The output beam of such devices is generally emitted in the direction normal to the 2D PC plane and has a small beam divergence angle of less than 1° owing to the large area of coherent oscillation [5, 6]. Furthermore, both the polarization and beam pattern [2, 6, 12, 13], as well as the beam direction [11], of the laser output can be controlled by appropriate design of the PC geometry. In addition, the lasing wavelengths of PC-SELs have spanned a wide range from near-infrared to blue-violet, mid-infrared, and terahertz regimes [610].

To further improve the device performance or exploit new functionalities of PC-SELs, it is important to develop an accurate first-principles model to describe the lasing properties in these structures. Realistic devices of PC-SELs require three-dimensional (3D) analysis, because both the surface emission in the vertical direction and the 2D optical feedback in the PC plane need to be taken into account. In addition, due to two important features of the PC laser cavities – a large area that is typically above 200 periods for one side length and increasingly complicated air-hole shape designs [2, 11, 21, 22] – it is difficult to model this type of laser using brute-force computer simulation methods such as the finite-difference time-domain (FDTD) [14, 15] or finite-element [10] methods.

One of the principal analytical tools that has been widely used to study laser devices incorporating periodic structures is the one-dimensional (1D) coupled-wave theory (CWT) originally proposed by Kogelnik and Shank in 1972 [16]. Recently, with the development of PC-SELs, this original theory has been extended, albeit specifically for a 2D model [4, 1720]. More recently, we proposed a 3D CWT model that affords an exact analytical treatment of the full 3D structure and achieves a very accurate and efficient analysis of the band-edge modes of square-lattice PC-SEL devices [2123]. Moreover, this 3D CWT framework has been extended in a fairly general fashion [24]. Now, it is not only applicable to the general centered-rectangular lattice PC structures [13] but also capable of treating the non-Γ-point modes in the vicinity of the second-order Γ point. However, all of our previous works on the 3D CWT have thus far been restricted to lattice structures for which the eigenstate of light can be described by four dominant wavevectors (i.e., basic waves). As a result, they are intrinsically applicable to lattice structure with C2v or C4v symmetry [25] and are unable to treat PCs with higher symmetry, e.g., triangular-lattice PCs with C6v symmetry [19,25]. In this study, we further extend the 3D CWT works [2124] to develop a coupled-wave model for triangular-lattice PC-SELs. We investigate various modal properties of interest, including the band structure, radiation constant, threshold gain, field intensity envelope profile, and far-field pattern (FFP). The calculated results of the band structure and FFP of the predicted lasing mode are compared with experimental data to validate our theory. Finally, the effect of the air-hole size on the mode selection of the laser device is studied and discussed.

The remainder of this paper is organized as follows. Section 2 describes the PC geometry and derivations of the coupled-wave equations. Section 3 presents numerical results, their comparison with experimental results, and discussions. Section 4 concludes this work by summarizing our findings.

2. Formulation

2.1. Device geometry and photonic-crystal design

A schematic of the PC-SEL device investigated in this study is shown in Fig. 1(a). The PC layer is embedded inside a multilayer structure that is assumed to support only a single fundamental guided mode. Details of the device are described in [26], and the structural parameters are listed in Table 1. The average refractive index of the PC layer is given by nav=fna2+(1f)nb2, where na is the refractive index of air, nb is the refractive index of the background dielectric material (GaAs), and f is the air-hole filling factor (i.e., the fraction of the area of a unit cell occupied by air holes). For circular-shape air holes, f=(2π/3)r2/a2, where r is the hole radius and a is the lattice constant. A typical band structure calculated by the 2D plane-wave expansion method [3] is shown in Fig. 1(b). We are particularly interested in the band-edge modes in the vicinity of the second-order Γ point (red circle), where the group velocity of light becomes zero. At these band edges, the PC lattice serves both to provide the 2D distributed feedback and to couple the in-plane guided waves into the surface normal to construct the output.

 

Fig. 1 (a) Schematic structure of a photonic-crystal surface-emitting laser device with a triangular lattice. (b) Band structure of triangular-lattice PCs calculated by the 2D plane-wave expansion method for transverse-electric (TE) mode [3]. The red circle indicates the second-order Γ point. The inset shows the high-symmetry points at the corners of the irreducible Brillouin zone (shaded light blue).

Download Full Size | PPT Slide | PDF

Tables Icon

Table 1. Structural parameters of the PC-SEL device.

Figure 2(a) shows a schematic of a triangular-lattice PC. The primitive translation vectors of the triangular lattice, a1 and a2, can be expressed as

a1=(3a/2,a/2),a2=(3a/2,a/2).
Then, the corresponding reciprocal lattice vector shown in Fig. 2(b) is given by
Gmn=mb1+nb2=m+n2β0x^+3(nm)2β0y^m0x(m,n)β0x^+n0y(m,n)β0y^,
where the primitive reciprocal lattice vectors are b1=(β0/2,3β0/2) and b2=(β0/2,3β0/2), β0=4π/3a, and m and n are integers.

 

Fig. 2 (a) Schematic of a triangular-lattice PC in real space. The blue arrows denote the primitive translation vectors a1 and a2, and a is the lattice constant. (b) Reciprocal lattice space of a triangular-lattice PC. The colored arrows indicate the six basic waves: R1, S1, R2, S2, R3, and S3 at the second-order Γ point, whose wavenumber is equal to β0=4π/3a.

Download Full Size | PPT Slide | PDF

Light propagating inside a PC must obey Bloch’s theorem, and the Bloch wave state u(r) can be expressed as

u(r)=Gmnamn(z)ei(Δk+Gmn)r,
where amn represents the field amplitude of a wave of the order of (m, n) and Δk = Δxβ0 + Δyβ0ŷ represents the wavenumber deviation from the Γ point [24]. Specifically, in the vicinity of the second-order Γ point (Δk ≃ 0), only the field amplitudes of the waves with |Gmn| = β0 are dominant [3, 21, 24]. We refer to these waves as basic waves: R1, S1, R2, S2, R3, and S3, as illustrated in Fig. 2(b).

2.2. Derivation of 3D coupled-wave equations

By following Refs. [21, 23, 24], we seek solutions to

××E(r)=k02n˜2(r)E(r),
where the time dependence of the electric field E(r) is eiωt, k0(= ω/c) is the free-space wavenumber, ω is the angular frequency, c is the velocity of light in free space, and ñ is the refractive index (a complex number) satisfying k02n˜2(r)k02n2(r)+2ik0n0(z)α˜(z)[23]; n0(z) represents the average refractive index of the material at position z, and α̃(z) represents the gain (α̃ > 0) or loss (α̃ < 0) in each region. For TE polarization, the electric field can be assumed as (Ex, Ey, 0) and expanded as
Ej(z)=m,nEj,m,n(z)eimx,mnβ0xiny,mnβ0y,j=x,y,
where mx,mnm0x(m, n) + Δx and ny,mn = n0y(m, n) + Δy. Similarly, the periodic function of n2(r) can be expanded as
n2(r)=n02(z)+m0,n0ξm,n(z)eim0xβ0xin0yβ0y,
where ξm,n(z) = 0 outside the PC layer, and air holes within the PC region are assumed to have perfectly vertical sidewalls [22]; thus, n02(z) and ξm,n are constant within every region.

By substituting Eqs. (5) and (6) into Eq. (4) and collecting all terms that are multiplied by the factor eimxβ0xinyβ0y, we obtain

[2z2+k02n02(z)+2ik0n0(z)α˜(z)ny2β02]Ex,m,n+mxnyβ02Ey,m,n2inyβ0Ex,m,ny+iβ0(mxEy,m,ny+nyEy,m,nx)=k02mm,nnξmm,nnEx,m,n,
[2z2+k02n02(z)+2ik0n0(z)α˜(z)mx2β02]Ey,m,n+mxnyβ02Ex,m,n2imxβ0Ey,m,nx+iβ0(mxEx,m,ny+nyEx,m,nx)=k02mm,nnξmm,nnEy,m,n,
z[Ex,m,nx+Ey,m,nyiβ0(mxEx,m,n+nyEy,m,n)]=0.
Here, we have neglected the second-order spatial derivatives [23].

By introducing a proper linear combination of Ex,m,n(z) and Ey,m,n(z), i.e., (nyEx,m,nmxEy,m,n) as described in [21, 24], we obtain

[2z2+k02n02+2ik0n0(z)α˜(z)(mx2+ny2)β02](nyEx,m,nmxEy,m,n)+iβ0[mxnyEy,m,ny(mx2+2ny2)Ex,m,ny]+iβ0[mxnyEx,m,nx+(ny2+2mx2)Ey,m,nx]=k02mm,nnξmm,nn(nyEx,m,nmxEy,m,n).

We are interested in the modes in the vicinity of the second-order Γ point, for which Δx, Δy ≪ 1. We assume that all of the six basic waves have an identical vertical field profile Θ0(z) [24], which satisfies

2Θ0z2+[k02n02β2]Θ0(z)=0,
where the propagation constant β(ββ0) and the vertical field profile Θ0(z) can be calculated by employing the transfer matrix method (TMM) [21], and Θ0 is normalized as Θ0(z)Θ0*(z)dz=1. We further assume that basic waves are transverse waves. Thus, for basic waves, we can express the Ex,m,n and Ey,m,n as
Ex,m,n=ρmnAm,n(x,y)Θ0(z),Ey,m,n=ηmnAm,n(x,y)Θ0(z),
with
ρmn=ny,mnmx,mn2+ny,mn2,ηmn=mx,mnmx,mn2+ny,mn2,
where (m,n) ∈ {(1, 0), (−1, 0), (0, 1), (0, −1), (1, 1), (−1, −1)} which correspond, in order, to {R1, S1, R2, S2, R3, S3} for Am,n(x, y), as shown in Fig. 2(b); ρmn and ηmn denote the perturbation of the basic waves’ polarizations [24]. By substituting Eqs. (12) and (13) in Eq. (10), we obtain
[2z2+k02n02+2ik0n0(z)α˜(z)(mx2+ny2)β02]Am,nΘ0(z)2iβ0(mxx+nyy)Am,nΘ0(z)=k02m,nξmm,nn(ρmnEx,m,n+ηmnEy,m,n).
By substituting Eq. (11) in Eq. (14), multiplying Θ0*(z) on both sides, and integrating over (−∞, ∞), we obtain a coupled-wave equation of the following form:
[δ+iα+(1mx2+ny2)β0]Am,ni(mxx+nyy)Am,n=|r2+s2|=1;rm,snκmn(rs)Ar,sk022β0|m2+n2|1ξmm,nnPC(ρmnEx,m,n+ηmnEy,m,n)Θ0*(z)dz,
where
κmn(rs)=k022β0PCξmr,ns(ρmnρrs+ηmnηrs)|Θ0(z)|2dz
represents direct couplings between basic waves, δ = ββ0 = neffω/cβ0 is the deviation from the Bragg condition, neff is the effective refractive index for the fundamental guided mode, α is the modal loss defined by α=k0β0n0(z)α˜(z)|Θ0(z)|2dz. Note that the integrals in Eqs. (15) and (16) extend only over the PC region because ξm,n = 0 outside that range. The field profiles [i.e., Ex,m′,n′(z) and Ey,m′,n′(z)] of radiative (m′ = n′ = 0) and high-order (|m′2 +n′2| > 1) waves can be expressed in terms of basic waves as described in Appendix A.

The 3D coupled-wave Eq. (15) derived above includes both the in-plane high-order coupling and vertical diffraction effects, neither of which was included or appropriately described in the previous 2D CWT models for triangular-lattice PC-SELs [4, 19, 27]. It is important to note that, Eq. (15) is a very general formula capable of treating any type of centered-rectangular lattice [13] by replacing the basic waves terms. Furthermore, the derivative terms of the basic waves, which are crucial factors for treating the finite-size laser structures but not considered in the previous 3D CWT model for centered-rectangular lattices [24], are now included.

3. Numerical results and discussion

In this section, we present numerical results of the coupled-wave equations derived above. By solving the coupled-wave Eq. (15), we can obtain various properties of interest, including the band structure, radiation constant, threshold gain, field intensity envelope profile within the device, and FFPs. To validate our CWT analysis, we also compare some of the calculated results with experimental observations. Unless otherwise noted in the following examples, the lattice constant a = 341 nm, the air-hole shape is circular with air-hole filling factor f = 0.15 (corresponding to an r/a ratio of 0.20, where r is the air-hole radius), and the structural parameters listed in Table 1 are used.

3.1. Band structure and radiation constant

By considering an infinitely periodic structure, that is, by neglecting the derivative terms in the coupled-wave Eq. (15), the band structure in the vicinity of the second-order Γ point can be obtained. The calculated band structure in the two characteristic directions of the triangular lattice (Γ–X and Γ–J directions) is shown in Fig. 3(a), where the six eigenmodes are referred to as A, B1, B2, C, D1, and D2, respectively, in the order of increasing frequency. The leakiness of these eigenmodes is quantified in terms of the radiation constant αr (αr = 2α: the modal power loss due to the surface emission) [21, 24]. Figure 3(b) shows the calculated radiation constant plotted against the in-plane wavevector. It can be seen from the figure that modes D1 and D2 lead to a substantial radiation loss at the band edge (i.e., at the Γ point), whereas modes A, B1, B2, and C are not leaky (α = 0). This contrast is closely related with the symmetry of their field patterns with respect to the air holes. Figure 3(c) shows the calculated field patterns of the band-edge modes AD2. We note that the electric fields of modes D1 and D2 are symmetric with respect to the air holes, and thus, they lead to significant loss. In contrast, the electric fields of modes A, B1, B2, and C are antisymmetric, leading to zero radiation loss [3, 21]. As a consequence, the antisymmetric (nonradiative) modes are favored for lasing. The lasing mode is determined, as described in the following section, by taking into account the total losses (both the surface emission and in-plane losses) in a finite-size laser structure.

 

Fig. 3 Calculated (a) band structure and (b) radiation constant of the eigenmodes in the vicinity of the second-order Γ point for a triangular-lattice PC-SEL with circular air holes. Six modes near the Γ point are referred to as modes: A, B1, B2, C, D1, and D2, in the order of increasing frequency. Modes B1 and B2, as well as D1 and D2, are doubly degenerate at the Γ point. (c) E-field vector distribution (arrows) and H-field patterns (in color) of the individual band-edge modes. The black circles indicate the air holes. Band-edge modes A and C are known as hexapole and monopole modes, respectively [3]. In the calculations, we use the structural parameters listed in Table 1, the lattice constant a = 341 nm, and the air-hole filling factor f = 0.15.

Download Full Size | PPT Slide | PDF

A noteworthy fact is that as the in-plane wavevector is slightly detuned from the band edge, most of the antisymmetric modes start to become significantly more lossy, regardless of the directions, thereby ensuring that band-edge modes keep lasing stably. However, the only exception is mode C, which has a slow-group-velocity region in the Γ–J directions. (We refer to modes in this region as flat-band modes.) Because the radiation constant becomes noticeably small within this region, the flat-band modes may also participate in lasing mode competition. This type of mode competition was reported in Ref. [28], where a multiple-mode lasing action occurring at both the band edge and the flat-band modes in the Γ–J direction was observed.

3.2. Threshold gain and field intensity envelope

In this section, we solve the coupled-wave Eq. (15) in a finite-size system to evaluate the threshold gain of the eigenmodes. We discretize Eq. (15) by using the staggered-grid finite-difference method [23]. It is important to note here that, for the triangular-lattice PC configuration, it is crucial to both discretize the computational domain on a hexagonal grid [29] instead of the commonly used square grid and use appropriate boundary conditions (Appendix B). We focus our analysis on the antisymmetric band-edge modes mentioned above (i.e., modes A, B1, B2, and C), which are most likely to be lasing modes. In the following calculations, we assume an absorbing boundary condition and the radius of the circular-shape computational region L is set to be L = 30 μm (Appendix B).

For finite-size structures, threshold gain and mode frequency correspond to the imaginary and real parts of the normalized eigenvalue, (δ + )L[23]. As an example, we show in Fig. 4(a) a plot of the normalized threshold gain (αL) as a function of deviation from the Bragg condition (δL). Though a large number of modes are found to exist as a result of the numerical calculation, we identify the fundamental modes (indicated by arrows) using the techniques described in Refs. [19, 23]. We note that there exists an additional mode that has a low threshold comparable to that of the fundamental band-edge modes. We refer to this mode as the W mode. The normalized threshold gains of the low-threshold modes indicated by arrows in Fig. 4(a) are listed in Table 2. As the lasing action occurs at the mode with the lowest threshold gain (loss), Table 2 indicates that mode C [indicated by red arrow in Fig. 4(a)] is favored for lasing with a large threshold-gain discrimination of over 20 cm−1 against the competing side modes (i.e., modes B1 and B2). This fact reveals that stable single-mode operation at mode C can be possibly achieved for triangular-lattice PC-SELs.

 

Fig. 4 (a) Normalized threshold gain (αL) as a function of normalized mode frequency deviation (δL). The fundamental band-edge modes (A, B1, B2, and C) and an additional mode W are indicated by arrows. (b) Field intensity envelopes of the modes indicated by arrows in (a). The data are calculated by using the same parameters as specified in the caption of Fig. 3. Note that the field intensity envelopes are plotted on hexagonal grids with a circular-shape computational domain (dashed circle). The radius of the circular domain, L = 30 μm, is discretized to span seven grid cells for which the eigenvalues converge well (see Appendix B for details).

Download Full Size | PPT Slide | PDF

Tables Icon

Table 2. Normalized threshold gain (αL) of the low-threshold modes indicated by arrows in Fig. 4(a).

To gain insight into the nature of the modes indicated in Fig. 4(a), it is useful to examine their field intensity envelopes throughout the entire laser cavity. The slowly-varying envelope profile (function) modulates the fast-varying Bloch waves in the PC lattice and can be determined by [20, 23]

I(x,y)=k=13|Rk(x,y)|2+|Sk(x,y)|2.
Figure 4(b) shows the field intensity envelopes of the individual modes. We note that energy of the fundamental band-edge modes (A, B1, B2, and C) is relatively confined to the central part of the laser cavity, whereas mode W exhibits a null intensity in the center. The exact physical interpretation of the origin of mode W is not yet clear. From the character of the doughnut-shape intensity envelope, we suggest that mode W might be a whispering-gallery-like mode that was discussed in Ref. [9]. We also note that the threshold-gain discrimination between the individual modes, otherwise equal to zero in infinitely periodic structures, may result from the differences in their in-plane losses [23]. A closer examination of the field intensities [Fig. 4(b)] distributed at the boundary (dashed circle) of the laser cavity reveals that the lowest-threshold mode C possesses the smallest amount of energy that escapes from the boundary (i.e., the in-plane losses). In contrast, mode A exhibits a much higher threshold because of the relatively larger in-plane losses. This intuitive explanation is confirmed by a quantitative study of the in-plane losses based on the formula derived in Ref. [23]. We find that, for the laser structure specified in Fig. 4, the in-plane losses dominate the threshold and that the in-plane losses of mode A are a factor of 2.7 larger than those of mode C.

3.3. Comparison of experimental and theoretical results and discussion

3.3.1. Comparison of band structures and FFPs

Figure 5(a) shows the measured band structure of a fabricated device with the same structural parameters as specified above. The band structure was mapped out around the Γ point by measuring the angle-dependent spontaneous emission spectra well below the lasing threshold [26]. The threshold current (Ith) at room temperature continuous wave (CW) operation was 25 mA. The measurements were performed in the Γ–X and Γ–J directions under a CW current level of 0.9Ith and in the measurements the sample was attached on a heat sink. The dashed curves in Fig. 5(a) replot the data shown in Fig. 3(a) in order to enable a direct comparison with the experimental results. It can be seen that the measured and calculated band structures are in satisfactory agreement. The slight discrepancy may be mainly attributed to the experimental details, such as the change in refractive index owing to carrier injection and thermal effects. We note that in the measured band structure, several additional bands (e.g., bands existing between modes D1 and C near the Γ point) were also observed. These bands are considered to be TM modes as discussed in Ref. [24]. As the laser structure described in Table 1 is not completely symmetric in the vertical direction, TM mode usually coexists with TE mode and may derive some optical gain from the multiple quantum-well active layer. However, since the optical gain is greater for the TE mode than for the TM mode [4], the existence of the TM bands rarely becomes an issue for our laser devices. It is also interesting to note that, especially near the Γ point, the intensity of spontaneous emission from the antisymmetric modes AC is extremely weak compared to that from modes D1 and D2. This is due to the difference in their radiative nature, as presented in Fig. 3(b). The lasing spectrum at the Γ point measured at 1.2Ith current level is shown in Fig. 5(b), which evidently suggests that single-mode lasing indeed occurs at band-edge mode C (yellow dashed line). This observation is consistent with the lasing mode predicted above.

 

Fig. 5 (a) Comparison of the measured (in color) and calculated (white dashed curves) band structures. (b) Lasing spectrum measured above the lasing threshold in the direction normal to the PC plane. The threshold current (Ith) at room temperature CW operation was 25 mA. Band structure and lasing spectrum were measured at CW current levels of 0.9Ith and 1.2Ith, respectively. The frequency of the lasing peak in (b) is 0.3434 (a/λ), indicating that the lasing mode is band-edge mode C (yellow dashed line).

Download Full Size | PPT Slide | PDF

Next, we examine the FFP and polarization profile of the lasing mode. (See Ref. [23] for details of the calculation method for FFP.) Figure 6(a) shows the calculated FFP of mode C. A doughnut-shape beam with a divergence angle of around 1° is obtained; this closely matches the measured FFP shown in Fig. 6(b). The polarization profiles shown in the insets indicate that mode C exhibits an azimuthal polarization, which stems from the monopolar nature of its field pattern, as depicted in Fig. 3(c).

 

Fig. 6 (a) Calculated and (b) measured FFPs of mode C. A scanning microscope image of the fabricated PC with a = 341 nm and f = 0.15 (where the r/a ratio is 0.20 and r is the air-hole radius) is shown in the left inset of (b). Ex (Ey) displayed in the right insets represents the x(y) component of the FFP. Parameters used for the calculations are the same as those shown in the caption of Fig. 4. The yellow arrows in (b) indicate the directions of the measured beam polarization. The beam divergence angle of the FFPs for both cases is around 1°, reflecting the large area of coherent oscillation.

Download Full Size | PPT Slide | PDF

3.3.2. Lasing mode control by tuning air-hole size

As has been demonstrated in the previous works on triangular-lattice PC-SELs [3, 12], single-mode lasing action may also be achieved at band-edge mode A (i.e., the hexapole mode). Based on our CWT analysis, we find that by tuning some structural parameters, e.g., by varying the air-hole size, mode A may be favored for lasing. As an example, we list in Table 3 the normalized threshold gain of the low-threshold modes AC for air holes with a larger filling factor of f = 0.26. The data given in Table 3 suggest that mode A is the lasing mode with a threshold-gain discrimination of around 20 cm−1 against the competing modes (i.e., modes B1 and B2). The calculated FFP and polarization profile of mode A are shown in Fig. 7(a). It is seen that mode A exhibits a doughnut-shape pattern with hexagonal symmetry and possesses a characteristic polarization: The electric-field vector is rotated 4π around the circumference of the doughnut beam (instead of 2π for the monopole mode C). The weak side lobes of the calculated FFP result from the radiation field being terminated abruptly in the numerical calculation.

 

Fig. 7 (a) Calculated and (b) measured FFPs of mode A. A scanning microscope image of the fabricated PC with a = 341 nm and f = 0.26 (where r/a = 0.27) is shown in the left inset of (b). Ex (Ey) displayed in the right insets represents the x (y) component of the FFP. Parameters used for the calculations are specified in the caption of Table 3. The yellow arrows in (b) indicate the directions of the measured beam polarization.

Download Full Size | PPT Slide | PDF

Tables Icon

Table 3. Normalized threshold gain (αL) of the low-threshold modes AC. Parameters used for the calculations are the same as those shown in Fig. 4 except that a larger air-hole filling factor f = 0.26 is used.

To further demonstrate the predictive power of our CWT analysis, we fabricated another device with f = 0.26. By measuring the band structure and the lasing spectrum in the same manner as described for Fig. 5, we confirmed that lasing indeed occurs at band-edge mode A. The measured FFP and polarization profile are shown in Fig. 7(b), and they are in good agreement with the calculated FFP shown in Fig. 7(a). The slight asymmetry in the measured FFP (e.g., the weak streak in the central part of the doughnut-shape beam) can be attributed to the nonuniform electrical pumping or the asymmetric geometry of the gain region determined by carrier diffusion in practical experimental situations.

4. Conclusions

A full 3D CWT was developed for triangular-lattice PC-SELs with TE polarization by extending our previous CWT model [24]. When compared to the previous model, the key feature of the extended theory is that six basic waves were introduced in the coupled-wave equation to describe the C6v symmetry of the eigenmodes in the vicinity of the second-order Γ point; additionally, both surface radiation and in-plane losses, which are critical for analyzing the finite-size laser structure, were included. By solving the extended coupled-wave equation, modal properties, including the band structure, radiation constant, threshold gain, field intensity profile, and FFP were studied. Comparison of the calculated band structure and FFP of the lasing mode with experimental results showed good agreement. Single-mode lasing at the monopole mode C predicted by our threshold gain analysis was confirmed by experimental data. Furthermore, it was shown that the hexapole mode A can be favored for lasing by using a larger air-hole size; this was also confirmed by the experiments.

Although the analyses performed in this work were restricted to circular air holes, the power of our theory lies in the fact that it can be applied to air holes of any arbitrary shape as it inherently incorporates a large number of high-order waves [21]. Moreover, the general coupled-wave equation derived in this work is applicable not only to triangular lattice but also to centered-rectangular lattices with C2v or C4v symmetry, thereby unifying the analysis of the modal properties of PC-SELs. Our future work will focus on a more systematic study of the modal properties of PC-SELs by varying structural parameters of PCs, such as the air-hole shape and the lattice angle, to obtain a comprehensive understanding the physics of PC-SELs and exploit their potential applications.

The presented 3D-CWT analysis is a linear theory, and it gives reliable solutions only at or near the threshold current. When the current level is far above the threshold, modal properties of the practical laser devices may deviate from predictions based on the linear analysis. We will leave the above-threshold 3D-CWT analysis that incorporates the optical nonlinearities such as hole burning, nonuniform gain distribution, and thermal effects to future work.

Appendix A: Solutions of radiative and high-order waves for triangular-lattice PCs

The radiative and high-order waves [i.e., Ex,m′,n′ and Ey,m′,n′ in Eq. (15)] can be solved in a manner similar to that described in our previous works [21,24]. The major difference is that the number of basic waves that serve as sources to excite radiative and high-order waves are six in this case, instead of four. Therefore, we shall only show the final analytical solutions.

Radiative waves

The radiative waves (for which m = n = 0) can be expressed as

(ΔExΔEy)=k02PCG(z,z)Θ0(z)dz(ξ1,0ρ1,0ξ1,0ρ1,0ξ0,1ρ0,1ξ0,1ρ0,1ξ1,1ρ1,1ξ1,1ρ1,1ξ1,0η1,0ξ1,0η1,0ξ0,1η0,1ξ0,1η01ξ1,1η1,1ξ1,1η1,1)V,
where G(z,z′) = −i/2βz · ez|zz′| is an approximated Green’s function [21,22], βz = k0n0(z), and V = (R1, S1, R2, S2, R3, S3)t.

High-order waves

The integrals of high-order waves (for which |m2 + n2| > 1) included in the second term of the right-hand side of the coupled-wave Eq. (15) can be expressed as

(PCEx,m,n(z)Θ0*(z)dzPCEy,m,n(z)Θ0*(z)dz)=1mx2+ny2(nymxmxny)(μm,n(1,0)μm,n(1,0)μm,n(0,1)μm,n(0,1)μm,n(1,1)μm,n(1,1)νm,n(1,0)νm,n(1,0)νm,n(0,1)νm,n(0,1)νm,n(1,1)νm,n(1,1))V(ςx,m,n(1,0)ςx,m,n(1,0)ςx,m,n(0,1)ςx,m,n(0,1)ςx,m,n(1,1)ςx,m,n(1,1)ςy,m,n(1,0)ςy,m,n(1,0)ςy,m,n(0,1)ςy,m,n(0,1)ςy,m,n(1,1)ςy,m,n(1,1))V,
where
μmn(rs)=k02PCξmr,ns(nyρrsmxηrs)Gm,n(z,z)Θ0(z)Θ0*(z)dzdz,
νmn(rs)=PC1n02ξmr,ns(mxρrs+nyηrs)|Θ0(z)|2dz,
Gm,n(z,z)=12βz,m,neβz,m,n|zz|,βz,m,n=(mx2+ny2)β02k02n02(z).

As a consequence, the overlap integrals appearing on the right-hand side of Eq. (15) can be replaced by terms only associated with basic waves. Finally, a coupled-wave equation for infinite periodic structures [for which the derivative terms in Eq. (15) are neglected] can be written in the matrix form as

(δ+iα)V=CV,
where C is a 6 × 6 matrix. The matrix elements of C can be written as
C=dk0,mn+Cb+Cr+Ch,
where
dk0,mn=(dk0;1,0000000dk0;1,0000000dk0;0,1000000dk0;0,1000000dk0;1,1000000dk0;1,1),
Cb=(0κ1,0(1,0)κ1,0(0,1)κ1,0(0,1)κ1,0(1,1)κ1,0(1,1)κ1,0(1,0)0κ1,0(0,1)κ1,0(0,1)κ1,0(1,1)κ1,0(1,1)κ0,1(1,0)κ0,1(1,0)0κ0,1(0,1)κ0,1(1,1)κ0,1(1,1)κ0,1(1,0)κ0,1(1,0)κ0,1(0,1)0κ0,1(1,1)κ0,1(1,1)κ1,1(1,0)κ1,1(1,0)κ1,1(0,1)κ1,1(0,1)0κ1,1(1,1)κ1,1(1,0)κ1,1(1,0)κ1,1(0,1)κ1,1(0,1)κ1,1(1,1)0),
Cr=(ζ1,0(1,0)ζ1,0(1,0)ζ1,0(0,1)ζ1,0(0,1)ζ1,0(1,1)ζ1,0(1,1)ζ1,0(1,0)ζ1,0(1,0)ζ1,0(0,1)ζ1,0(0,1)ζ1,0(1,1)ζ1,0(1,1)ζ0,1(1,0)ζ0,1(1,0)ζ0,1(0,1)ζ0,1(0,1)ζ0,1(1,1)ζ0,1(1,1)ζ0,1(1,0)ζ0,1(1,0)ζ0,1(0,1)ζ0,1(0,1)ζ0,1(1,1)ζ0,1(1,1)ζ1,1(1,0)ζ1,1(1,0)ζ1,1(0,1)ζ1,1(0,1)ζ1,1(1,1)ζ1,1(1,1)ζ1,1(1,0)ζ1,1(1,0)ζ1,1(0,1)ζ1,1(0,1)ζ1,1(1,1)ζ1,1(1,1)),
Ch=(χ1,0(1,0)χ1,0(1,0)χ1,0(0,1)χ1,0(0,1)χ1,0(1,1)χ1,0(1,1)χ1,0(1,0)χ1,0(1,0)χ1,0(0,1)χ1,0(0,1)χ1,0(1,1)χ1,0(1,1)χ0,1(1,0)χ0,1(1,0)χ0,1(0,1)χ0,1(0,1)χ0,1(1,1)χ0,1(1,1)χ0,1(1,0)χ0,1(1,0)χ0,1(0,1)χ0,1(0,1)χ0,1(1,1)χ0,1(1,1)χ1,1(1,0)χ1,1(1,0)χ1,1(0,1)χ1,1(0,1)χ1,1(1,1)χ1,1(1,1)χ1,1(1,0)χ1,1(1,0)χ1,1(0,1)χ1,1(0,1)χ1,1(1,1)χ1,1(1,1)),
and
dk0;m,n=(mx(m,n)2+ny(m,n)21)β0,
ζmn(rs)=k042β0PCξm,nξr,s(ρmnρrs+ηmnηrs)G(z,z)Θ0(z)Θ0*(z)dzdz,
χmn(rs)=k022β0m2+n2>1ξmm,nn(ρmnςx,m,n(rs)+ηmnςy,m,n(rs)).
Here, a large number of high-order waves are included by truncating the summation terms in Eq. (A14) at |m′,n′| ≤ 10 [21]. Note that the matrix dk0,mn represents the variation of the in-plane wavenumbers resulting from deviation of the individual basic wave vectors from the second-order Γ point, and Cb, Cr, and Ch correspond to the coupling effects of basic, radiative, and high-order waves, respectively.

Appendix B: Staggered-grid finite-difference method on hexagonal grids

The coupled-wave equation for finite systems can be solved by using a staggered-grid finite-difference method. We rewrite the coupled-wave Eq. (15) for the band-edge modes (for which Δx = Δy = 0) in the following form:

(δ+iα)[R1S1R2S2R3S3]=[C][R1S1R2S2R3S3]+i[12R1x32R1y12S1x+32S1y12R2x+32R2y12S2x32S2yR3xS3x],
where C is a 6 × 6 matrix composed of coupling coefficients defined by Eqs. (A7)(A11).

To model a triangular-lattice PC having a sixfold symmetry, we discretize the computational domain on a hexagonal grid (see, e.g., Sec. 3.7.2 of Ref. [29]). The location of the six basic-wave components and the coupling coefficients in the vicinity of the (j,k)th hexagonal cell of the grid is illustrated in Fig. 8(a). Then, Eq. (B1) becomes

12(δ+iα)[R1j+1,k+R1j,kS1j+1,k+S1j,kR2j,k+1+R2j,kS2j,k+1+S2j,kR3j+1,k+1+R3j,kS3j+1,k+1+S3j,k]=12[C][R1j+1,k+R1j,kS1j+1,k+S1j,kR2j,k+1+R2j,kS2j,k+1+S2j,kR3j+1,k+1+R3j,kS3j+1,k+1+S3j,k]+i1h[R1j+1,kR1j,kS1j+1,k+S1j,kR2j,k+1R2j,kS2j,k+1+S2j,kR3j+1,k+1R3j,kS3j+1,k+1+S3i,k],
where h is the distance between two adjacent grid points and the coupling coefficients Cmn (1 ≤ m,n ≤ 6) are assumed to be constant throughout the laser cavity. Figure 8(b) shows a schematic of the circular-shape computational domain (yellow shaded region) considered in this work. The computational domain physically corresponds to the gain region determined by the electrode geometry. Because the shape of the electrode at the bottom of the fabricated device is circular (see Ref. [13] for details), we consider a circular computational domain. We define the radius of the circular-shape gain region as L = Nh (N: integer). Schematic depicted in Fig. 8(b) corresponds to the case when N = 2. By multiplying 2L on both sides of Eq. (B2), we obtain
L(δ+iα)[R1j+1,k+R1j,kS1j+1,k+S1j,kR2j,k+1+R2j,kS2j,k+1+S2j,kR3j+1,k+1+R3j,kS3j+1,k+1+S3j,k]=L[C][R1j+1,k+R1j,kS1j+1,k+S1j,kR2j,k+1+R2j,kS2j,k+1+S2j,kR3j+1,k+1+R3j,kS3j+1,k+1+S3j,k]+i2N[R1j+1,kR1j,kS1j+1,k+S1j,kR2j,k+1R2j,kS2j,k+1+S2j,kR3j+1,k+1R3j,kS3j+1,k+1+S3j,k],
which can be treated as a generalized eigenvalue problem with spatially dependent eigenvectors (R1j,k,S1j,k,R2j,k,S2j,k,R3j,k,S3j,k,R1j+1,k,S1j+1,k,R2j+1,k,S2j+1,k,R3j+1,k,S3j+1,k,)t and normalized eigenvalues (δ + )L. In this work, we assume the following absorbing boundary condition: The field amplitude of the basic waves always starts from zero at the boundaries. This generalizes a similar concept for the 1D distributed feedback laser structure originally proposed by Kogelnik and Shank [16]. The boundary consists of the set of points lying immediately outside but closest to one half of the circular domain on the side opposite to the wave’s propagation. As an example, the boundary of R3 is indicated by black squares in Fig. 8(b). We investigated the impact of discretization on the solutions of Eq. (B3) by varying N and confirmed that well-converged solutions can be obtained at N = 7 with a short calculation time (∼1 minute with a personal computer).

 

Fig. 8 (a) Location of field components (basic waves) and coupling coefficients in the vicinity of the (j,k)th hexagonal cell of the grid. Positions of the unknown field components (colored hollow dots) are staggered from the positions of the known coupling coefficients (solid dots). The colored hollow dots correspond to the points that are updated using the finite-difference scheme, while the solid points are points that are not solved. (b) Schematic of a circular computational domain (yellow shaded region) discretized on the hexagonal grids with L = 2h (L is the radius of the circular shape and h is the distance between two adjacent grid points). The green hollow dots inside the black squares indicate the boundary of R3; these are set to be zero in the calculations.

Download Full Size | PPT Slide | PDF

Acknowledgments

The authors are grateful to Dr. T. Asano, Dr. A. Oskooi, Dr. Y. Kurosaka, T. Nakamura, and T. Okino for their helpful discussions and valuable suggestions. This work was partly supported by the Core Research for Evolution Science and Technology (CREST) and Consortium for Photon Science and Technology (C-PhoST) programs of the Japan Science and Technology Agency. One of the authors (Y. Liang) is supported by research fellowships from the Japan Society for the Promotion of Science (JSPS).

References and links

1. M. Imada, 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. 75, 316–318 (1999). [CrossRef]  

2. S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science 293, 1123–1125 (2001). [CrossRef]   [PubMed]  

3. M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B 65, 195306 (2002). [CrossRef]  

4. I. Vurgaftman and J. R. Meyer, “Design optimization for high-brightness surface-emitting photonic-crystal distributed-feedback lasers,” IEEE J. Quantum Electron. 39, 689–700 (2003). [CrossRef]  

5. D. Ohnishi, T. Okano, M. Imada, and S. Noda, “Room temperature continuous wave operation of a surface-emitting two-dimensional photonic crystal diode laser,” Opt. Express 12, 1562–1568 (2004). [CrossRef]   [PubMed]  

6. E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London) 441, 946 (2006). [CrossRef]  

7. H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science 319, 445–447 (2008). [CrossRef]  

8. M. Kim, 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. 88, 191105 (2006). [CrossRef]  

9. Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London) 457, 174–178 (2009). [CrossRef]  

10. L. Mahler and A. Tredicucci, “Photonic engineering of surface-emitting terahertz quantum cascade lasers,” Laser Photonics Rev. 5, 647–658 (2011).

11. Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics 4, 447–450 (2010). [CrossRef]  

12. S. Iwahashi, Y. Kurosaka, K. Sakai, K. Kitamura, N. Takayama, and S. Noda, “Higher-order vector beams produced by photonic-crystal lasers,” Opt. Express 19, 11963–11968 (2011). [CrossRef]   [PubMed]  

13. S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B 85, 035304 (2012). [CrossRef]  

14. H. Y. Ryu, M. Notomi, and Y. H. Lee, “Finite-difference time-domain investigation of band-edge resonant modes in finite-size two-dimensional photonic crystal slab,” Phys. Rev. B 68, 045209 (2003). [CrossRef]  

15. M. Yokoyama and S. Noda, “Finite-difference time-domain simulation of two-dimensional photonic crystal surface-emitting laser,” Opt. Express 13, 2869–2880 (2005). [CrossRef]   [PubMed]  

16. H. Kogelnik and C. V. Shank, “Coupled-wave theory of distributed feedback lasers,” J. Appl. Phys. 43, 2327–2335 (1972). [CrossRef]  

17. M. Toda, “Proposed cross grating single-mode DFB laser,” IEEE J. Quantum Electron. 28, 1653–1662, (1992). [CrossRef]  

18. K. Sakai, E. Miyai, and S. Noda, “Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse-electric-like mode,” Appl. Phys. Lett. 89, 021101 (2006). [CrossRef]  

19. K. Sakai, J. Yue, and S. Noda, “Coupled-wave model for triangular-lattice photonic crystal with transverse electric polarization,” Opt. Express 16, 6033–6040 (2008). [CrossRef]   [PubMed]  

20. K. Sakai, E. Miyai, and S. Noda, “Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization,” IEEE J. Quantum Electron. 46, 788–795 (2010). [CrossRef]  

21. Y. Liang, 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 84, 195119 (2011). [CrossRef]  

22. C. Peng, 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 19, 24672–24686 (2011). [CrossRef]   [PubMed]  

23. Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express 20, 15945–15961 (2012). [CrossRef]   [PubMed]  

24. C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B 86, 035108 (2012). [CrossRef]  

25. K. Sakoda, “Symmetry, degeneracy, and uncoupled modes in two-dimensional photonic lattices,” Phys. Rev. B 52, 7982 (1995). [CrossRef]  

26. K. Sakai, 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. 23, 1335–1340 (2005). [CrossRef]  

27. M. Koba and P. Szczepanski, “The threshold mode structure analysis of the two-dimensional photonic crystal lasers,” Prog. Electromagn. Res. 125, 365–389 (2012). [CrossRef]  

28. K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B 82, 539–541 (2006). [CrossRef]  

29. A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, Norwood, 2005).

References

  • View by:
  • |
  • |
  • |

  1. M. Imada, 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.75, 316–318 (1999).
    [CrossRef]
  2. S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
    [CrossRef] [PubMed]
  3. M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B65, 195306 (2002).
    [CrossRef]
  4. I. Vurgaftman and J. R. Meyer, “Design optimization for high-brightness surface-emitting photonic-crystal distributed-feedback lasers,” IEEE J. Quantum Electron.39, 689–700 (2003).
    [CrossRef]
  5. D. Ohnishi, T. Okano, M. Imada, and S. Noda, “Room temperature continuous wave operation of a surface-emitting two-dimensional photonic crystal diode laser,” Opt. Express12, 1562–1568 (2004).
    [CrossRef] [PubMed]
  6. E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
    [CrossRef]
  7. H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
    [CrossRef]
  8. M. Kim, 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.88, 191105 (2006).
    [CrossRef]
  9. Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
    [CrossRef]
  10. L. Mahler and A. Tredicucci, “Photonic engineering of surface-emitting terahertz quantum cascade lasers,” Laser Photonics Rev.5, 647–658 (2011).
  11. Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
    [CrossRef]
  12. S. Iwahashi, Y. Kurosaka, K. Sakai, K. Kitamura, N. Takayama, and S. Noda, “Higher-order vector beams produced by photonic-crystal lasers,” Opt. Express19, 11963–11968 (2011).
    [CrossRef] [PubMed]
  13. S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B85, 035304 (2012).
    [CrossRef]
  14. H. Y. Ryu, M. Notomi, and Y. H. Lee, “Finite-difference time-domain investigation of band-edge resonant modes in finite-size two-dimensional photonic crystal slab,” Phys. Rev. B68, 045209 (2003).
    [CrossRef]
  15. M. Yokoyama and S. Noda, “Finite-difference time-domain simulation of two-dimensional photonic crystal surface-emitting laser,” Opt. Express13, 2869–2880 (2005).
    [CrossRef] [PubMed]
  16. H. Kogelnik and C. V. Shank, “Coupled-wave theory of distributed feedback lasers,” J. Appl. Phys.43, 2327–2335 (1972).
    [CrossRef]
  17. M. Toda, “Proposed cross grating single-mode DFB laser,” IEEE J. Quantum Electron.28, 1653–1662, (1992).
    [CrossRef]
  18. K. Sakai, E. Miyai, and S. Noda, “Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse-electric-like mode,” Appl. Phys. Lett.89, 021101 (2006).
    [CrossRef]
  19. K. Sakai, J. Yue, and S. Noda, “Coupled-wave model for triangular-lattice photonic crystal with transverse electric polarization,” Opt. Express16, 6033–6040 (2008).
    [CrossRef] [PubMed]
  20. K. Sakai, E. Miyai, and S. Noda, “Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization,” IEEE J. Quantum Electron.46, 788–795 (2010).
    [CrossRef]
  21. Y. Liang, 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. B84, 195119 (2011).
    [CrossRef]
  22. C. Peng, 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. Express19, 24672–24686 (2011).
    [CrossRef] [PubMed]
  23. Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express20, 15945–15961 (2012).
    [CrossRef] [PubMed]
  24. C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
    [CrossRef]
  25. K. Sakoda, “Symmetry, degeneracy, and uncoupled modes in two-dimensional photonic lattices,” Phys. Rev. B52, 7982 (1995).
    [CrossRef]
  26. K. Sakai, 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.23, 1335–1340 (2005).
    [CrossRef]
  27. M. Koba and P. Szczepanski, “The threshold mode structure analysis of the two-dimensional photonic crystal lasers,” Prog. Electromagn. Res.125, 365–389 (2012).
    [CrossRef]
  28. K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
    [CrossRef]
  29. A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, Norwood, 2005).

2012

C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
[CrossRef]

M. Koba and P. Szczepanski, “The threshold mode structure analysis of the two-dimensional photonic crystal lasers,” Prog. Electromagn. Res.125, 365–389 (2012).
[CrossRef]

S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B85, 035304 (2012).
[CrossRef]

Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express20, 15945–15961 (2012).
[CrossRef] [PubMed]

2011

S. Iwahashi, Y. Kurosaka, K. Sakai, K. Kitamura, N. Takayama, and S. Noda, “Higher-order vector beams produced by photonic-crystal lasers,” Opt. Express19, 11963–11968 (2011).
[CrossRef] [PubMed]

C. Peng, 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. Express19, 24672–24686 (2011).
[CrossRef] [PubMed]

Y. Liang, 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. B84, 195119 (2011).
[CrossRef]

L. Mahler and A. Tredicucci, “Photonic engineering of surface-emitting terahertz quantum cascade lasers,” Laser Photonics Rev.5, 647–658 (2011).

2010

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization,” IEEE J. Quantum Electron.46, 788–795 (2010).
[CrossRef]

2009

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

2008

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

K. Sakai, J. Yue, and S. Noda, “Coupled-wave model for triangular-lattice photonic crystal with transverse electric polarization,” Opt. Express16, 6033–6040 (2008).
[CrossRef] [PubMed]

2006

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

M. Kim, 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.88, 191105 (2006).
[CrossRef]

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse-electric-like mode,” Appl. Phys. Lett.89, 021101 (2006).
[CrossRef]

2005

M. Yokoyama and S. Noda, “Finite-difference time-domain simulation of two-dimensional photonic crystal surface-emitting laser,” Opt. Express13, 2869–2880 (2005).
[CrossRef] [PubMed]

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

2004

2003

H. Y. Ryu, M. Notomi, and Y. H. Lee, “Finite-difference time-domain investigation of band-edge resonant modes in finite-size two-dimensional photonic crystal slab,” Phys. Rev. B68, 045209 (2003).
[CrossRef]

I. Vurgaftman and J. R. Meyer, “Design optimization for high-brightness surface-emitting photonic-crystal distributed-feedback lasers,” IEEE J. Quantum Electron.39, 689–700 (2003).
[CrossRef]

2002

M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B65, 195306 (2002).
[CrossRef]

2001

S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
[CrossRef] [PubMed]

1999

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

1995

K. Sakoda, “Symmetry, degeneracy, and uncoupled modes in two-dimensional photonic lattices,” Phys. Rev. B52, 7982 (1995).
[CrossRef]

1992

M. Toda, “Proposed cross grating single-mode DFB laser,” IEEE J. Quantum Electron.28, 1653–1662, (1992).
[CrossRef]

1972

H. Kogelnik and C. V. Shank, “Coupled-wave theory of distributed feedback lasers,” J. Appl. Phys.43, 2327–2335 (1972).
[CrossRef]

Barbieri, S.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Beere, H. E.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Bewley, W. W.

M. Kim, 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.88, 191105 (2006).
[CrossRef]

Busch, K.

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

Canedy, C. L.

M. Kim, 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.88, 191105 (2006).
[CrossRef]

Chassagneux, Y.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Chutinan, A.

M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B65, 195306 (2002).
[CrossRef]

S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
[CrossRef] [PubMed]

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

Colombelli, R.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Crewett, J.

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

Davies, A. G.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Diem, M.

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

Forberich, K.

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

Gombert, A.

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

Hagness, S. C.

A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, Norwood, 2005).

Imada, M.

D. Ohnishi, T. Okano, M. Imada, and S. Noda, “Room temperature continuous wave operation of a surface-emitting two-dimensional photonic crystal diode laser,” Opt. Express12, 1562–1568 (2004).
[CrossRef] [PubMed]

M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B65, 195306 (2002).
[CrossRef]

S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
[CrossRef] [PubMed]

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

Iwahashi, S.

C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
[CrossRef]

S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B85, 035304 (2012).
[CrossRef]

Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express20, 15945–15961 (2012).
[CrossRef] [PubMed]

Y. Liang, 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. B84, 195119 (2011).
[CrossRef]

S. Iwahashi, Y. Kurosaka, K. Sakai, K. Kitamura, N. Takayama, and S. Noda, “Higher-order vector beams produced by photonic-crystal lasers,” Opt. Express19, 11963–11968 (2011).
[CrossRef] [PubMed]

C. Peng, 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. Express19, 24672–24686 (2011).
[CrossRef] [PubMed]

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

Jianglin, Y.

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

Khanna, S. P.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Kim, C. S.

M. Kim, 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.88, 191105 (2006).
[CrossRef]

Kim, M.

M. Kim, 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.88, 191105 (2006).
[CrossRef]

Kitamura, K.

Koba, M.

M. Koba and P. Szczepanski, “The threshold mode structure analysis of the two-dimensional photonic crystal lasers,” Prog. Electromagn. Res.125, 365–389 (2012).
[CrossRef]

Kogelnik, H.

H. Kogelnik and C. V. Shank, “Coupled-wave theory of distributed feedback lasers,” J. Appl. Phys.43, 2327–2335 (1972).
[CrossRef]

Kunishi, W.

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

Kurosaka, Y.

S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B85, 035304 (2012).
[CrossRef]

S. Iwahashi, Y. Kurosaka, K. Sakai, K. Kitamura, N. Takayama, and S. Noda, “Higher-order vector beams produced by photonic-crystal lasers,” Opt. Express19, 11963–11968 (2011).
[CrossRef] [PubMed]

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

Lee, Y. H.

H. Y. Ryu, M. Notomi, and Y. H. Lee, “Finite-difference time-domain investigation of band-edge resonant modes in finite-size two-dimensional photonic crystal slab,” Phys. Rev. B68, 045209 (2003).
[CrossRef]

Lemmer, U.

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

Liang, Y.

C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
[CrossRef]

Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express20, 15945–15961 (2012).
[CrossRef] [PubMed]

Y. Liang, 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. B84, 195119 (2011).
[CrossRef]

C. Peng, 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. Express19, 24672–24686 (2011).
[CrossRef] [PubMed]

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

Lindle, J. R.

M. Kim, 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.88, 191105 (2006).
[CrossRef]

Linfield, E. H.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Mahler, L.

L. Mahler and A. Tredicucci, “Photonic engineering of surface-emitting terahertz quantum cascade lasers,” Laser Photonics Rev.5, 647–658 (2011).

Maineult, W.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Matsubara, H.

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

Meyer, J. R.

M. Kim, 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.88, 191105 (2006).
[CrossRef]

I. Vurgaftman and J. R. Meyer, “Design optimization for high-brightness surface-emitting photonic-crystal distributed-feedback lasers,” IEEE J. Quantum Electron.39, 689–700 (2003).
[CrossRef]

Miyai, E.

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization,” IEEE J. Quantum Electron.46, 788–795 (2010).
[CrossRef]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse-electric-like mode,” Appl. Phys. Lett.89, 021101 (2006).
[CrossRef]

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

Mochizuki, M.

M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B65, 195306 (2002).
[CrossRef]

S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
[CrossRef] [PubMed]

Murata, M.

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

Noda, S.

C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
[CrossRef]

S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B85, 035304 (2012).
[CrossRef]

Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express20, 15945–15961 (2012).
[CrossRef] [PubMed]

Y. Liang, 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. B84, 195119 (2011).
[CrossRef]

S. Iwahashi, Y. Kurosaka, K. Sakai, K. Kitamura, N. Takayama, and S. Noda, “Higher-order vector beams produced by photonic-crystal lasers,” Opt. Express19, 11963–11968 (2011).
[CrossRef] [PubMed]

C. Peng, 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. Express19, 24672–24686 (2011).
[CrossRef] [PubMed]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization,” IEEE J. Quantum Electron.46, 788–795 (2010).
[CrossRef]

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

K. Sakai, J. Yue, and S. Noda, “Coupled-wave model for triangular-lattice photonic crystal with transverse electric polarization,” Opt. Express16, 6033–6040 (2008).
[CrossRef] [PubMed]

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse-electric-like mode,” Appl. Phys. Lett.89, 021101 (2006).
[CrossRef]

M. Yokoyama and S. Noda, “Finite-difference time-domain simulation of two-dimensional photonic crystal surface-emitting laser,” Opt. Express13, 2869–2880 (2005).
[CrossRef] [PubMed]

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

D. Ohnishi, T. Okano, M. Imada, and S. Noda, “Room temperature continuous wave operation of a surface-emitting two-dimensional photonic crystal diode laser,” Opt. Express12, 1562–1568 (2004).
[CrossRef] [PubMed]

M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B65, 195306 (2002).
[CrossRef]

S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
[CrossRef] [PubMed]

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

Notomi, M.

H. Y. Ryu, M. Notomi, and Y. H. Lee, “Finite-difference time-domain investigation of band-edge resonant modes in finite-size two-dimensional photonic crystal slab,” Phys. Rev. B68, 045209 (2003).
[CrossRef]

Ohnishi, D.

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

D. Ohnishi, T. Okano, M. Imada, and S. Noda, “Room temperature continuous wave operation of a surface-emitting two-dimensional photonic crystal diode laser,” Opt. Express12, 1562–1568 (2004).
[CrossRef] [PubMed]

Okano, T.

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

D. Ohnishi, T. Okano, M. Imada, and S. Noda, “Room temperature continuous wave operation of a surface-emitting two-dimensional photonic crystal diode laser,” Opt. Express12, 1562–1568 (2004).
[CrossRef] [PubMed]

Peng, C.

C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
[CrossRef]

Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express20, 15945–15961 (2012).
[CrossRef] [PubMed]

Y. Liang, 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. B84, 195119 (2011).
[CrossRef]

C. Peng, 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. Express19, 24672–24686 (2011).
[CrossRef] [PubMed]

Ritchie, D. A.

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

Ryu, H. Y.

H. Y. Ryu, M. Notomi, and Y. H. Lee, “Finite-difference time-domain investigation of band-edge resonant modes in finite-size two-dimensional photonic crystal slab,” Phys. Rev. B68, 045209 (2003).
[CrossRef]

Saito, H.

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

Sakaguchi, T.

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

Sakai, K.

C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
[CrossRef]

S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B85, 035304 (2012).
[CrossRef]

Y. Liang, C. Peng, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave analysis for square-lattice photonic-crystal lasers with transverse electric polarization: Finite-size effects,” Opt. Express20, 15945–15961 (2012).
[CrossRef] [PubMed]

Y. Liang, 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. B84, 195119 (2011).
[CrossRef]

S. Iwahashi, Y. Kurosaka, K. Sakai, K. Kitamura, N. Takayama, and S. Noda, “Higher-order vector beams produced by photonic-crystal lasers,” Opt. Express19, 11963–11968 (2011).
[CrossRef] [PubMed]

C. Peng, 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. Express19, 24672–24686 (2011).
[CrossRef] [PubMed]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization,” IEEE J. Quantum Electron.46, 788–795 (2010).
[CrossRef]

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

K. Sakai, J. Yue, and S. Noda, “Coupled-wave model for triangular-lattice photonic crystal with transverse electric polarization,” Opt. Express16, 6033–6040 (2008).
[CrossRef] [PubMed]

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse-electric-like mode,” Appl. Phys. Lett.89, 021101 (2006).
[CrossRef]

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

Sakoda, K.

K. Sakoda, “Symmetry, degeneracy, and uncoupled modes in two-dimensional photonic lattices,” Phys. Rev. B52, 7982 (1995).
[CrossRef]

Sasaki, G.

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

Shank, C. V.

H. Kogelnik and C. V. Shank, “Coupled-wave theory of distributed feedback lasers,” J. Appl. Phys.43, 2327–2335 (1972).
[CrossRef]

Szczepanski, P.

M. Koba and P. Szczepanski, “The threshold mode structure analysis of the two-dimensional photonic crystal lasers,” Prog. Electromagn. Res.125, 365–389 (2012).
[CrossRef]

Taflove, A.

A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, Norwood, 2005).

Takayama, N.

Tanaka, Y.

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

Toda, M.

M. Toda, “Proposed cross grating single-mode DFB laser,” IEEE J. Quantum Electron.28, 1653–1662, (1992).
[CrossRef]

Tokuda, T.

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

Tredicucci, A.

L. Mahler and A. Tredicucci, “Photonic engineering of surface-emitting terahertz quantum cascade lasers,” Laser Photonics Rev.5, 647–658 (2011).

Vurgaftman, I.

M. Kim, 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.88, 191105 (2006).
[CrossRef]

I. Vurgaftman and J. R. Meyer, “Design optimization for high-brightness surface-emitting photonic-crystal distributed-feedback lasers,” IEEE J. Quantum Electron.39, 689–700 (2003).
[CrossRef]

Yokoyama, M.

M. Yokoyama and S. Noda, “Finite-difference time-domain simulation of two-dimensional photonic crystal surface-emitting laser,” Opt. Express13, 2869–2880 (2005).
[CrossRef] [PubMed]

S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
[CrossRef] [PubMed]

Yoshimoto, S.

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

Yue, J.

Appl. Phys. B

K. Forberich, M. Diem, J. Crewett, U. Lemmer, A. Gombert, and K. Busch, “Lasing action in two-dimensional organic photonic crystal lasers with hexagonal symmetry,” Appl. Phys. B82, 539–541 (2006).
[CrossRef]

Appl. Phys. Lett.

M. Imada, 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.75, 316–318 (1999).
[CrossRef]

M. Kim, 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.88, 191105 (2006).
[CrossRef]

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave model for square-lattice two-dimensional photonic crystal with transverse-electric-like mode,” Appl. Phys. Lett.89, 021101 (2006).
[CrossRef]

IEEE J. Quantum Electron.

K. Sakai, E. Miyai, and S. Noda, “Coupled-wave theory for square-lattice photonic crystal lasers with TE polarization,” IEEE J. Quantum Electron.46, 788–795 (2010).
[CrossRef]

I. Vurgaftman and J. R. Meyer, “Design optimization for high-brightness surface-emitting photonic-crystal distributed-feedback lasers,” IEEE J. Quantum Electron.39, 689–700 (2003).
[CrossRef]

M. Toda, “Proposed cross grating single-mode DFB laser,” IEEE J. Quantum Electron.28, 1653–1662, (1992).
[CrossRef]

IEEE J. Sel. Areas Commun.

K. Sakai, 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.23, 1335–1340 (2005).
[CrossRef]

J. Appl. Phys.

H. Kogelnik and C. V. Shank, “Coupled-wave theory of distributed feedback lasers,” J. Appl. Phys.43, 2327–2335 (1972).
[CrossRef]

Laser Photonics Rev.

L. Mahler and A. Tredicucci, “Photonic engineering of surface-emitting terahertz quantum cascade lasers,” Laser Photonics Rev.5, 647–658 (2011).

Nat. Photonics

Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics4, 447–450 (2010).
[CrossRef]

Nature (London)

Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield, and A. G. Davies, “Electrically pumped photonic-crystal terahertz lasers controlled by boundary conditions,” Nature (London)457, 174–178 (2009).
[CrossRef]

E. Miyai, K. Sakai, T. Okano, W. Kunishi, D. Ohnishi, and S. Noda, “Lasers producing tailored beams,” Nature (London)441, 946 (2006).
[CrossRef]

Opt. Express

Phys. Rev. B

S. Iwahashi, K. Sakai, Y. Kurosaka, and S. Noda, “Centered-rectangular lattice photonic-crystal surface-emitting lasers,” Phys. Rev. B85, 035304 (2012).
[CrossRef]

H. Y. Ryu, M. Notomi, and Y. H. Lee, “Finite-difference time-domain investigation of band-edge resonant modes in finite-size two-dimensional photonic crystal slab,” Phys. Rev. B68, 045209 (2003).
[CrossRef]

M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B65, 195306 (2002).
[CrossRef]

Y. Liang, 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. B84, 195119 (2011).
[CrossRef]

C. Peng, Y. Liang, K. Sakai, S. Iwahashi, and S. Noda, “Three-dimensional coupled-wave theory analysis of centered-rectangular lattice photonic crystal with transverse-electric-like mode,” Phys. Rev. B86, 035108 (2012).
[CrossRef]

K. Sakoda, “Symmetry, degeneracy, and uncoupled modes in two-dimensional photonic lattices,” Phys. Rev. B52, 7982 (1995).
[CrossRef]

Prog. Electromagn. Res.

M. Koba and P. Szczepanski, “The threshold mode structure analysis of the two-dimensional photonic crystal lasers,” Prog. Electromagn. Res.125, 365–389 (2012).
[CrossRef]

Science

H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic-crystal surface-emitting laser at blue-violet wavelengths,” Science319, 445–447 (2008).
[CrossRef]

S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science293, 1123–1125 (2001).
[CrossRef] [PubMed]

Other

A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. (Artech House, Norwood, 2005).

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1
Fig. 1

(a) Schematic structure of a photonic-crystal surface-emitting laser device with a triangular lattice. (b) Band structure of triangular-lattice PCs calculated by the 2D plane-wave expansion method for transverse-electric (TE) mode [3]. The red circle indicates the second-order Γ point. The inset shows the high-symmetry points at the corners of the irreducible Brillouin zone (shaded light blue).

Fig. 2
Fig. 2

(a) Schematic of a triangular-lattice PC in real space. The blue arrows denote the primitive translation vectors a1 and a2, and a is the lattice constant. (b) Reciprocal lattice space of a triangular-lattice PC. The colored arrows indicate the six basic waves: R1, S1, R2, S2, R3, and S3 at the second-order Γ point, whose wavenumber is equal to β 0 = 4 π / 3 a.

Fig. 3
Fig. 3

Calculated (a) band structure and (b) radiation constant of the eigenmodes in the vicinity of the second-order Γ point for a triangular-lattice PC-SEL with circular air holes. Six modes near the Γ point are referred to as modes: A, B1, B2, C, D1, and D2, in the order of increasing frequency. Modes B1 and B2, as well as D1 and D2, are doubly degenerate at the Γ point. (c) E-field vector distribution (arrows) and H-field patterns (in color) of the individual band-edge modes. The black circles indicate the air holes. Band-edge modes A and C are known as hexapole and monopole modes, respectively [3]. In the calculations, we use the structural parameters listed in Table 1, the lattice constant a = 341 nm, and the air-hole filling factor f = 0.15.

Fig. 4
Fig. 4

(a) Normalized threshold gain (αL) as a function of normalized mode frequency deviation (δL). The fundamental band-edge modes (A, B1, B2, and C) and an additional mode W are indicated by arrows. (b) Field intensity envelopes of the modes indicated by arrows in (a). The data are calculated by using the same parameters as specified in the caption of Fig. 3. Note that the field intensity envelopes are plotted on hexagonal grids with a circular-shape computational domain (dashed circle). The radius of the circular domain, L = 30 μm, is discretized to span seven grid cells for which the eigenvalues converge well (see Appendix B for details).

Fig. 5
Fig. 5

(a) Comparison of the measured (in color) and calculated (white dashed curves) band structures. (b) Lasing spectrum measured above the lasing threshold in the direction normal to the PC plane. The threshold current (Ith) at room temperature CW operation was 25 mA. Band structure and lasing spectrum were measured at CW current levels of 0.9Ith and 1.2Ith, respectively. The frequency of the lasing peak in (b) is 0.3434 (a/λ), indicating that the lasing mode is band-edge mode C (yellow dashed line).

Fig. 6
Fig. 6

(a) Calculated and (b) measured FFPs of mode C. A scanning microscope image of the fabricated PC with a = 341 nm and f = 0.15 (where the r/a ratio is 0.20 and r is the air-hole radius) is shown in the left inset of (b). Ex (Ey) displayed in the right insets represents the x(y) component of the FFP. Parameters used for the calculations are the same as those shown in the caption of Fig. 4. The yellow arrows in (b) indicate the directions of the measured beam polarization. The beam divergence angle of the FFPs for both cases is around 1°, reflecting the large area of coherent oscillation.

Fig. 7
Fig. 7

(a) Calculated and (b) measured FFPs of mode A. A scanning microscope image of the fabricated PC with a = 341 nm and f = 0.26 (where r/a = 0.27) is shown in the left inset of (b). Ex (Ey) displayed in the right insets represents the x (y) component of the FFP. Parameters used for the calculations are specified in the caption of Table 3. The yellow arrows in (b) indicate the directions of the measured beam polarization.

Fig. 8
Fig. 8

(a) Location of field components (basic waves) and coupling coefficients in the vicinity of the (j,k)th hexagonal cell of the grid. Positions of the unknown field components (colored hollow dots) are staggered from the positions of the known coupling coefficients (solid dots). The colored hollow dots correspond to the points that are updated using the finite-difference scheme, while the solid points are points that are not solved. (b) Schematic of a circular computational domain (yellow shaded region) discretized on the hexagonal grids with L = 2h (L is the radius of the circular shape and h is the distance between two adjacent grid points). The green hollow dots inside the black squares indicate the boundary of R3; these are set to be zero in the calculations.

Tables (3)

Tables Icon

Table 1 Structural parameters of the PC-SEL device.

Tables Icon

Table 2 Normalized threshold gain (αL) of the low-threshold modes indicated by arrows in Fig. 4(a).

Tables Icon

Table 3 Normalized threshold gain (αL) of the low-threshold modes AC. Parameters used for the calculations are the same as those shown in Fig. 4 except that a larger air-hole filling factor f = 0.26 is used.

Equations (34)

Equations on this page are rendered with MathJax. Learn more.

a 1 = ( 3 a / 2 , a / 2 ) , a 2 = ( 3 a / 2 , a / 2 ) .
G m n = m b 1 + n b 2 = m + n 2 β 0 x ^ + 3 ( n m ) 2 β 0 y ^ m 0 x ( m , n ) β 0 x ^ + n 0 y ( m , n ) β 0 y ^ ,
u ( r ) = G m n a m n ( z ) e i ( Δ k + G m n ) r ,
× × E ( r ) = k 0 2 n ˜ 2 ( r ) E ( r ) ,
E j ( z ) = m , n E j , m , n ( z ) e i m x , m n β 0 x i n y , m n β 0 y , j = x , y ,
n 2 ( r ) = n 0 2 ( z ) + m 0 , n 0 ξ m , n ( z ) e i m 0 x β 0 x i n 0 y β 0 y ,
[ 2 z 2 + k 0 2 n 0 2 ( z ) + 2 i k 0 n 0 ( z ) α ˜ ( z ) n y 2 β 0 2 ] E x , m , n + m x n y β 0 2 E y , m , n 2 i n y β 0 E x , m , n y + i β 0 ( m x E y , m , n y + n y E y , m , n x ) = k 0 2 m m , n n ξ m m , n n E x , m , n ,
[ 2 z 2 + k 0 2 n 0 2 ( z ) + 2 i k 0 n 0 ( z ) α ˜ ( z ) m x 2 β 0 2 ] E y , m , n + m x n y β 0 2 E x , m , n 2 i m x β 0 E y , m , n x + i β 0 ( m x E x , m , n y + n y E x , m , n x ) = k 0 2 m m , n n ξ m m , n n E y , m , n ,
z [ E x , m , n x + E y , m , n y i β 0 ( m x E x , m , n + n y E y , m , n ) ] = 0 .
[ 2 z 2 + k 0 2 n 0 2 + 2 i k 0 n 0 ( z ) α ˜ ( z ) ( m x 2 + n y 2 ) β 0 2 ] ( n y E x , m , n m x E y , m , n ) + i β 0 [ m x n y E y , m , n y ( m x 2 + 2 n y 2 ) E x , m , n y ] + i β 0 [ m x n y E x , m , n x + ( n y 2 + 2 m x 2 ) E y , m , n x ] = k 0 2 m m , n n ξ m m , n n ( n y E x , m , n m x E y , m , n ) .
2 Θ 0 z 2 + [ k 0 2 n 0 2 β 2 ] Θ 0 ( z ) = 0 ,
E x , m , n = ρ m n A m , n ( x , y ) Θ 0 ( z ) , E y , m , n = η m n A m , n ( x , y ) Θ 0 ( z ) ,
ρ m n = n y , m n m x , m n 2 + n y , m n 2 , η m n = m x , m n m x , m n 2 + n y , m n 2 ,
[ 2 z 2 + k 0 2 n 0 2 + 2 i k 0 n 0 ( z ) α ˜ ( z ) ( m x 2 + n y 2 ) β 0 2 ] A m , n Θ 0 ( z ) 2 i β 0 ( m x x + n y y ) A m , n Θ 0 ( z ) = k 0 2 m , n ξ m m , n n ( ρ m n E x , m , n + η m n E y , m , n ) .
[ δ + i α + ( 1 m x 2 + n y 2 ) β 0 ] A m , n i ( m x x + n y y ) A m , n = | r 2 + s 2 | = 1 ; r m , s n κ m n ( r s ) A r , s k 0 2 2 β 0 | m 2 + n 2 | 1 ξ m m , n n P C ( ρ m n E x , m , n + η m n E y , m , n ) Θ 0 * ( z ) d z ,
κ m n ( r s ) = k 0 2 2 β 0 P C ξ m r , n s ( ρ m n ρ r s + η m n η r s ) | Θ 0 ( z ) | 2 d z
I ( x , y ) = k = 1 3 | R k ( x , y ) | 2 + | S k ( x , y ) | 2 .
( Δ E x Δ E y ) = k 0 2 P C G ( z , z ) Θ 0 ( z ) d z ( ξ 1 , 0 ρ 1 , 0 ξ 1 , 0 ρ 1 , 0 ξ 0 , 1 ρ 0 , 1 ξ 0 , 1 ρ 0 , 1 ξ 1 , 1 ρ 1 , 1 ξ 1 , 1 ρ 1 , 1 ξ 1 , 0 η 1 , 0 ξ 1 , 0 η 1 , 0 ξ 0 , 1 η 0 , 1 ξ 0 , 1 η 0 1 ξ 1 , 1 η 1 , 1 ξ 1 , 1 η 1 , 1 ) V ,
( P C E x , m , n ( z ) Θ 0 * ( z ) d z P C E y , m , n ( z ) Θ 0 * ( z ) d z ) = 1 m x 2 + n y 2 ( n y m x m x n y ) ( μ m , n ( 1 , 0 ) μ m , n ( 1 , 0 ) μ m , n ( 0 , 1 ) μ m , n ( 0 , 1 ) μ m , n ( 1 , 1 ) μ m , n ( 1 , 1 ) ν m , n ( 1 , 0 ) ν m , n ( 1 , 0 ) ν m , n ( 0 , 1 ) ν m , n ( 0 , 1 ) ν m , n ( 1 , 1 ) ν m , n ( 1 , 1 ) ) V ( ς x , m , n ( 1 , 0 ) ς x , m , n ( 1 , 0 ) ς x , m , n ( 0 , 1 ) ς x , m , n ( 0 , 1 ) ς x , m , n ( 1 , 1 ) ς x , m , n ( 1 , 1 ) ς y , m , n ( 1 , 0 ) ς y , m , n ( 1 , 0 ) ς y , m , n ( 0 , 1 ) ς y , m , n ( 0 , 1 ) ς y , m , n ( 1 , 1 ) ς y , m , n ( 1 , 1 ) ) V ,
μ m n ( r s ) = k 0 2 P C ξ m r , n s ( n y ρ r s m x η r s ) G m , n ( z , z ) Θ 0 ( z ) Θ 0 * ( z ) d z d z ,
ν m n ( r s ) = P C 1 n 0 2 ξ m r , n s ( m x ρ r s + n y η r s ) | Θ 0 ( z ) | 2 d z ,
G m , n ( z , z ) = 1 2 β z , m , n e β z , m , n | z z | , β z , m , n = ( m x 2 + n y 2 ) β 0 2 k 0 2 n 0 2 ( z ) .
( δ + i α ) V = CV ,
C = d k 0 , m n + C b + C r + C h ,
d k 0 , m n = ( d k 0 ; 1 , 0 0 0 0 0 0 0 d k 0 ; 1 , 0 0 0 0 0 0 0 d k 0 ; 0 , 1 0 0 0 0 0 0 d k 0 ; 0 , 1 0 0 0 0 0 0 d k 0 ; 1 , 1 0 0 0 0 0 0 d k 0 ; 1 , 1 ) ,
C b = ( 0 κ 1 , 0 ( 1 , 0 ) κ 1 , 0 ( 0 , 1 ) κ 1 , 0 ( 0 , 1 ) κ 1 , 0 ( 1 , 1 ) κ 1 , 0 ( 1 , 1 ) κ 1 , 0 ( 1 , 0 ) 0 κ 1 , 0 ( 0 , 1 ) κ 1 , 0 ( 0 , 1 ) κ 1 , 0 ( 1 , 1 ) κ 1 , 0 ( 1 , 1 ) κ 0 , 1 ( 1 , 0 ) κ 0 , 1 ( 1 , 0 ) 0 κ 0 , 1 ( 0 , 1 ) κ 0 , 1 ( 1 , 1 ) κ 0 , 1 ( 1 , 1 ) κ 0 , 1 ( 1 , 0 ) κ 0 , 1 ( 1 , 0 ) κ 0 , 1 ( 0 , 1 ) 0 κ 0 , 1 ( 1 , 1 ) κ 0 , 1 ( 1 , 1 ) κ 1 , 1 ( 1 , 0 ) κ 1 , 1 ( 1 , 0 ) κ 1 , 1 ( 0 , 1 ) κ 1 , 1 ( 0 , 1 ) 0 κ 1 , 1 ( 1 , 1 ) κ 1 , 1 ( 1 , 0 ) κ 1 , 1 ( 1 , 0 ) κ 1 , 1 ( 0 , 1 ) κ 1 , 1 ( 0 , 1 ) κ 1 , 1 ( 1 , 1 ) 0 ) ,
C r = ( ζ 1 , 0 ( 1 , 0 ) ζ 1 , 0 ( 1 , 0 ) ζ 1 , 0 ( 0 , 1 ) ζ 1 , 0 ( 0 , 1 ) ζ 1 , 0 ( 1 , 1 ) ζ 1 , 0 ( 1 , 1 ) ζ 1 , 0 ( 1 , 0 ) ζ 1 , 0 ( 1 , 0 ) ζ 1 , 0 ( 0 , 1 ) ζ 1 , 0 ( 0 , 1 ) ζ 1 , 0 ( 1 , 1 ) ζ 1 , 0 ( 1 , 1 ) ζ 0 , 1 ( 1 , 0 ) ζ 0 , 1 ( 1 , 0 ) ζ 0 , 1 ( 0 , 1 ) ζ 0 , 1 ( 0 , 1 ) ζ 0 , 1 ( 1 , 1 ) ζ 0 , 1 ( 1 , 1 ) ζ 0 , 1 ( 1 , 0 ) ζ 0 , 1 ( 1 , 0 ) ζ 0 , 1 ( 0 , 1 ) ζ 0 , 1 ( 0 , 1 ) ζ 0 , 1 ( 1 , 1 ) ζ 0 , 1 ( 1 , 1 ) ζ 1 , 1 ( 1 , 0 ) ζ 1 , 1 ( 1 , 0 ) ζ 1 , 1 ( 0 , 1 ) ζ 1 , 1 ( 0 , 1 ) ζ 1 , 1 ( 1 , 1 ) ζ 1 , 1 ( 1 , 1 ) ζ 1 , 1 ( 1 , 0 ) ζ 1 , 1 ( 1 , 0 ) ζ 1 , 1 ( 0 , 1 ) ζ 1 , 1 ( 0 , 1 ) ζ 1 , 1 ( 1 , 1 ) ζ 1 , 1 ( 1 , 1 ) ) ,
C h = ( χ 1 , 0 ( 1 , 0 ) χ 1 , 0 ( 1 , 0 ) χ 1 , 0 ( 0 , 1 ) χ 1 , 0 ( 0 , 1 ) χ 1 , 0 ( 1 , 1 ) χ 1 , 0 ( 1 , 1 ) χ 1 , 0 ( 1 , 0 ) χ 1 , 0 ( 1 , 0 ) χ 1 , 0 ( 0 , 1 ) χ 1 , 0 ( 0 , 1 ) χ 1 , 0 ( 1 , 1 ) χ 1 , 0 ( 1 , 1 ) χ 0 , 1 ( 1 , 0 ) χ 0 , 1 ( 1 , 0 ) χ 0 , 1 ( 0 , 1 ) χ 0 , 1 ( 0 , 1 ) χ 0 , 1 ( 1 , 1 ) χ 0 , 1 ( 1 , 1 ) χ 0 , 1 ( 1 , 0 ) χ 0 , 1 ( 1 , 0 ) χ 0 , 1 ( 0 , 1 ) χ 0 , 1 ( 0 , 1 ) χ 0 , 1 ( 1 , 1 ) χ 0 , 1 ( 1 , 1 ) χ 1 , 1 ( 1 , 0 ) χ 1 , 1 ( 1 , 0 ) χ 1 , 1 ( 0 , 1 ) χ 1 , 1 ( 0 , 1 ) χ 1 , 1 ( 1 , 1 ) χ 1 , 1 ( 1 , 1 ) χ 1 , 1 ( 1 , 0 ) χ 1 , 1 ( 1 , 0 ) χ 1 , 1 ( 0 , 1 ) χ 1 , 1 ( 0 , 1 ) χ 1 , 1 ( 1 , 1 ) χ 1 , 1 ( 1 , 1 ) ) ,
d k 0 ; m , n = ( m x ( m , n ) 2 + n y ( m , n ) 2 1 ) β 0 ,
ζ m n ( r s ) = k 0 4 2 β 0 P C ξ m , n ξ r , s ( ρ m n ρ r s + η m n η r s ) G ( z , z ) Θ 0 ( z ) Θ 0 * ( z ) d z d z ,
χ m n ( r s ) = k 0 2 2 β 0 m 2 + n 2 > 1 ξ m m , n n ( ρ m n ς x , m , n ( r s ) + η m n ς y , m , n ( r s ) ) .
( δ + i α ) [ R 1 S 1 R 2 S 2 R 3 S 3 ] = [ C ] [ R 1 S 1 R 2 S 2 R 3 S 3 ] + i [ 1 2 R 1 x 3 2 R 1 y 1 2 S 1 x + 3 2 S 1 y 1 2 R 2 x + 3 2 R 2 y 1 2 S 2 x 3 2 S 2 y R 3 x S 3 x ] ,
1 2 ( δ + i α ) [ R 1 j + 1 , k + R 1 j , k S 1 j + 1 , k + S 1 j , k R 2 j , k + 1 + R 2 j , k S 2 j , k + 1 + S 2 j , k R 3 j + 1 , k + 1 + R 3 j , k S 3 j + 1 , k + 1 + S 3 j , k ] = 1 2 [ C ] [ R 1 j + 1 , k + R 1 j , k S 1 j + 1 , k + S 1 j , k R 2 j , k + 1 + R 2 j , k S 2 j , k + 1 + S 2 j , k R 3 j + 1 , k + 1 + R 3 j , k S 3 j + 1 , k + 1 + S 3 j , k ] + i 1 h [ R 1 j + 1 , k R 1 j , k S 1 j + 1 , k + S 1 j , k R 2 j , k + 1 R 2 j , k S 2 j , k + 1 + S 2 j , k R 3 j + 1 , k + 1 R 3 j , k S 3 j + 1 , k + 1 + S 3 i , k ] ,
L ( δ + i α ) [ R 1 j + 1 , k + R 1 j , k S 1 j + 1 , k + S 1 j , k R 2 j , k + 1 + R 2 j , k S 2 j , k + 1 + S 2 j , k R 3 j + 1 , k + 1 + R 3 j , k S 3 j + 1 , k + 1 + S 3 j , k ] = L [ C ] [ R 1 j + 1 , k + R 1 j , k S 1 j + 1 , k + S 1 j , k R 2 j , k + 1 + R 2 j , k S 2 j , k + 1 + S 2 j , k R 3 j + 1 , k + 1 + R 3 j , k S 3 j + 1 , k + 1 + S 3 j , k ] + i 2 N [ R 1 j + 1 , k R 1 j , k S 1 j + 1 , k + S 1 j , k R 2 j , k + 1 R 2 j , k S 2 j , k + 1 + S 2 j , k R 3 j + 1 , k + 1 R 3 j , k S 3 j + 1 , k + 1 + S 3 j , k ] ,

Metrics