Figures of Merit for 2d Surface Plasmon Waveguides and Application to Metal Stripes

Three figures of merit, useful as quality measures for 2D surface plasmon waveguides, are discussed and applied to help trade-off mode confinement against attenuation for the symmetric mode propagating along metal stripes. Different stripe geometries are considered, and Au, Ag and Al are compared as the stripe metal over the wavelength range from 200 to 2000 nm. Depending on which figure of merit is used, and on how mode confinement is measured, different preferred designs emerge. For instance, given a mode area, narrow thick stripes are better than wide thin ones, but given a distance from the light line, the opposite is true. Each of the metals analyzed show wavelength regions where their performance is best. The figures of merit are generally applicable and should be useful to help compare, assess and optimize designs in other 2D surface plasmon waveguides or in other absorbing waveguides. Plasmonics: Localization and guiding of electromagnetic energy in metal/dielectric structures, " J. Surface-polariton-like waves guided by thin, lossy metal films, " Phys. Rev. B 33, 5186-5201 (1986). 9. P. Berini, " Plasmon-polariton modes guided by a metal film of a finite width, " Opt. Plasmon polariton waves guided by thin lossy metal films of finite width: bound modes of symmetric structures, " Phys. Demonstration of integrated optics elements based on long-ranging surface plasmon polaritons, " Opt. Passive integrated optics elements based on long-range surface plasmon polaritons, " J. Demonstration of Bragg gratings based on long-ranging surface plasmon polariton waveguides, " Opt. " Surface plasmon propagation in microscale metal stripes, " Appl. observation of surface plasmon polariton propagation on thin metal stripes, " Phys. Plasmon polaritons of metallic nanowires for controlling submicron propagation of light, " Phys. Geometries and materials for subwavelength surface plasmon modes, " J. Plasmon slot waveguides: Towards chip-scale propagation with subwavelength-scale localization " Phys. Channel plasmon-polariton guiding by subwavelength metal grooves, " Phys. Extended long range plasmon waves in finite thickness metal film and layered dielectric materials, " Opt. Simulations of hybrid long-range plasmon modes with application to 90° bends, " Opt. Efficient and accurate numerical analysis of multilayer planar optical waveguides in lossy anisotropic media, " Opt.range surface plasmon-polariton mode cutoff and radiation in embedded strip waveguides, " J. Characterization of long-range surface-plasmon-polariton waveguides, " J.


Introduction
The optical properties of metals have been studied and are known to exhibit a negative real part of permittivity at optical frequencies [1].This property allows the metal-dielectric interface to support a surface plasmon-polariton (SPP) mode which is bound to the interface through the coupling of electromagnetic waves to oscillations in conduction electrons in the metal.The SPP has been studied extensively, and has been reviewed by Raether [2], Welford [3], Barnes [4,5], Maier and Atwater [6] among others.The high loss of the single interface SPP makes it, in general, incapable of long range transmission.
A thin metal slab bounded by dielectric supports bound supermodes, labeled s b and a b [7,8].As the thickness of the slab approaches zero the attenuation of the s b mode, or longrange SPP (LRSPP), decreases and the mode becomes less confined eventually evolving into the TEM wave of the background dielectric.When the width of the thin metal slab is made finite, defining a stripe, and the metal stripe is embedded in a homogeneous background dielectric, then a more complicated set of modes emerge [9,10,11] since the stripe provides 2D confinement in the plane transverse to the direction of propagation.For a sufficiently thin or narrow metal stripe, the ss b 0 mode becomes long-range (it is the LRSPP).Numerous LRSPP passive devices such as Mach-Zehnder interferometers, Y-junctions, directional couplers [12,13] and Bragg gratings [14], excited by butt-coupling with an optical fiber, have been experimentally tested and shown to closely agree with theory.Similar devices were also reported in [15,16].The metal stripe exposed to air has also been studied [17][18][19].
The metal slab and stripe exhibit a trade-off with regards to the confinement and attenuation of the LRSPP supported therein: The confinement and attenuation of the LRSPP rise and fall together [10,20].This trade-off also applies to other SPP waveguides such as the single-interface, the metal cladded dielectric slab [21,22] and the channel waveguide [23][24][25], for example.The attenuation-confinement trade-off in surface plasmon waveguides was recently discussed in [26], where three FoMs (figures of merit), denoted M 1 1D , M 2 and M 3 , were proposed to assist with trade-off analyses.The FoMs were then used to study three 1D waveguides: the single-interface, the metal slab and the metal cladded dielectric slab [26].
In this paper, we extend the M 1 1D FoM proposed in [26] to 2D waveguides, defining the new variant M 1 2D (Section 2), and we use the three FoMs (M 1 2D , M 2 and M 3 ) to study the LRSPP in the metal stripe as a function of geometry (Section 3) and metal choice (Au, Ag and Al -Section 4).Section 5 gives a brief summary and conclusions.

