OSA's Digital Library

Optics Express

Optics Express

  • Editor: Andrew M. Weiner
  • Vol. 21, Iss. 2 — Jan. 28, 2013
  • pp: 1430–1439
« Show journal navigation

Nonlocal effects in a hybrid plasmonic waveguide for nanoscale confinement

Qiangsheng Huang, Fanglin Bao, and Sailing He  »View Author Affiliations


Optics Express, Vol. 21, Issue 2, pp. 1430-1439 (2013)
http://dx.doi.org/10.1364/OE.21.001430


View Full Text Article

Acrobat PDF (1250 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

The effect of nonlocal optical response is studied for a novel silicon hybrid plasmonic waveguide (HPW). Finite element method is used to implement the hydrodynamic model and the propagation mode is analyzed for a hybrid plasmonic waveguide of arbitrary cross section. The waveguide has an inverted metal nano-rib over a silicon-on-insulator (SOI) structure. An extremely small mode area of~10−6λ2 is achieved together with several microns long propagation distance at the telecom wavelength of 1.55μm. The figure of merit (FoM) is also improved in the same time, compared to the pervious hybrid plasmonic waveguide. We demonstrate the validity of our method by comparing our simulating results with some analytical results for a metal cylindrical waveguide and a metal slab waveguide in a wide wavelength range. For the HPW, we find that the nonlocal effects can give less loss and better confinement. In particular, we explore the influence of the radius of the rib’s tip on the loss and the confinement. We show that the nonlocal effects give some new fundamental limitation on the confinement, leaving the mode area finite even for geometries with infinitely sharp tips.

© 2013 OSA

1. Introduction

In order to realize electronic photonic integrated circuits, we need to develop optical waveguides capable of guiding light with deep subwavelength confinement. It is well known that surface plasmon polariton (SPP) waveguides can provide subwavelength confinement in the scale of 100nm or less, by coupling the electromagnetic waves to electron oscillations at the metallic surface. In the past years, people have demonstrated several novel structures, such as metal slot waveguides [1

1. L. Liu, Z. H. Han, and S. L. He, “Novel surface plasmon waveguide for high integration,” Opt. Express 13(17), 6645–6650 (2005). [CrossRef] [PubMed]

, 2

2. Z. Han, A. Y. Elezzabi, and V. Van, “Experimental realization of subwavelength plasmonic slot waveguides on a silicon platform,” Opt. Lett. 35(4), 502–504 (2010). [CrossRef] [PubMed]

], channel SPP (CPP) waveguides [3

3. D. F. Pile and D. K. Gramotnev, “Channel plasmon-polariton in a triangular groove on a metal surface,” Opt. Lett. 29(10), 1069–1071 (2004). [CrossRef] [PubMed]

, 4

4. S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, J. Y. Laluet, and T. W. Ebbesen, “Channel plasmon subwavelength waveguide components including interferometers and ring resonators,” Nature 440(7083), 508–511 (2006). [CrossRef] [PubMed]

], and wedge plasmon polariton (WPP) waveguides [5

5. D. Pile, T. Ogawa, D. K. Gramotnev, T. Okamoto, M. Haraguchi, M. Fukui, and S. Matsuo, “Theoretical and experimental investigation of strongly localized plasmons on triangular metal wedges for subwavelength waveguiding,” Appl. Phys. Lett. 87(6), 061106–061103 (2005). [CrossRef]

8

8. A. Boltasseva, V. S. Volkov, R. B. Nielsen, E. Moreno, S. G. Rodrigo, and S. I. Bozhevolnyi, “Triangular metal wedges for subwavelength plasmon-polariton guiding at telecom wavelengths,” Opt. Express 16(8), 5252–5260 (2008). [CrossRef] [PubMed]

]. However, they lead to high propagation loss. Recently, several hybrid plasmonic waveguides have been proposed to achieve both subwavelength mode confinement and relatively low loss [9

9. R. F. Oulton, V. J. Sorger, D. A. Genov, D. Pile, and X. Zhang, “A hybrid plasmonic waveguide for subwavelength confinement and long-range propagation,” Nat. Photonics 2(8), 496–500 (2008). [CrossRef]

13

13. D. X. Dai, Y. C. Shi, S. L. He, L. Wosinski, and L. Thylen, “Gain enhancement in a hybrid plasmonic nano-waveguide with a low-index or high-index gain medium,” Opt. Express 19(14), 12925–12936 (2011). [CrossRef] [PubMed]

], which makes them attractive as potential candidates for the realization of nano-photonic circuits. As enhancing light confinement while maintaining the transmission loss is a goal for most of the hybrid waveguide designs, the trade-off between the mode confinement and loss still exists. Currently, the miniaturization trend in the design of plasmonic devices is reaching the limit where the local solutions of macroscopic Maxwell’s equation can no longer accurately describe their electromagnetic properties due to a longitude mode near the metal surface. To overcome this, a sophisticated nonlocal material model is required, such as the hydrodynamic model of the electron gas as discussed by Boardman [14

14. A. D. Boardman, Electromagnetic Surface Modes (Wiley, 1982).

]. Recently, a number of nonlocal solutions of Maxwell’s equations in two dimensions (2D) or 2.5D or 3D (Cartesian coordinates or cylindrical coordinates) have been proposed based on numerical approaches such as FDM [15

15. J. M. McMahon, S. K. Gray, and G. C. Schatz, “Nonlocal optical response of metal nanostructures with arbitrary shape,” Phys. Rev. Lett. 103(9), 097403 (2009). [CrossRef] [PubMed]

], FDTD [16

16. J. M. McMahon, S. K. Gray, and G. C. Schatz, “Calculating nonlocal optical properties of structures with arbitrary shape,” Phys. Rev. B 82(3), 035423 (2010). [CrossRef]

] or FEM [17

17. G. Toscano, S. Raza, A. P. Jauho, N. A. Mortensen, and M. Wubs, “Modified field enhancement and extinction by plasmonic nanowire dimers due to nonlocal response,” Opt. Express 20(4), 4176–4188 (2012). [CrossRef] [PubMed]

22

22. C. Ciracì, R. T. Hill, J. J. Mock, Y. Urzhumov, A. I. Fernández-Domínguez, S. A. Maier, J. B. Pendry, A. Chilkoti, and D. R. Smith, “Probing the ultimate limits of plasmonic enhancement,” Science 337(6098), 1072–1074 (2012). [CrossRef] [PubMed]

] for frequency domain stationary application or transient propagation application. An insightful transformation-optics method [23

23. A. I. Fernández-Domínguez, A. Wiener, F. J. García-Vidal, S. A. Maier, and J. B. Pendry, “Transformation-optics description of nonlocal effects in plasmonic nanostructures,” Phys. Rev. Lett. 108(10), 106802 (2012). [CrossRef] [PubMed]

] has also been proposed to study the optical properties of plasmonic nanostructures and has revealed the nonlocal effect to the reduced field enhancement for e.g. metallic dimer of touching metal nanowires. An experiment has also shown that the intrinsic nonlocality of the metal limits the field enhancement [22

22. C. Ciracì, R. T. Hill, J. J. Mock, Y. Urzhumov, A. I. Fernández-Domínguez, S. A. Maier, J. B. Pendry, A. Chilkoti, and D. R. Smith, “Probing the ultimate limits of plasmonic enhancement,” Science 337(6098), 1072–1074 (2012). [CrossRef] [PubMed]

]. In the analysis of propagation mode for some plasmonic waveguide, early formulation treatments have been given for a metal slab waveguide [24

24. G. C. Aers, B. V. Paranjape, and A. D. Boardman, “Non-radiative surface plasma-polariton modes of inhomogeneous metal circular cylinders,” J. Phys. F 10(1), 53–65 (1980). [CrossRef]

, 25

25. R. Ruppin, “Effect of non-locality on nanofocusing of surface plasmon field intensity in a conical tip,” Phys. Lett. A 340(1-4), 299–302 (2005). [CrossRef]

], a metal cylindrical waveguide [26

26. R. Ruppin, “Non-local optics of the near field lens,” J. Phys. Condens. Matter 17(12), 1803–1810 (2005). [CrossRef]

] or aggregates thereof [27

27. F. J. García de Abajo, “Nonlocal effects in the plasmons of strongly interacting nanoparticles, dimers, and waveguides,” J. Phys. Chem. C 112(46), 17983–17987 (2008). [CrossRef]

]. For a plasmonic waveguide of large size with apex feature, the analytical method in [27

27. F. J. García de Abajo, “Nonlocal effects in the plasmons of strongly interacting nanoparticles, dimers, and waveguides,” J. Phys. Chem. C 112(46), 17983–17987 (2008). [CrossRef]

] is invalid. Although there has been an enormous advance in the theoretical analysis of nonlocal effects in plasmonic structures, to the best of our knowledge, there is no report on a propagation mode analysis of nonlocal effect in a plasmonic waveguide of arbitrary cross-section.

In this letter, we introduce a novel highly confining hybrid plasmonic waveguide with an inverse triangular metal rib separated from a high-index material by a low-index material. The mode area of the hybrid mode can be reduced to about 10−6λ2, while still maintaining propagation distances exceeding several microns at telecom wavelengths. We also investigate how nonlocal effects alter the confinement of the hybrid plasmonic waveguide. However, we don’t consider the quantum effect at the metal surface [28

28. D. C. Marinica, A. K. Kazansky, P. Nordlander, J. Aizpurua, and A. G. Borisov, “Quantum plasmonics: Nonlinear effects in the field enhancement of a plasmonic nanoparticle dimer,” Nano Lett. 12(3), 1333–1339 (2012). [CrossRef] [PubMed]

, 29

29. R. Esteban, A. G. Borisov, P. Nordlander, and J. Aizpurua, “Bridging quantum and classical plasmonics with a quantum-corrected model,” Nat Commun 3, 825 (2012). [CrossRef] [PubMed]

]. This is because the effect of electron tunneling between the metal and the dielectric is weaker than that between the metal nanoparticles. First, we present a FEM nonlocal solution of Maxwell’s equation using the weak formulation. Then we validate the simulation method by comparing the simulation results with the analytical results for a metal cylindrical waveguide and a metal slab waveguide. Moreover, we compare the results simulated by the FEM nonlocal method and the traditional FEM Drude method. At last, we investigate how the nonlocal effects change the sensitivity of the electromagnetic field to the radius of the tip. We demonstrate that the nonlocal effects drastically change the nanoscale confinement performance of plasmonic tips when the radius of the tip is comparable to the scale of the longitudinal mode penetration depth.

2. Theoretical formalism

In the bulk metal, the dispersion relation of the longitude mode can be described by the following equation:
ω(ω+iγ)=ωp2+β2k2.
(3)
Thus, we can write the longitudinal mode penetration depth δL in the bulk metal
δL=1Im(k)=βIm(ω(ω+iγ)ωp2).
(4)
Note that, when ω is higher or β is larger, δL becomes larger. Thus, the nonlocal effect becomes important at high frequencies. The penetration depth is in the order of 0.1nm in noble metals, such as gold or silver [20

20. A. Wiener, A. I. Fernández-Domínguez, A. P. Horsfield, J. B. Pendry, and S. A. Maier, “Nonlocal effects in the nanofocusing performance of plasmonic tips,” Nano Lett. 12(6), 3308–3314 (2012). [CrossRef] [PubMed]

]. Therefore, the theoretical investigation of plasmonic phenomena in this sub-nanometer structure requires the hydrodynamic model.

We follow the approach outlined in Ref [17

17. G. Toscano, S. Raza, A. P. Jauho, N. A. Mortensen, and M. Wubs, “Modified field enhancement and extinction by plasmonic nanowire dimers due to nonlocal response,” Opt. Express 20(4), 4176–4188 (2012). [CrossRef] [PubMed]

]. to solve the above equations and use a highly adaptive mesh through a commercial FEM package, Comsol Multiphysics. We take advantage of the fact that J and E have the same propagation constant. We can write Eq. (2) into a weak form resulting in an integral expression [30

30. P. Monk, Finite Element Methods for Maxwell's Equations (Oxford University Press, 2003).

]
Ω(β2((×J)·(×J˜)(J)·(J˜))+ω(ω+iγ)J·J˜iωωp2ε0E·J˜)dΩ=0,
(5)
where J˜ denotes the vector-valued test function of J. Note that one may obtain a more simplied (but less accurate) integral expression as [21

21. K. R. Hiremath, L. Zschiedrich, and F. Schmidt, “Numerical solution of nonlocal hydrodynamic Drude model for arbitrary shaped nano-plasmonic structures using Nédélec finite elements,” J. Comput. Phys. 231(17), 5890–5896 (2012). [CrossRef]

] Ω(β2(·J)(·J˜)+ω(ω+iγ)J·J˜iωωp2ε0E·J˜)dΩ=0. The additional boundary condition is the continuity of the normal component of E, and the continuity of the normal component of J at the surface due to the finite values of charge and current densities [31

31. A. R. Melnyk and M. J. Harrison, “Theory of optical excitation of plasmons in metals,” Phys. Rev. B 2(4), 835–850 (1970). [CrossRef]

].

3. Validation for benchmark problems

For a validation of the present numerical approach, we simulate a metal cylindrical waveguide and a metal slab waveguide, where the cylindrical symmetry allows for an analytical solution in terms of modified Bessel and Hankel functions [24

24. G. C. Aers, B. V. Paranjape, and A. D. Boardman, “Non-radiative surface plasma-polariton modes of inhomogeneous metal circular cylinders,” J. Phys. F 10(1), 53–65 (1980). [CrossRef]

, 25

25. R. Ruppin, “Effect of non-locality on nanofocusing of surface plasmon field intensity in a conical tip,” Phys. Lett. A 340(1-4), 299–302 (2005). [CrossRef]

], and the slab waveguide allows for an analytical solution in terms of exponential functions [26

26. R. Ruppin, “Non-local optics of the near field lens,” J. Phys. Condens. Matter 17(12), 1803–1810 (2005). [CrossRef]

], both for Drude and nonlocal responses. Figure 1
Fig. 1 The real part of the effective index, the propagation distance and electric field distribution for the metal cylindrical waveguide and the metal slab waveguide: (a1) - (a4) the radius of cylinder a = 2nm, (b1) – (b4) the thickness of slab d = 2nm. All panels show comparisons of numerical simulations to analytical results both for local response (β = 0) and for nonlocal response (β = 0.0036c0). All numerical curves overlap the corresponding analytical curves.
summarizes the results for our benchmark problems, and illustrates the strong dependence of the nonlocal effects on the subwavelength size of the nanostructures. Figure 1 shows the effective index neff, the propagation distance Lprop and the mode profile for the metal cylindrical waveguide and the metal slab waveguide. The propagation distance is given by Lprop = 1/(2Im(neff)k0), where k0 is the wave number in vacuum. We take parameters for silver, namely, the plasma frequency ηωp = 8.59eV, Drude damping coefficient ηγ = 0.075eV, and nonlocal parameter β = 0.0036c0 [20

20. A. Wiener, A. I. Fernández-Domínguez, A. P. Horsfield, J. B. Pendry, and S. A. Maier, “Nonlocal effects in the nanofocusing performance of plasmonic tips,” Nano Lett. 12(6), 3308–3314 (2012). [CrossRef] [PubMed]

] in Fig. 1. However, the interband transition is not considered here, nor the dependence of β on the frequency due to transition absorption.

Figures 1(a1)-(a4) compare the numerical results for nanowires of radius R = 2nm with the exact analytical solution, both for local (β = 0) and nonlocal response (β = 0.0036c0). For these tiny nanowires, the nonlocal effect can be considerable. To give two examples, the relative difference of Lprop between the nonlocal and local responses in the figure can be up to 4.6%. The relative difference of neff between the nonlocal and local responses in the figure can be up to 2.5%. The mode profile in the metal is different due to the fact that the nonlocal parameter β will increase the longitudinal mode penetration depth δL. Effectively, the electric field penetration into the metal is increased at the same time.

The analytical and numerical curves overlap almost completely for the Drude model, and likewise for the nonlocal model. Hence only two of the four curves are visible. The relative errors of the numerically computed Lprop and neff for the nonlocal response are always less than 0.2% and 0.4%, respectively, while for the local response the relative errors are always less than 10−5.

Figures 1(b1)-(b4) show the corresponding four curves as Fig. 1(a1)-(a4), but now for a metal slab waveguide with thickness d = 2nm. As one can see, the maximum relative differences in Lprop and neff mode profile between the nonlocal and local models are 1.8% and 0.9%, respectively. This indicates that the infinity along one dimension of the slab waveguide (which has only one-dimensional confinement) reduces the difference between the local and nonlocal models as compared to the cylinder waveguide, which has a two-dimensional confinement. The maximum relative errors in the numerically computed Lprop and neff for the nonlocal response in Fig. 1(b) are 0.16% and 0.07%, respectively, while for the local response it is always smaller than 10−5.

In summary, our numerical implementation for the metal cylindrical waveguide and the metal slab waveguide is accurate for the distribution of the electromagnetic field in a wide range of wavelength and length scales.

4. Novel silicon hybrid plasmonic waveguide

Having addressed the metal cylindrical waveguide and the metal slab waveguide, we now consider a silicon-based hybrid plasmonic waveguide with an inverted metal nano-rib, which can reduce the mode area to the order of 10−6λ2, while the propagation distance is still several microns simultaneously at telecom wavelengths. The cross section of the present silicon hybrid plasmonic waveguide is shown in Fig. 2(a)
Fig. 2 (a). The cross section of the present hybrid plasmonic waveguide with an inverted metal nano-rib. (b) The distribution of electromagnetic energy density for the present hybrid plasmonic waveguide. (c) The distribution of electromagnetic energy density for the hybrid plasmonic waveguide without the metal nano-rib. (d) The distribution of electromagnetic energy density for the only inverted metal nano-rib. The structure parameters are wco = 22nm, hmetal = 10nm, hrib = 10nm, hSi = 50nm, θ = 10°, g = 0.5nm and R = 0.5nm.
. It consists of a SOI rib with an inverted metal nano-rib on the top. There is a low-index material (e.g., air) between the SOI rib and the metal structure. This novel hybrid plasmonic waveguide combines the advantage of the silicon hybrid plasmonic waveguide structure (which has low propagation loss and strong field confinement in the vertical direction) and the advantage of the WPP waveguide (low propagation loss and strong field confinement in the horizontal direction). When the metal rib height is zero (i.e., a metal slab) in Fig. 2(c), or the SOI rib is absent in Fig. 2(d), the fundamental mode field is not well confined in any of the two directions and the mode area is about 10−4λ2. However, when there is an inverted metal nano-rib, the field distribution of the guide mode is significantly modified and the field is strongly confined at the tip of the metal rib (see Fig. 2(b)) without much decrease in the propagation distance.

Figure 3
Fig. 3 The effective index, mode area, propagation distance and the field distributions. (a) The effective index, (b) Propagation distance and the mode area of HWP. Green (Blue) line shows the Drude (nonlocal) results. (c) and (d): electromagnetic energy density distribution for the Drude and nonlocal models, respectively.
displays the behavior of HPW modes when wco = 22nm, hmetal = 10nm, hrib = 10nm, hSi = 50nm, θ = 10°, g = 0.5nm and R = 0.5nm. The frequency-dependent permittivity of the metal, SiO2 and Si are the same as those in Ref [32

32. D. X. Dai, Y. C. Shi, S. L. He, L. Wosinski, and L. Thylen, “Silicon hybrid plasmonic submicron-donut resonator with pure dielectric access waveguides,” Opt. Express 19(24), 23671–23682 (2011). [CrossRef] [PubMed]

]. We do not take interband transition into consideration here, nor the dependence of β on the frequency due to transition absorption. Panel (a) shows the dispersion relation of the fundamental mode using the Drude method (line) or the nonlocal method (line). The mode has no cutoff wavelength. Panel (b) shows the propagation distance Lprop and the mode area Am. The mode area Am is defined as the ratio of the total mode energy to the peak energy density [9

9. R. F. Oulton, V. J. Sorger, D. A. Genov, D. Pile, and X. Zhang, “A hybrid plasmonic waveguide for subwavelength confinement and long-range propagation,” Nat. Photonics 2(8), 496–500 (2008). [CrossRef]

]
Am=Wmmax{W}=1max{W}WdS,
(6)
where W is the time averaged electromagnetic energy density. In the Drude model, the average energy density in the metal can be expressed as follows [33

33. R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A 299(2-3), 309–312 (2002). [CrossRef]

]
W=14(Re(d(εω)dω)|E|2+μ0|H|2).
(7)
In the nonlocal model the energy conservation law reads [34

34. F. Forstmann and H. Stenschke, “Electrodynamics at metal boundaries with inclusion of plasma waves,” Phys. Rev. Lett. 38(23), 1365–1368 (1977). [CrossRef]

]
t(12ε0E2+12μ0H2+12ωp2ε0(J2+β2ω2(·J)2))=·(E×Hiβωp2ε0ω(·J)·J)+γωp2ε0J2.
(8)
The left-hand term is the time derivative of the time-dependent energy density in the metal structure with a negative sign. Taking an average over time, we have the energy density
W=14(ε0|E|2+1ωp2ε0|J|2+μ0|H|2)+β24ωp2ω2ε0|·J|2.
(9)
While in the absence of the nonlocal term, i.e., β = 0 in the Drude model, Eq. (9) can be converted to Eq. (7) [33

33. R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A 299(2-3), 309–312 (2002). [CrossRef]

]. On the right side, the second term represents the dissipation caused by the Joule heating, with which we can have a deeper insight of the relative energy loss.
p=Smetalγ2ωp2ε0ω|J|2dSWm,
(10)
where p is the ratio of the total dissipation energy to the total mode energy. When λ = 1.55μm, we have p = 2.17% in the nonlocal model and p = 2.18% in the Drude model for this HPW. Here we can see that a relatively lower loss is obtained in the nonlocal model. This agrees with the simulation result that the propagation distance is longer in the nonlocal model than that in the Drude model. An amazing point is that in the meantime the mode area is smaller in the nonlocal model than that in the Drude model because the continuity of the normal component of E in the nonlocal model enables a peak energy density located in the metal side at the interface between the metal and dielectric. When the propagation distance is 2.5μm, the mode area is 4.4 × 10−6λ2 in the traditional Drude model at 1.55μm. The mode area is 2.8 × 10−6λ2 and the propagation distance is 2.6μm in the hydrodynamic model at 1.55μm. To the best of our knowledge, this is the best result that has ever been reported about the mode area (note that the energy density distribution given in Ref [35

35. R. Hao, E. Li, and X. Wei, “Two-dimensional light confinement in cross-index-modulation plasmonic waveguides,” Opt. Lett. 37(14), 2934–2936 (2012). [CrossRef] [PubMed]

]. is not correct because they used the Drude model at the sharp corner and did not use the correct definition of Eq. (7) for the energy density).

Yet we cannot say that the trade-off between the mode confinement and the propagation length has been broken through since not that much improvement lies here. From the Panels 3(a), (b), we can also see, the propagation distance is decreased as λ decreases, while the ratio of the mode area to the λ2 is increased as λ decreases. Panels 3(c), (d) show the electromagnetic energy density distribution at a wavelength of 1.55μm for both the nonlocal model and the Drude model. The difference between the energy density distributions of the nonlocal model and the Drude model is associated with the boundary condition and the changes in the induced-charge distribution. In the Drude model, the charge is strictly a surface charge, while in the nonlocal model the charge density is finite and the longitudinal mode penetration depth δL tends to spatially smear out the charge distribution. Effectively, this smearing increases the electric field penetration into the metal, which causes the energy density strong in the thin layer near the metal surface.

We further explore the influence of the gap on the confinement and loss of this hybrid plasmonic waveguide. We sweep the gap from 0.5nm to 7nm at a fixed wavelength of 1.55μm. In Figs. 4(a)
Fig. 4 (a),(b) The mode area Am and the propagation distance Lprop as g varies from 0.5nm to 7nm at λ = 1.55μm.
and 4(b), when the gap height decreases from 7nm to 3nm, the mode area keeps decreasing while the propagation distance does not change notably. When the gap height is reduced further, the mode area decreases to the order of 10−6λ2. The results of the nonlocal hydrodynamic model and the traditional Drude model have the same tendency. However, the difference between the mode areas simulated by the two models increases as the gap increases.

We also use a dimensionless figure of merit (FoM) to evaluate the confinement and propagation distance. FoM is defined as the ratio of the propagation distance (Lprop) to the effective mode radius (Rm) [36

36. R. Buckley and P. Berini, “Figures of merit for 2D surface plasmon waveguides and application to metal stripes,” Opt. Express 15(19), 12174–12182 (2007). [CrossRef] [PubMed]

]
FoM=2LpropRm=2LpropπAm,
(11)
where Rm is defined as the radius of the effective mode area (Am). The FoM for our novel hybrid plasmonic waveguide can be as high as 2700 (or 4000) in the Drude (or hydrodynamic) model. While the FoMs of previously studied wedge plasmonic waveguide, cylinder plasmonic waveguide, channel plasmonic waveguide and the hybrid plasmonic waveguide in [37

37. R. F. Oulton, G. Bartal, D. F. Pile, and X. Zhang, “Confinement and propagation characteristics of subwavelength plasmonic modes,” New J. Phys. 10(10), 105018 (2008). [CrossRef]

] are just around 1500 at the smallest mode area of 2.5 × 10−3 λ2 at the telecom wavelength of 1.55μm in the Drude model. Obviously, our novel HWP can reduce the mode area and in the mean time improve the FoM.

The finite value of the charge density and the longitude plasmonic mode also allow resolving the field in the proximity of very sharp corners and tips, which is similar to the finite SER enhancement factor for an infinitely sharp tip, as discussed in Ref [19

19. G. Toscano, S. Raza, S. Xiao, M. Wubs, A. P. Jauho, S. I. Bozhevolnyi, and N. A. Mortensen, “Surface-enhanced Raman spectroscopy: nonlocal limitations,” Opt. Lett. 37(13), 2538–2540 (2012). [CrossRef] [PubMed]

]. In Fig. 5(a)
Fig. 5 (a) The mode area Am as R varies from 1 nm to 0.005nm at λ = 1.55μm. (b) Normalized energy density along x = 0 for both the Drude model (green) and the nonlocal model (blue) at R = 0.005nm. The shaded grey and brown areas represent the silicon and metal regions, respectively. It shows that the energy density distribution has a sharp peak at the rib’s tip in the Drude model, while such a phenomenon does not exist in the nonlocal model.
, we reduce R from 1nm down to 0.005nm and see the difference between the nonlocal model and the Drude model. As expected, the difference between the two models increase as the tip radius R approaches the length scale of the longitudinal mode penetration depth δL (~0.1nm). In the Drude model, when R is approaching zero, the mode area monotonously decreases. When r is 0.005nm, the mode area approaches 10−8λ and a sharp peak appears in the energy density distribution as shown in Fig. 5(b). However, there is a fundamental saturation of the mode area (Am = 3 × 10−7λ2 at r = 0.005nm) in the nonlocal model, as the nonlocal effect relaxes the sharpness of the tip and the energy density distribution is smeared out as shown in Fig. 5(b).

5. Conclusions

In conclusion, we use finite element method for implementing the hydrodynamic model to analyze the propagation mode of a plasmonic waveguide of arbitrary cross section. To the best of our knowledge, this is the first time that the nonlocal effect is numerically implemented for the mode analysis of a plasmonic waveguide. We have proposed a SOI hybrid plasmonic waveguide with an inverted metal nano-rib, which enables nanoscale light confinement. Our theoretical investigation has shown that the mode area could be of 10−6λ2 size in a wide wavelength range from 1.3μm to 2μm, and the propagation distance still remains several microns. We have studied theoretically how the nonlocal mode affects this new hybrid plasmonic waveguide. We have found that the nonlocal effects can reduce the loss and improve the confinement. We have shown that when the radius of the rib’s tip is approaching zero, the mode area gets saturated in the nonlocal model, instead of monotonous decrease in the Drude model.

Acknowledgment

We thank Peipeng Xu and Yingchen Wu for valuable discussions. This work is partially supported by the National Natural Science Foundation of China (Nos. 60990322 and 61178062), the National High Technology Research and Development Program (863) of China (No. 2012AA012201) and Swedish VR (No. 621-2011-4620).

References and links

1.

L. Liu, Z. H. Han, and S. L. He, “Novel surface plasmon waveguide for high integration,” Opt. Express 13(17), 6645–6650 (2005). [CrossRef] [PubMed]

2.

Z. Han, A. Y. Elezzabi, and V. Van, “Experimental realization of subwavelength plasmonic slot waveguides on a silicon platform,” Opt. Lett. 35(4), 502–504 (2010). [CrossRef] [PubMed]

3.

D. F. Pile and D. K. Gramotnev, “Channel plasmon-polariton in a triangular groove on a metal surface,” Opt. Lett. 29(10), 1069–1071 (2004). [CrossRef] [PubMed]

4.

S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, J. Y. Laluet, and T. W. Ebbesen, “Channel plasmon subwavelength waveguide components including interferometers and ring resonators,” Nature 440(7083), 508–511 (2006). [CrossRef] [PubMed]

5.

D. Pile, T. Ogawa, D. K. Gramotnev, T. Okamoto, M. Haraguchi, M. Fukui, and S. Matsuo, “Theoretical and experimental investigation of strongly localized plasmons on triangular metal wedges for subwavelength waveguiding,” Appl. Phys. Lett. 87(6), 061106–061103 (2005). [CrossRef]

6.

T. Ogawa, D. Pile, T. Okamoto, M. Haraguchi, M. Fukui, and D. K. Gramotnev, “Numerical and experimental investigation of wedge tip radius effect on wedge plasmons,” J. Appl. Phys. 104(3), 033102–033106 (2008). [CrossRef]

7.

E. Moreno, S. G. Rodrigo, S. I. Bozhevolnyi, L. Martín-Moreno, and F. J. García-Vidal, “Guiding and focusing of electromagnetic fields with wedge plasmon polaritons,” Phys. Rev. Lett. 100(2), 023901 (2008). [CrossRef] [PubMed]

8.

A. Boltasseva, V. S. Volkov, R. B. Nielsen, E. Moreno, S. G. Rodrigo, and S. I. Bozhevolnyi, “Triangular metal wedges for subwavelength plasmon-polariton guiding at telecom wavelengths,” Opt. Express 16(8), 5252–5260 (2008). [CrossRef] [PubMed]

9.

R. F. Oulton, V. J. Sorger, D. A. Genov, D. Pile, and X. Zhang, “A hybrid plasmonic waveguide for subwavelength confinement and long-range propagation,” Nat. Photonics 2(8), 496–500 (2008). [CrossRef]

10.

D. X. Dai and S. L. He, “A silicon-based hybrid plasmonic waveguide with a metal cap for a nano-scale light confinement,” Opt. Express 17(19), 16646–16653 (2009). [CrossRef] [PubMed]

11.

D. X. Dai and S. L. He, “Low-loss hybrid plasmonic waveguide with double low-index nano-slots,” Opt. Express 18(17), 17958–17966 (2010). [CrossRef] [PubMed]

12.

Y. S. Bian, Z. Zheng, Y. Liu, J. S. Liu, J. S. Zhu, and T. Zhou, “Hybrid wedge plasmon polariton waveguide with good fabrication-error-tolerance for ultra-deep-subwavelength mode confinement,” Opt. Express 19(23), 22417–22422 (2011). [CrossRef] [PubMed]

13.

D. X. Dai, Y. C. Shi, S. L. He, L. Wosinski, and L. Thylen, “Gain enhancement in a hybrid plasmonic nano-waveguide with a low-index or high-index gain medium,” Opt. Express 19(14), 12925–12936 (2011). [CrossRef] [PubMed]

14.

A. D. Boardman, Electromagnetic Surface Modes (Wiley, 1982).

15.

J. M. McMahon, S. K. Gray, and G. C. Schatz, “Nonlocal optical response of metal nanostructures with arbitrary shape,” Phys. Rev. Lett. 103(9), 097403 (2009). [CrossRef] [PubMed]

16.

J. M. McMahon, S. K. Gray, and G. C. Schatz, “Calculating nonlocal optical properties of structures with arbitrary shape,” Phys. Rev. B 82(3), 035423 (2010). [CrossRef]

17.

G. Toscano, S. Raza, A. P. Jauho, N. A. Mortensen, and M. Wubs, “Modified field enhancement and extinction by plasmonic nanowire dimers due to nonlocal response,” Opt. Express 20(4), 4176–4188 (2012). [CrossRef] [PubMed]

18.

S. Raza, G. Toscano, A. P. Jauho, M. Wubs, and N. A. Mortensen, “Unusual resonances in nanoplasmonic structures due to nonlocal response,” Phys. Rev. B 84(12), 121412 (2011). [CrossRef]

19.

G. Toscano, S. Raza, S. Xiao, M. Wubs, A. P. Jauho, S. I. Bozhevolnyi, and N. A. Mortensen, “Surface-enhanced Raman spectroscopy: nonlocal limitations,” Opt. Lett. 37(13), 2538–2540 (2012). [CrossRef] [PubMed]

20.

A. Wiener, A. I. Fernández-Domínguez, A. P. Horsfield, J. B. Pendry, and S. A. Maier, “Nonlocal effects in the nanofocusing performance of plasmonic tips,” Nano Lett. 12(6), 3308–3314 (2012). [CrossRef] [PubMed]

21.

K. R. Hiremath, L. Zschiedrich, and F. Schmidt, “Numerical solution of nonlocal hydrodynamic Drude model for arbitrary shaped nano-plasmonic structures using Nédélec finite elements,” J. Comput. Phys. 231(17), 5890–5896 (2012). [CrossRef]

22.

C. Ciracì, R. T. Hill, J. J. Mock, Y. Urzhumov, A. I. Fernández-Domínguez, S. A. Maier, J. B. Pendry, A. Chilkoti, and D. R. Smith, “Probing the ultimate limits of plasmonic enhancement,” Science 337(6098), 1072–1074 (2012). [CrossRef] [PubMed]

23.

A. I. Fernández-Domínguez, A. Wiener, F. J. García-Vidal, S. A. Maier, and J. B. Pendry, “Transformation-optics description of nonlocal effects in plasmonic nanostructures,” Phys. Rev. Lett. 108(10), 106802 (2012). [CrossRef] [PubMed]

24.

G. C. Aers, B. V. Paranjape, and A. D. Boardman, “Non-radiative surface plasma-polariton modes of inhomogeneous metal circular cylinders,” J. Phys. F 10(1), 53–65 (1980). [CrossRef]

25.

R. Ruppin, “Effect of non-locality on nanofocusing of surface plasmon field intensity in a conical tip,” Phys. Lett. A 340(1-4), 299–302 (2005). [CrossRef]

26.

R. Ruppin, “Non-local optics of the near field lens,” J. Phys. Condens. Matter 17(12), 1803–1810 (2005). [CrossRef]

27.

F. J. García de Abajo, “Nonlocal effects in the plasmons of strongly interacting nanoparticles, dimers, and waveguides,” J. Phys. Chem. C 112(46), 17983–17987 (2008). [CrossRef]

28.

D. C. Marinica, A. K. Kazansky, P. Nordlander, J. Aizpurua, and A. G. Borisov, “Quantum plasmonics: Nonlinear effects in the field enhancement of a plasmonic nanoparticle dimer,” Nano Lett. 12(3), 1333–1339 (2012). [CrossRef] [PubMed]

29.

R. Esteban, A. G. Borisov, P. Nordlander, and J. Aizpurua, “Bridging quantum and classical plasmonics with a quantum-corrected model,” Nat Commun 3, 825 (2012). [CrossRef] [PubMed]

30.

P. Monk, Finite Element Methods for Maxwell's Equations (Oxford University Press, 2003).

31.

A. R. Melnyk and M. J. Harrison, “Theory of optical excitation of plasmons in metals,” Phys. Rev. B 2(4), 835–850 (1970). [CrossRef]

32.

D. X. Dai, Y. C. Shi, S. L. He, L. Wosinski, and L. Thylen, “Silicon hybrid plasmonic submicron-donut resonator with pure dielectric access waveguides,” Opt. Express 19(24), 23671–23682 (2011). [CrossRef] [PubMed]

33.

R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A 299(2-3), 309–312 (2002). [CrossRef]

34.

F. Forstmann and H. Stenschke, “Electrodynamics at metal boundaries with inclusion of plasma waves,” Phys. Rev. Lett. 38(23), 1365–1368 (1977). [CrossRef]

35.

R. Hao, E. Li, and X. Wei, “Two-dimensional light confinement in cross-index-modulation plasmonic waveguides,” Opt. Lett. 37(14), 2934–2936 (2012). [CrossRef] [PubMed]

36.

R. Buckley and P. Berini, “Figures of merit for 2D surface plasmon waveguides and application to metal stripes,” Opt. Express 15(19), 12174–12182 (2007). [CrossRef] [PubMed]

37.

R. F. Oulton, G. Bartal, D. F. Pile, and X. Zhang, “Confinement and propagation characteristics of subwavelength plasmonic modes,” New J. Phys. 10(10), 105018 (2008). [CrossRef]

OCIS Codes
(130.2790) Integrated optics : Guided waves
(240.6680) Optics at surfaces : Surface plasmons
(250.5300) Optoelectronics : Photonic integrated circuits
(260.3910) Physical optics : Metal optics
(160.4236) Materials : Nanomaterials

ToC Category:
Optics at Surfaces

History
Original Manuscript: November 1, 2012
Revised Manuscript: December 16, 2012
Manuscript Accepted: January 2, 2013
Published: January 14, 2013

Citation
Qiangsheng Huang, Fanglin Bao, and Sailing He, "Nonlocal effects in a hybrid plasmonic waveguide for nanoscale confinement," Opt. Express 21, 1430-1439 (2013)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-21-2-1430


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. L. Liu, Z. H. Han, and S. L. He, “Novel surface plasmon waveguide for high integration,” Opt. Express13(17), 6645–6650 (2005). [CrossRef] [PubMed]
  2. Z. Han, A. Y. Elezzabi, and V. Van, “Experimental realization of subwavelength plasmonic slot waveguides on a silicon platform,” Opt. Lett.35(4), 502–504 (2010). [CrossRef] [PubMed]
  3. D. F. Pile and D. K. Gramotnev, “Channel plasmon-polariton in a triangular groove on a metal surface,” Opt. Lett.29(10), 1069–1071 (2004). [CrossRef] [PubMed]
  4. S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, J. Y. Laluet, and T. W. Ebbesen, “Channel plasmon subwavelength waveguide components including interferometers and ring resonators,” Nature440(7083), 508–511 (2006). [CrossRef] [PubMed]
  5. D. Pile, T. Ogawa, D. K. Gramotnev, T. Okamoto, M. Haraguchi, M. Fukui, and S. Matsuo, “Theoretical and experimental investigation of strongly localized plasmons on triangular metal wedges for subwavelength waveguiding,” Appl. Phys. Lett.87(6), 061106–061103 (2005). [CrossRef]
  6. T. Ogawa, D. Pile, T. Okamoto, M. Haraguchi, M. Fukui, and D. K. Gramotnev, “Numerical and experimental investigation of wedge tip radius effect on wedge plasmons,” J. Appl. Phys.104(3), 033102–033106 (2008). [CrossRef]
  7. E. Moreno, S. G. Rodrigo, S. I. Bozhevolnyi, L. Martín-Moreno, and F. J. García-Vidal, “Guiding and focusing of electromagnetic fields with wedge plasmon polaritons,” Phys. Rev. Lett.100(2), 023901 (2008). [CrossRef] [PubMed]
  8. A. Boltasseva, V. S. Volkov, R. B. Nielsen, E. Moreno, S. G. Rodrigo, and S. I. Bozhevolnyi, “Triangular metal wedges for subwavelength plasmon-polariton guiding at telecom wavelengths,” Opt. Express16(8), 5252–5260 (2008). [CrossRef] [PubMed]
  9. R. F. Oulton, V. J. Sorger, D. A. Genov, D. Pile, and X. Zhang, “A hybrid plasmonic waveguide for subwavelength confinement and long-range propagation,” Nat. Photonics2(8), 496–500 (2008). [CrossRef]
  10. D. X. Dai and S. L. He, “A silicon-based hybrid plasmonic waveguide with a metal cap for a nano-scale light confinement,” Opt. Express17(19), 16646–16653 (2009). [CrossRef] [PubMed]
  11. D. X. Dai and S. L. He, “Low-loss hybrid plasmonic waveguide with double low-index nano-slots,” Opt. Express18(17), 17958–17966 (2010). [CrossRef] [PubMed]
  12. Y. S. Bian, Z. Zheng, Y. Liu, J. S. Liu, J. S. Zhu, and T. Zhou, “Hybrid wedge plasmon polariton waveguide with good fabrication-error-tolerance for ultra-deep-subwavelength mode confinement,” Opt. Express19(23), 22417–22422 (2011). [CrossRef] [PubMed]
  13. D. X. Dai, Y. C. Shi, S. L. He, L. Wosinski, and L. Thylen, “Gain enhancement in a hybrid plasmonic nano-waveguide with a low-index or high-index gain medium,” Opt. Express19(14), 12925–12936 (2011). [CrossRef] [PubMed]
  14. A. D. Boardman, Electromagnetic Surface Modes (Wiley, 1982).
  15. J. M. McMahon, S. K. Gray, and G. C. Schatz, “Nonlocal optical response of metal nanostructures with arbitrary shape,” Phys. Rev. Lett.103(9), 097403 (2009). [CrossRef] [PubMed]
  16. J. M. McMahon, S. K. Gray, and G. C. Schatz, “Calculating nonlocal optical properties of structures with arbitrary shape,” Phys. Rev. B82(3), 035423 (2010). [CrossRef]
  17. G. Toscano, S. Raza, A. P. Jauho, N. A. Mortensen, and M. Wubs, “Modified field enhancement and extinction by plasmonic nanowire dimers due to nonlocal response,” Opt. Express20(4), 4176–4188 (2012). [CrossRef] [PubMed]
  18. S. Raza, G. Toscano, A. P. Jauho, M. Wubs, and N. A. Mortensen, “Unusual resonances in nanoplasmonic structures due to nonlocal response,” Phys. Rev. B84(12), 121412 (2011). [CrossRef]
  19. G. Toscano, S. Raza, S. Xiao, M. Wubs, A. P. Jauho, S. I. Bozhevolnyi, and N. A. Mortensen, “Surface-enhanced Raman spectroscopy: nonlocal limitations,” Opt. Lett.37(13), 2538–2540 (2012). [CrossRef] [PubMed]
  20. A. Wiener, A. I. Fernández-Domínguez, A. P. Horsfield, J. B. Pendry, and S. A. Maier, “Nonlocal effects in the nanofocusing performance of plasmonic tips,” Nano Lett.12(6), 3308–3314 (2012). [CrossRef] [PubMed]
  21. K. R. Hiremath, L. Zschiedrich, and F. Schmidt, “Numerical solution of nonlocal hydrodynamic Drude model for arbitrary shaped nano-plasmonic structures using Nédélec finite elements,” J. Comput. Phys.231(17), 5890–5896 (2012). [CrossRef]
  22. C. Ciracì, R. T. Hill, J. J. Mock, Y. Urzhumov, A. I. Fernández-Domínguez, S. A. Maier, J. B. Pendry, A. Chilkoti, and D. R. Smith, “Probing the ultimate limits of plasmonic enhancement,” Science337(6098), 1072–1074 (2012). [CrossRef] [PubMed]
  23. A. I. Fernández-Domínguez, A. Wiener, F. J. García-Vidal, S. A. Maier, and J. B. Pendry, “Transformation-optics description of nonlocal effects in plasmonic nanostructures,” Phys. Rev. Lett.108(10), 106802 (2012). [CrossRef] [PubMed]
  24. G. C. Aers, B. V. Paranjape, and A. D. Boardman, “Non-radiative surface plasma-polariton modes of inhomogeneous metal circular cylinders,” J. Phys. F10(1), 53–65 (1980). [CrossRef]
  25. R. Ruppin, “Effect of non-locality on nanofocusing of surface plasmon field intensity in a conical tip,” Phys. Lett. A340(1-4), 299–302 (2005). [CrossRef]
  26. R. Ruppin, “Non-local optics of the near field lens,” J. Phys. Condens. Matter17(12), 1803–1810 (2005). [CrossRef]
  27. F. J. García de Abajo, “Nonlocal effects in the plasmons of strongly interacting nanoparticles, dimers, and waveguides,” J. Phys. Chem. C112(46), 17983–17987 (2008). [CrossRef]
  28. D. C. Marinica, A. K. Kazansky, P. Nordlander, J. Aizpurua, and A. G. Borisov, “Quantum plasmonics: Nonlinear effects in the field enhancement of a plasmonic nanoparticle dimer,” Nano Lett.12(3), 1333–1339 (2012). [CrossRef] [PubMed]
  29. R. Esteban, A. G. Borisov, P. Nordlander, and J. Aizpurua, “Bridging quantum and classical plasmonics with a quantum-corrected model,” Nat Commun3, 825 (2012). [CrossRef] [PubMed]
  30. P. Monk, Finite Element Methods for Maxwell's Equations (Oxford University Press, 2003).
  31. A. R. Melnyk and M. J. Harrison, “Theory of optical excitation of plasmons in metals,” Phys. Rev. B2(4), 835–850 (1970). [CrossRef]
  32. D. X. Dai, Y. C. Shi, S. L. He, L. Wosinski, and L. Thylen, “Silicon hybrid plasmonic submicron-donut resonator with pure dielectric access waveguides,” Opt. Express19(24), 23671–23682 (2011). [CrossRef] [PubMed]
  33. R. Ruppin, “Electromagnetic energy density in a dispersive and absorptive material,” Phys. Lett. A299(2-3), 309–312 (2002). [CrossRef]
  34. F. Forstmann and H. Stenschke, “Electrodynamics at metal boundaries with inclusion of plasma waves,” Phys. Rev. Lett.38(23), 1365–1368 (1977). [CrossRef]
  35. R. Hao, E. Li, and X. Wei, “Two-dimensional light confinement in cross-index-modulation plasmonic waveguides,” Opt. Lett.37(14), 2934–2936 (2012). [CrossRef] [PubMed]
  36. R. Buckley and P. Berini, “Figures of merit for 2D surface plasmon waveguides and application to metal stripes,” Opt. Express15(19), 12174–12182 (2007). [CrossRef] [PubMed]
  37. R. F. Oulton, G. Bartal, D. F. Pile, and X. Zhang, “Confinement and propagation characteristics of subwavelength plasmonic modes,” New J. Phys.10(10), 105018 (2008). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1 Fig. 2 Fig. 3
 
Fig. 4 Fig. 5
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited