## Abstract

Analysis of photonic crystal coupled resonator optical waveguide (CROW) structures with a highly dispersive background medium is presented. A finite-difference time-domain algorithm was employed which contains an exact representation of the permittivity of a three-level atomic system which exhibits electromagnetically induced transparency (EIT). We find that the coupling strength between nearest-neighbor cavities in the CROW decreases with increasing steepness of the background dispersion, which is continuously tunable as it is directly related to the control field Rabi frequency. The weaker coupling decreases the speed of pulse propagation through the waveguide. In addition, due to the dispersive nature of the EIT background, the CROW band shape is tuned around a fixed k-point. Thus, the EIT background enables dynamic tunability of the CROW band shape and the group velocity in the structure at a fixed operating point in momentum space.

© 2007 Optical Society of America

## 1. Introduction

Signal routing on integrated optics platforms can be achieved through total internal reflection [1], Bragg reflection [2], and coupled optical resonators [3]. In the latter, a series of high-*Q* micro-resonators are aligned within proximity to one another such that photons can tunnel from one resonator to its nearest-neighbor and thus propagate as if in a waveguide. Figure 1 shows a schematic of such a coupled-resonator optical waveguide (CROW) formed by an infinite straight line of defect cavities in a two-dimensional (2D) photonic crystal (PhC) [4, 5, 6]. Yariv et al. showed through a tight-binding (TB) analysis that the dispersion relation and propagation velocity of a pulse in such a structure are determined solely by a coupling factor, *κ* :

and

respectively, where *ω _{res}* is the mode frequency of a single, isolated resonator of the waveguide,

*R*is the spacing between resonators, and

*k*is the wave vector. Therefore, the group velocity, Eq. (2), of the pulse can be decreased by decreasing

*κ*. However,

*κ*is typically a fixed quantity based on the design of the structure, e.g. the intrinsic

*Q*of a the individual resonators and/or the spacing

*R*. Yet, for realizing active components in an optical circuit, a dynamic control over

*v*

*is desirable. The problem then becomes how one can manipulate κ in a seemingly fixed system. We propose that a tunable, highly dispersive process such as electromagnetically induced transparency (EIT) [7] will facilitate control over the propagation velocity in the CROW structure by modifying the coupling factor between cavities.*

_{g}EIT has been demonstrated in multilevel atomic systems [7, 8]. In a three-level Λ system, three atomic energy levels are situated such that the transition probability amplitudes from the ground levels *b* and *c* to the upper level *a* show destructive quantum interference under the two-photon resonance condition. Phenomenologically, a transparency window opens between two absorption peaks in the susceptibility, which allows a probe pulse to propagate through the otherwise optically opaque medium. At the same time, this transparency window occurs at a highly dispersive region of the susceptibility. It has been reported that such a lossless, highly dispersive EIT medium within an optical resonator cavity has the effect of dramatically narrowing the cavity linewidth and increasing the cavity *Q*-factor [9, 10, 11]. In the following work, we consider a PhC CROW structure (Figure 1) which possesses a background exhibiting EIT. Thus, the system consists of a series of optical resonators embedded in a highly dispersive medium whose dispersive properties are linked to the Rabi frequency of the control field excitation. We show that control over the dispersion relation of the waveguide is gained through the control field Rabi frequency, which in turn enables control of *v*
* _{g}* in the waveguide.

We studied the proposed CROW system through two-dimensional (2D) finite-difference time-domain (FDTD) calculations capable of modeling EIT phenomena [12]. The main advantage of our implementation is the direct representation of the EIT permittivity [13, 14] as a complex multipole expansion which is readily incorporated into the FDTD algorithm via the auxiliary differential equation (ADE) method [15, 16, 17]. The ADE method uses the link between the frequency dependent permittivity *ε* and the electric flux density, *D*⃗, via the polarization, *P*⃗* _{l}*. This link introduces a phasor polarization current,

*J*⃗

*(*

_{l}*t*)=

*∂P*⃗