Figures of merit for 2D waveguides
An exp(jωt) time dependence is assumed with mode propagation occurring along the +z axis with an exp(-γ z z) dependence, where ω is the angular frequency and γ z = α z + jβ z is the complex propagation constant with α z and β z the attenuation and phase constants, respectively.The complex effective index N eff is given by N eff = γ z /β 0 = k eff + jn eff where β 0 = 2π/λ 0 is the free-space phase constant and λ 0 is the free-space operating wavelength.

Figures of merit: M 1 1D , M 2 and M 3
Three FoMs were defined in [26] to provide objective measures of comparison for purely bound surface plasmon modes in SPP waveguides.Their definition is based on forming benefit-to-cost ratios where the benefit is confinement and the cost is attenuation.Different ways of measuring confinement led to different definitions.The first FoM, M 1 1D = 1/δ w α z , uses the inverse mode size (δ w ) as its confinement measure, where δ w is the distance between the 1/e field magnitude points of the main transverse electric field component relative to the global field maximum.The second FoM, M 2 = (β z -β 1 )/α z = (n eff -n 1 )/k eff measures confinement as the mode's distance from the light line in the dielectric.The third FoM, M 3 = 1/λ g α z = β z /2πα z = n eff /2πk eff , uses the inverse guided wavelength λ g as its confinement measure.M 3 is proportional to the quality factor (Q) when dispersion is negligible.

Definition of the figure of merit M 1 2D for 2D waveguides
The definition of M 2 , and M 3 holds for modes in 2D waveguides, but the definition of M 1 , which is based on mode size, depends on the dimensionality of the structure as emphasized in [26].For 1D waveguides, the mode size is the width δ w , leading to M 1 1D .For 2D waveguides, the mode size is an area, leading to a new definition for M 1 , denoted M 1 2D .The mode size is taken as the area A e bounded by the closed 1/e field magnitude contour relative to the global field maximum.The 2D spatial distribution in the transverse plane of the main transverse electric field component is used to find this contour and the area A e .The contour and A e are determined numerically, so highly deformed modes, sometimes supported by SPP waveguides, are easily and unambiguously handled.
Taking the confinement measure as (π/A e ) 1/2 , instead of simply 1/A e , seems preferable since this measure tends to the inverse mode radius as the mode becomes circular and it leads to a unit-less FoM.This confinement measure also tends to zero as the mode expands or as the waveguide evolves into a 1D structure (as it should since confinement is lost along one of the transverse dimensions).Based on these considerations, the M 1 2D FoM is defined as:

