Coupled resonator optical waveguide structures with highly dispersive media Curtis W. Neff, L. Mauritz Andersson and Min Qiu Laboratory of Optics, Photonics and Quantum Electronics, Department of Microelectronics and Applied Physics, Royal Institute of Technology, 164 40 Kista, Sweden
[email protected],
[email protected]
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 nearestneighbor 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 OCIS codes: (230.4555) Coupled resonators; (130.5296) Photonic crystal waveguides; (190.5530) Pulse propagation and solitons
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). 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). 4. E. Yablonovitch, “Inhibited Spontaneous Emission in Solid-State Physics and Electronics,” Phys. Rev. Lett. 58, 2059–2062 (1987). 5. S. L. McCall, P. M. Platzman, R. Dalichaouch, D. Smith, S. Schultz, “Microwave propagation in two-dimensional dielectric lattices,” Phys. Rev. Lett. 67, 2017–2020 (1991). 6. M. Bayindir, B. Temelkuran, and E. Ozbay, “Tight-binding description of the coupled defect modes in threedimensional photonic crystal,” Phys. Rev. Lett. 84, 2140–2143 (2000). 7. K.-J. Boller, A. Imamoglu, and S. E. Harris, “Observation of electromagnetically induced transparency,” Phys. Rev. Lett. 66, 2593–2596 (1991). 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). 9. M. D. Lukin, M. Fleischhauer, M. O. Scully, and V. L. Velichansky, “Intracavity electromagnetically induced transparency,” Opt. Lett. 23, 295–297 (1998). 10. G. M¨uller, M. M¨uller, A. Wicht, R.-H. Rinkleff, and K. Danzmann, “Optical resonator with steep internal dispersion,” Phys. Rev. A 56, 2385–2389 (1997). 11. M. Soljaˇci´c, E. Lidorikis, L. V. Hau, and J. D. Joannopoulos, “Enhancement of microcavity lifetimes using highly dispersive materials,” Phys. Rev. E 71, 026602 (2005). 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). 13. M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University Press, Cambridge, England, 1997).
#84643 - $15.00 USD
(C) 2007 OSA
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10362
14. P. J¨anes, J. Tidstr¨om, and L. Thyl´en, “Limits on optical pulse compression and delay bandwidth product in electromagnetically induced transparency media,” J. Lightwave Technol. 23, 3893–3899 (2005). 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). 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 complexconjugate pole-residue pairs,” Microwave and Wireless 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). 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´e 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´e 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). 23. J. Tidstr¨om, P. J¨anes, and L. M. Andersson, “Delay bandwidth product of electromagnetically induced transparency media,” Phys. Rev. A 75, 53803 (2007). 24. M. F. Yanik and S. Fan, “Stopping light all optically,” Phys. Rev. Lett. 92, 083901 (2004).
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
ω (k) = ωres [1 + κ cos(kR)]
(1)
vg (k) = −ωres Rκ sin(kR),
(2)
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 g 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. 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
#84643 - $15.00 USD
(C) 2007 OSA
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10363
Fig. 1. Cross-section schematic of a CROW structure formed in a 2D PhC square lattice of high index rods where single rods were removed to create the resonant cavities.
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, Pl . This link introduces a phasor polarization current, Jl (t) = ∂ Pl /∂ 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]:
χ (Δ) = −C
Δ + iγbc , (Δ + iγbc )(Δ + iγab ) − Ω2c /4
(3)
where Δ = ω − ωab is the probe beam detuning from the resonance frequency ω ab of the abtransition, γab is the loss due to the ab-transition, Ω c is the Rabi frequency of the control field, γbc 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, ε = ε 0 (1 + χ ), of the EIT medium using an complex multipole expansion of Eq. (3):
ε = εb +
#84643 - $15.00 USD
(C) 2007 OSA
A1 A2 + , B 1 + iω B 2 + iω
(4)
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10364
EIT Permittivity 4 3
Re{ε}
2 1 0
a
−1 −2 0.25
Εc
Εp
0.3 0.35 0.4 0.45 Normalized frequency (a/λ) Ωc=0.08
Ωc=0.03
6 5
c (a)
4 Im{ε}
b
3 2 1 0 0.25
0.3 0.35 0.4 Normalized frequency (a/λ)
0.45
(b) Fig. 2. (a) Schematic of the energy levels in a three level Λ system. Representative plot of the real and imaginary parts of the permittivity for an EIT medium at two Ωc values of the control field: 0.08 (solid) and 0.03 (dashed)
where εb is the background dielectric constant and A 1,2 and B1,2 are the coefficients: iC γbc − γab −1∓ A1,2 = 1/2 , 2 (γbc − γab )2 − Ω2
(5)
c
and
1/2 1 . (6) γbc + γab ∓ (γbc − γab )2 − Ω2c 2 Figure 2(b) is a representative plot of Eq. (4) assuming ε b = 1.0, ωab = 0.35, γab = 0.002, C = 0.01s−1 , γbc = 0, and the two Rabi frequency values 0.08 (solid line) and 0.03 (dashed line) (all units are normalized by the lattice constant, a, of the PhC). Adjusting Ω c 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, ℜ(ε ), in the transparency region is inversely related to Ω c . Thus, vg of a pulse propagating in the EIT medium decreases as Ω c decreases since vg = dω /dk, ie. the steeper the dispersion, the slower the group velocity. B1,2 = −iωab +
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 = 3a (see Fig. 1). The cavities are formed in 2D square lattice of #84643 - $15.00 USD
(C) 2007 OSA
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10365
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 3a × 14a and contains and array of circular rods with diameter d = 0.4a 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 = 3a. 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.025a 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´e 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.1025v g/c0 of the band occurs near the zone center. We can calculate the coupling factor from κ = Δω /(2ω res ) 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 κ 0 = 0.0144. For the EIT background isolated resonator, the value of ω res changes for each value of Ω c because of a frequency pulling effect [21]. 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 Ω c 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 ω 0 = ωab around which the bands pivot as Ω c is tuned. When Ωc = 0.01, slightly less than Δω of the vacuum CROW, the band is very flat with a maximum group velocity less than 2.0 × 10 −3vg /c0 and a Δω < 2.0 × 10 −4a/λ . 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 Ωc = 0.003, the group velocity is below 1.8 × 10 −4vg /c0 , approximately 570 times slower than that of the vacuum CROW. For values of Ω c ≤ 0.002, the absorption of the EIT media (recall γab = 0.002 in our calculations) would reduce the Q-factor of the cavities, thus setting a lower limit for the group velocity. (However, we note that this limit occurs only for a steadystate control field. If Ω c 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 #84643 - $15.00 USD
(C) 2007 OSA
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10366
0.386
0.11
(a)
c
g 0
Group velocity (v /c )
Ω =0.05 c
0.382 0.38
Ωc=0.03 Ωc=0.02 Ω =0.01
0.378
c
(b)
Vacuum CROW
0.1
Ω =0.1 c Ω =0.08
0.384 Normalized Frequency ( a/λ)
Vacuum CROW
0.08 Ωc=0.1 0.06 Ω =0.08 c 0.04
Ω =0.05 c
0.376 0.02
Ω =0.03 c Ω =0.02
0.374
c
0
0.1 0.2 Wavevector (π/a)
0 0
0.3
0.1 0.2 Wavevector (π/a)
0.3
Fig. 3. Evolution of EIT CROW bands (a) and group velocities (b) for case 1: Ωc ≥ Δω .
−3
0.38
8 Ωc=0.02
7
0.3798
Ω =0.01
0.3797
0.3795
0.3794
c
Ω =0.008 c Ω =0.005 c Ωc=0.003
Group velocity (vg/c0)
Normalized Frequency ( a/ λ)
(a)
x 10
c
6 5 4 Ωc=0.003
3 2
0.3792
(b)
Ω =0.02
Ωc=0.01
Ω =0.008 c
1
Ω =0.005 c
0.3791 0
0.1 0.2 Wavevector (π/a)
0.3
0 0
0.1 0.2 Wavevector (π/a)
0.3
Fig. 4. Evolution of EIT CROW bands (a) and group velocities (b) for case 2: Ωc ≤ Δω .
#84643 - $15.00 USD
(C) 2007 OSA
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10367
1
0.8
CROW EIT Fit EIT group velocity (vg/c0)
Norm. coupling coeff. (κ/κ0)
0.9
0.7 0.6 0.5 0.4 0.3 0.2 0.1 0
−2
−1
10
10
log Ω
c
Fig. 5. Coupling coefficient (normalized by the vacuum CROW κ0 ) as a function of control field Rabi frequency. The points and crosses are taken from the FDTD and TB analysis data, respectively, and the solid line is a plot of the EIT group velocity relationship (right axis) with a fitted α parameter.
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. On the right axis in Figure 5, we plot the v g vs. Ωc relationship of the EIT background [23]. Working from the group velocity at the ab-transition frequency we can derive an expression for vg (Ωc ): c0 c0 c0 α Ω2c vg = = = , (7) ω 4C 1 + α Ω2c 1 + 2ab 2 1 + ωab ∂ χ 2 ∂Δ
Ω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. Ωc 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 α = 1.940 × 10 2. The shape of κ (Ωc ) can be explained as follows: Since the coupling is dependent upon the overlap of the field between adjacent resonators, κ ∝ 1/Q. At the same time, as Soljaˇci´c et al. have shown, Q increases as the steepness of the dispersion increases, thus v g ∝ 1/Q. It follows that κ ∝ v g , meaning that in the EIT CROW system, manipulation of the coupling coefficient is directly coupled to a controllable parameter, namely Ω 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 #84643 - $15.00 USD
(C) 2007 OSA
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10368
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, κ (Ω c ) follows that of the EIT group velocity, vg (Ωc ). 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 CROW implementation could be applied in optical buffer or optical memory schemes. 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).
#84643 - $15.00 USD
(C) 2007 OSA
Received 28 Jun 2007; revised 24 Jul 2007; accepted 27 Jul 2007; published 1 Aug 2007
6 August 2007 / Vol. 15, No. 16 / OPTICS EXPRESS 10369