## Abstract

The geometry of photonic crystal waveguides with ring-shaped holes is optimized to minimize dispersion in the slow light regime. We found geometries with a nearly constant group index in excess of 20 over a wavelength range of 8 nm. The origin of the low dispersion is related to the widening of the propagating mode close to the lower band gap edge.

©2007 Optical Society of America

## 1. Introduction

Photonic crystal waveguides (PhCW’s) have very different dispersion properties compared to conventional waveguides due to their periodic boundaries. In particular, PhCW’s exhibit slow group velocity near the Brillouin zone edge [1]–[5]. Reduced group velocity increases light-matter interaction and nonlinearities, which can be utilized to realize more compact integrated optics devices. However, the slow modes usually have a very high group velocity dispersion (GVD), which leads to optical signal degradation in telecommunications systems. Therefore it is important to realize slow-light structures with tailored dispersion properties. A number of approaches has been used to achieve this, for example using W2 waveguides [6] or modification of the hole radius [7] or period [8] in the first one or two hole row(s) next to the waveguide channel.

Our approach is to use photonic crystal waveguides with ring-shaped holes (RPhCW’s) [9]–[12], for which we have measured a group index up to 20[11]. Due to an extra free parameter in the lattice design, PhCs with ring-shaped holes allow fine-tuning of the waveguide dispersion properties. In this paper we optimize the ring dimensions so that W1 waveguides show minimal dispersion with a relatively high group index.

## 2. Optimization of the ring geometry

Figure 1 shows a RPhCW fabricated into silicon-on-insulator substrate with a 240 nm thick top silicon layer [11]. The waveguide is a single line defect of missing holes in a triangular lattice of rings with outer and inner radii *R _{out}* = 0.344

*a*and

*R*= 0.203

_{in}*a*, respectively.

2D simulations with the plane wave method [13] are used. The background refractive index is the effective refractive index of the guided TE polarized mode in the SOI slab. For the 240 nm thick silicon slab on silicon dioxide, which supports only a single TE mode at 1550 nm, we use an effective refractive index of 2.84.

The plane wave simulation gives the dispersion relation *u*(β), where *u* is the normalized frequency $\left(u=\frac{a}{\lambda}\right)$ and β is the propagation constant. The group velocity *v _{g}* is defined as

It can be obtained using the built-in function in the MIT Photonic Bands package, where *v _{g}* is calculated for the mode for each value of β by operating the magnetic field using the Hellmann-Feynman theorem [14]. The group index is

Figure 2 shows the calculated dispersion relation of the waveguide in Fig. 1. For comparison, the dispersion relation of the conventional PhC waveguide with the same air fill factor is shown as a scatter plot. The group index of the RPhCW is larger than that of the PhCW between 1610 nm and 1620 nm. The region between 1615 nm and 1620 nm is the most interesting; here the RPhCW has *n _{g}* > 25 and a smaller group velocity dispersion parameter

*D*than the conventional PhCW.

*D*is defined as

Figure 3 shows the shift of the band gap edge and the even waveguide mode frequency as a function of the ring outer radius, while keeping the ring width *R _{out}* –

*R*constant. The cut-off frequency increases with

_{in}*R*, but not as much as the band gap edge frequency, therefore the band gap edge comes closer to the mode. This affects the shape of the dispersion relation of the waveguide mode, particularly near the transition from the index-guided mode to band gap guided mode around $\beta =0.35\frac{2\pi}{a}.$ The dispersion diagram of the RPhCW with

_{out}*R*= 0.38

_{out}*a*shows a nearly constant slope between $\beta =0.35\frac{2\pi}{a}.$ and $\beta =0.45\frac{2\pi}{a}.$

Slow light in the corrugated waveguides is explained as an interaction between the forward and backward propagating modes. The interaction is at its strongest at the edge of the Brillouin zone, where the group velocity vanishes. When moving away from the zone boundary, the interaction gets weaker due to phase mismatch between the forward and backward propagating modes. However, as the mode of the RPhCW is at its closest to the bandgap edge near $\beta =0.35\frac{2\pi}{a}.$ it penetrates deeper into the surrounding photonic crystal lattice (Fig. 4). This increases the effect of the periodicity and compensates for the phase mismatch. This explains why the group index is nearly constant over a relatively large wavelength range.

Due to the large mode area, a supercell width of 35*a* is needed in the plane wave simulations in order to prevent the effect of the periodic boundary conditions. In order to verify the results, the group index was also deduced from the Fabry-Pérot oscillations in the RPhCW transmission spectrum simulated by the finite-difference time-domain (FDTD) method [15].

Figure 5 shows the dispersion relation and the group index of a W1 RPhCW with *R _{out}*=0.38

*a*and

*R*=0.24