Geometric study of metal stripe waveguides
The waveguides considered are shown in Fig. 1, and consist of (a) the metal stripe (w < ∞) and slab (w = ∞), (b) and (c) 2 and 3 symmetrically coupled (SC) metal stripes, and (d) the cladded metal stripe.The cladded metal slab (w = ∞) was analyzed in [27], and the cladded metal stripe was recently reported in [28] and added to this paper during revisions.The metal stripe on a thin dielectric layer or membrane (not shown) was recently introduced in [29] and bears points of similarity to the cladded metal stripe [28].The metal slab was analyzed using the transfer matrix method [30].A commercial software package based on the finite element ), with the surrounding dielectric being SiO 2 (ε r,1 = n 1 2 = 2.085 [32]) and, additionally, vacuum (ε 0 ) in the case of the cladded stripe (Fig. 1(d)).Figure 2(c) shows that as t decreases, the mode size increases for all structures considered, as expected.From Fig. 2(d), it is noted that M 1 1D and M 1 2D increase for all structures as t and/or w decrease, indicating that α z decreases more rapidly than the modes' expansion.Given a mode size (A e /π) 1/2 , single narrow thick stripes are better than wide thin ones or coupled ones, since they generate less attenuation yielding a larger M 1 2D .For a specific t narrower stripes produce a larger M 1 2D than wider ones.The modes' distance from the light line, plotted in Fig. 2(e), decreases with t and w as expected since the modes evolve into the TEM wave of the background as the metal vanishes.From Fig. 2(f) it is noted that M 2 increases sharply with decreasing t, reaching a peak beyond which it tends to 0 as t 0. These peaks are located at t = 18 nm for the w = 2 μm stripe, at t = 10 nm for the w = 8 μm stripe, at t = 9.5 nm for the wswsw = 22222 μm SC stripes, and at t ~ 0 for the s b mode in the slab.On the thicker side of the peaks, k eff decreases more rapidly than the confinement (n eff -n 1 ) as t is reduced, but the opposite holds true on the thinner side of the peak.Given a distance from the light line (n eff -n 1 ), single wide thin stripes perform better than narrow thick ones or coupled ones, since they generate less attenuation, yielding a larger M 2 .For a specific t wider waveguides produce a larger M 2 than narrower ones.These trends are opposite to those observed from M 1 2D .
λ g plotted in Fig. 2(g) increases with decreasing t and w.M 3 , plotted in Fig. 2(h), shows a similar trend to the other FoMs in that M 3 increases with decreasing t.This implies that α z decreases more rapidly than λ g increases as t is reduced.Given a t narrower waveguides are better than wider ones leading to a larger M 3 , as observed for M 1 2D .The cladded metal stripe (Fig. 1(d)) [28] is similar to the metal stripe on a thin dielectric membrane [29], in that as the thickness of the dielectric changes, the ss b 0 mode may become more confined and attenuated.Another point of similarity rests with the conditions for ss b 0 confinement, which are that its n eff must be larger than n 1 and than the n eff of the TM 0 mode in the dielectric slab present to the left and right of the metal stripe.If TM-TE mode conversion is expected, say due to discontinuities, then it should also be larger than the n eff of the TE 0 mode of the dielectric slab, but this condition is essentially otherwise irrelevant as was demonstrated experimentally in [29] since the ss b 0 and TE 0 modes are substantially orthogonal.In a symmetric slab, the TM 0 and TE 0 modes are guided for all dielectric thicknesses.Another point of similarity rests with the excitation of the waveguides, in that if the source and ss b 0 mode fields are not well matched then light becomes trapped in the dielectric slab and may interfere with the ss b 0 mode.This may be problematic in structures such as couplers and Mach-Zehnder interferometers.It is noted that M 1 2D , M 2 and M 3 of the cladded metal stripe (Fig. 1(d)) are larger than those of the metal stripe (Fig. 1(a), d = ∞) over a good range of dimensions, indicating that the former can provide a better trade-off between confinement and attenuation, as noted in [28].It must be borne in mind, however, that the metal stripe must have a smaller t or w or both as d decreases in order to maintain the same α z and thus the same range as the d = ∞ case.
Producing high quality metal stripes can be challenging for t < 20 nm [33].3(c).There is also an inner contour(s) that exists near the surface of the stripe(s), which depends on the penetration of |E y | into the metal.In the case of Figs.3(a) and (c), the area A e is the area within the outer closed contour, not including the area within the inner closed contour since the field enclosed by the latter has a magnitude below the 1/e level.The movies associated with Figs.3(a) and (b) show how the mode evolves with t.Both waveguides show the field and 1/e contour extending further into the dielectric as t is reduced.This decrease in confinement also causes the mode field to become more circular, more closely matching that of a single mode fiber.This can be seen in the movie in Fig. 3(a) with the mode field almost circular for t ≤ 20 nm.From the movie in Fig. 3(b), weaker coupling between the three SC stripes is observed from the shape of the mode field as t and the confinement increase.The movie in Fig. 3(c) shows the mode size decreasing with d until d ~ 900 nm where the 1/e field begins to extend further into the vacuum region.

