## Abstract

We present a detailed study of oscillating modes in a slab waveguide with air core and anisotropic metamaterial cladding. It is shown that, under specific dielectric configurations, slow and even stopped electromagnetic wave can be supported by such an air waveguide. We propose a linearly tapped waveguide structure that could lead the propagating light to a complete standstill. Both the theoretical analysis and the proposed waveguide have been validated by full-wave simulation based on finite-difference time-domain method.

© 2009 Optical Society of America

## 1. Introduction

Devices that can buffer or delay the light pulse are key components for both optical communication networks and quantum information processing systems. Large pulse delay (or ‘slow light’) is achievable with many methods/structures including electromagnetically induced transparency [1], quantum-dot semiconductor optical amplifiers [2], photonic crystal waveguide [3], direct coupled resonators [4], coherent population oscillations [5], stimulated Brillouin scattering [6], stimulated Raman scattering [7] and surface plasmon polaritons [8–10]. Recently, several groups proposed the slow light structures based on metamaterials [11–15]. K. L. Tsakmakidis *et al*. bring together the realms of metamaterials and slow light by studying a three-layer slab heterostructure composed of negative-index (‘left-handed material’) core and positive-index (‘right-handed material’) cladding [13]. The physics of slowing down the group velocity is been explained by negative Goos-Hänchen lateral displacements on the media interfaces. They discover that for different excitation frequencies the guided oscillatory fields stop at correspondingly different waveguide thicknesses in a linearly tapered waveguide, which is the so-called ‘trapped rainbow’ effect. However, fabricating an isotropic left-handed material (LHM) at optical frequencies is a significant challenge, while several anisotropic metamaterials have been successfully demonstrated by experiments [16, 17]. Moreover, the metamaterials used in the waveguide are always lossy and the attenuation of the slow light will be obvious due to the interaction between the guided light and the lossy core materials. In our previous work [18], we have proposed a planar air waveguide with anisotropic metamaterial (AMM) cladding to slow down and even trap the light. Similar structure based on negative-refractive-index photonic crystal was also proposed [19].

In this letter, we report our intensive investigation on a symmetric slab waveguide where an air core is surrounded by an AMM cladding. It is shown that the air waveguide could support different oscillating modes. By deriving closed-form expressions for the cutoff conditions of the oscillating modes, we prove that, with careful choice of the material parameters, there exist guiding modes for which extremely low or even zero group velocity can be attained. We also demonstrate by finite-difference time-domain (FDTD) method that a linearly tapped waveguide can be used to efficiently and coherently bring light pulses to a complete standstill.

## 2. Air waveguide design and mode analysis

The geometry of the waveguide is shown in Fig. 1, where an air slab of thickness *d* is surrounded by the AMM with partially negative permittivity and/or permeability components. To simplify the analysis, we assume the permittivity and permeability tensors that are
simultaneously diagonal as *ε̿* = *diag*[*ε _{x}*,

*ε*,

_{y}*ε*],

_{z}*μ̿*=

*diag*[

*μ*,

_{x}*μ*,

_{y}*μ*], with their principal axis paralleling the coordinate axis used in Fig. 1. We consider primarily

_{z}*p*-polarized (transverse magnetic, TM) waves, in which the magnetic field is along the

*y*axis. Similar analysis can be carried out for the transverse electric modes. From the Maxwell equations and the boundary conditions, we could obtain the dispersion equation that defines a set of allowed oscillating modes:

where the following ratios are introduced: *η* = *ε*
_{0}
*p*/*ε _{z}κ*,

*ip*and

*κ*standing for the transverse component of the wave vector in the AMM and the air core. Complex guided modes, i.e. modes with a complex propagation constant, have been observed in some completely lossless wave guiding structures [20, 21]. To giving a complete mode spectrum, we define the complex propagation constant along z direction

*k*=

_{z}*β*+

*iα*, and the effective index

*n*= (

_{eff}*β*+

*iα*)/

*k*

_{0}, here

*k*

_{0}representing the free-space wave number. Therefore

*p*and

*κ*can be obtained from the dispersion relations of the TM waves in air and AMM, i.e.,

*p*=

*k*

_{0}(

*n*

_{eff}^{2}

*ε*/

_{z}*ε*-

_{x}*ε*

_{z}*μ*/

_{y}*ε*

_{0}

*μ*

_{0})

^{1/2},

*κ*=

*k*

_{0}(1 -

*n*

_{eff}^{2})

^{1/2}, respectively. If

