Surface Plasmon-polaritons on thin metal films - IMI (insulator-metal-insulator) structure -

Transcription

1 Surface Plasmon-polaritons on thin metal films - IMI (insulator-metal-insulator) structure - Dielectric 3 Metal 2 Dielectric 1

2 References Surface plasmons in thin films, E.N. Economou, Phy. Rev. Vol.182, (1969) Surface-polariton-like waves guided by thin, lossy metal films, J.J. Burke, G. I. Stegeman, T. Tamir, Phy. Rev. B, Vol.33, (1986) Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of symmetric structures, P. Berini, Phy. Rev. B, Vol.61, (2000) Geometries and materials for subwavelength surface plasmon modes, Rashid Zia, Mark D. Selker, Peter B. Catrysse, and Mark L. Brongersma, J. Opt. Soc. Am. A, Vol. 21, (2004) Plasmon slot waveguides: Towards chip-scale propagation with subwavelength-scale localization, J. A. Dionne,* L. A. Sweatlock, and H. A. Atwater, Phy. Rev. B, Vol.73, (2006)

3 Introduction: When the film thickness becomes finite. mode overlap

4 Introduction: Possibility of Propagation Range Extension Long-Range SP: weak surface confinement, low loss frequency Short-Range SP: strong surface confinement, high loss in-plane wavevector

5 Introduction: Extremely long-range SP? Symmetrically coupled LRSP frequency Anti-symmetrically coupled LRSP in-plane wavevector

6 Introduction: Dependence of dispersion on film thickness h = 10 60nm practically forbidden

7 Surface-polariton-like waves guided by thin, lossy metal films J.J. Burke, G. I. Stegeman, T. Tamir, Phy. Rev. B, Vol.33, (1986). Dispersion relations for waves guided by a thin, lossy metal film surrounded by dielectric media Characteristic of "spatial transients" : Usual symmetric and antisymmetric branches each split into a pair of waves one radiative (leaky waves) and the other nonradiative (bound waves). Symmetric modes : the transverse electric field does not exhibit a zero inside the metal film Antisymmetric modes : the transverse electric field has a zero inside the film. h 1 m = - R i I 3 z x 1 > 3 Burke, PRB 1986

8 Dispersion relation for thin metal films (3 layers) obtained from the Maxwell equations Hi( x, z, t) = e yh0 fi( z)exp i( βx ωt) [ 3 ] ( ) [ 2 ] 0 [ 2 ] ( ) [ sz] ( z ) Bexp s ( z h) in medium 3 z h fi( z) = Ahexp s ( z h) + A exp s z in medium 2 0 z h exp 1 in medium 1 0 s = β k j j j 0 i H y i df j ( z) Ex = = H0 exp i x ω z ω dz i E= H Ey = 0 ω β β Ez = Hy = H0 f j ( z) exp[ iβ x] ω ω [ β ] s3 B exp [ s 3( z h )] ( z h ) 3 i s2 Ex = H0exp[ iβ x] A exp s2( z h) A0exp s2z 0 z h ω 2 s1 exp[ sz 1 ] ( z 0) 1 { h [ ] [ ]} ( )

9 From the boundary conditions, ( ) ( ) z = 0: Hx 1 = Hx2 exp[ s2h] Ah + A0 = 1 z = h: Hx2 = Hx3 Ah + exp[ s2h] A0 = B 2s1 z = 0 : Ex 1 = Ex2 exp[ s2h] Ah A0 = 1s2 s = = = 2 3 z h: Ex2 Ex3 Ah exp[ s2h] A0 B 3s2 From the equations at z = 0, A h, A o, and B can be determined by, 2s 1 cosh[ sh 2 ] + sinh[ sh 2 ] exp [ s3( z h) ] ( j= 3: z h) 1s2 s f j ( z) = cosh s2z + sinh s2z j = 2 : 0 z h 1s2 exp[ sz 1 ] ( j= 1: z 0) 2 1 [ ] [ ] ( ) Therefore, anther equations at z = h gives the dispersion relation, 2 2 ( 1 3s2 2s1s3) [ s2h] 2s2( 3s1 1s3) + tanh + + = 0 s = β k j j j 0

