## Abstract

We investigate slow light propagation in monomode photonic crystal waveguides with different spectral features such as constant group index, high bandwidth and low group velocity dispersion. The form of the waveguide mode alters dramatically and spans three different spectral intervals by tuning the size of the boundary holes. Namely, slope of the band gap guided mode changes sign from negative to positive toward the Brillouin zone edge. In between there is a transition region where modes have nearly zero slopes. Maximum group index occurs at these turning points at the expense of high dispersion and narrow bandwidth. The apparent trade-off relationship between group index and bandwidth is revealed systematically. We show that as the radius of the innermost hole is increased above a certain value, the former one decreases and the latter one increases both exponentially but with a different ratio. The product of average group index and bandwidth is defined as a figure of merit which reaches up to a value of approximately 0.30 after a detailed parametric search. The findings of the frequency domain analysis obtained by plane wave expansion method are confirmed via finite-difference time-domain study.

©2010 Optical Society of America

## 1. Introduction

Endeavoring to reduce the speed of light may seem to be against the high speed optical communication occurring at close proximity to the speed of light. However, in some circumstances it is imperative to manipulate the speed of light in a controlled manner. For example, nonlinear optics, optical memories, buffers, delay lines and electro- and magneto-optic effects may gain benefit from slow light phenomena due to enhanced light-matter interaction and spatial compression of optical pulses [1–6]. Consequently, the concept of tailoring the speed of light has been an intense research area. The approaches in terms of underlying fundamentals can be mainly grouped into two categories. The first group utilizes the material dispersion in atomic media or solid-state systems. The second category belongs to the structural dispersion of specially designed photonic devices [7–9]. The pioneering studies in the field have been exploited the first type of slow light mechanism which achieved a record breaking speed of 17 m/s [10]. The two major problems associated with these approaches are the limited bandwidth nature and the practical constraints due to the requirement of special environments [11]. Carefully engineered dielectric structures show great promise to tame and manipulate the speed of light by lowering speed many orders of magnitude even though that value is still far beyond the achievement of the approaches belonging the first group.

The beauty of reducing speed of light based upon structural dispersion comes from the relatively easy engineering of the complex structures at will. Photonic crystals (PC) are periodic dielectric structures that may exhibit band gap phenomena so that we can modify them by inducing defects of certain type to obtain cavities and waveguides [12]. As a result, light can be confined in the desired regions or guided through certain symmetry directions. Photonic crystal waveguides (PCW) show dispersive characteristics so that the group velocity varies enormously within a certain *k*-space. The two regions which have low group velocities are near the Brillouin zone edges, namely, $({k}_{x},\text{\hspace{0.17em}}{k}_{y})$ = $(0,\text{\hspace{0.17em}}0)2\pi /a$ or $(0,\text{\hspace{0.17em}}0.5)2\pi /a$. But near these points, band shape resembles quadratic behavior which naturally accompanies large group velocity dispersion (GVD) and small bandwidth. As a result, the optical signal gets distorted after propagating a certain distance [13]. One solution to tackle with the GVD lies behind the idea of obtaining flat bands. So having high group index $({n}_{g})$ values with small GVD and large bandwidth are important research objectives [14–29] because of the diverse applications of slow light as mentioned above. One crucial aspect that should be kept in mind while obtaining flat bands with constant slopes is the relatively simple structures that can be easily realized. In this work, we obtain different families of dispersion diagrams which ensure large and constant group index in addition to high bandwidth operation. This aspect of the work consists of one part of the study.

The second part as explained later in the paper is to explore systematically the trade-off relation between bandwidth and group index. We emphasized the importance of side rows of holes to reduce the speed of light. The band shape changes dramatically when the innermost holes’ radii reach a certain value as the rest of the background holes’ radii are kept constant. An inverse relation persists between the bandwidth and average group index. One may have either a large group index value of 590 or a small value of 14 with the corresponding normalized bandwidths of $5.35\times {10}^{-5}$ and 0.011, respectively for a certain filling factor. As we can see later, while one parameter diminishes by more than one order of magnitude the other parameter increases by more than two orders of magnitude. In the previous studies, attention was given to the flat-band creation and a figure of merit is defined as the product of “average group index and bandwidth”. The trade-off relation between the two parameters was mentioned in previous studies but none of them reported an exponential dependency. Flat-band occurs at two different spectral regions, one is close to dielectric band and the second region is close to air band. In between the two regions, there is a slow light waveguide mode whose band slope is almost zero. This in turn yields the most dispersive mode but at the same time the largest group index. The optimization of group index and bandwidth product was carried out in one of the flat bands region and a value of approximately 0.30 is obtained. The finite-difference time-domain (FDTD) method is used to verify propagation of slow light pico-second optical pulse [30]. An appropriate bandwidth of optical pulse (12 nm) is at 1550 nm selected and relevant parameters are extracted by the time-of-flight method. Both index and band gap guided modes are excited and the pulse separation in time domain was observed.

