Plasmonic Modes of Metamaterial-Based Slot Waveguides

Most metamaterials exhibit pronounced anisotropic properties that are crucial for the understanding of their superior optical behavior, especially when they are integrated into the structure of a plasmonic waveguide. In this paper, we analytically solve the dispersion relation for a slot plasmonic waveguide filled with an anisotropic-stratified metamaterial and reveal that it supports two modes featuring relatively long propagation lengths in the limit of vanishing slot thickness. We classify these modes according to their physical origin and study the variation of their dispersion properties with material parameters.


Introduction
The ultimate goal of optics is to enable a perfect control of the interaction between light and matter.This goal has been brought closer by the recent advances in nanotechnology that have made possible the fabrication of optical metamaterials [1,2].The unusual electromagnetic properties of metamaterials are expected to enable a new generation of optical devices.In developing design strategies and new concepts for such devices, it is paramount that anisotropic properties of metamaterials are considered along with their other material features.Moreover, even the ways in which common devices operate require revisions when ordinary materials in their design are replaced by anisotropic metamaterials.A considerable amount of theoretical effort has been recently devoted to the analysis of optical propagation through different types of metamaterial structures, including uniaxial dielectrics [3] and indefinite media [4,5], metal-dielectric heterostructures [6] and superlattices [7], and strongly anisotropic waveguides [8].In this paper, we reexamine the guiding properties of slot plasmonic waveguides filled with an anisotropic medium.Our work is intended to demonstrate that integration of plasmonic waveguides with anisotropic optical metamaterials not only brings additional freedom to their design, but can also lead to new physical phenomena that may benefit the waveguide performance.
The plasmonic waveguide discussed here consists of an anisotropic medium of thickness 2h embedded between two metals of permittivity ε m .We assume that the medium's permittivity is described by a constant, diagonal tensor ε = diag(ε xx , ε yy , ε zz ) with its principle axes parallel to the waveguide's edges.Even though the permeability has similar anisotropic properties, only one component of its tensor affects the transverse magnetic (TM) modes, which are of primary interest for plasmonic waveguides.This allows us to describe the permeability using a single parameter μ.As is well known, the evolution of the electric field E = (E x , 0, E z ) and the magnetic field H = (0, H y , 0) of a TM mode is governed by the propagation constant β.In the case of an anisotropic core layer, β obeys the dispersion relation where , ω is the frequency of the surface plasmon polariton (SPP), c is the speed of light in vacuum, and the ± signs correspond to the symmetric and antisymmetric modes, respectively.Equation ( 1) is applicable to a broad range of metamaterial-based plasmonic waveguides and requires the specification of metamaterial design for further analysis.We restrict ourselves to the simplest scenario in which the metamaterial is created via stacking layers of isotropic materials into a periodic heterostructure.If the layers' thicknesses are much smaller than the optical wavelength, then the permittivity tensor of the heterostructure has only two different diagonal components given by [9] where the subscripts and ⊥ designate directions parallel and perpendicular to the layers, n is the number of layers within one heterostructure period, and f j is the filling factor for the jth layer of permittivity ε j (all filling factors must add up to unity).Such a stratified metamaterial offers fabrication of the three types of plasmonic waveguides shown in Figure 1.Their permittivities are Although all of them can, in principle, be created using modern fabrication techniques, we focus here on the last structure that is relatively easy to fabricate.Apart from the ease of fabrication, the third type of plasmonic waveguide in Figure 1 is the only one that supports symmetric TM modes with complex β values in the limit of vanishing ε , which can be realized by properly matching metamaterial constituents and compensating for the absorption losses inside them (note that the condition ε = 0 is not equivalent to an epsilon-near-zero (ENZ) regime [10], because the permittivity is a complex tensor).Within the slot region, these modes are characterized by the components E x and H y that are independent of x, and by E z ∝ x.Solving (1) in the limit ε → 0 yields the following four values of β: where K = ε m (kh) 2 and r = ε m /ε ⊥ is the ratio of the dielectric constants.Since both of these parameters are generally complex, β = β + iβ is also complex, and its real and imaginary parts provide the phase velocity and energy loss of various SPP modes.Because β in (3) satisfies (1) approximately in the limit of small h (when |q|h 1) even when ε / = 0, the following analysis allows us to elucidate some general features of the SPP behavior in metamaterialbased plasmonic waveguides.
For simplicity, we focus on the stratified heterostructure composed of two different nonmagnetic (μ = 1) materials, one of which is the metal used for waveguide cladding.The permittivity of the second material is then fixed by the condition ε = 0 (this material should provide gain, since metals are essentially lossy).It then follows that the ratio r depends solely on the filling factor f of the metal; r = 2 − 1/ f .We use silver as a metal in our simulations, take its permittivity in the form of a seven-pole Drude-Lorentz formula given in the work by Pannipitiya et al. [11], and introduce an absorption parameter γ, such that ε m (γ) = ε + iγε .For definiteness, we also choose β > 0 and refer to β as the "damping factor," thus implying that our waveguides do not amplify SPPs.
According to (3), the waveguide in Figure 1(c) supports a maximum of two SPP modes regardless of slot's thickness, provided that ε = 0.The dispersion curves of these two modes are plotted in Figure 2 for three values of γ and two sets of material parameters f and h.In the lossless case (γ = 0, green curves), the modes can be grouped into three distinct classes: (i) propagating modes with β = 0 and β / = 0; (ii) complex modes with β / = 0 and β / = 0; (iii) evanescent modes with β / = 0 and β = 0.The first two types of modes may travel either forward or backward, except for certain "degenerate" frequencies for which the group velocity and total energy flow of SPPs vanish.In our example, the degenerate frequencies correspond to the points A, B, and C of the energy spectrum.The degeneracy of waves traveling in the opposite directions is removed by absorption (see pink curves for γ = 0.2) since the two modes develop a slightly different lateral confinement.
The complex modes of a lossy waveguide arise through quasimixing of the three types of modes as γ is gradually varied from 0 to 1.Such transformation of the propagating, complex, and/or evanescent modes is accompanied by their degeneracy removal, which is also the scenario followed by the quadrupole modes of plasmonic nanowires [12].Indeed, the twofold degenerate propagating mode at point B and the twofold degenerate evanescent mode at point C split into two nondegenerate complex modes for γ = 1 (see Figure 2).The process of modes transformation allows one to track their  origin and better understand waveguide performance.For example, because of the strong damping intrinsic to evanescent modes, one can expect that the segments of spectra β (ω) emerging from them in the case of real losses will be lesser affected by the dispersion properties of ε (ω) than the segments that stem from undamped propagating modes.This conclusion is confirmed by Figures 2(b) and 2(d), where the damping factors are separately shown for the forward (right curves) and backward (left curves) propagating SPP modes.Notice that the spectra in Figure 2 are congregate of the results valid for different plasmonic waveguides, rather than spectra corresponding to a specific structure.As we have seen, dispersions of the real and imaginary parts of β drastically depend on the filling factor of metal and slot thickness.To study this dependance, we set ω = 0.2ω p , where ω p is the plasma frequency of silver, and plot in Figure 3 β and β as functions of slot thickness for five values of f .By looking at this figure from left to right, one can draw up the following picture of SPP behavior.When the filling factor of metal is small, two SPP modes of different Advances in OptoElectronics origins are supported by the waveguide.One of themmode 1 shown by red curves-springs from the propagating mode of a lossless waveguide and, therefore, exhibits small damping.The other-mode 2 plotted in blue-arises from the evanescent mode when h < h 0 = k −1 |ε m | −1/2 ≈ 25 nm and from the propagating mode when h > h 0 (h 0 is the skin depth of a metal in the limit f → 0).The first mode propagates backward regardless of h, whereas the second one changes its traveling direction from forward to backward around h = h 0 /2.Because modes 1 and 2 originate from the propagating modes for h > h 0 and f < 0.5, their damping rates are almost the same for this range of parameters h and f .When f = 0.5, the ratio r = 0 since the transverse component of the permittivity tensor becomes infinite, and only mode 2 with β = ±k ε m (1 + K) survives.For f > 0.5, mode 1 becomes forward propagating, while mode 2 gets strongly attenuated for all h.The overall dynamics of such changes with f is illustrated in the supplementary animation.
It should be emphasized that these results are based on the effective-medium theory and approximate.For a more accurate treatment, the inclusion of nonlocality effects is required [13,14].An important feature of the modes given in (3) is that they can exhibit relatively small damping depending on the values of the parameters f and h.For example, low damping occurs for both modes for f < 0.2 and h > h 0 .Moreover, the damping of mode 1 vanishes in the limit of small h for all values of f / = 0.5.This behavior is opposite to that of SPPs in ordinary (isotropic) metal-dielectric-metal waveguides, whose attenuation rate diverges as 1/h [15].Also noteworthy is that the confinement of mode 1, determined by the evanescent decay length scaling like δ 1 ∝ h Re √ r 2 when h → 0, improves with the reduction of the waveguide thickness, whereas the confinement of mode 2 worsens in the limit of small h, due to the divergence of the decay length δ 2 ∝ (hk 2 0 Re ε 2 m (1 − r) 2 ) −1 . Hence, we expect that metamaterialbased plasmonic waveguides may enable strong localization of optical energy in an SPP mode that can propagate over relatively long lengths.
In conclusion, we have found that a maximum of two symmetric TM modes with essentially different propagation properties are supported by plasmonic slot waveguides filled with a stratified metamaterial exhibiting a relatively small permittivity in the plane of the layers.

Figure 1 :
Figure 1: (Color online) Three types of plasmonic waveguides made using a stratified metamaterial.In all cases, SPP modes propagate in the z direction and the waveguide extends to infinity in the y direction.

Figure 2 :
Figure 2: (Color online) Transformation of [(a) and (c)] SPP energy spectrum and [(b) and (d)] damping-factor dispersion with changing absorption parameter γ.The SPP frequency is normalized to the plasma frequency ω p of silver.Solid green curves show propagating modes, whereas dashed curves show complex and evanescent modes.

Figure 3 :
Figure 3: (Color and multimedia online) Transformation of the real (upper panels) and imaginary (lower panels) parts of β at a specific frequency ω = 0.2ω p as the filling factor f of the metal is increased from 0.2 to 0.8.Insets show a magnified view of mode 2.