Photonic crystal nanobeam cavities are versatile platforms of interest for optical communications, optomechanics, optofluidics, cavity QED, etc. In a previous work [Appl. Phys. Lett. 96, 203102 (2010)], we proposed a deterministic method to achieve ultrahigh Q cavities. This follow-up work provides systematic analysis and verifications of the deterministic design recipe and further extends the discussion to air-mode cavities. We demonstrate designs of dielectric-mode and air-mode cavities with Q > 109, as well as dielectric-mode nanobeam cavities with both ultrahigh-Q (> 107) and ultrahigh on-resonance transmissions (T > 95%).
©2011 Optical Society of America
High quality factor (Q), small mode volume (V)  optical cavities provide powerful means for modifying the interactions between light and matter , and have many exciting applications including quantum information processing , nonlinear optics , optomechanics [5, 6], optical trapping  and optofluidics . Photonic crystal cavities (PhC) [9, 10] have demonstrated numerous advantages over other cavity geometries due to their wavelength-scale mode volumes and over-million Q-factors [11–31]. Although small mode volumes of PhC cavities can be easily achieved by design, ultrahigh Q factors are typically obtained using extensive parameter search and optimization. In a previous work , we proposed a deterministic method to design an ultrahigh Q PhC nanobeam cavity and verified our designs experimentally. The proposed method does not rely on any trial-and-error based parameter search and does not require any hole shifting, re-sizing and overall cavity re-scaling. The key design rules we proposed, that result in ultrahigh Q cavities, are (i) zero cavity length (L = 0), (ii) constant length of each mirror segment (‘period’=a) and (iii) a Gaussian-like field attenuation profile, provided by linear increase in the mirror strength.
In this follow-up work, we provide numerical proof of the proposed principles, and systematically optimize the design recipe to realize a radiation limited cavity and waveguide coupled cavity. Furthermore, we extend the recipe to the design of air-mode cavities, whose optical energies are concentrated in the low-index region of the structure.
Nanobeam cavities have recently emerged as a powerful alternative to the slab-based 2-D PhC cavities [17–22]. Nanobeams can achieve Qs on par with those found in slab-based geometries, but in much smaller footprints, and are the most natural geometries for integration with waveguides [23–31]. Our deterministically designed cavities have similar structures to the mode-gap cavities proposed by Notomi et al. , and later demonstrated experimentally by Kuramochi et al. , as well as our own work . We note that the same design principle discussed here could be directly applied to realize ultra-high Q cavities based on dielectric stacks that are of interest for realization of vertical-cavity surface emitting lasers (VCSELs) and sharp filters. Finally, it is important to emphasize that while our method is based on the framework of Fourier space analysis [35–37], alternative approach, based on phase-matching between different mirror segments, could also be used to guide the design, as well as to explain the origin of deterministic ultra-high Q-factors in our devices [38, 39].
2. Numerical verification of the deterministic design approach
In this section, we use finite-difference time-domain (FDTD) simulations to systematically study the design principles proposed in . Figure 1(a) shows the schematics of the proposed cavity structure . It consists of an array of air-holes in decreasing radii, etched into a ridge waveguide. The hole-to-hole distances (“periodicity”) are constant. The structure is symmetric with respect to the dashed line in Fig. 1(a). In contrast to the majority of other cavity designs, current structure has no additional cavity length inserted between the two mirrors (L=0), that is the hole-to-hole spacing between the two central holes is the same as the rest of the structure (a). This minimizes the cavity loss and the mode volume simultaneously. The cavity loss is composed of the radiation loss into the free space (characterized by Q rad) and the coupling loss to the feeding waveguide (Q wg). Q wg can be increased simply by adding more gratings along the waveguide. Q rad can be increased by minimizing the spatial Fourier harmonics of the cavity mode inside the lightcone, achieved by creating a Gaussian-like attenuation profile . The optical energy is concentrated in the dielectric region in the middle of the cavity (Fig. 1(b)). In order to achieve the Gaussian-like attenuation, we proposed to use a linearly increasing mirror strength along the waveguide , which was achieved by tapering the hole radii.
First, we analyze the ideal tapering profile using plane wave expansion method and verify the results with 3D FDTD simulations. The dielectric profile of the structure in the middle plane of the cavity can be expressed as40] in the beam direction (x̂), 40]
The dispersion relation can be obtained by solving the master equation :32], the frequency can be written as (δ is the detuning from the mid-gap frequency) and the wavevector as k = (1 + iγ)π/a. Substituting this into the master equation, we obtain . The cavity resonance asymptotes to the dielectric band-edge of the center mirror segment: (j represents the j th mirror segment counted from the center), at which point the mirror strength γj =1 = 0. γ increases with j. With ε air = 1 and ε Si = 3.462, we calculate in Fig. 2(a) the γ – j relation for different tapering profiles. It can be seen that quadratic tapering profile results in linearly increasing mirror strengths, needed for Gaussian field attenuation . To verify this, we perform FDTD simulation and obtain the band diagram (Fig. 2(b)) and γ – f relation (Fig. 2(c)). As shown in Fig. 2(d), linearly increasing mirror strength is indeed achieved after quadratic tapering.
Next, with the optimized tapering profile, the cavity is formed by putting two such mirrors back to back, leaving a cavity length L in between (Fig. 1(a)). Figure 2(e) shows the simulated Q-factors for various Ls. Highest Q rad is achieved at zero cavity length (L=0), which supports the prediction in  based on 1D model.
Third, we verify that the cavity mode has a Gaussian-like attenuation profile. Figure 2(g) shows the Hz-field distribution in the plane right above the cavity, obtained from 3D FDTD simulation. As shown in Fig. 2(h), this field distribution can be ideally fitted with Hz = sin(πx/a)exp(−σx 2)exp(−ξy 2), with a = 0.33, σ = 0.14 and ξ = 14. The fitted value a agrees with the “period”, and σ agrees with that extracted value from Fig. 2(d): . Figure 2(i) shows Hz distribution along the dashed line in Figs. 2(g) and 2(h). Therefore, we conclude that zero cavity length, fixed periodicity and a quadratic tapering of the filling fraction results in a Gaussian field profile, which leads to a high-Q cavity .
Finally, as we have pointed out in , current method results in a cavity whose resonance is asymptotically approaching the dielectric band-edge frequency of the central mirror segment (circled in Fig. 2(b)). The deviation from the band-edge frequency can be calculated using perturbation theory [41, 42]:Eq. (6), the frequency offset δλ/λ v.s N can be obtained. Figure 2(f) shows the frequency offset for different total number of mirror pairs (N), calculated from the perturbation theory, as well as using FDTD simulations. It can be seen that the deviation decreases as the number of modulated mirror segments increases, and is below 1% for N > 15.
Therefore, we verify that an ultrahigh-Q, dielectric-mode cavity resonant at a target frequency can be designed using the following algorithm:
- Determine a target frequency. For example in our case we want f target = 200THz. Since the cavity resonant frequency is typically 1% smaller than the dielectric band-edge of the central segment, estimated using the perturbation theory, we shift-up the target frequency by 1%, i.e. f adjusted = 202THz.
- Pick the thickness of the nanobeam - this is often pre-determined by the choice of the wafer. For example, in our case, the thickness of the nanobeam is 220nm, determined by the thickness of the device layer of our silicon-on-insulator (SOI) wafer.
- Choose periodicity according to a = λ 0/2n eff, where n eff is effective mode index of the cavity and can be estimated by numerical modeling of a strip waveguide that nanobeam cavity is based on. However, we found that the absolute value of the periodicity is not crucial in our design, as long as there exists a bandgap. Therefore, we pick n eff = 2.23, which is a median value of possible effective indices in the case of free standing silicon nanobeam (n eff ∈ (1, 3.46)). This results in a = 330nm.
- Set nanobeam width. Large width increases the effective index of the cavity mode, pulls the mode away from the light line, thus reducing the in-plane radiation loss. On the other hand, a large beam width will allow for higher order modes with the same symmetry as the fundamental mode of interest. Using band diagram simulations, we found that the width of 700nm is good trade-off between these two conditions (Fig. 2(b)).
- Set the filling fraction of the first mirror section such that its dielectric band-edge is at the adjusted frequency: 202THz in our case. Band diagram calculations based on unit cells are sufficient for this analysis. We found that an optimal filling fraction in our case is f start = 0.2 (Fig. 2(b)).
- Find the filling fraction that produces the maximum mirror strength for the target frequency. This involves calculating the mirror strength for several filling fractions (Fig. 2(c)), each of which takes one or two minutes on a laptop computer. In our case we found that f end = 0.1.
- Pick the number of mirror segments (N) to construct the Gaussian mirror: we found that N ≥ 15 (on each side) are generally good to achieve high radiation-Qs.
- Create the Gaussian mirror by tapering the filling fractions quadratically from f start (=0.2 in our case) to f end (=0.1) over the period of N segments. From the above analysis, the mirror strengths can be linearized through quadric tapering (Fig. 2(d)).
- Finally, the cavity is formed by putting two Gaussian mirrors back to back, with no additional cavity length in between (L = 0). To achieve a radiation-limited cavity (Q wg >> Q rad), 10 additional mirrors with the maximum mirror strength are placed on both ends of the Gaussian mirror. We will show in the next section, no additional mirrors are needed to achieve a waveguide-coupled cavity (Q rad >> Q wg).
Besides the structure that were proposed in  (Fig. 1), the alternative structure which has the air-hole in the symmetry plane, as shown in Fig. 3(a), also satisfies (i)–(ix). Both structures result in dielectric-mode cavities, since the bandgap of each mirror segments red-shifts away from the center of the cavity, and thus a potential well is created for the dielectric band-edge mode of the central segment. The difference is that the energy maximum in the air-hole centered cavity is no longer located in the middle of the structure, but instead in the dielectric region next to the central hole (Fig. 3(b)). Figure 3(c) shows the Hz field profile in the plane right above the cavity, obtained from FDTD simulation. Figure 3(d) shows the fitted field profile using the same parameters that are used in the original structure shown in Figs. 2(g)–2(i), but with sine function replaced by cosine function. Figure 3(e) shows the Hz distribution along the dashed line in Figs. 3(c) and 3(d).
Armed with the analytical field profile of the cavities: (Fig. 2) and (Fig. 3), we can obtain the radiation losses and far fields of the cavities using the Fourier space analysis . The Fourier transforms can be analytically obtained and . Under σa 2 << 1, both distributions have their Fourier components strongly localized at k = ±π/a, as is verified by FDTD simulations in Fig. 4(a) and 4(b). Since is an odd function, it always has a zero Fourier component at k = 0. Therefore, dielectric-centered cavities (Fig. 1) should have higher Q-factors. However, in high-Q cavity designs, σa 2 << 1 is satisfied and thus both dielectric-centered and air-centered cavities have comparable Q-factors. FDTD simulation shows that the above and cavities have Q tot = 3.8 × 108 and Q tot = 3.5 × 108 respectively. The mode volume of the cavity is 0.67(λ res/n Si)3, smaller than the cavity (V = 0.76(λ res/n Si)3).
The far field radiation patterns (obtained using FDTD simulations) of the two cavities are shown in Figs. 4(c) and 4(d). The powers, in both cases, are radiated at shallow angles (> 70° zenith angle) to the direction of the waveguide. The Hodd cavity has even less radiated power at small zenith angles, consistent with the above analysis. By integrating the zenith and azimuth angle dependent far field emission, we found that 32% and 63% of the power emitted to +ẑ direction can be collected by a NA=0.95 lens, respectively for Hodd cavity and Heven cavity.
3. Ultra-high Q, dielectric-mode photonic crystal nanobeam cavities
3.1. Radiation-Q limited and waveguide-coupled cavities
Since the dielectric-centered cavity has smaller V than the one, we focus our discussion in the case. Using the above design algorithm, we design the Gaussian mirror and put 10 additional mirrors with the maximum mirror strength on both ends of the Gaussian mirror to obtain the radiation-limited cavity (Q wg >> Q tot). We find in Fig. 5 that Q tot increases exponentially and V increases linearly as the total number of mirror pairs in the Gaussian mirror (N) increases. A record ultra-high Q of 5.0 × 109 is achieved while maintaining the small mode volume of 0.9 × (λ res/n Si)3 at N = 30.
Our design strategy has an additional important advantage over other types of photonic crystal cavities [17–31], that is: the cavity naturally couples to the feeding waveguide, as the hole radii decrease away from the center of the cavity. High-Q and high transmissions (T) cavities are possible with the above design steps (i)–(ix), with no additional “coupling sections” needed. We study T and Q total dependence on the total number of mirror pair segments in the Gaussian mirror (N) in Fig. 5(b). Partial Q-factors (Q rad, Q wg) are obtained from FDTD simulations, and T is obtained using . As shown in Fig. 5(b), we achieve a nanobeam cavity with Q = 1.3 × 107, T = 97% at N = 25.
3.2. Higher order modes of the dielectric-mode cavity
The ultra-high Q mode that we deterministically designed is the fundamental mode of the cavity. Meanwhile, higher order cavity modes also exist. The number of higher order modes depends on the width of the photonic band gap and total number of mirror segments in the Gaussian mirror. To reduce the simulation time, we study the higher order modes of a waveguide-coupled cavity, that has a total number of 12 mirror pair segments, possessing a moderate Q-factor. Figure 6(a) shows the transmission spectrum obtained from FDTD simulation, by exciting the input waveguide with a waveguide mode, and monitoring the transmission through the cavity at the output waveguide. The band-edge modes are observed at wavelengths longer than 1.6μm and shorter than 1.3μm. Figures 6(b)–6(d) shows the major field-component (Ey) distribution of the three cavity modes. As expected, the eigenmodes alternate between symmetric and anti-symmetric modes. Symmetry plane is defined perpendicular to the beam direction, in the middle of the cavity (dashed line in Fig. 6). The total Q-factors of modes I–III are 10,210, 1,077 and 286 respectively. Effective mode volumes of them are 0.55, 0.85 and 1.06 respectively. We note that transversely odd modes are well separated from the transversely symmetric cavity modes, hence were not considered in Fig. 6.
4. Ultra-high Q, air-mode photonic crystal nanobeam cavities
4.1. Radiation-Q limited cavity
An air-mode cavity concentrates the optical energy in the low index region of the cavity. Therefore, these cavities are of interest for applications where strong interactions between light and material placed in the low index region of the cavity is required, including nonlinear optics , optical trapping , biochemical sensing  and light-atom interaction . The ultra-high Q air-mode nanobeam cavity is realized by pulling the air-band mode of photonic crystal into its bandgap, which can also be designed using the same design principles that we developed for dielectric-mode cavities. In contrast to the dielectric-mode case, the resonant frequency of the air-mode cavity is determined by the air band-edge frequency of the central mirror segment. Then, to create the Gaussian confinement, the bandgaps of the mirror segments should shift to higher frequencies as their distances from the center of the cavity increase. This can be achieved by progressively increasing the filling fractions of the mirror segments away from the center of the structure (instead of decreasing in the dielectric-mode cavity case). One way to accomplish this is to increase the size of the holes away from the center of the cavity. While this may be suitable for non-waveguide coupled (radiation-Q limited) cavities, it is not ideal for a waveguide-coupled cavity. For this reason, we employ the design that relies on tapering of the waveguide width instead of the hole size. Similar geometry was recently proposed by Ahn et. al.  for the design of a dielectric-mode photonic crystal laser.
The same design steps can be followed as in the dielectric-mode cavity case, with the following changes: First, the adjusted frequency (198THz) is 1% lower than the target frequency (200THz). (The thickness of the nanobeam is 220nm and period is 330nm, same as previous case.) Second, the nanobeam width at the center of the cavity is w start = 1μm (Fig. 7(a)), with the hole radii kept constant at 100nm. Third, to create the Gaussian mirror, the beam widths are quadratically tapered from w start = 1μm to w end = 0.7μm, which produces the maximum mirror strength (band diagrams shown in Fig. 7(a)). This procedure involves calculating the mirror strength for several beam widths (Fig. 7(b)), each takes one or two minutes on a laptop computer. As shown in Fig. 7(c), the mirror strengths are linearized after the quadratic tapering. In order to achieve a radiation-Q limited cavity, 10 additional mirror segments are placed at both ends of the Gaussian mirror that has beam width w end = 0.7μm.
Similar in the dielectric-mode cavity cases, and air-mode cavities can be formed by placing the air and dielectric in the central symmetric plane of the cavity, respectively. Again, we will focus on , air-mode cavities and the conclusions will be valid to the cavities as well. Figure 8(a) shows the total Q of nanobeam cavities that have different total number of mirror pair segments in the Gaussian mirrors. We have achieved a record ultra-high Q of 1.4 × 109, air-mode nanobeam cavity. As shown in Fig. 8(a), the effective mode volumes of the air-mode cavities are much larger than the dielectric-mode cavities.
4.2. Cavity strongly coupled to the feeding waveguide
As we have pointed out, the tapering-width approach (as compared to taping hole radii) offers a natural way of coupling the nanobeam air-mode cavity to the feeding waveguide. Since the width of the beam is decreasing, the cavity naturally couples to the feeding waveguide. We study T and Q total dependence on the total number of mirror pair segments in the Gaussian mirror (N) using FDTD simulations. As shown in Fig. 8(b), we achieve nanobeam cavity with Q = 3.0 × 106, T = 96% at N = 25.
4.3. Higher order modes of the air-mode cavity
The ultra-high Q cavity that we were able to design is the fundamental mode of the cavity. Higher order modes coexist with the fundamental modes inside the band gap. Fig. 9(a) shows the transmission spectrum of a waveguide-coupled air-mode nanobeam cavity, that has 15 mirror pair segments in the Gaussian mirror. The band-edge modes are observed at wavelengths longer than 1.6μm. The modes in the range of 1.2μm to 1.35μm are formed by the higher order band modes in Fig. 7(a). Figures 9(b)–9(c) show the major field-component distribution (Ey) of the two cavity modes inside the bandgap. The total Q-factors of these two modes are 23,935 and 5,525 respectively. The effective mode volumes are 2.32 and 3.01 respectively.
We have presented a detailed analysis and a deterministic design of the ultra-high Q photonic crystal nanobeam cavities. With this method, Q > 109 radiation-limited cavity, and Q > 107, T > 95% waveguide-coupled cavity are deterministically designed. These Q-factors are comparable with those found in whispering gallery mode (WGM) cavities [45–47]. Meanwhile, the mode volumes are typically two or three orders of magnitude smaller than WGM ones. Furthermore, energy maximum can be localized in either the dielectric region or air region with this method. Although we demonstrate designs for TE-like, transversely symmetric cavity modes, the design method is universal, and can be applied to realize nanobeam cavities that support TM-polarized modes, as well as line-defect 2D photonic crystal cavities. We believe that the proposed method will greatly ease the processes of high Q nanobeam cavity design, and thus enable both fundamental studies in strong light and matter interactions, and practical applications in novel light sources, functional optical components (filters, delay lines, sensors) and densely integrated photonic circuits.
We acknowledge numerous fruitful discussions with M.W. MuCutcheon and P. B. Deotare. This work is supported by NSF Grant No. ECCS-0701417 and NSF CAREER grant.
References and links
1. Quality factor is defined as , and mode volume is defined as V = ∫ dVε|E|2/[ε|E|2]max.
3. J. L. O’Brien, A. Furusawa, and J. Vuckovic, “Photonic quantum technologies,” Nat. Photonics 3, 687–695 (2009). [CrossRef]
4. J. Leuthold, C. Koos, and W. Freude, “Nonlinear silicon photonics,” Nat. Photonics 4, 535–544 (2010). [CrossRef]
6. D. Van Thourhout and J. Roels, “Optomechanical device actuation through the optical gradient force,” Nat. Photonics 4, 211–217 (2010). [CrossRef]
7. D. G. Grier, “A revolution in optical manipulationm,” Nature 424, 21–27 (2003). [CrossRef]
11. J. S. Foresi, P. R. Villeneuve, J. Ferrera, E. R Thoen, G. Steinmeyer, S. Fan, J. D. Joannopoulos, L. C. Kimerling, H. I. Smith, and E. P. Ippen, “Photonic-bandgap microcavities in optical waveguides,” Nature 390, 143–145 (1997). [CrossRef]
13. J. Ctyroky, “Photonic bandgap structures in planar waveguides,” J. Opt. Soc. Am. A 18, 435–441 (2001). [CrossRef]
14. M. R. Watts, S. G. Johnson, H. A. Haus, and J. D. Joannopoulos, “Electromagnetic cavity with arbitrary Q and small modal volume without a complete photonic bandgap,” Opt. Lett. 27, 1785–1787 (2002). [CrossRef]
15. J. M. Geremia, J. Williams, and H. Mabuchi, “Inverse-problem approach to designing photonic crystals for cavity QED experiments,” Phys. Rev. E 66, 066606 (2002). [CrossRef]
16. M. Burger, S. J Osher, and E. Yablonovitch, “Inverse problem techniques for the design of photonic crystals,” IEICE Trans. Electron. E87C, 258–265 (2004).
18. B. S. Song, S. Noda, T. Asano, and Y. Akahane, “Ultra-high-Q photonic double-heterostructure nanocavity,” Nat. Mater. 4, 207–210 (2005). [CrossRef]
19. S. Tomljenovic-Hanic, C. M. de 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] [PubMed]
20. 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 (2006). [CrossRef]
21. K. Nozaki, S. Kita, and T. Baba, “Room temperature continuous wave operation and controlled spontaneous emission in ultrasmall photonic crystal nanolaser,” Opt. Express 15, 7506–7514 (2007). [CrossRef] [PubMed]
22. Y. Tanaka, T. Asano, and S. Noda, “Design of photonic crystal nanocavity with Q-factor of ∼109,” J. Lightwave Technol. 26, 1532 (2008). [CrossRef]
24. P. Velha, E. Picard, T. Charvolin, E. Hadji, J. C. Rodier, P. Lalanne, and E. Peyrage, “Ultra-high Q/V Fabry-Perot microcavity on SOI substrate,” Opt. Express 15, 16090–16096 (2007). [CrossRef] [PubMed]
25. S. Reitzenstein, C. Hofmann, A. Gorbunov, M Strauß, S. H. Kwon, C. Schneider, A. Loffler, S. Hofling, M. Kamp, and A. Forchel, “AlAs/GaAs micropillar cavities with quality factors exceeding 150000,” Appl. Phys. Lett. 90, 251109 (2007). [CrossRef]
26. A. R. Md Zain, N. P. Johnson, M. Sorel, and R. M. De La Rue, “Ultra high quality factor one dimensional photonic crystal/photonic wire microcavities in silicon-on-insulator (SOI),” Opt. Express 16, 12084 (2008). [CrossRef]
28. M. W. McCutcheon and M. Loncar, “Design of a silicon nitride photonic crystal nanocavity with a Quality factor of one million for coupling to a diamond nanocrystal,” Opt. Express 16, 19136–19145 (2008). [CrossRef]
29. L. D. Haret, T. Tanabe, E. Kuramochi, and M. Notomi, “Extremely low power optical bistability in silicon demonstrated using 1D photonic crystal nanocavity,” Opt. Express 17, 21008–21117 (2009). [CrossRef]
31. P. B. Deotare, M. W. McCutcheon, I. W. Frank, M. Khan, and M. Loncar, “High quality factor photonic crystal nanobeam cavities,” Appl. Phys. Lett. 94, 121106 (2009). [CrossRef]
32. Q. Quan, P. B. Deotare, and M. Loncar, “Photonic crystal nanobeam cavity strongly coupled to the feeding waveguide,” Appl. Phys. Lett. 96, 203102 (2010). [CrossRef]
33. E. Kuraamochi, H. Taniyama, T. Tanabe, K. Kawasaki, Y-G. Roh, and M. Notomi, “Ultrahigh-Q one-dimensional photonic crystal nanocavities with modulated mode-gap barriers on SiO2 claddings and on air claddings,” Opt. Express 18, 15859–15869 (2010). [CrossRef]
34. Q. Quan, I. B. Burgess, S. K. Y. Tang, D. L. Floyd, and M. Loncar, “High-Q/V photonic crystal nanobeam cavities in an ultra-low index-contrast polymeric optofluidic platform,” arXiv:1108.2669 (2010).
35. J. Vuckovic, M. Loncar, H. Mabuchi, and A. Scherer, “Optimization of the Q factor in photonic crystal microcavities,” IEEE J. Quantum Electron. 38, 850–856 (2002). [CrossRef]
36. K. Srinivasan and O. Painter, “Momentum space design of high-Q photonic crystal optical cavities,” Opt. Express 10, 670–684 (2002). [PubMed]
38. M. Palamaru and P. Lalanne, “Photonic crystal waveguides: Out-of-plane losses and adiabatic modal conversion,” Appl. Phys. Lett. 78, 1466–1468 (2001). [CrossRef]
39. P. Lalanne, S. Mias, and J. P. Hugonin, “Two physical mechanisms for boosting the quality factor to cavity volume ratio of photonic crystal microcavities,” Opt. Express 12, 458–467 (2004). [CrossRef] [PubMed]
40. K. Sakoda, Optical Properties of Photonic Crystals, 2nd Ed. (Springer, 2005).
41. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Molding the Flow of Light, 2nd Ed. (Cambridge University Press, 2007).
42. S. G. Johnson, M. Ibanescu, M. A. Skorobogatiy, O. Weisbergs, J. D. Joannopoulos, and Y. Fink, “Perturbation theory for Maxwell’s equations with shifting material boundaries,” Phys. Rev. E 65, 066611 (2002). [CrossRef]
43. J. Vuckovic, M. Loncar, H. Mabuchi, and A. Scherer, “Optimization of three-dimensional micropost microcavities for cavity quantum electrodynamics,” Phys. Rev. E 65, 016608 (2001).
45. D. W. Vernooy, A. Furusawa, N. P. Georgiades, V. S. Ilchenko, and H. J. Kimble, “Cavity QED with high-Q whispering gallery modes,” Phys. Rev. A 57, R2293–R2296 (1998). [CrossRef]