_{in}*a*. The group index has a quasi constant value of 25 over a wavelength range of 8 nm.

Modes below the SiO_{2} light line are theoretically lossless. Therefore it is desirable to extend the part of the slow mode that is below the light line by decreasing its frequency. This can be realized by increasing the effective refractive index of the silicon slab, i.e., by increasing the silicon layer thickness *h*. With *h* = 400 nm (*n _{eff}*=3.18), we find that the slow light region is almost entirely below the light line. Figure 6 shows the band diagram and group index of a RPhCW with

*R*=0.385

_{out}*a*,

*R*=0.235

_{in}*a*and

*h*= 400 nm. One can see that both the bandgap and the mode frequency have dropped compared to Fig. 5. The group index is 37 ± 3 over a wavelength range of 8 nm.

As can be seen from the definition of the group velocity dispersion parameter *D*(Eq. 3), nearly constant group index regions exhibit small dispersion. Figure 7 shows the GVD curves calculated from the group index for the RPhCW’s in Figs. 5 and 6. Group velocity dispersion minima can be seen in the nearly constant group index region with *D* < 1ps/(mm ∙ nm) over 8 nm for the RPhCW of Fig. 5 and *D* < 1ps/(mm ∙ nm) over 3 nm for the RPhCW of Fig. 6. In a conventional PhCW, *D* increases monotonically when coming closer to cut-off wavelength.

## 3. Conclusion

We showed that photonic crystal waveguides with ring-shaped holes enable dispersion tailored slow light structures. We show two RPhCW’s exhibiting nearly constant group index regimes with *n _{g}* ≈ 25 and

*n*≈ 37, respectively, both over a wavelength range of 8 nm, with group velocity dispersion parameters

_{g}*D*< 1ps/(mm∙nm) over a few nanometers. Small dispersion is crucial in high bit rate telecommunications applications. In a wavelength division multiplexed telecommunications system with a channel spacing of 100 GHz (≈0.8 nm), the wavelength range of the useful slow light region achieved in this work is equivalent to about 10 channels.

## Acknowledgments

This work was financed by the Academy of Finland under the PHC-OPTICS project and the European Union under project PHAT.

## References and links

**01. **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**, 253902 (2001). [CrossRef] [PubMed]

**02. **X. Letartre, C. Seassal, C. Grillet, P. Rojo-Romeo, P. Viktorovitch, M. Le Vassor d’Yerville, D. Cassagne, and C. Jouanin, “Group velocity and propagation losses measurement in a single-line photonic-crystal waveguide on InP membranes,” Appl. Phys. Lett. **79**, 2312–2314 (2001). [CrossRef]

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

**04. **H. Gersen, T. J. Karle, R. J. P. Engelen, W. Bogaerts, J. P. Korterik, N. F. van Hulst, T. F. Krauss, and L. Kuipers, “Real-Space Observation of Ultraslow Light in Photonic Crystal Waveguides,” Phys. Rev. Lett. **94**, 073903 (2005). [CrossRef] [PubMed]

**05. **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,” Nature **438**, 65–69 (2005). [CrossRef] [PubMed]

**06. **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**, 219–226 (2007). [CrossRef] [PubMed]

**07. **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**, 9444–9450 (2006). [CrossRef] [PubMed]

**08. **A. Jafarpour, A. Adibi, Y. Xu, and R. K. Lee, “Mode dispersion in biperiodic photonic crystal waveguides,” Phys. Rev. B **68**, 233102 (2003). [CrossRef]

**09. **M. Mulot, A. Säynätjoki, S. Arpiainen, H. Lipsanen, and J. Ahopelto, “Photonic crystal slabs with ring-shaped holes in a triangular lattice,” proceedings of the 3rd European Symposium on Photonic Crystals (ESPC 2005), Barcelona, Spain, 3–7 July, 2005.

**10. **H. Kurt and D. S. Citrin, “Annular Photonic Crystals,” Opt. Express **13**, 10316–10326 (2005). [CrossRef] [PubMed]

**11. **M. Mulot, A. Säynätjoki, S. Arpiainen, H. Lipsanen, and J. Ahopelto, “Slow light propagation in photonic crystal waveguides with ring-shaped holes,” accepted for publication in J. Opt. A : Pure and Appl. Opt

**12. **A. Säynätjoki, M. Mulot, S. Arpiainen, J. Ahopelto, and H. Lipsanen, “Photonic crystals with ring-shaped holes,” PECS VII conference, Monterey, USA, 8-11 April, 2007.

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

**14. **R. P. Feynman, “Forces in molecules,” Phys. Rev. **56**, 340–343 (1939). [CrossRef]

**15. **M. Qiu, F2P software, http://www.imit.kth.se/info/FOFU/PC/F2P/