*/*

_{l}*∂t*, from which time-step update-equations can be derived. Since it adds only one extra step to the FDTD update algorithm, the ADE method does not increase the computation load significantly.

#### 1.1. Electromagnetically induced transparency

Figure 2(a) is a schematic of the three-level EIT system considered in this work (the Λ scheme). The susceptibility such a system can be expressed by [13, 14]:

where Δ=*ω*-*ω _{ab}* is the probe beam detuning from the resonance frequency

*ω*of the

_{ab}*ab*transition, γ

*is the loss due to the*

_{ab}*ab*-transition, Ω

*is the Rabi frequency of the control field, γ*

_{c}*is the decoherence rate, and C is a materials constant related to the density of excitation centers and the dipole moment. We express the permittivity,*

_{bc}*ε*=

*ε*

_{0}(1+

*χ*), of the EIT medium using an complex multipole expansion of Eq. (3):

where *ε _{b}* is the background dielectric constant and

*A*

_{1,2}and

*B*

_{1,2}are the coefficients:

and

Figure 2(b) is a representative plot of Eq. (4) assuming *ε _{b}*=1.0,

*ω*=0.35,

_{ab}*γ*=0.002, C=0.01s

_{ab}^{-1},

*γ*=0, and the two Rabi frequency values 0.08 (solid line) and 0.03 (dashed line) (all units are normalized by the lattice constant,

_{bc}*a*, of the PhC). Adjusting Ω

*has two effects on the optical properties of the EIT medium. First, the frequency width of the transparency window increases with Ω*

_{c}*since the separation between the absorption peaks is equal to Ω*

_{c}*. Secondly, the slope of the dispersion, ℜ(*

_{c}*ε*), in the transparency region is inversely related to Ω

*. Thus,*

_{c}*v*of a pulse propagating in the EIT medium decreases as Ω

_{g}*decreases since*

_{c}*v*=d

_{g}*ω*/d

*k*, ie. the steeper the dispersion, the slower the group velocity.

#### 1.2. Coupled resonator optical waveguide structure

The CROW of this study consists of an infinitely long line of high-*Q*, point defect PhC cavities separated by a distance of *R*=3*a* (see Fig. 1). The cavities are formed in 2D square lattice of high-dielectric rods, *ε*=11.56, embedded in an EIT medium. The EIT medium has the same physical parameters as those given in section 1.1 except that *ω _{ab}* is chosen so that the EIT resonance frequency matches the mid-band frequency of the vacuum CROW,

*ω*

_{0}=0.3796, the subscript denoting the vacuum background structure. The size of the computational domain is 3

*a*× 14

*a*and contains and array of circular rods with diameter

*d*=0.4

*a*arranged in a square lattice. The point defect cavities are formed by removing a rod in the center of the array. Periodic boundaries are used along the long sides of the domain to form an infinitely long series of defect cavities separated by the distance

*R*=3

*a*. Perfectly matched layer [18] absorbing boundaries on the short sides of the domain truncate the structure in the transverse direction. A spatial resolution of Δ

*x*=Δ

*y*=0.025

*a*is used in all of the calculations. A transverse-magnetic polarized field (magnetic field components normal to the rod axis) is used in the calculations since a photonic band gap exists for this polarization, and the fields are randomized throughout the domain at the beginning of the calculation. Data collected from a point detector positioned near the center of the point defect is used to calculate the resonant frequency of the CROW using the Padé approximation and Baker’s algorithm [19, 20].

## 2. Results and discussion

We investigate the effects of the EIT medium on the optical response of the CROW by calculating the dispersion relationship of the CROW for different values for Ω* _{c}*. For the vacuum background CROW, the frequency width Δ

*ω*of the waveguide band is 0.0109, and the maximum group velocity

*v*

_{g,0}=0.1025

*v*/

_{g}*c*

_{0}of the band occurs near the zone center. We can calculate the coupling factor from

*κ*=Δ

*ω*/(2

*ω*) if the resonant frequency of the isolated resonator

_{res}*ω*is known [6]. A separate calculation reveals that for the vacuum background single-defect cavity

_{res}*ω*=0.3793, thus giving a coupling factor of

_{res}*κ*

_{0}=0.0144. For the EIT background isolated resonator, the value of

*ω*changes for each value of Ω

_{res}*because of a frequency pulling effect [21].*

_{c}For convenience, we divide the response of the EIT CROW into two regimes: case 1) when the frequency width of the EIT transparency exceeds the Δ*ω* of the CROW and case 2) when the opposite is true.

Figures 3(a) and (b) show the CROW bands and group velocities as a function of wavevector for case 1. The vacuum background band is included as a comparison to the EIT-CROW bands under different values of Ω* _{c}*. The points are the results from FDTD calculations while the solid lines are the results given by Eqs. (1) and (2), after calculating

*κ*from the FDTD data, as done previously for the vacuum CROW. As shown, the bands become flatter as Ω

*is decreased (e.g. as the dispersion of the EIT media increases). In addition, the bands share a common point at the zone center at*

_{c}*ω*

_{0}=

*ω*around which the bands pivot as Ω

_{ab}*is tuned. When Ω*

_{c}*=0.01, slightly less than Δ*

_{c}*ω*of the vacuum CROW, the band is very flat with a maximum group velocity less than 2.0×10

^{-3}

*v*/

_{g}*c*

_{0}and a Δ

*ω*<2.0×10