## 2. Slow light photonic crystal waveguide structure and frequency domain analysis

The intact photonic crystal (PC) used in this paper is a two dimensional structure with triangular lattice air holes in dielectric background as shown in Fig. 1 . The air hole radius is ${r}_{b}=0.30a$, where a is the lattice constant and background refractive index is 3.46. The symmetry directions are also shown in the figure. The waveguide is constructed by filling the center holes with the same background dielectric material along $\mathrm{\Gamma}\mathrm{{\rm K}}$ direction. The dispersion diagram of the PCW as shown in Fig. 2 is calculated by plane wave expansion method [31] using MIT’s freely available software MPB. The waveguide mode of even type is represented by a solid line. The odd mode is also included in the figure as a dashed-line. The group index value that corresponds to the waveguide mode is demonstrated in Fig. 3 . The group index is defined as ${n}_{g}=c/{v}_{g}$, where c is the speed of light. The group velocity as represented by ${v}_{g}$ is the first order derivative of waveguide mode, i.e., ${v}_{g}=\left|\partial \omega /\partial k\right|$. Due to the polarization sensitivity of 2D PCs only TE modes (electric fields are in the plane) are considered in the study.

Considering the out-of-plane losses, we focus on only region of interest in the reciprocal space where mode stays under the light line. This part of the region is apparent if we only present dispersion diagrams within the interval from $k=0.25(2\pi /a)$ to band-edge. Figure 3 shows group index variation of a regular PCW. As evident from this plot, it diverges very fast at the close proximity of cut-off frequency. The two main problems of such a waveguide mode are the narrow bandwidth operation and high GVD. The way of eliminating these problems lies behind the idea of obtaining flat-bands (non-zero, constant slope). Flat bands with small slopes provide constant group index and have very low GVD.

As the field distribution of the guided mode is highest in the center of the waveguide and decreases along transverse to propagation direction, the first rows of holes affect the field distribution and the dispersion diagram more than the second and third rows. Therefore, in this study we modified the side rows holes’ size and changed symmetrically their radii along the waveguide.

The dispersion diagram of PCW contains various curves as shown in Fig. 4
when we change the first row of holes’ radii from ${r}_{d}=0.20a$ to $0.4750a$ with an incremental step of 0.0125*a*. The background holes radii are kept constant at ${r}_{b}=0.30a$. Within the resolution of our parameter scan, we realized that when hole radius is around 0.36*a* and beyond that first flat-band region emerges in the dispersion diagram. For the values of hole radii between 0.30*a* and 0.36*a*, the waveguide mode is very dispersive and there is no signature of flat-band region. The second flat-band regime occurs when ${r}_{d}$ takes values around 0.20*a*. Most of the previous studies targeted to tune geometric parameters within a small interval. Hence, some of the general trends that may be captured were unrevealed. As a result, large structural variations reveal important observations. The first observation in the dispersion diagram which is also a trivial one is the bands movement to higher/lower frequencies as we increase/decrease the side row of holes’ radii. The band flattens as the waveguide mode comes close to continuum of dielectric band. The flatness disappears and mode becomes highly dispersive at intermediate values of${r}_{d}$ due to hyperbolic forms of bands. We will demonstrate later in the paper that, there is always a turning point (${r}_{d}=0.35a$, designated with a marker in the figure) as we change the filling factor of the PC by altering ${r}_{d}$. This is the second important observation that can be drawn from the dispersion diagram. The slope of the band is almost zero at the turning point and below/above this value the bands have negative/positive slopes. The upward movement speed increases even though the amount of increment is kept the same. The three spectral intervals are roughly separated by double-arrows in Fig. 4. While the lower and upper parts contain flat-band regions suitable for slow light application the middle part is highly dispersive.