10 Dispersion relation when h When h >> c/ω p (classical skin depth), tanh(s 2 h) 1, ( ) ω i m βi= 1,3 h = c + i m The solutions consist of decoupled surface-plasmon polaritons (SPP) : SP 1 : propagating along the 1 - m interface SP 3 : propagating along the 3 - m interface If we assume that β >> β and = + i R I m R I m 2 ω p = 1 2 ω Burke, PRB 1986

11 If β > ω/ c (light line), i S = β k > i i o [ 3 ] ( ) [ 2 ] 0 [ 2 ] ( ) [ sz] ( z ) Bexp s ( z h) in medium 3 z h fi( z) = Ahexp s ( z h) + A exp s z in medium 2 0 z h exp 1 in medium 1 0 Hence S i is real. But, there are two types of solutions for the semi-infinite media SPPs : (1) S i > 0, SPPs are nonradiative, or bound (2) S i < 0, SPPs are grow exponentially with distance from the interface, which are physically rejected because of their non-guiding property Therefore, SPPs, are bound at the semi-infinite media only when S i are real and positive. For a finite film thickness, the two allowed semi-infinite SPPs are coupled. One of the SPPs could become leaky (radiative) in the 3 medium when 1 > 3 If β < ω/ c (light line), i S < 0 2 i Hence S i is imaginary. The field is a plane wave radiating away form the metal boundary. Burke, PRB 1986

12 Waves guided by symmetric structures ( = ) 1 3 There are four types of solutions satisfying : S = β k > i i o We can estimate the properties of the solutions as follows: h h 0 ± S1 =± S3 (1) : S1 > 0 & S3 > 0 S1 S3 (2): < 0 & < 0 ± S1 = S3 (3) : S1 > 0 & S3 < 0 S1 S3 (4): < 0 & > 0 (1) : S > 0 & S > The fields decay exponentially into both 1 and 3, and the wave fronts are tilted in towards the metal film. -> nonradiative waves There are two types of solutions for S i > 0 : symmetric (s), antisymmetric (a) (2): S < 0 & S < SP solutions No solution The fields grow exponentially with wave front tilted to carry energy away from the metal 1 3 m H y Leaky ( 1 ) Leaky ( 3 ) 1 m 3 H y Bound (s) 1 m 3 H y Bound (a) Two nonradiative, Fano modes (3) : S > 0 & S < SP solution The field in 1 and the metal is guided by the interface, the filed in 3 grows exponentially (leaky). 1 m 3 Bound ( 1 ) Leaky ( 3 ) Four SPP modes (4): S < 0 & S > SP solution The field in 3 and the metal is guided by the interface, the filed in 1 grows exponentially (leaky). 1 m 3 H y Leaky ( 1 ) Bound ( 3 ) H y Two radiative, leaky modes Burke, PRB 1986

13 Two nonradiative, Fano modes 1 H y 1 H y (1) : S > 0 & S > m 3 Bound (s) m 3 Bound (a) Symmetric bound (s b ) Asymmetric bound (a b )

14 Surface plasmon dispersion for thin films Two modes appear Drude model 1 (ω)=1-(ω p /ω) 2 L - (asymmetric) L + (symmetric) Thinner film: Shorter SP wavelength Example: λ HeNe = 633 nm λ SP = 60 nm L - Propagation lengths: cm!!! (infrared)

15 Ex) λ = 1.55 μm, = i11.58, = = ( 2 2 ) 1 3s2 2s1s3 [ s2h] 2s2( 3s1 1s3) + tanh + + = 0 s b a b s = β k j j j 0 1 m H y 1 m H y 3 Bound (s) 3 Bound (a) a b 0.1 a b β r /k s b β i /k E-3 1E-4 s b thickness (h : nm) 1E-5 1E-6 1E thickness (h : nm)

