We propose an ultrahigh quality factor (Q) photonic crystal slab cavity created by the local modulation of the air hole radius in the waveguide. In the cavity, photons are confined between two mirror regions with larger air holes and the lifetime of photons is greatly enhanced by introducing tapered regions with linearly increasing air hole radii. Q and mode volume are investigated for the cavities with various lengths and air hole size offsets of the tapered region with linearly increasing air hole radii by three-dimensional finite-difference time-domain method. The behaviors are analyzed by the mode patterns in real space and wavevector space. We obtain a numerical Q up to 8.8×107 for a mode volume of 1.6 (λ/n)3. Concerning the waveguide coupling, the cavity shows 80% coupling efficiency while keeping Q higher than 106.
© 2008 Optical Society of America
Since E. Yablonovitch and S. John firstly proposed the concept of photonic bandgap (PBG) in 1987 [1,2], the photonic crystal (PhC) cavity has attracted particular attention owing to its ability to confine photons within wavelength dimensions by the PBG. Especially two-dimensional (2D) PhC slab cavities have been widely studied by utilizing well established fabrication processes in the semiconductor industry [3–18]. In such 2D PhC slab structures, photons can be confined horizontally by the 2D PBG and vertically by total internal reflection (TIR). However, the quality factor (Q) in the 2D PhC cavity is mainly limited by the leaky components due to insufficient TIR.
An important design rule has been proposed to obtain high Qs in 2D PhC slab cavities [4,5]: The electric field in a cavity should slowly vary like a Gaussian function to eliminate leaky components out of the slab. A photonic heterostructure cavity with different lattice constants along the PhC waveguide was proposed to satisfy this design rule and showed the remarkable high experimental Qs of 2.5×106 in Si [6–8] and 105 in GaAs . Kuramochi et al. proposed a different heterostructure by locally modulating the waveguide width and also demonstrated a high Q of 1.3×106 experimentally [10,11]. Recently, other types of photonic heterostructure cavities have been studied by a local index modulation of a waveguide using fluidic infiltration or photosensitive materials [12–14]. On the other hand, it was reported that PhCs with different air hole radii can confine photons by shifting PBG .
In this paper, we propose a new heterostructure PhC cavity by the local modulations of the radii of the air holes along a line defect waveguide, which completes the series of the heterostructure cavities realizable by lithography. For the proposed heterostructure cavity, we investigate the mode patterns in the real and the wavevector space and the mode behaviors such as Q and mode volume (Vm) by three-dimensional (3D) finite-difference time-domain (FDTD) method. We discuss the effects of the length and the air hole size offset of the tapered regions on Q and Vm.
2. Design of heterostructure cavity using air hole modulation
Our cavity is based on a single line defect waveguide (W1) in a triangular lattice of air holes drilled in a slab surrounded by air. In our cavity, the photon is mainly concentrated between two mirror regions in which the radius of air holes, rwg, is enlarged compared to the background PhC, rpc. In order to obtain more gently varying electric field, the tapered regions with linearly increasing air holes are introduced between the cavity and the mirror regions, as shown as red boxes in Fig. 1(a) and (b). The air hole radii in the tapered region are linearly interpolated between rPC and rwg. Here the air hole size offset is the differences between rPC and rwg. The heterostructure cavity in this paper is described by two structure parameters, the length of tapered regions and rwg of the mirror region. For example of our nomenclature, S5 and S3 cavities have 5 and 3 layers from the cavity center to the mirror region with rwg~0.28a, including 4 and 2 tapered layers, respectively, as shown in Fig. 1. The photons are confined in the cavities by mode-gap along the waveguide, PBG in the photonic crystal, and TIR in the vertical direction of the slab. The refractive index of a slab is 3.4, which corresponds to indexes of Si and GaAs at λ=1.55 µm. We assumed a lattice constant a of 400 nm, a slab thickness of 250 nm, and a hole radius in the PhC region, rPC, of 0.25a, respectively. The size of surrounding PhC is sufficiently large, 48(x axis)×24(y axis) lattice periods. Under these conditions, the total Q is almost equal to the vertical Q in most of the cases.
Dispersion characteristics of the PhC W1 waveguide mode are plotted in Fig. 2(a), which is calculated by two-dimensional (2D) plane wave expansion method. The mode in the white area is guided inside the PhC slab by TIR. As rwg is enlarged, the effective waveguide width decreases so that the frequency of W1 bandedge blue-shifts. When a waveguide with rwg~0.25a is connected to a waveguide with rwg~0.28a, the mode near the bandedge of the waveguide with rwg~0.25a is placed within the mode-gap of the waveguide with rwg~0.28a, as shown in Fig. 2(a). Therefore, the waveguide with enlarged rwg can work as a mirror for the waveguide mode with smaller rwg near the bandedge. The magnetic field pattern of the bandedge mode, which is the origin of the heterostructure cavity mode, is illustrated in Fig. 2(d). As expected by the wavevector 0.5(2π/a) of the bandedge, the field oscillates with a period with 2a along the Γ-K.
The typical mode patterns and the far-field patterns of the heterostructure cavity mode are illustrated in Fig. 3, which are obtained for the S5 cavity with the largest vertical Q~8.8×107 in the case of rwg~0.27a. This value is a very high theoretical Q value for PhC nanocavities, which is comparable with the numerical Q obtained by Kuramochi et al. . Here the Vm is 1.62 (λ/n)3 and the normalized frequency is 0.255 (a/λ). Because all heterostructure cavities are based on the idea that the near-bandedge mode of the W1 waveguide with lower frequency is confined by the mode-gap of a W1 waveguide with higher frequency, the mode patterns shown in Fig. 3 have almost the same symmetries and distributions for all investigated cavities in this paper as well as other heterocavities reported by other groups [6–14]. Hence the mode behaviors for the length and the air hole size offset of tapered regions discussed in the paper can be commonly applied for all other heterostructure cavities. Of course the air hole size offset corresponds to the lattice constant, the waveguide width, and the refractive index changes in other cavities.
Figure 3(a) shows the calculated electric field intensity profile. It shows that the mode is strongly confined near the center, even for the cavity in the image where the change of the air hole radii along the waveguide is not recognizable, which prove that the mode-gap due to the modulation of the radii of the air holes works as a mirror. The vertical magnetic field profile (Hz) of Fig. 3(b) is almost same with the one of the W1 bandedge mode shown in Fig. 2(d).
In general, the electric field symmetry of the PhC cavity mode explains how the mode can be collected in the vertical direction of the slab by an objective lens, which is the most common way to investigate the photoluminescence of PhC cavities. Figures 3(d) and (e) show Ex and Ey field distributions. In the Ex field distribution, the vertical emission is prohibited because of the odd mirror symmetries for both x=0 and y=0. However, in the Ey field distribution, since both mirror symmetries are even, the cavity mode can emit y-polarized photons vertically. Actually, the collected electric field within the NA(Numerical Aperture)=0.6 in the far-field simulation shows a strong y-directional polarized emission, as shown in Fig. 3(f). The polarization of the mode can be used to identify the heterostructure cavity mode in an experiment.
The calculated far-field pattern is plotted in Fig. 3(c). Here the far-field radiation pattern, which can be observed sufficiently far (r≫λ) from a cavity, is calculated by Fourier-transforming the in-plane field components obtained by 3D FDTD . The cavity mode emits most of photons with ±75° emission angle from the vertical direction of the slab and the amounts of the emitted photons within ±30° (NA=0.6) is negligible. Because the heterostructure cavity mode is designed to minimize the components within the light cone so that the mode components near at kx=ky=0 in wavevector space ought to be negligible. Therefore, as long as the ultrahigh Q is maintained in the cavity, the PhC heterostructure cavity is not appropriate as a photon source with vertical emission. However, since the mode in the heterostructure cavity is based on the W1 waveguide mode, the mode is expected to be able to couple to the W1 waveguide very well, which will be illustrated later in this paper.
3. Effect of the length and the air hole size offset of the tapered region
For a slowly changing cavity structure, the tapered region can be introduced for any heterostructure cavity. It was reported that the cavity with one or two tapered layers showed a higher Q . However, in our knowledge, the optimal length and the taper offset in the tapered region have not yet been investigated for a high Q cavity. In this section, we investigate the effect of the length and the air hole size offset of the tapered region on Q and Vm, which can be utilize to enhance any heterostructure cavity using tapered regions.
3.1 Effect of the length of the tapered region
The Q and the Vm of S5, S3, S2 and S1 cavities are investigated over a broad range of rwg. Here, S5, S3, S2, and S1 cavities have 4, 2, 1, and 0 tapered layers from the center to the mirror regions. The Vm is defined in terms of the electric field density in the cavity . Cavities with smaller Vm have broader spreading in the distribution of the wavevector space leading to higher losses into the vertical direction. Consequently, such cavities with smaller Vm tend to have a lower Q [17,18]. If we compare the cavities with the same Vm to exclude the effect of Vm among various cavities, we can focus on the effect of the length of tapered regions.
The normalized frequency, Q, and Vm for the cavities are shown in table 1. The Qv (vertical Q) is inversely proportional to the optical losses out of photonic crystal slab, which limits total Q in most of the cases. The Vm of the cavities is 1.45~1.46 (λ/n)3 and the normalized frequency is around 0.255~0.257 (a/λ). The normalized frequency does not show any significant change for all cavities with various rwg because the mode is originated from the bandedge of W1 waveguide with rwg~0.25a so that the frequency of the cavity is always close to it of the waveguide. The Qv of the cavities increases dramatically with the length of the tapered region. For example, the Qv of the S5 cavity is 3.0×107, which is 100 times larger than 2.8×105 of the S1 cavity.
In order to inspect the behavior of the larger Qv for the cavity with longer tapered regions, the mode profiles of the heterostructure cavities are investigated for the different lengths of the tapered regions. As the amount of the mode within the light cone determines the vertical loss in PhC slabs, it is easier to analyze the behaviors of Q through the mode distribution in the wavevector space. The mode distributions of the cavities S1 and S5 in the wavevector space are illustrated in Fig. 4(a) and (b), which are obtained through 2D Fourier transforms of Ex and Ey fields calculated by the 3D FDTD method. The modes of the investigated cavities are localized with some spreading at the four points, indicated by red color in Fig. 4(a) and (b). The fact that kx of the mode is near ±0.5(2π/a) means that the mode is very close to the bandedge of W1 waveguide mode. The spreading of ky is broader than that of kx because the mode is confined in a shorter distance along the y-direction compared with the confinement along the x-direction, as shown in Fig. 3(a).
On the other hand, in contrast to the distribution of the S5 cavity in Fig. 4(b), fringes over the wavevector space are observed in the distribution of the S1 cavity in Fig. 4(a). These fringes result in non-negligible components within the light cone (indicated by a yellow circle) in the mode distribution of S1 cavity. The components inside the light cone correspond to the vertical loss which cannot satisfy TIR condition in PhC slabs. Therefore the S1 cavity with the shortest tapered regions has a smaller Qv than the S5 cavity with the longest tapered regions. Since the leaky components within the light cone are generated along the x-axis in real space , we investigated mode distribution of S1, S2, S3, and S5 cavities along kx direction in Fig. 4(c).
S2, S3, and S5 cavities have similar spreading near kx=±0.5(2π/a) but not S1 since the mode volumes are almost same. The spreading of the S1 cavity is not clearly observable due to large fringes. Although the magnitudes of the fringes are different for each cavity, there always exist some fringes between kx=±0.5(2π/a) points. As the length of tapered regions increases (S1→S5), the fringes inside the light cone (indicated a yellow box) become smaller which results in a higher Qv. The fringes corresponding to higher order wavevector components are originated from the abrupt change of the electric field envelop . Therefore, the smoothly changing cavity structure building by the longer tapered regions can make a more gently varying electric field envelop, which contributes to a higher Qv.
In Fig. 4(d), Qv increases exponentially as the length of the tapered region, but increases only 1.5 times from the S3 to S5 cavities. According to this graph, it is expected that the Qv would be saturated as the tapered region becomes longer than that of S5 cavity. Because the change of the air hole radii in the tapered region of the S5 cavity is gentle enough to make slowly varying electric field envelop for high Q cavity. On the other hand, in the real fabrication, less modification of air holes makes it simpler to fabricate a cavity. Therefore S3 cavity which has comparable Q with S5 cavity and less modulation of radii can be one of good candidates for an ultrahigh Q cavity.
3.2 Effect of the air hole size offset of the tapered region
In order to investigate the effect of the air hole size offset in tapered regions, we have simulated Q, Vm, and normalized frequency of the S3 cavity for a range of rwg from 0.26a to 0.36a. As the air holes in the tapered region increases linearly from rpc to rwg and rpc is fixed to 0.25a, larger rwg gives a larger air hole size offset in the tapered region. The total Q factor (Qt) decreases exponentially from 2.0×107 to 1.0×105 with larger rwg after a maximum at 0.27a, as shown in Fig. 5(a). Therefore, smaller air hole size offset in the tapered regions can give a higher Qv. It will be discussed in Fig. 6 why Qv always decreases with larger rwg. The total Q and Vm of the S5 cavity show similar behaviors with the S3 cavity, as shown in Fig. 5(d). But Q and Vm of the S5 cavity are always larger than S3 cavity due to more gently changing air hole radii in the tapered region.
On the other hand, Vm decreases linearly with larger rwg from 2.0 to 0.9 (λ/n)3. It should be noted that Qt is always higher than 106 for a broad range of rwg from 0.27a to 0.32a while Vm is kept on the order of a cubic wavelength. It means the cavity has a large tolerance with respect to fabrication errors causing overall air hole size changes. Qv, Qx, and Qy are estimated from the energy flow out of the slab surface, along the waveguide, and through the PhC region, respectively, and are shown in Fig. 5(b). In most of the cases, Qv limits Qt over various rwg except for the cavity with rwg~0.26a where Qx limits Qt.
At the cavity with rwg~0.26a, Qx drops to 50,000. Even though we put many PhC layers along the x-direction, Qx was still small enough to limit Qt because the bandedge frequency of the mirror region with rwg~0.26a is very close to the frequency of the cavity mode resulting in a small mode-gap and weak photon confinement. In contrast to that the frequency of the bandedge mode of the mirror region increases rapidly with rwg, the frequency of the heterostructure cavity mode increases slowly with rwg since only a part of the mode experiences the mirror region with larger rwg, as shown in Fig. 5(c). Thus the mode-gap corresponding to the region below the bandedge mode frequency of the mirror waveguide with smaller rwg becomes smaller so that photon confinement is getting weaker. Therefore, a cavity consisting of two kinds of waveguides with too small difference of refractive index, lattice constant, waveguide width, or air hole radius is not suited for a high Q. Since the fabrication error can broaden the dispersions of PhCs, the small mode-gap would disappear in that case. On the other hand, Vm becomes smaller in the cavity with larger rwg due to stronger confinement by the larger mode-gap.
In Fig. 6, we compare the mode patterns of the S3 cavities with rwg~0.27a and rwg~0.36a in the real space and the wavevector space to analyze the reason for the higher Qv in the cavity with smaller rwg. Qv and Vm are 2.0×107 and 1.46 (λ/n)3 for the cavity with rwg~0.27a, 1.0×105 and 0.93 (λ/n)3 for the cavity with rwg~0.36a. The Qv of the cavity rwg~0.27a is 200 times larger than one in the case of rwg~0.36a. The real space mode pattern for the cavity with smaller rwg in Fig. 6(a) spreads more along the waveguide direction than one for the cavity with larger rwg in Fig. 6(b), which gives larger mode volume. As mentioned previously, the mirror region with larger rwg has a larger mode-gap which results in stronger confinement in the real space.
The mode patterns in the wavevector space show two differences in Fig. 6(c) and (d), which becomes even clearer in the mode distributions along kx direction in Fig. 6(e). First, wavevector spreading along the kx at the kx=±0.5(2π/a) for the cavity with rwg~0.27a is apparently smaller than it is for the cavity with rwg~0.36a. This smaller spreading due to larger Vm induces less leaky components within light cone, which contributes partly to larger Qv in the cavity with smaller rwg. Second, there are relatively larger fringes inside light cone in the cavity mode with rwg~0.36a in Fig. 6(e). It is explained by that larger air hole size offset in the cavity with larger rwg makes rapid change of the cavity structure so that higher order wavevector components appear within light cone. It also partly results in larger Qv in the cavity with smaller rwg.
4. Light extraction using waveguides
As shown in the far-field of Fig. 3(c), since a heterostructure cavity emits photons with large emission angles, only small portion of the emitted photons in the cavity can be extracted to the vertical direction of the slab. On the other hand, the heterostructure cavity would be expected to be a good candidate as an in-plane light source since the cavity is realized by modifying a waveguide. Qt and extraction efficiency, η, are calculated in the S3 cavity with rwg~0.28a as the length of mirror layer connected to output waveguide increases. The extraction efficiency is defined as the fraction of the energy coupled into the output waveguide of total emitted energy. The simulated structure is a S3 cavity with finite number of mirror layers connected with the output waveguide regions with rwg~0.25a at the both sides.
Qy and Qv are kept to be almost constant, however, Qx, which is inversely proportional to the energy portion coupling to the output waveguides, decreases with a smaller number of mirror layers, as shown in Fig. 7(a). This means that the output waveguide works as a pathway of photons which does not induce any vertical scattering or loss through PhC region. This coupling behavior is due to the fact that the mode of the output waveguide with rwg~0.25a basically has the same mode pattern of the heterostructure cavity mode so that the most of photons couple to the output waveguide in the case of a small number of mirror layers, the extraction efficiency becomes larger, as shown in Fig. 7(b). When the S3 cavity is connected to the output waveguide with 5 mirror layers, Qt is over 106 and 80% of the photons emitted from the cavity can couple to the output waveguide.
In conclusion, we demonstrate a new type of heterostructure PhC cavity by the local modulation of air hole radii along a W1 waveguide. Introducing the tapered regions in the heterostructure cavity increases Qv greatly. Owing to the smoothly changing cavity structure, Qv of the heterostructure cavity increases by a factor of 100 as the tapered regions increases from 0 to 4 lattice periods. The cavity with a two lattice periods long tapered regions has an ultrahigh Q of 2.0×107. The smaller air hole size offset of the tapered regions is also able to give a higher Qv due to a gently varying envelop of the electric field. In contrast to the small vertical emission of the heterostructure cavity, the cavity mode can couple into the output waveguide with 80% coupling efficiency while maintaining Qt~106 in the in-plane waveguide coupling. Thus the heterostructure cavity can be a good candidate for an in-plane photon source. The mode behaviors of the cavity with air hole modulation discussed in this paper can give clues to design any heterostructure cavity with same mode distributions based on W1 waveguide and mode-gap.
The authors would like to thank Min-Kyo Seo and Yong-Hee Lee from Nanolaser laboratory, KAIST for supporting their FDTD code and cluster computer system.
References and links
6. B. S. Song, S. Noda, T. Asano, and Y. Akahane, “Ultra-high-Q photonic double-heterostructure nanocavity,” Nature Mater. 4, 207–210 (2005). [CrossRef]
7. S. Noda, M. Fujita, and T. Asano, “Spontaneous-emission control by photonic crystals and nanocavities,” Nature Photon. 1, 449–458 (2007). [CrossRef]
9. R. Herrmann, T. Sünner, T. Hein, A. Löffler, M. Kamp, and A. Forchel, “Ultrahigh-quality photonic crystal cavity in GaAs,” Opt. Lett. 31, 1299–1301(2006). [CrossRef]
10. E. Kuramochi, M. Notomi, S. Mitsugi, A. Shinya, T. Tanabe, and T. Watanabe, “Ultrahigh-Q photonic crystal nanocavities realized by the local width modulation of a line defect,” Appl. Phys. Lett. 88, 041112 1-3 (2006). [CrossRef]
11. T. Tanabe, M. Notomi, E. Kuramochi, A. Shinya, and H. Taniyama, “Trapping and delaying photons for one nanosecond in an ultrasmall high-Q photonic-crystal nanocavity,” Nature Photon. 1, 49–52 (2007). [CrossRef]
12. S. T. Hanic, M. J. Steel, C. M. Sterke, and D. J. Moss, “High-Q cavities in photosensitive photonic crystals,” Opt. Lett. 32, 542–544 (2007). [CrossRef]
13. S. T. Hanic, C. M. Sterke, and M. J. Steel, “Design of high-Q cavities in photonic crystal slab heterostructures by air-holes infiltration,” Opt. Express 14, 12451–12456 (2006). [CrossRef]
14. C. L. C. Smith, D. K. C. Wu, M. W. Lee, C. Monat, S. T. Hanic, C. Grillet, B. J. Eggleton, D. Freeman, Y. Ruan, S. Madden, B. L. Davies, H. Giessen, and Y. H. Lee, “Microfluidic photonic crystal double heterostructures,” Appl. Phys. Lett. 91, 121103-1-3 (2007). [CrossRef]
16. S. H. Kim, S. Kim, and Y. H. Lee, “Vertical beaming of wavelength-scale photonic crystal resonantors,” Phys. Rev. B.73, 235117 1-13 (2006). [CrossRef]
17. H. Y. Ryu, M. Notomi, and Y. H. Lee, “High-quality-factor and small-mode-volume hexapole modes in photonic-crystal slab nanocavities,” Appl. Phys. Lett. 83, 4294–4296 (2003). [CrossRef]