^{-4}

*a*/

*λ*. This represents over a 50 times slowing of the group velocity in comparison to the vacuum background CROW.

Figures 4 (a) and (b) present the bands and group velocities for the structure under case 2, respectively. A single band from just above the transition between case 1 and 2 is included for comparison (Ω* _{c}*=0.02). As shown, the bands in case 2 are significantly flatter than in case 1. For Ω

*=0.003, the group velocity is below 1.8×10*

_{c}^{-4}

*v*/

_{g}*c*

_{0}, approximately 570 times slower than that of the vacuum CROW. For values of Ω

*≤ 0.002, the absorption of the EIT media (recall*

_{c}*γ*=0.002 in our calculations) would reduce the

_{ab}*Q*-factor of the cavities, thus setting a lower limit for the group velocity. (However, we note that this limit occurs only for a steady-state control field. If Ω

*were ramped down adiabatically to zero as the pulse propagates in a finite CROW structure, the pulse would be halted and stored in the EIT [22].) In addition to the slower group velocity, the rate of change in the slope of the bands with change in Ω*

_{c}*decreases in case 2 in comparison to case 1. This can be seen more clearly in a plot of the group velocity in the EIT CROW as a function of Ω*

_{c}*. Since this quantity is directly related to the coupling factor by Eq. (2), we present*

_{c}*κ*(Ω

*) in Figure 5.*

_{c}On the right axis in Figure 5, we plot the *v _{g}* vs. Ω

*relationship of the EIT background [23]. Working from the group velocity at the*

_{c}*ab*-transition frequency we can derive an expression for

*v*(Ω

_{g}*):*

_{c}where *α* is (2C*ω _{ab}*)

^{-1}. From the values given earlier, we calculate

*α*=1.317×10

^{2}for the EIT background medium and plot Eq. (7) in Figure 5. We find that the relationship of the EIT CROW

*κ*vs. Ω

*follows that of the EIT background. We used a minimization algorithm to fit the FDTD data points to an equation of the form of Eq. (7) to find*

_{c}*α*=1.940×10

^{2}. The shape of

*κ*(Ω

*) can be explained as follows: Since the coupling is dependent upon the overlap of the field between adjacent resonators,*

_{c}*κ*∝1/

*Q*. At the same time, as Soljačić et al. have shown,

*Q*increases as the steepness of the dispersion increases, thus

*v*∝1/

_{g}*Q*. It follows that

*κ*∝

*v*, meaning that in the EIT CROW system, manipulation of the coupling coefficient is directly coupled to a controllable parameter, namely Ω

_{g}*.*

_{c}## 3. Conclusion

In conclusion, we have successfully modeled a PC CROW structure embedded in a three-level EIT medium with the FDTD method. The ADE method was employed to incorporate the frequency dependence of *ε* of the highly dispersive EIT background into the time-stepping algorithm. We found that by adjusting the control field strength, dynamic control the shape of the CROW band was enabled, whereas the shape of the CROW band is typically fixed since it is controlled by the coupling strength of the cavities, eg. distance *R* between the defects. We classified two regimes of the dynamic system by the relative frequency band width of the vacuum CROW structure verses the band width of the EIT transparency window. In the two regimes, the rate of change in the coupling factor (or the group velocity) as a function of the control field Rabi frequency differed, being larger in the case where Ω* _{c}* ≥ Δ

*ω*. We also found that the relationship of the EIT-CROW coupling factor,

*κ*(Ω

*) follows that of the EIT group velocity,*

_{c}*v*(Ω

_{g}*). In the spirit of the TB analysis, one can draw a parallel between the atomic separation distance and the resonator cavity separation in the CROW. Therefore, one can also think of the dispersion of the media or Ω*

_{c}*in the cavity as the ‘atomic potential’. As Ω*

_{c}*decreases, the atomic potential increases, reducing the coupling between modes. To achieve a usable slow light effect, a large bandwidth is required to accept the input pulse into the system before the band is flattened and the pulse stopped [24]. Through our investigation, we have shown an implementation of the CROW which would allow dynamic control of the propagation speed of a pulse through the system. We chose the 2D square lattice PhC system because it has a monopole point defect mode and simple band structure, making it a ‘clean’ system for this study; however, these effects are not limited to this configuration. More importantly, the unique dispersion characteristics of the EIT in the CROW fixes a pivot point in the CROW bands, around which the bands change shape. This pivot offers a fixed operation point in momentum space which could useful in coupling in and out of the structure. The EIT and CROWimplementation could be applied in optical buffer or optical memory schemes.*

_{c}## Acknowledgments