16 Waves guided by asymmetric structures ( ) 1 3 There are also four types of solutions satisfying ± S =± S & ± S = S One antisymmetric mode is always obtained. The "symmetric" solutions are of two types, nonradiative (bound) and nonradiative (leaky): S = β k > i i o Bound and Leaky modes when Si β iko 0 1 m H y Bound (s) (1) : S > 0 & S > m 3 = > Growing modes H y Bound (a) (2): S < 0 & S < 0 (3) : S1 > 0 & S3 < 0 (4): S1 < 0 & S3 > 0 For example, S < 0 & β > 0 1 m 3 1R 3I when S = S - is and β = β -iβ i ir ii i ir ii H y Growing (z) z x Growing (x) 1 m H y Leaky ( 1 ) 1 m Bound ( 1 ) 1 m H y Leaky ( 1 ) 3 Leaky ( 3 ) 3 Leaky ( 3 ) H y 3 Bound ( 3 ) Burke, PRB 1986

17

18 Waves guided by asymmetric structures ( ) 1 3 Nonradiative mode the fields in the dielectric decay exponentially away from the film and the wave fronts are tilted into the metal film, in order to remove energy from the dielectric media (for dissipation in the metal) as the wave attenuates with propagation distance. Leaky (radiative) mode The wave energy is localized in one of the dielectrics, say 1, at that dielectric-metal interface. The wave amplitude decays exponentially across the film, and then grows exponentially into the other dielectric medium, 3 in this case. In the 1 medium, the wave fronts are tilted towards the film to supply energy from 1 for both dissipation in the metal, and radiation into 3. Leaky (radiative) mode The analogous case of localization in 3, and radiation into 1. Growing mode The field amplitude grows both with propagation distance as exp(β I x), and into one of the dielectrics as exp(s R z). Since the wave-front tilt is into the film (as opposed to away from it for leaky waves), these waves are dependent on externally incident fields supplying energy to make the total wave amplitude grow. Burke, PRB 1986

19 Leaky waves They only have meaning in a limited region of space above the film and require some transverse plane (say x=o) containing an effective source that launches a localized wave in one dielectric near its metal-dielectric boundary. The field decays across the metal and couples to radiation fields in the opposite dielectric. The 1 field amplitude grows exponentially for only a finite distance z, In this sense, the solutions do not violate boundary conditions as z - infinite. For fields radiated at an angle θ relative to the surface, the angular spectrum of the radiated plane waves is Thus, the radiated power is The wave attenuation due to radiation loss can be estimated from the solutions by calculating the Poynting vector for energy leaving, for example, the 1 - m boundary, per meter of wave front, θ Burke, PRB 1986

20 Surface plasmons in thin, multilayer films E.N. Economou, Phy. Rev. Vol.182, (1969) LOCAL THEORY FOR MULTIPLE-FILM SYSTEM Maxwell equations with a local current-field relation as follows: Local approximation to the current-field relation with The local approximation satisfies the dielectric function in the metal, For electric (or TM) waves that H z = H x = E y = 0 in all of the media, the solution for any component of the fields can thus be represented in the form