Wavelength response of metal stripe waveguides
The FoMs are now used to assess the wavelength response of metal stripe waveguides for three different metals embedded in SiO 2 : Ag, Au and Al.The optical properties of these materials were taken from the experimentally determined data compiled in Palik [32].A cubic spline was used to interpolate the permittivity at the desired wavelengths.The waveguides used in the analysis were chosen based on their similar proximity to the light line at λ 0 = 1550 nm (all have n eff -n 1 ~ 5×10 -4 , Fig. 2(e)): A single stripe with w = 2 μm, t = 40 nm, a single stripe with w = 4 μm, t = 20 nm, and a pair of SC stripes with wsw = 222 μm, t = 30 nm.
Figure 4(a) shows n eff versus λ 0 of the ss b 0 mode, along with the light line in SiO 2 ; Fig. 4(b) shows the corresponding α z .An interesting feature arose in the case of Au: n eff of the single interface SPP crosses the light line at λ 0 = 473 nm (slightly beyond the energy asymptote), while n eff of the slab and of the stripes do not cross the light line at this wavelength.The modes appear to remain bound for λ 0 ≥ 250 nm.This trend does not show up in the case of Ag.The modes supported by the Al waveguides did not reach their energy asymptote for the range of wavelengths used, which makes it unclear as to whether Al would show this feature.Fig. 4(b) shows that lower loss can be achieved when operating deeper into the infrared for all the waveguides and metals analyzed.For all metals, as the wavelength decreases, the ss b 0 mode becomes more localized to the center of the stripe along the metal/dielectric interfaces and the effects of the corners and the dielectric regions on either side diminish.
Figure 4(c) shows that the mode size increases when operating further into infrared.At short wavelengths the curves exhibit features and small peaks which are explained by inspecting the mode fields, revealing that they increase in strength near the corners at specific wavelengths.This behavior also perturbs the 1/e contour associated with the modes.
Figure 4(d) shows that M 1 2D increases with increasing wavelength, whereas M 1 1D remained flat in the Drude region as observed from Fig. 3 of [26].The modes supported by the Ag waveguides experience a small drop in M 1 2D near λ 0 = 1630 nm while the Al waveguides show a drop for λ 0 between 700 and 900 nm.Looking at Fig. 4(b), an increase α z is observed in these regions, which is not matched by an equivalent increase in confinement, causing the dip in M 1 2D .In the wavelength range around 1630 nm, where the Ag structures show a decrease in M 1 2D , Au out-performs.Al is a good choice for λ 0 ≤ 400 nm.
From Fig. 4(e) it is noted that the distance from the light line decreases sharply at long wavelengths.Interestingly, M 2 plotted in Fig. 4(f), exhibits maxima for Au and Ag.The Ag structures exhibit a peak with its location depending slightly on the geometry: the w = 4 μm t = 20 nm, wsw = 222 μm t = 30 nm, and w = 2 μm t = 40 nm structures have peaks at λ 0 = 845, 835 and 825 nm, respectively.The corresponding Au structures also exhibit peaks in M 2 near these wavelengths, but they are broader and lower than those of the Ag structures.Ag outperforms Au over a broad wavelength range but beyond 1630 nm Au seems better.M 2 decreases sharply on the long wavelength side of the peak, in keeping with the rapid loss of confinement noted earlier from Fig. 4(e).Based on M 2 , it is preferable to operate Au and Ag waveguides near the short wavelength limit of the Drude region, as was already observed for Ag in [26].The modes supported by the Al waveguides show a decrease in M 2 for λ 0 in the range of 400 to 900 nm.The Al waveguides again perform better than the Au and Ag ones for λ 0 ≤ 400 nm, while being considerably outperformed throughout the rest of the spectrum.Figure 4(f) shows the quality factor Q. (In the calculation of Q, the derivative in Eq. ( 18) of [26] was estimated using finite difference approximants.)M 3 is shown in Fig. 4(h).When the modes experience minimal dispersion, as they do in the Drude region, then Q ~ πM 3 , as is observed by comparing these Figs.In the short wavelength region the modes experience a large change in their n eff as they approach the energy asymptote, causing a greater difference between Q and M 3 .Both sets of curves show a similar trend to M 1 2D , where the Au and Ag