*α*≠ 0 ,

*p*and

*κ*are assumed to be complex numbers. For each mode and different material parameters, the dispersion diagrams,

*n*versus the reduced slab thickness

_{eff}*k*

_{0}

*d*, may be computed by numerical method.

For wave propagation in the air waveguide, the *z* component of the Poynting vector is:

Integrating *S _{z}* over the waveguide’s cross section, we obtain the power flow

*P*in different layers and the total power flow

_{i}*P*= Σ

_{tot}^{3}

_{i=1}

*P*. It is inferred from Eq. (2) that the net power flow can become negative in the AMM cladding when

_{i}*ε*< 0, but remains positive in the air core. If

_{x}*P*< 0, the mode is a backward mode since the total energy flow is opposite to the phase velocity. Otherwise, the mode is a forward mode. Under specific conditions the energy flow in the AMM cancels out that in the air core, i.e.,

_{tot}*P*= 0. One can prove that the group velocity reduce to zero simultaneously. Although the above equation of the total power flow is more complicated when complex guided mode appears, the result turns out to be simply zero [21].

_{tot}It is clear from the dispersion relations that TM mode solutions are closely related with the values of *ε _{x}*,

*ε*and

_{z}*μ*. For oscillating modes with real propagation constant,

_{y}*p*and

*κ*should be real, so that the parameters should be restricted by 0 <

*n*< 1 and

_{eff}*n*

_{eff}^{2}

*ε*/

_{z}*ε*>

_{x}*ε*

_{z}*μ*/

_{y}*ε*

_{0}

*μ*

_{0}. Table 1 summarizes the six possible sets of material parameters with which TM oscillating modes could exist. The complete dispersion diagrams belong to different sets, including ordinary and complex guided modes, are illustrated in Fig. 2. It shows that the waveguide could present various wave propagating properties under different cases. The notation “TM” in Fig. 2 identifies the polarization, followed by a subscript to track the number of field nodes in the core region and a superscript “

*f*” or “

*b*” to designate whether the mode is forward or backward, i.e., TM

_{m+1}

^{f}for the (

*m*+1)

^{th}-order forward mode. The insets in Fig. 2 show the corresponding field patterns for different modes.

When *ε _{x}* < 0 (Fig. 2(c), 2(e), 2(f)), one may notice that some of the modes have forward wave branch, as well as backward wave branch with antiparallel phase (

*v*) and group (

_{p}*v*) velocities. The two branches merge at a mode-degeneracy point, corresponding to a ‘critical’ reduced slab thickness with zero group velocity. Due to the difficult to achieve magnetic response at optical range, we restrict on non-magnetic AMM with permeability being 1. Therefore, we focus on Set E. Corresponding to the parameters used in Fig. 2(e), Fig. 3(a) and 3(b) denote the variation of the group velocity and the normalized power flow

_{g}*P*=

*P*/(|

_{tot}*P*

_{1}|+|

*P*

^{2}|+|

*P*

_{3}|) with reduced slab thickness. It is clear that the total power

*P*vanishes and the group velocity reduces to zero at the corresponding degeneracy points. For complex guided modes, the total power flow and the group velocity are both zero, as we predicted previously.

_{tot}## 3. Slow light waveguide analysis and simulation

The unique properties of the propagating modes can be used to slow down and even trap light. In Fig. 3(b), the group velocity decreases while the slab thickness reduces to the critical thickness for TM_{3} or TM_{4} mode at certain excitation frequency. Therefore, we design a linearly tapered waveguide to stop the wave, as illustrated in Fig. 4(a), with the same material parameters used in Fig. 3. Here we choose TM_{3}
^{f} mode as the operating mode because single-mode propagating can be excited while the air core thickness is correctly designed. FDTD method is employed to simulate the light propagating in the structure. Figure 4(b) demonstrates snapshots of the propagation of a monochromatic Gaussian modulated pulse with p-polarized magnetic-field component (with wavelength of 600nm) in the structure at different moments. The light pulse enters the air waveguide from the wide end (the width of air core is 1.4*μ*m). The width of the other end is 200nm, which is less than the ‘critical’ thickness (435nm). The AMM in our simulation is modeled by a lossless, dispersive material obeying Drude model with *ε _{x}* = -0.1,

*ε*= -1 and