Funding for this work was through the Swedish Foundation for Strategic Research (SSF) under the INGVAR program, the SSF Strategic Research Center in Photonics, and the Swedish Research Council (VR).

## References and links

**1. **D. Marcuse, *Theory of Dielectric Optical Waveguides*,(Academic Press, New York, 1974).

**2. **P. Yeh and A. Yariv, “Bragg reflection waveguides,” Opt. Commun. **19**, 427–430 (1976). [CrossRef]

**3. **A. Yariv, Y. Xu, R. K. Lee, and A. Scherer, “Coupled-resonator optical waveguide: a proposal and analysis,” Opt. Lett. **24**, 711–713 (1999). [CrossRef]

**4. **E. Yablonovitch, “Inhibited Spontaneous Emission in Solid-State Physics and Electronics,” Phys. Rev. Lett. **58**, 2059–2062 (1987). [CrossRef] [PubMed]

**5. **S. L. McCall, P. M. Platzman, R. Dalichaouch, D. Smith, and S. Schultz, “Microwave propagation in two-dimensional dielectric lattices,” Phys. Rev. Lett. **67**, 2017–2020 (1991). [CrossRef] [PubMed]

**6. **M. Bayindir, B. Temelkuran, and E. Ozbay, “Tight-binding description of the coupled defect modes in three-dimensional photonic crystal,” Phys. Rev. Lett. **84**, 2140–2143 (2000). [CrossRef] [PubMed]

**7. **K.-J. Boller, A. Imamoglu, and S. E. Harris, “Observation of electromagnetically induced transparency,” Phys. Rev. Lett. **66**, 2593–2596 (1991). [CrossRef] [PubMed]

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

**9. **M. D. Lukin, M. Fleischhauer, M. O. Scully, and V. L. Velichansky, “Intracavity electromagnetically induced transparency,” Opt. Lett. **23**, 295–297 (1998). [CrossRef]

**10. **G. Müller, M. Müller, A. Wicht, R.-H. Rinkleff, and K. Danzmann, “Optical resonator with steep internal dispersion,” Phys. Rev. A **56**, 2385–2389 (1997). [CrossRef]

**11. **M. Soljačić, E. Lidorikis, J. D. Hau, and Joannopoulos, “Enhancement of microcavity lifetimes using highly dispersive materials,” Phys. Rev. E **71**, 026602 (2005). [CrossRef]

**12. **C. W. Neff, L. M. Andersson, and M. Qiu, “Modelling electromagnetically induced transparency media using the finite-difference time-domain method,” New J. Phys. **9**, 48 (2007). [CrossRef]

**13. **M. O. Scully and M. S. Zubairy, *Quantum Optics* (Cambridge University Press, Cambridge, England, 1997).

**14. **P. Jänes, J. Tidström, and L. Thylén, “Limits on optical pulse compression and delay bandwidth product in electromagnetically induced transparency media,” J. Lightwave Technol. **23**, 3893–3899 (2005). [CrossRef]

**15. **M. Okoniewski, M. Mrozowski, and M. Stuchly, “Simple treatment of multi-term dispersion in FDTD,” IEEE Micro. Guided Wave Lett. **7**, 121–123 (1997). [CrossRef]

**16. **A. Taflove and S. C. Hagness, *Computational electrodynamics: the finite-difference time-domain method*, 2nd ed. (Artech House, Boston, 2000).

**17. **M. Han, R. Dutton, and S. Fan, “Model dispersive media in finite-difference time-domain method with complex-conjugate pole-residue pairs,” Microwave andWireless Components Letters, IEEE [see also IEEE Micro. Guided Wave Lett.] **16**, 119–121 (2006).

**18. **J.-P. Berenger, “A perfectly matched layer for the absorption of electromagnetic waves,” J. Comp. Phys. **114**, 185–200 (1994). [CrossRef]

**19. **W.-H. Guo, W.-J. Li, and Y.-Z. Huang, “Computation of resonant frequencies and quality factors of cavities by FDTD technique and Padé approximation,” Microwave and Wireless Components Letters, IEEE [see also IEEE Micro. Guided Wave Lett.] **11**, 223–225 (2001).

**20. **G. A. Baker and J. L. Gammel, *The Padé Approximant in Theoretical Physics*, (Academic, New York, 1970).

**21. **
manuscript in preparation

**22. **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 (London) **409**, 490–493 (2001). [CrossRef]

**23. **J. Tidström, P. Jönes, and L. M. Andersson, “Delay bandwidth product of electromagnetically induced transparency media,” Phys. Rev. A **75**, 53803 (2007). [CrossRef]

**24. **M. F. Yanik and S. Fan, “Stopping light all optically,” Phys. Rev. Lett. **92**, 083901 (2004). [CrossRef] [PubMed]