Fig. 1 .
Fig. 1.Cross sectional view of surface plasmon waveguides.(a) Single stripe, (b) pair of SC stripes, (c) three SC stripes, (d) cladded stripe.The metal stripes (Figs.1(a)-(c)) are discussed first.Figs.2(a) and (b) give the effective index (n eff ) and attenuation (α z , k eff ) of the s b and ss b 0 modes as a function of t.The usual trends of n eff n 1 (vanishing confinement) and α z , k ef f 0 (vanishing attenuation) as t 0 are noted.

Fig. 2 .
Fig. 2. (a) n eff and (b) α z and k eff of the ss b 0 and s b modes.(c) Mode size for the s b mode in the slab (right axis) and the ss b 0 mode in the stripes (left axis).(d) M 1 1D (right axis) and M 1 2D (left axis).(e) Distance from the light line and (f) M 2 .(g) Guided wavelength and (h) M 3 .The gray curves are for the cladded stripe (Fig. 1(d)).

Figure 2 (
Figure 2(a) shows n eff of the ss b 0 mode in the cladded metal stripe for w = 4 μm and t = 20 nm, as well as n eff of the TM 0 and TE 0 modes in the dielectric slab alone, as a function of d.From Fig. 2(a) we see that n eff of the ss b 0 mode is always larger than that of the TM 0 mode, while it falls below that of the TE 0 mode for d ~ 1.9 μm.n eff of the ss b 0 mode decreases as d decreases, which is expected since the mode fields extend further into the vacuum.In Fig. 2(b) α z increases with d until a maximum, from which it drops quickly.The mode size plotted in Fig. 2(c) decreases with decreasing d until a minimum is reached around d = 900 nm.For smaller d the mode extends deeper into the vacuum, explaining the decreasing α z and n eff .

Fig. 3 .
Fig. 3. Spatial distribution of |E y | associated with the ss b 0 mode in various waveguides.Quarter symmetry is used with the origin, (x = y = 0), being the center of the mode.The 1/e field contour is also plotted as the thin black curve.In all cases, the fields are normalized such that max(|E y |) = 1. (a) Single stripe with w = 8 μm and t = 70 nm; the associated movie shows a sweep over t.(b) Three SC stripes with wswsws = 22222 μm and t = 70 nm; the associated movie shows a sweep over t.(c) Cladded stripe with w = 4 μm, t = 20 nm and d = 0.8 μm; the associated movie shows a sweep over d.(d) Pair of SC Au stripes, wsw = 222 μm, t = 30 nm and λ 0 = 1000 nm; the associated movie shows a sweep overλ 0 .Figures 3(a), (b) and (c) give the spatial distribution of |E y | associated with the ss b 0 mode in a single stripe having w = 8 μm and t = 80 nm, in three SC stripes having wswsw = 22222 μm and t = 70 nm and in a cladded stripe having w = 4 μm and t = 20 nm, respectively.The 1/e field contour is plotted as the thin black curve and the SiO 2 -vaccuum interface is shown as the dashed white line in Fig.3(c).There is also an inner contour(s) that exists near the surface of the stripe(s), which depends on the penetration of |E y | into the metal.In the case of Figs.3(a) and (c), the area A e is the area within the outer closed contour, not including the area within the inner closed contour since the field enclosed by the latter has a magnitude below the 1/e level.The movies associated with Figs.3(a) and (b) show how the mode evolves with t.Both waveguides show the field and 1/e contour extending further into the dielectric as t is reduced.This decrease in confinement also causes the mode field to become more circular, more closely matching that of a single mode fiber.This can be seen in the movie in Fig.3(a) with the mode field almost circular for t ≤ 20 nm.From the movie in Fig. 3(b), weaker coupling

Fig. 4 .
Fig. 4. (a) n eff and (b) α z versus λ 0 .(c) Mode size and (d) M 1 2D .(e) Distance from the light line and (f) M 2 .(g) Quality factor Q and (h) M 3 .(green: Au, blue: Ag, magenta: Al) Figure 3(d) gives the spatial distribution of |E y | associated with the ss b 0 mode in two SC Au stripes having wsw = 222 μm and t = 30 nm.The associated movie shows how the mode evolves with decreasing λ 0 .As the mode approaches its energy asymptote, it becomes more confined and localized to the center of the stripe along the metal/dielectric interfaces.At short λ 0 , the coupling between the stripes weakens, eventually disappearing altogether.