We demonstrate a novel type of slow light photonic crystal waveguide which can produce unusual “U” type group index - frequency curves with constant group index over large bandwidth. By shifting the boundaries of this waveguide, flexible control of (10 << 210) with large bandwidth (1nm<Δλ<43nm centered at 1550nm) and normalized Delay-Bandwidth Product (0.1363<DBP<0.3143) are achieved. Additionally, depending on the chosen waveguide geometry, extremely low group velocity dispersion (GVD<0.5 ps.nm−1.mm−1), with controllable group velocity dispersion of both signs is obtained.
©2010 Optical Society of America
Slow light technology shows an amazing future for optical buffering, optical logical gates and all-optical signal processing in the next-generation photonic networks and optical integrated circuits [1–5]. Up to now, there are mainly three methods to exploit slow light effects: electromagnetically induced transparency , coherent population oscillation , and optical coupled resonators or photonic crystal waveguides (PCW) [4,5]. Slow light in PCW has received more attention because it works at room temperature, provides larger bandwidths if compared with the other two solutions, and can be fabricated using standard semiconductor planar technologies.
Among all PCW devices, W1 waveguides formed by removing the central row of holes/rods in a perfect two-dimensional planar photonic crystal (PC) have been studied thoroughly [8,9]. These waveguides sustain guided mode inside photonic band gap and produce slow light at the band edge where the dispersion relation has a parabolic form. However, large group index in W1 waveguides is only obtained in extremely narrow bandwidth and large group velocity dispersion (GVD) effects accompany the slow light regime, with as a result causes large signal distortion during pulse propagation .
So far, two main solutions in PCW have been proposed to eliminate the drawbacks due to limited bandwidth and obtain dispersion free slow light. The first one is based on nearly-zero dispersion by adjusting the PC geometry, including changing the hole size  or hole position  of the first two rows, and introducing annular holes for the whole lattices  or only the first row [14,15]. The second one is based on dispersion compensation, achieved by adding two opposite dispersion regions in a chirped waveguide [16,17].
We explore in this paper a much simpler solution to provide slow waves in the largest possible bandwidth. The proposed solution is based on a single parameter change with respect to W1 waveguide but brings advantages to achievable bandwidth for large group indices and low GVD values.
This paper is organized as follows. Frequency band engineering of the proposed structure is studied firstly. The physical explanations of the related phenomena are then given by comparing the electric-field distributions of several modified PC waveguides. The ultra-low dispersion properties of the proposed structure are studied, showing that it is well suited to the versatile control of slow waves with chosen group index and bandwidth.
The new waveguides could be useful for future integrated photonic applications such as optical delay lines and dispersion compensating devices.
2. The new waveguide geometry
The proposed structure starts from a W1 waveguide obtained by removing one row of holes in a triangular PC lattice. Then, the two rows of air holes that border the waveguide axis are deliberately shifted by a δx distance. Standard silicon on insulator (SOI) is considered as the typical reference technology. The considered indices for air (top layer), silicon layer (middle layer), and silica (bottom layer) are 1.0, 3.45 and 1.45, respectively. The silicon layer of SOI wafer has a thickness of 283 nm, and the radius of the air hole is 0.286a, with a the PC lattice constant. The whole structure is depicted in Fig. 1 (a).
Obviously, PCW geometry can be treated as a one-dimensional waveguide. Yet, shifting the two bordering of holes introduces another one-dimensional corrugation, which drastically changes the slow light nature. We found the group index – normalized frequency curves turn to the form of letter “U” when δx is increased to appropriate values. Figure 1(b) shows the situation for δx = 0.27a as a typical example. In this very first situation, it is observed that a nearly flat group index with values above 100 is obtained in a fairly large bandwidth.
3. Physical analysis
To understand this surprising phenomenon, we studied in detail the influence of the δx shift quantity on the waveguide dispersion curves using the plane wave expansion method (PWE). The investigated guide modes are supposed to transverse-electric (TE)-like polarization, i.e. with magnetic field parallel to the vertical direction. Due to the computational cost of three-dimensional simulations, two-dimensional simulations with a slab effective index 2.98 were considered here.
Figure 2(a) shows the calculated dispersion diagrams of the proposed structure when δx is gradually varied from 0 to 0.5a. Results for δx in the range between 0.5a and a would be the same as those obtained for δx between 0 and 0.5a due to the waveguide symmetry in x direction, and are thus not depicted. When there is no shift (δx = 0, regular PCW) the waveguide modes show two distinct regions below the light line. These are the fast light region (index-guiding regime) and slow light region (photonic gap guiding regime) . The two regions have a turning point, i.e. the slope of the curve changes from linear one to a quadratic form. When δx starts to increase, the waveguide mode shifts to upper frequencies. The band nature then changes drastically due to the different variations at the two areas around k = 0.36(2π/a) and k = 0.50(2π/a). For the area near the first Brillouin zone (BZ) edge (k = (0.5)2π/a), band moves upwards particularly fast, while for the area near the light line, the band moves less pronounced. As a result, the tailored band curve becomes obviously curl up. When δx is increased from 0a to 0.2a, bands shift little to the higher frequencies and became more and more flat, which means the group index () gradually increases. When δx is about 0.25a, the flattest band curve is obtained, corresponding to the largest. After 0.25a, with continuous increase of δx, band curves are no longer flat. However, with continuous increase of δx, quasi-linear evolution curves are observed between the high moving speed part and the low moving speed part. This predicts a constant over a large frequency range [11,12,14,15].
The evolution of the PCW dispersion diagram with δx can be interpreted from the analysis of field profiles within the structure for different wave vectors and different values of δx. Figsures 2(b) (c) (d) give the electric field distributions for seven δx values at three specific wave numbers: K = 0.36881(2π/a), K = 0.44059(2π/a), K = 0.495(2π/a), respectively.
As shown in Fig. 2(b), when the wave vector is near the light line, field leak through sideways is limited. This is due to the fact that the confinement of E-field is then dominated by the index-like effect . However, when the wave vector is close to BZ edge as in Fig. 2(d), field is less confined within the central waveguide region because confinement of waves turns to be dominated by band gap effect . As a result, the waveguide mode at the BZ edge is more sensitive to the shift of the air holes if compared with those near the light line. This explains why the dispersion diagram evolution when δx varies is small around the light line near K = 0.36881(2π/a). Moreover, the degree of the localization of the field within the air holes increases with increasing δx values. As waves at the BZ edge are more and more strongly confined inside air holes, the effective index seen by the field decreases, which means that bands are pushed towards higher values.
To better illustrate the waveguide properties, band diagram and group index curves are simultaneously plotted in Fig. 3 for three values of δx. The black line in Fig. 3 represents un-shifted structure (a normal PCW), the brown circled line and the blue triangle line represent the shifted structure with δx = 0.28a and 0.5a, respectively. Two inflection points “a” and “c” should be noticed because their group velocities both go down to zero but have different origins. Point “c” is the well-known slow light band edge formed by the interference of forward and backward waves . However, point “a” is believed to be formed by shifting the boundary of the waveguide in our proposed structure. It is clearly shown in Fig. 3(b) that inflection point “a” corresponds to a new sharp peak in the middle of K-path (in the region K = 0.36(2π/a)~0.38(2π/a)). This means we now have a new slow light center in the middle of the K-path. Regarding to the slow light at the band edge, we have two slow light peaks in the K-path, which form this U type.
Somehow, the cost to pay to obtain these original dispersion curves lies in the fact that the new PhC waveguides are not strictly single-mode for all δx values. It is clear from Fig. 3 that two modes with the same even symmetry with respect to the waveguide axis are present at the same normalized wavelength: a/λ = 0.255 and δx = 0.5a can be typically considered to see this. It is yet emphasized here that one of these two modes is situated above the light line, making it intrinsically leaky. As shown in , such modes strongly decay when light propagation lengths are more than 100µm.
All PhC waveguide structures proposed in this work operate in single mode situation below the light line. This makes them suitable for the practical realization of slow wave devices.
Moreover, due to the same geometry as W1 waveguide, light coupling into the proposed slow light waveguides could be optimized using previous coupling approaches to W1 waveguides, e.g.: gradual lattice tapers [19,20] or efficient injectors smartly exploiting evanescent modes . As one step in that direction, we will not focus hereafter on the issue of light coupling into the new proposed waveguides, but rather analyse their slow light performances in terms of group index, bandwidth, and GVD.
4. Slow light performance
4.1 Group index characteristics
The related curves with normalized frequency under different values of δx are plotted in Fig. 4 . To clearly illustrate all situations, larger curves are plotted in Fig. 4(a) and smaller curves are plotted in Fig. 4(b), respectively.
As shown in Fig. 4, all calculated curves turn to an unusual form of “decrease—constant — increase”, hereafter called the “U-type” curves. As a consequence, nearly constant over a wide bandwidth have been achieved due to the flatness of the U-type curve center.
Furthermore, a general trend can be drawn here: when δx increases, the average group index () value in the constant region decreases, at the same time the corresponding bandwidth increases. This trend predicts there exists a trade-off between values and bandwidth values in our device, like in all other slow light devices.
4.2 Delay-bandwidth product
The concept of group index describes how slow the light is. However, for slow light device, only talking about is meaningless , because we must care about the bandwidth at the same time. The delay-bandwidth product (DBP) is a good indication of the highest slow light capacity that the device potentially provides . We will focus on the normalized delay-bandwidth product (NDBP) value in this paper, because it is universal for comparison between devices having different lengths and different operating frequencies. The NDBP is defined by:
Up to now, there are many papers studying about flat band slow light [11–17,22–25]. However, it is difficult to make comparison between those papers, because different papers have given NDBP according to different flat judgments. For example, ref  has pointed out NDBP value of 0.24 within 10% variation with respect to their mean value, and ref  has pointed out value of 0.6 within variation much larger than 100% with respect to their mean value. It is not correct to say ref  has better NDBP performance than ref .
Before we start to investigate the NDBP, the concept of “flat ratio” is introduced for the first time to clearly define the flatness within the related bandwidth.
The flat ratio (FR) is defined as:
Table 1 summaries the values of average group index and NDBP values. In accordance with refs [11,12,23], the flat ratio considered here is 0.2 corresponding to a maximum of 10% variation with respect to the mean value.
Results gathered in Table 1 show that NDBP values increase when constant decreases. This tendency is due to the fact that the corresponding bandwidth values increase more rapidly. In comparison with the previous works, it is observed that comparable NDPB values are obtained with a much simpler PC tuning geometry [11–14].
To the best of our knowledge, the largest constant with a bandwidth reported in PCs is 105 with a bandwidth 2.3 nm , but it requires a complicated structure. Our results provide even larger constant with a similar delay-bandwidth product performance, for example: = 211 with ∆ω/ω = 0.000645 at δx = 0.25a, = 140 with ∆ω/ω = 0.001 at δx = 0.26a. It is underlined here that these results reveal the possibility to a versatile control of and bandwidth with only one simple structural parameter.
5. Dispersion performance
The GVD parameter Dλ is calculated as:
Figure 6 shows the calculated values of Dλ of the proposed structure. It is clearly shown from Fig. 6 that there exists a bandwidth of “zero” dispersion in each situation, and the bandwidth of “zero” dispersion increases rapidly when δx increases. These “zero” dispersion areas correspond to the constant areas as shown in Fig. 4.
For low situation, the corresponding GVD values are also low, and flat bandwidth with “zero” dispersion is thus observed. For high situations (0.25a ≤ δx ≤ 0.27a), GVD values are separately shown in the inset picture of Fig. 6, for clarity.
Interestingly, each GVD curve also presents a negative GVD frequency range and a positive GVD frequency range. This reveals the possibility of dispersion compensating applications  using our proposed structure. The slow light bandwidth could be effectively enlarged by cascading two sections of waveguides with opposite GVD values for the same wavelength. These aspects will be explored in a future work.
To show the advantage of our structure in dispersion, the bandwidth corresponding to absolute GVD value limited within 2 ps.nm−1.mm−1 are checked. We should note that this value is extremely small if compared with other papers [13–15,25,26]. However, even in this strict consideration, our proposed structure shows very good low dispersion properties. In the situation 0.35a ≤ δx ≤ 0.5a, it can be seen that most of the GVD values are smaller than 0.5 ps.nm−1.mm−1. The obtained bandwidth for low GVD is much larger than previous study [11,14,15,25], e.g.: in the situation δx = 0.5a, our structure provides as large as more than 40 nm bandwidth (centered at 1550nm).
In conclusion, a new photonic crystal waveguide geometry is proposed to obtain slow waves in the largest possible bandwidth and lowest group velocity dispersion. The proposed geometry is based only on the shift of the waveguide bordering rows of holes. Thus, it is characterized by its simplicity when compared to previously proposed solutions [11,12,14–16,24]. As fast waves dominated by the index guiding mechanism are almost unaffected by the introduced shift while slow waves mainly dominated by gap guided mechanisms are, the structure produces unusual “U” type -frequency curves. Additionally, control of normalized delay-bandwidth product has been demonstrated (0.1363 for = 210 < NDPB < 0.3141 for = 10), and the obtained values are favorably compared with previous works. The GVD properties of the proposed structure have been also analyzed. It turns out that nearly-zero dispersion can be achieved within large bandwidth. Additionally, the new waveguide allows a possible control of positive or negative GVD values, opening opportunities for dispersion compensation devices.
References and links
1. Y. A. Vlasov, M. O’Boyle, H. F. Hamann, and S. J. McNab, “Active control of slow light on a chip with photonic crystal waveguides,” Nat. 438(7064), 65–69 (2005). [CrossRef]
2. R. Won, “Slow light now and then,” Nat. Photonics 2(8), 454–455 (2008). [CrossRef]
3. T. F. Krauss, “Why do we need slow light?” Nat. Photon. 2(8), 448–450 (2008). [CrossRef]
4. T. Baba, “Slow light in photonic crystals,” Nat. Photon. 2(8), 465–473 (2008). [CrossRef]
5. T. F. Krauss, “Slow light in photonic crystal waveguides,” J. Phys. D Appl. Phys. 40(9), 2666–2670 (2007). [CrossRef]
8. 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. 113903, 87 (2001). [CrossRef] [PubMed]
9. M. Loncar, T. Doll, J. Vuckovic, and A. Scherer, “Design and fabrication of silicon photonic crystal optical waveguides,” J. Lightwave Technol. 18(10), 1402–1411 (2000). [CrossRef]
10. E. Dulkeith, F. N. Xia, L. Schares, W. M. J. Green, and Y. A. Vlasov, “Group index and group velocity dispersion in silicon-on-insulator photonic wires,” Opt. Express 14(9), 3853–3863 (2006). [CrossRef] [PubMed]
11. 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]
12. 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]
14. J. Hou, D. Gao, H. Wu, and R. Hao, “Flat Band Slow Light in Symmetric Line Defect Photonic Crystals Waveguide,” IEEE Photon. Tech. Lett.21, 1571-1573 (2009). [CrossRef]
15. L. Dai and C. Jiang, “Photonic crystal slow light waveguides with large delay-bandwidth product,” Appl. Phys. B 95(1), 105–111 (2009). [CrossRef]
16. T. Baba, T. Kawaaski, H. Sasaki, J. Adachi, and D. Mori, “Large delay-bandwidth product and tuning of slow light pulse in photonic crystal coupled waveguide,” Opt. Express 16(12), 9245–9253 (2008). [CrossRef] [PubMed]
19. D. Marris-Morini, E. Cassan, D. Bernier, G. Maire, and L. Vivien, “Ultracompact tapers for light coupling into two-dimensional slab photonic-crystal waveguides in the slow light regime,” Opt. Eng. 47(1), 014602 (2008). [CrossRef]
21. C. Martijn de Sterke, K. B. Dossou, T. P. White, L. C. Botten, and R. C. McPhedran, “Efficient coupling into slow light photonic crystal waveguide without transition region: role of evanescent modes,” Opt. Express 17(20), 17338–17343 (2009). [CrossRef] [PubMed]
22. D. Mori, S. Kubo, H. Sasaki, and T. Baba, “Experimental demonstration of wideband dispersion-compensated slow light by a chirped photonic crystal directional coupler,” Opt. Express 15(9), 5264–5270 (2007). [CrossRef] [PubMed]
23. J. Jágerská, N. Le Thomas, V. Zabelin, R. Houdré, W. Bogaerts, P. Dumon, and R. Baets, “Experimental observation of slow mode dispersion in photonic crystal coupled-cavity waveguides,” Opt. Lett. 34(3), 359–361 (2009). [CrossRef] [PubMed]
24. 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]
25. 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]
26. R. J. P. Engelen, Y. Sugimoto, Y. Watanabe, J. P. Korterik, N. Ikeda, N. F. van Hulst, K. Asakawa, and L. Kuipers, “The effect of higher-order dispersion on slow light propagation in photonic crystal waveguides,” Opt. Express 14(4), 1658–1672 (2006). [CrossRef] [PubMed]