To gather group index variation of the bands in the three regions, we prepared Fig. 5 . The band which has the smallest slope is also indicated with a marker. There are two interesting regimes in the figure in terms of slow light perspective. Bands with constant group indices appear both at lower and upper parts of the spectrum. The first one is close to dielectric band continuum states. On the contrary, the second part is close to continuum states of the air band. The same types of lines correspond to same modes in the plots of $\omega (k)$ and ${n}_{g}(\omega )$ in Figs. 4 and 5, respectively. Hence, the two figures complement each other. The useful low dispersion area in Fig. 4 can be easily detected by inspecting Fig. 5. The constant group index regions plotted with red-dotted lines imitate a U-shaped behavior which is the region of low GVD. We selected three representative points (shown with arrows) in the group index spectrum where the corresponding group indices are 27.9, 43.4 and 10.4, respectively to investigate e-field distributions. The result is shown in Fig. 6 . The e-field penetrates more towards transverse to propagation direction in Fig. 6(a) than the rest. The field components reach up to the 4th row. This is because of the relevant mode closeness to dielectric continuum states. On the other hand, it decays very quickly over the 1st and 2nd rows for the other two cases in Figs. 6(b) and 6(c). The unique behavior of e-field distributions for different bands tells us that one may play with the rows up to 3rd or 4th row to engineer slow light in lower part of the spectrum. However, 1st and 2nd row have the most pronouncing role to manipulate mode shape for slow light in the upper part of the group index spectrum.

The nice behavioral change between group index and bandwidth in the upper part of the group index spectrum (Fig. 5) directed us to extract quantitative values. We re-produced group index variation for ${r}_{d}$ changing from 0.3625*a* to 0.45*a* keeping ${r}_{b}$ fixed at 0.30*a*. Different *y*-scales are used in Figs. 7(a)
and 7(b). From these plots, constant-group index and the corresponding bandwidth values are extracted. We can see that there are constant group index regions with different bandwidths. The ${n}_{g}$ values are U-shaped and the bottom of this U shape is very appropriate for obtaining constant group index with small GVD and high bandwidth. We extracted the value of minimum group index and bandwidth from Fig. 7. The bandwidth is defined as $\Delta \omega /{\omega}_{0}$, where $\Delta \omega $ is the difference between maximum and minimum frequencies and ${\omega}_{0}$ is the center frequency of the U-shape. The maximum and minimum frequencies are determined so that corresponding ${n}_{g}$ values are within the 10% of the minimum group index. The corresponding variations of bandwidth and group index are shown in Figs. 8(a)
and 8(b), respectively. As we can see, while bandwidth increases exponentially, the group index decreases also exponentially but with a different decay rate.

Since the group index and bandwidth values are inversely proportional, we had to determine a figure of merit (FOM). We first evaluated average group index and bandwidth values and then defined the product of them as a FOM. The result is shown in Fig. 9
. The FOM increases almost five times reaching a value of 0.16 as we increase air holes radii of the first rows from 0.3625*a* to 0.45*a*.

The GVD graph calculated by evaluating ${\partial}^{\text{2}}k/\partial {\omega}^{2}$ is depicted in Figs. 10(a)
and 10(b) for radii of holes 0.3875*a* and 0.45*a*, respectively. The low GVD reduces down to the order of 10^{5}
${\text{ps}}^{\text{2}}\text{/km}$. We should note that GVD changes symmetrically with positive and negative parts which can be used for dispersion compensation application.

## 3. Time domain analysis of slow light propagation in photonic crystal waveguide

In this section, the slow light results are reinforced by the time domain study using FDTD method. It allows us to visualize slow light propagation inside PCW. The structure length is taken to be 300*a*. The computational domain is surrounded by perfectly matched layer boundary condition. Different detection points are placed throughout the structure and the distance between each detection point is 32*a*. We made simulations on the structure with side row radii equal to 0.4375*a*. The central frequency of Gaussian source is fixed at $\omega a/2\pi c=0.2452$with a normalized bandwidth of 0.001. We excited the slow light mode by sending a tilted input source. After a certain distance which is proportional to the difference between the slopes of the waveguide mode at two different *k* points, the separation of the input pulse into fast (index-guided) and slow modes (band-gap guided) is easily observable.