21 A. Single Metal-Dielectric Interface The solution for E x that remains finite at infinity is The continuity of E and H fields across the boundary gives the dispersion relation as If the dielectric constant for the insulator i =1, ( ) ( ) / 1 / metal insulator K R K = m d x m d k c ω = ) (1 ) ( p d d p sp x c k k ω ω ω ω ω + = =

22 Retardation (radiative loss) It is obvious that retardation effects are important for q = (k/k p ) < 1 and that they do not play any role for q>>1. Therefore, we are interested in region IV. Dispersion of bulk plasmon: 2 c c ω p ω = ω =±, ω =± ω + if m = 1 2 ω k k p c k m m

23 B. Insulating Film between Two Semi-Infinite Metals Branches I and II are adequately described by the longitudinal electrostatic theory when k >> k p, If i =1, the low-frequency k < k p, part of I, If k p d i << l, as is usually the case, branch III is If k p d i << l, the k << k p portion of II is and, at k = 0

24 C. Metal Film in Vacuum R ( K / ) ( K / ) metal insulator K 2 2 ω im, im, = k 2 c * Burke, PRB 1986 : 2 2 S = k k 2, n= 1, 2, 3. n n o When k < k p, for branch II (antisymmetric oscillation) When k < k p, for branch I (symmetric oscillation) When k > k p, both I and II approach

25 D. Swihart's Geometry When k << k p, for branch I Branch II just R=1, corresponds to oscillations on the external interface (insulator/metal). Branch III, The intersection of branch III with the line ω=ck occurs, when Branch IV behaves for small k as

26 E. Two Metal Films of Different Thicknesses Where, When k << k p, branch II, III

27 The branch II and III Where, Branch II Branch III Branch II Branch III

28 F. Two Dielectric Films For k<<k p branches I and II are given by In the special case when d 1 =d 2 =d i they correspond to symmetric and antisymmetric solutions, And, in the limit when Branch III starts, when k= 0,

29 Branch III starts, when k= 0, For k<k p Branch IV is, for small k, It corresponds to an oscillation which couples the two junctions, and, when one of the thicknesses of the dielectric films becomes large, the coupling is broken and the oscillation is confined to one of them.

30 G. Three Metal Films When k<<k p branches I and II are given by When k<<k p branches III and IV: When k<<k p branches V and VI

31 H. Periodic structure of alternating metal and insulating films Periodicity of alternating thicknesses d m and d i, implies that the eigensolutions obey the Floquet-Bloch theorem; namely, The secular equation is Solutions are found in the shaded regions. When k<<k p curves III and IV can be taken as straight lines with phase velocities, Any intermediate solution has phase velocity

32 The upper region within which solutions lie is bounded by I, a portion of IIa, and IIb. If (k p 2 d i d m )<< 1, Curve I is while curve IIa is Solutions near I are given by When k>>k p the secular equation is

33 Conclusions of Economou Modes of SPO in multiple-film systems can be classified into two main groups. One group contains those modes whose dispersion curves start from zero frequency at k=o, increase as k increases, but remain below the line ω= ck. The other group starts at k= 0 from ω= ω p, or a value slightly less than ω p, and remains close to the line ω= ω p. For very large k, all the dispersion curves of both groups converge asymptotically to the classical surface plasmon frequency ω p /root(2). In addition to these two groups, some uninteresting modes may appear with dispersion curves that lie just below the curve ω 2 = ω p2 + c 2 k 2 which corresponds to the trivial solution of zero fields. For normal metals this description is valid only for high enough frequencies so that oscillation damping is negligible. On the other hand, for superconducting metals, the picture is valid not only for the high-frequency region but also for low frequencies. In Multiple-film structures radiative SPO exists, which should have observable effects in the radiation properties of these structures. In particular, there seems to be a possibility of obtaining intense radiation as the number of the films increases.

34 Dispersion relation for thin metal strips with finite widths metal strip dielectric

35 Finite film thickness and width Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of symmetric structures, P. Berini, Phy. Rev. B, Vol.61, (2000)

36 E( xyzt,,, ) = E ( xye, ) H( xyzt,,, ) = H ( xye, ) 0 0 i( βz ωt) i( β z ωt) E o, H o : polarization direction z-axis : propagation direction From the Maxwell equations ( ) = k 1 2 r H 0 H Assume that all media be isotropic. The magnetic field on x-y (transverse plane) satisfies ( H ) ( H ) ( k β ) H = t r t t r t t t 0 r t where = t i + j x y H = ( H i + H j) e i β z ωt t x y ( ) This eigenvalue problem can be solved by a numerical method with proper boundary conditions, such as one of FDM, FEM, MoL, Here, we use the FDM (finite difference method).

37 y

38 Ex) λ = 1.55 μm, = i11.58, = = 2.25, w = 5μm FDM z y x x y sa b o ss b o E-3 β r /k o β i /k 0 1E-4 1E-5 1E-6 sa b o ss b o thickness of metal (nm) 1E thickness of metal (nm)

39

40

41

42