_{z}*μ*= 1 at 600nm wavelength. While the guided oscillating fields propagate along the tapered waveguide, the group velocity progressively decreases. At a certain moment, the front part of the light pulse always travels slower than the rear part; therefore the pulse exhibits a spatial compression due to the reduction in group velocity while propagating along the structure, as demonstrated in Fig. 4(b). Eventually, the group velocity becomes zero at the ‘critical’, pre-determined, waveguide thickness (around 80T in Fig. 4(b)). At that position, the pulse is extremely compressed into a small region and stays there. Meanwhile, the amplitudes of the field components increase to about 4 times of the input values. The stopped light pulse should stay there for a long time. However, due to the variation of the core thickness

_{y}*d*with distance

*z*is not adiabatic, we can not expect to observe the splitting of the pulse in its frequency constituents and “stopping” of each frequency at a different location [21]. The amplitude of the trapped light will progressively decrease since it excites a backward-propagating forward-mode after staying for a while and the energy gradually spread out to the input end and fade out. If we try to totally stop a light, a very long adiabatically tapering structure should be used, which means that it is not easy to realize a perfect trap of the light. It should be noticed that this is totally different from the reflection by a cut-off ordinary waveguide since the light spreads out not in a pulse form as what is reflected in a cut-off case.

From previous analysis, we believe that for different excitation frequencies the guided lights will stop correspondingly at different guide thicknesses in the trapper. FDTD simulations are performed to verify this analytical result (Fig. 5). It shows that the full-wave simulation agrees with the analytical result very well. The propagating light having higher frequency stops at the position with smaller waveguide thickness. If the incident light is a white light, it can be altogether trapped within a fixed area, spanning a continuous range of waveguide thicknesses.

Considering practical realization, we propose a layered medium (shown in Fig. 6), consisting of alternating silicon dioxide layer (*ε*
_{1} = 3.9) and silver layer (*ε*
_{2} < 0 at optical frequency), which can be effectively mimicked as an anisotropic metamaterial [15, 23]. The effective dielectric permittivity tensor of such structure (with *N* as the volume fraction of silver) is given by [24]

$${\epsilon}_{z}=\genfrac{}{}{0.1ex}{}{{\epsilon}_{1}{\epsilon}_{2}}{N{\epsilon}_{1}\text{}+\left(1-N\right){\epsilon}_{2}}.$$

At optical frequency silver behaves as a plasmonic material with *ε*
_{2} < 0. According to the parameter data in Ref. 25, we have for example, *ε _{x}* = -0.125,

*ε*= -3.31 with the silver volume fraction

_{z}*N*= 70% at 351nm. Thus by fabricating a light trapper with such multilayered composite material, we could probably stop the 351nm light at the position when the waveguide thickness reaches 292nm.

## 4. Conclusion

In summary, we have investigated the guided modes in a planar air waveguide with AMM cladding. It is shown that both forward and backward waves can propagate under certain conditions, differing considerably from those of a conventional dielectric waveguide. The propagating electromagnetic waves can be stopped at the critical point where the forward and backward modes degenerate together, as a result of the power flows in the two media totally cancelling out. A linearly tapered air waveguide is introduced to bring the propagating wave to a complete standstill. We believe that such air waveguide have two distinguished advantages, i.e., less influence of material losses on the guiding light, and easier fabrication with layered metal/dielectric structures. Therefore, the proposed air waveguide is more suitable for most conceived applications.

## Acknowledgments

This work is supported by the National Basic Research Program of China (No. 2004CB719800) and the National Nature Science Foundations (No. 60671002 and No. 60801001).

## References and links

**1. **M. D. Lukin and A. Imamoglu, “Controlling photons using electromagnetically induced transparency,” Nature **413**, 273–4276 (2001). [CrossRef]

**2. **E. Gehrig, M. van der Poel, J. Mørk, and O. Hess, “Dynamic Spatiotemporal Speed Control of Ultrashort Pulses in Quantum-Dot SOAs,” IEEE J. Quant. Electron. **42**, 1047–1054 (2006). [CrossRef]

**3. **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]

**4. **A. Melloni, F. Morichetti, and M. Martinelli, “Optical slow wave structures,” Opt. Photon. News **14**, 44 (2003). [CrossRef]

**5. **M. S. Bigelow, N. N. Lepeshkin, and R. W. Boyd, “Superluminal and Slow Light Propagation in a Room-Temperature Solid,” Science **301**, 200–202 (2003). [CrossRef]

**6. **Y. Okawachi, M. S. Bigelow, J. E. Sharping, Z. Zhu, A. Schweinsberg, D. J. Gauthier, R. W. Boyd, and A. L. Gaeta, “Tunable All-Optical Delays via Brillouin Slow Light in an Optical Fiber,” Phys. Rev. Lett. **94**, 153902 (2005). [CrossRef]