The group index extraction is performed by tracing arrival times of the central peaks of the slow light. The lattice constant *a* is fixed at 380 nm targeting the optical regime of 1550 nm. The bandwidths of the input pulse and constant group index region are approximately 6.3 nm and 12.6 nm, respectively. The pulse is allowed to propagate for more than 10000 *a*/*c* time. By measuring delay time of the pulses at different locations, we obtain Fig. 11(a)
. The slope of the curve which approximates FDTD data gives group index value of ${n}_{g}=22.85$. The delay time information of the first detection point is scaled to zero time to take it as a reference. This number is very close to the frequency domain result of ${n}_{g}=20.59$. It takes 8.69 ps for the input pulse to travel 300*a* length PCW. The delay time can be increased either by using longer PCW or using different rod radii value. By using different rod radii for innermost holes, larger group index value can be attained at the expense of narrowing the bandwidth.

Figure 11(a) also indicates time delay information for fast light. As can be seen from Fig. 7(b), fast light has a smaller group index (lower tail part of dispersion line). As a result, the slope of the curve corresponding to fast light is almost three times smaller than slow light. The temporal evolutions of slow light optical pulses while propagating down the waveguide are presented in Figs. 11(b)–11(e) at four selected detection points, 32*a*, 96*a*, 160*a* and 224*a*, respectively. The absence of broadening and distortion in the pulse waveforms are an indication of negligible GVD. The fast light reaches the end of the waveguide earlier than slow light. Hence, the reflected portion of it interferes with the traveling slow light. As a result, a small broadening in the pulse form appears only in Figs. 11(d) and 11(e).

The pulse propagation through PCW is presented in Fig. 12
. There are two detection windows used for the observation of light propagation, one is placed in the middle of the waveguide while the other one is located at the end of the waveguide. First fast light appears at the observation plane as expected due to small group index as compared with the slow light which has almost three times larger group index. As can be seen from Fig. 12(a), e-field distribution resembles well the index guided mode pattern. After certain time delay, the slow light arrives with distinct field patterns such as diffusing more towards the transverse direction and appearing beating patterns along the longitudinal direction as indicated in Fig. 12(b). To prove further the dispersion free propagation of the slow light, we provided another snap-shot of the e-field of slow light at the end of the waveguide. That means the field in Fig. 12(b) is allowed to travel almost a distance of 300*a* to reach the end of the waveguide. It can be seen from the plot in Fig. 12(c) that the pulse propagation does not suffer from any group velocity dispersion. Media 1 shows the optical pulse propagation along the waveguide.

## 4. Optimization of figure of merit and loss analysis

The previous analyses were performed by keeping background holes’ radii at ${r}_{b}=0.30a$. It is worth considering the other cases where ${r}_{b}$ takes values smaller and larger than 0.30a. As a result, we consider the change in the figure of merit when we scan different combinations of ${r}_{b}$ and ${r}_{d}$.

Figure 13(a)
shows the color-map plot of average group index and bandwidth product with respect to ${r}_{b}$ and ${r}_{d}$. The FOM value reaches up to 0.30 at the upper right corner of the figure where ${r}_{b}$ and ${r}_{d}$ take larger values, 0.42*a* and 0.48*a*, respectively. To underpin the constituent parameters variation clearly, we show group index and bandwidth variation along the dotted line which corresponds to ${r}_{b}=0.36a$in Figs. 13(b) and c. The insets show the enlarged view of the selected portions. When group index takes the largest value for certain ${r}_{d}$, the corresponding bandwidth is at the smallest value. Hence, defined FOM is a signature designating the inverse relation between the two parameters. The optimizations were performed only considering the first two rows of holes. It may be possible to enhance FOM further if one considers also the role of second row of holes. If higher FOM is sought for, one should use larger ${r}_{b}$ and ${r}_{d}$ values. On the other hand, larger group index with small bandwidth is available for various ${r}_{b}$ values when ${r}_{d}$ is kept relatively smaller. FOM values in the previous studies can be as low as on the order of 0.01 [28] and may increase up to 0.85 by judicious waveguide design [26]. The FOM of the present work as compared to them stays within an intermediate value.The parametric scan of radii allowed us to identify one additional feature of slow light PCW. For each radius of background holes, there is a certain ${r}_{d}$ value of first row of holes where the bands have the smallest slope (close to zero within the numerical resolution). The U-shaped group index region starts to appear above this ${r}_{d}$ value. Table 1
shows the specific ${r}_{d}$ values which are transition regions between lower flat bands and upper flat bands. The modes are highly dispersive at and around these ${r}_{d}$ values. When we increase ${r}_{b}$ from 0.28*a* to 0.38*a*, the turning points first decrease and then start to increase. Table 1 informs us which region to search for slow light with flat band if we fix the filing factor of PCW. The turning point is reminiscent of the existence of lower flat band. As we increase ${r}_{b}$ above 0.38*a*, continuum states of dielectric band intervene strongly with the waveguide mode and lower flat band disappears. Hence the trend in Table 1 can be presented up to a certain ${r}_{b}$ value.

A detailed loss analysis is out of the scope of the present work. However, to have an idea on the performance of the proposed slow light structure under the loss, we performed some additional work. Among the various loss mechanisms that may occur in PCW, we focused on the loss due to backward scattering. For this purpose, we introduced randomly varying hole diameters within an interval of 0.4375*a*±0.01*a*. Such randomness may occur as a consequence of an imperfect structure fabrication. The background radii of holes are kept constant at ${r}_{b}=0.30a$. The expected group index value of the waveguide without inclusion of loss parameter is ${n}_{g}=22.85$. However, uniformly distributed and randomly varying hole diameters’ should modify the dispersion curve. As a result, group index may increase or decrease correspondingly. The length of the perturbed part within the PCW is varied up to 40*a* while keeping the rest of the structure unchanged. The result is presented in Fig. 14
. The estimated loss coefficient extracted within the limited resolution of FDTD discretization is approximately 0.05dB/*a*. At first glance, such a loss value may seem to be high if a comparison is made with the literature [14]. However, careful examination of the way of introducing fabrication loss component into the current analysis provides a reasonable explanation. The holes’ radii of first rows are artificially perturbed from the original value, 0.437*a* to 0.4275*a* and 0.4475*a* (extreme cases). In these situations, dispersion analyses yield constant group indices of the following values 26.92 and 14.56 with different band widths. As compared to the loss free case, group index change occurs within a large interval of 20.59 ± (~6). Besides, slow light waveguide mode moves to higher/lower frequencies and in some cases there may be no overlap between the constant group index regions. As a result, slow light mode excited for ${n}_{g}=20.59$ faces with a huge scattering loss due to group index mismatch and band movement. Uniformly distributed perturbation parameters also exaggerate the fabrication disorder loss. Instead, a more realistic distribution such as Gaussian distribution can be adapted. Finally, field distribution of the slow light mode plays a crucial role for the loss amount. Strongly concentrated fields especially around the perturbed area are more amenable to high loss. It is possible to tailor the disorder induced losses in slow light waveguides by manipulating appropriately the locations of the e-field concentration [32]. The electric field concentration should be distributed at regions where disordered places do not overlap with the high intensity of the field. The out of plane losses in PCWs will not be a major concern as long as the operating spectral region is selected to be below the light line [33,34]. In the present work, slow light region appears below the light as can be seen from Fig. 2.

## 5. Conclusion

In conclusion, we have studied tuning of group index and bandwidth values by means of obtaining flat dispersion curves in a simplistic manner. We performed detailed analysis to engineer the dispersion and show that complete scan of parameters is very beneficial for slow light concept to identify and categorize different spectral trends that appear in dispersion diagrams. The frequency domain calculations confirm that when the radii of the side rows of holes increase above certain value for each ${r}_{b}$ then unique properties in the waveguide mode appear due to the interaction with continuum states (from below/above: dielectric/air band). The inverse exponential relation between group index and bandwidth is explored. The defined figure of merit shows increment with respect to increase in the radii of the side rows of holes reaching a value of approximately 0.30. The group velocity of 0.0017*c* is achieved in the modified photonic crystal waveguide. The transmissions of pico-second optical pulses were realized by time domain study. The presented results obtained by simple geometrical modifications via detailed parameter scan guide us how to obtain large bandwidth, constant group index and low GVD for slow light applications.

## Acknowledgements

The authors gratefully acknowledge the financial support of the Scientific and Technological Research Council of Turkey (TUBITAK), Project no: 108T717. HK also acknowledges support from the Turkish Academy of Sciences Distinguished Young Scientist Award (TUBA GEBIP).

## References and links

**1. **R. W. Boyd, D. J. Gauthier, and A. L. Gaeta, “Applications of slow light in telecommunications,” Opt. Photon. News **17**(4), 18–23 (2006). [CrossRef]

**2. **T. F. Krauss, “Why do we need slow light,” Nat. Photonics **2**(8), 448–450 (2008). [CrossRef]

**3. **T. Baba, “Slow light in photonic crystals,” Nat. Photonics **2**(8), 465–473 (2008). [CrossRef]

**4. **F. Xia, L. Sekaric, and Y. Vlasov, “Ultracompact optical buffers on a silicon chip,” Nat. Photonics **1**(1), 65–71 (2007). [CrossRef]

**5. **T. Tanabe, M. Notomi, S. Mitsugi, A. Shinya, and E. Kuramochi, “Fast bistable all-optical switch and memory on a silicon photonic crystal on-chip,” Opt. Lett. **30**(19), 2575–2577 (2005). [CrossRef] [PubMed]

**6. **C. Liu, Z. Dutton, C. H. Behroozi, and L. V. Hau, “Observation of coherent optical information storage in an atomic medium using halted light pulses,” Nature **409**(6819), 490–493 (2001). [CrossRef] [PubMed]

**7. **T. Baba and D. Mori, “Slow light engineering in photonic crystals,” J. Phys. D Appl. Phys. **40**(9), 2659–2665 (2007). [CrossRef]

**8. **T. F. Krauss, “Slow light in photonic crystal waveguides,” J. Phys. D Appl. Phys. **40**(9), 2666–2670 (2007). [CrossRef]

**9. **J. B. Khurgin, “Slow light in various media: a tutorial,” Adv. Opt. Photon. **2**(3), 287–318 (2010). [CrossRef]

**10. **L. Hau, S. Harris, Z. Dutton, and C. Behroozi, “Light speed reduction to 17 metres per second in an ultracold atomic gas,” Nature **397**(6720), 594–598 (1999). [CrossRef]

**11. **R. S. Tucker, P.-C. Ku, and C. J. Chang-Hasnain, “Slow-light optical buffers-capabilities and fundamental limitations,” J. Lightwave Technol. **23**(12), 4046–4066 (2005). [CrossRef]

**12. **J. D. Joannopoulas, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of the Light. Princeton, NJ: Princeton Univ. Press, (1995).

**13. **M. Notomi, K. Yamada, A. Shinya, J. Takahashi, C. Takahashi, and I. Yokohama, “Extremely large group-velocity dispersion of line-defect waveguides in photonic crystal slabs,” Phys. Rev. Lett. **87**(25), 253902 (2001). [CrossRef] [PubMed]

**14. **L. H. Frandsen, A. V. Lavrinenko, J. Fage-Pedersen, and P. I. Borel, “Photonic crystal waveguides with semi-slow light and tailored dispersion properties,” Opt. Express **14**(20), 9444–9450 (2006). [CrossRef] [PubMed]

**15. **Y. Hamachi, S. Kubo, and T. Baba, “Slow light with low dispersion and nonlinear enhancement in a lattice-shifted photonic crystal waveguide,” Opt. Lett. **34**(7), 1072–1074 (2009). [CrossRef] [PubMed]

**16. **M. D. Settle, R. J. P. Engelen, M. Salib, A. Michaeli, L. Kuipers, and T. F. Krauss, “Flatband slow light in photonic crystals featuring spatial pulse compression and terahertz bandwidth,” Opt. Express **15**(1), 219–226 (2007). [CrossRef] [PubMed]

**17. **A. Yu. Petrove and M. Eich, “Zero dispersion at small group velocities in photonic crystal waveguides,” Appl. Phys. Lett. **85**(21), 4866–4868 (2004). [CrossRef]

**18. **H. Kurt, H. Benisty, T. Melo, O. Khayam, and C. Cambournac, “Slow-light regime and critical coupling in highly multimode corrugated waveguides,” J. Opt. Soc. Am. B **25**(12), C1–C14 (2008). [CrossRef]

**19. **J. Ma and C. Jiang, “Demonstration of Ultraslow Modes in Asymmetric Line-Defect Photonic Crystal Waveguides,” IEEE Photon. Technol. Lett. **20**(14), 1237–1239 (2008). [CrossRef]

**20. **J. Hou, D. Gao, H. Wu, R. Hao, and Z. Zhou, “Flat Band Slow Light in Symmetric Line Defect Photonic Crystal Waveguides,” IEEE Photon. Technol. Lett. **21**(20), 1571–1573 (2009). [CrossRef]

**21. **J. Li, T. P. White, L. O’Faolain, A. Gomez-Iglesias, and T. F. Krauss, “Systematic design of flat band slow light in photonic crystal waveguides,” Opt. Express **16**(9), 6227–6232 (2008). [CrossRef] [PubMed]

**22. **S. Rawal, R. K. Sinha, and R. M. De La Rue, “Slow light miniature devices with ultra-flattened dispersion in silicon-on-insulator photonic crystal,” Opt. Express **17**(16), 13315–13325 (2009). [CrossRef] [PubMed]

**23. **R. Hao, E. Cassan, H. Kurt, X. Le Roux, D. Marris-Morini, L. Vivien, H. Wu, Z. Zhou, and X. Zhang, “Novel slow light waveguide with controllable delay-bandwidth product and utra-low dispersion,” Opt. Express **18**(6), 5942–5950 (2010). [CrossRef] [PubMed]

**24. **O. Khayam and H. Benisty, “General recipe for flatbands in photonic crystal waveguides,” Opt. Express **17**(17), 14634–14648 (2009). [CrossRef] [PubMed]

**25. **R. Hao, E. Cassan, H. Kurt, J. Hou, X. Le Roux, D. Marris-Morini, L. Vivien, G. Dingshan, Z. Zhou, and X. Zhang, “Novel Kind of Semislow Light Photonic Crystal Waveguides With Large Delay-Bandwidth Product,” IEEE Photon. Technol. Lett. **22**(11), 844–846 (2010). [CrossRef]

**26. **L. Dai and C. Jiang, “Ultrawideband Low Dispersion Slow Light Waveguides,” J. Lightwave Technol. **27**(14), 2862–2868 (2009). [CrossRef]

**27. **L. Dai and C. Jiang, “Photonic crystal slow light waveguides with large delay-bandwidth product,” Appl. Phys. B **95**(1), 105–111 (2009). [CrossRef]

**28. **J. Ma and C. Jiang, “Flatband slow light in asymmetric line-defect photonic crystal waveguide featuring low group velocity and dispersion,” IEEE J. Quantum Electron. **44**(8), 763–769 (2008). [CrossRef]

**29. **F. Wang, J. Ma, and C. Jiang, “Dispersionless slow wave in novel 2-D photonic crystal line defect waveguides,” J. Lightwave Technol. **26**(11), 1381–1386 (2008). [CrossRef]

**30. **A. Taflove, Computational Electrodynamics – The Finite-Difference Time Domain Method. Norwood, MA: Artech House, 2000.

**31. **S. G. Johnson and J. D. Joannopoulos, “Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis,” Opt. Express **8**(3), 173–190 (2001). [CrossRef] [PubMed]

**32. **A. Petrov, M. Krause, and M. Eich, “Backscattering and disorder limits in slow light photonic crystal waveguides,” Opt. Express **17**(10), 8676–8684 (2009). [CrossRef] [PubMed]

**33. **E. Kuramochi, M. Notomi, S. Hughes, A. Shinya, T. Watanabe, and L. Ramunno, “Disorder-induced scattering loss of line-defect waveguides in photonic crystal slabs,” Phys. Rev. B **72**(16), 161318 (2005). [CrossRef]

**34. **M. Notomi, A. Shinya, S. Mitsugi, E. Kuramochi, and H. Ryu, “Waveguides, resonators and their coupled elements in photonic crystal slabs,” Opt. Express **12**(8), 1551–1561 (2004). [CrossRef] [PubMed]