**7. **K. Lee and N. M. Lawandy, “Optically induced pulse delay in a solid-state Raman amplifier,” Appl. Phys. Lett. **78**, 703–705 (2001). [CrossRef]

**8. **M. I. Stockman, “Nanofocusing of optical energy in tapered plasmonic waveguides,” Phys. Rev. Lett. **93**, 137404 (2004). [CrossRef]

**9. **Z. Fu, Q. Gan, Y. J. Ding, and F. J. Bartoli, “From waveguiding to spatial localization of THz waves within a plasmonic metallic grating,” IEEE J. Sel. Top. Quantum Electron. **14**, 486–490 (2008). [CrossRef]

**10. **Q. Gan, Z. Fu, Y. J. Ding, and F. J. Bartoli, “Ultrawide-bandwidth slow-light system based on THz plasmonic graded metallic grating structures,” Phys. Rev. Lett. **100**, 256803 (2008). [CrossRef]

**11. **J. He and S. He, “Slow propagation of electromagnetic waves in a dielectric slab waveguide with a left-handed material substrate,” IEEE Microw. Wirel. Compon. Lett. **16**, 96–98 (2005). [CrossRef]

**12. **K. L. Tsakmakidis, A. Klaedtke, D. P. Aryal, C. Jamois, and O. Hess, “Single-mode operation in the slow-light regime using oscillatory waves in generalized left-handed heterostructures,” Appl. Phys. Lett. **89**, 201103 (2006). [CrossRef]

**13. **K. L. Tsakmakidis, A. D. Boardman, and O. Hess, “‘Trapped rainbow’ storage of light in metamaterials,” Nature **450**, 397–401 (2007). [CrossRef]

**14. **K. Y. Kim, “Tunnelling-induced temporary light trapping in negative-index-clad slab waveguide,” Jpn. J. Appl. Phys. **47**, 4843–4845 (2008). [CrossRef]

**15. **Y. J. Huang, W. T. Lu, and S. Sridhar, “Nanowire waveguide made from extremely anisotropic metamaterials,” Phys. Rev. A **77**, 063836 (2008). [CrossRef]

**16. **G. Dolling, M. Wegener, C. M. Soukoulis, and S. Linden, “Negative-index metamaterial at 780 nm wavelength,” Opt. Lett. **32**, 53–55 (2007). [CrossRef]

**17. **J. Yao, Z. Liu, Y. Liu, Y. Wang, C. Sun, G. Bartal, A. M. Stacy, and X. Zhang, “Optical Negative Refraction in Bulk Metamaterials of Nanowires,” Science **321**, 930 (2008). [CrossRef]

**18. **T. Jiang and Y. J. Feng, “Slow and frozen waves in a planar air waveguide with anisotropic metamaterial cladding,” in *Proceedings of the 2008 International Conference on Microwave and Millimeter Wave Technology* (Nanjing, China, Apr. 21-24, 2008), pp. 1627–1630. [CrossRef]

**19. **J. He, Y. Jin, Z. Hong, and S. He, “Slow light in a dielectric waveguide with negative-refractive-index photonic crystal cladding,” Opt. Express **16**, 11077(2008). [CrossRef]

**20. **C. Li, Q. Sui, and F. Li, “Complex guided wave solutions of grounded dielectric slab made of metamaterials,” PIER **51**, 187–195 (2005). [CrossRef]

**21. **W. Shu and J. M. Song, “Complete mode spectrum of a grounded dielectric slab with double negative metamaterials,” PIER **65**, 103–123 (2006). [CrossRef]

**22. **A. W. Snyder and J. D. Love, *Optical Waveguide Theory* (Chapman and Hall, New York, 1983).

**23. **B. Wood, J. B. Pendry, and D. P. Tsai, “Directed subwavelength imaging using a layered metal-dielectric system,” Phys. Rev. B **74**, 115116(2006). [CrossRef]

**24. **R. Wangberg, J. Elser, E. E. Narimanov, and V. A. Podolskiy, “Nonmagnetic nanocomposites for optical and infrared negative-refractive-index media,” J. Opt. Soc. Am. B **23**, 498–505 (2006). [CrossRef]

**25. **D. W. Lynch and W. R. Hunter, *Handbook of Optical Constants of Solids*,
E.D. Palik, ed. (Academic Press, New York, 1985).