Extended long range plasmon waves in finite thickness metal film and layered dielectric materials

In this paper, we show that the propagation distance of the long range plasmon wave mode guided by a finite thickness gold metal film can be extended several orders of magnitude longer if we place intermediate dielectric layers on both sides of the metal film and choose the layer thickness properly. The propagation distance goes to infinite if the intermediate layer thickness approaches a critical thickness. © 2006 Optical Society of America OCIS codes: (240.6680) Surface plasmons; (130.2790) Guided waves References and links 1. H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings (Springer-Verlag, Berlin Heidelberg, 1988). 2. E. N. Economou, “Surface plasmons in thin films,” Phys. Rev. 182, 539-554 (1969). 3. J. J. Burke, G. I. Stegeman, and T. Tamir, “Surface-polariton-like waves guided by thin, lossy metal films,” Phys. Rev. B 33, 5186-5201 (1986). 4. D. Sarid, “Long-range surface-plasma waves on very thin metal films,” Phys. Rev. Lett. 47, 1927-1930 (1981). 5. P. Berini, “Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of symmetric structures,” Phys. Rev. B 61, 10484-10503 (2000). 6. P. Berini, “Plasmon-polariton waves guided by thin lossy metal films of finite width: Bound modes of symmetric structures,” Phys. Rev. B 63, 125417 (2001). 7. P. Berini, R. Charbonneau, N. Lahoud, and G. Mattiussi, “Characterization of long-range surface-plasmonpolariton waveguides,” J. Appl. Phys. 98, 043109 (2005). 8. S. J. Al-Bader, “Optical transmission on metallic wires fundamental modes,” IEEE J. Quantum Electron. 40, 325-329 (2004). 9. B. Lamprecht, J.R. Krenn, G. Schider, H. Ditlbacher, M. Salerno, N. Felidj, A. Leitner, and F.R. Aussenegg, "Surface plasmon propagation in microscale metal stripes," Appl. Phys. Lett. 79, 51-53 (2001). 10. A. Degiron and D. Smith, “Numerical simulations of long-range plasmons,” Opt. Express 14, 1611-1625 (2006), http://www.opticsinfobase.org/abstract.cfm?URI=oe-14-4-1611. 11. R. Charbonneau, P. Berini, E. Berolo, and E. Lisicka-Shrzek, "Experimental observation of plasmon polariton waves supported by a thin metal film of finite width ," Opt. Lett. 25, 844-846 (2000). 12. R. Charbonneau, N. Lahoud, G. Mattiussi, and P. Berini, "Demonstration of integrated optics elements based on long-ranging surface plasmon polaritons," Opt. Express 13, 977-984 (2005). 13. K. Leosson, T. Nikolajsen, A. Boltasseva, and S. I. Bozhevolnyi, "Long-range surface plasmon polariton nanowire waveguides for device applications," Opt. Express 14, 314-319 (2006). 14. L. Holland, Vacuum Deposition of Thin Films (Chapman and Hall, London, 1966). 15. G. I Stegeman and J. J. Burke, “Long-range surface plasmons in electrode structures,” Appl. Phys. Lett. 43, 221-223, (1983). 16. F. Y. Kou and T. Tamir, "Range extension of surface plasmons by dielectric layers," Opt. Lett. 12, 367(1987). 17. L. Wendler and R. Haupt, “Long-range surface plasmon-polaritons in asymmetric layer structures,” J. Appl. Phys. 59, pp. 3289-3291 (1986). 18. J. Xia, A. K. Jordan, and J. A. Kong, “Electromagnetic inverse-scattering theory for inhomogeneous #74676 $15.00 USD Received 7 September 2006; revised 3 November 2006; accepted 17 November 2006 (C) 2006 OSA 11 December 2006 / Vol. 14, No. 25 / OPTICS EXPRESS 12409 dielectrics: the local reflection model,” J. Opt. Soc. Am. A. 9, 740-748 (1992). 19. E. Anemogiannis, E. N Glytsis, and T. K Gaylord, “Determination of guided and leaky modes in lossless and lossy planar multilayer optical waveguides: reflection pole method and wavevector density method,” J. Lightwave Technol. 17, 929-941 (1999). 20. A. Papoulis, Circuits and Systems (Holt, Rinehart and Winston, Inc., New York, 1980). 21. P. B. Johnson and R. W. Christy, “Optical Constants of the Noble Metals,” Phys. Rev. B 6, 4370–4379 (1972). 22. G. I. Stegeman, J. J. Burke, and D. G. Hall, “Surface-polariton like waves guided by thin, lossy metal films,” Opt. Lett. 8, 383-386 (1983). 23. J. J. Burke, G. I. Stegeman, and T. Tamir, “Surface-polariton-like waves guided by thin, lossy metal films,” Phys. Rev. B 33, 5186-5201 (1986). 24. J. A. Dionne, L. A. Sweatlock, H. A. Atwater, and A. Polman “Planar metal plasmon waveguides: frequency-dependent dispersion, propagation, localization, and loss beyond the free electron model,” Phys. Rev. B 72, 075405 (2005).


Introduction
A surface plasmon wave is the propagation of the free electron density oscillation on the surface of a metal in contact with a dielectric medium.The propagation of the free electron density oscillation creates a localized electromagnetic wave.Surface plasmon waves can only be supported along the boundary of two optical materials with opposite signs of the real parts of the dielectric constants, such as a dielectric and a metal [1,2,3].Surface plasmon waves attenuate rapidly during the propagation due to the intrinsic electron oscillation damping loss in the metal.
A thin metal film embedded in a homogeneous dielectric medium with a symmetric structure supports two guided plasmon wave modes.One mode, which has the symmetric field profile with respect to the center of the metal layer, is called the symmetric mode.Another mode, which has the anti-symmetric field profile with respect to the center of the metal layer, is called the anti-symmetric mode.The symmetric mode has a relatively long propagation distance compared with the surface plasmon wave along the surface of the bulk metal in a dielectric medium.This symmetric surface plasmon wave mode is also called the long range surface plasmon (LRSP) mode.The anti-symmetric mode has a zero in the center of the metal layer and attenuates rapidly during propagation.The anti-symmetric mode is also called the short range surface plasmon mode.
Long range surface plasmon waves have been studied by many authors [4,5,6,7,8,9,10,11,12,13]. The fundamental reasons that thin metal films (or strips) in symmetric structures can support long range plasmon waves are: (1) most energy is located in the surrounding dielectric material for thin metal films; (2) the surface plasmon waves associated with two boundaries travel with the same velocity due to the symmetry of the guide structure.
Although long range plasmon wave modes can be supported by thin metal films in homogeneous dielectric media with symmetric structures, the long range modes still attenuate due to the electron oscillation damping loss in the metal.To increase the travel distance of the long range plasmon mode, a simple strategy is to reduce the metal layer thickness.But it is experimentally difficult to deposit homogeneous metal film of less than 15 nm thickness because metals, such as gold (Au), typically form nanoscale islands in the initial deposition process [14].
Surface plasmon waves guided by multiple parallel thin metal films placed closely together have been studied [15].Range extension of plasmon waves has been studied with multiple dielectric layers in asymmetric waveguide structures [16].The extended long range modes in [16] were leaky modes.Long range modes of thin metal films between two different semi-infinite uniform dielectric materials have been studied in [17].The extended long range modes in the asymmetric waveguides were also leaky modes.In this work, we studied the guided plasmon wave modes along a finite thickness thin metal (Au) film in inhomogeneous layered dielectric materials with symmetric waveguide structures.We found that if we place intermediate dielectric layers with different dielectric constant on both sides of the metal film, the propagation distance of the symmetric long range mode can be extended several orders of magnitude.Also we found that when the intermediate dielectric layer thickness approaches certain critical thickness, the propagation distance of one long range goes to infinite.

The reflection pole method
The plasmon waveguide structure we studied is shown in the Fig. 1.Fig. 1 is the cross section of a thin metal film in inhomogeneous layered dielectric materials.A semi-infinite dielectric medium (ε 0 ) is in the bottom.A thin metal film (Au) layer (ε 2 = ε m ) of thickness d 2 is in the middle.One intermediate dielectric layer (ε 1 ) of thickness d 1 is place between the metal layer and the dielectric medium in the bottom cladding.Another intermediate dielectric layer (ε 2 ) of thickness d 2 is place between the metal layer and the dielectric medium in the top cladding.Another semi-infinite dielectric medium (ε 4 ) is on the top of the top intermediate dielectric layer.The two intermediate dielectric layers have the same dielectric constant (ε 1 = ε 3 ) and the same thickness (d 1 = d 3 ).The top dielectric cladding is the same material as the bottom dielectric cladding (ε 0 = ε 4 ).Therefore, the plasmon waveguide structure is symmetrical with respect to the center of the metal layer.
We are only interested in the guided plasmon modes in this waveguide structure.We use the reflection pole method [18,19] to calculate the plasmon wave modes.Since the plasmon wave modes are TM mode, we write the magnetic fields in the different regions as In the above equation, are the x components of complex wave vectors in different regions, γ is the propagation constant (complex) in the z direction.i A and i B are coefficients corresponding to the propagating waves in x and -x directions in the i-th region.o k is the wave number in the free space.
We assume that an incident transverse magnetic (TM) electromagnetic wave is incident to the boundary of ε 0 and ε 1 from the bottom of the waveguide structure.At each boundary, the tangential components of the magnetic field yi H and the electrical field    1) exp( 1) exp( 1) exp( 1 The relationship between the incident, the reflection, and the transmission is given in The reflection coefficient is The reflection coefficient (R) is a function of the complex propagation constant . The guided plasmon modes are the modes which propagate in the z direction and exponentially decay away from the metal layer in the x and -x directions t versus the real part of the propagation constant index.From the Fig. 2, we see two rapid changes of the phase of 22 t .These two rapid phase changes correspond to two guided plasmon wave modes.Once we know the approximate values of the mode index from the plot, we can numerically solve the equation t 22 =0 with the initial guess of its roots.
The complex mode indices we found are n 1 =1.45229459-0.00003345585j,n 2 =1.525501873-0.013489423j.Both modes are guided modes.The first mode is the symmetric mode (s b ) which has a propagation distance (1/e of the intensity) of 3686.8 micron.The second mode is the anti-symmetric (a b ) mode, which has a propagation distance of 9.144 micron.The mode size (full width at 1/e of the maximum intensity) of the symmetric mode is 3.0429 micron.The mode size of the anti-symmetric mode is 0.54047 micron.
The surface plasmon mode properties of thin metal films embedded in infinite homogeneous dielectric materials have been studied by many authors [4,22,23,24] through solving the dispersion equation.To verify the modes obtained with the reflection pole method, we derived the dispersion equation by applying the boundary conditions of the Maxwell equations and numerically search the roots of the dispersion equation.We found two solutions: 1.452293989-0.00000334539jand 1.52550212-0.01348915j.These two solutions agree well with the results obtained by using the reflection pole method.This validates the RPM method.  of 1000 nm, we plot the phase (φ) of t 22 as the function of the real part mode index in Fig. 3 (a).We see the stair case phase π changes at certain values of the real part of the mode index.

Higher dielectric constant intermediate layer
According to Bode plot theory [20], each π phase change corresponds to the real part of the root of t 22 .The imaginary part of the root of t 22 corresponds to the width of the phase derivative curve (dφ/dn) of t 22 .Fig. 3 (b) shows the phase derivative versus the real part of the mode index.Once we know the approximate values of mode indices from the plot, we numerically solved the equation t 22 =0 with the initial guesses of its roots.For the intermediate dielectric layer thickness of 1000 nm, we found four complex solutions: n 1 =1.732403005-0.0001470494j,n 2 =1.80394551-0.00522046879j,n 3 =1.989189193-0.000206846j,n 4 =2.194229366-0.0341893452j.To verify these solutions are indeed the poles of the reflection coefficient, we input these complex numbers to the t 22 , and find |Re[t 22 ]| < 10 -21 and |Im[t 22 ]| < 10 -21 .These four solutions are all guided plasmon wave modes.The properties of these guided modes are listed in the Table 1.The longest range plasmon wave mode is the mode (a) which has the smallest imaginary part of the mode index.The propagation distance of this mode is 838.8 micron.The propagation distance is defined as the travel distance where the intensity decreases a factor of 1/e.The mode (c) is the second long range plasmon mode.The propagation distance of the second long range mode is 596.314micron.The modes (b) and (d) are the two short range modes.The propagation distances of the short range modes are one order of magnitude less than the long range modes.
We plotted the magnetic field mode profiles (H y ) of these four guided plasmon modes in Fig. 4. The modes (c) and (d) are the fundamental modes.The modes (a) and (b) are the higher order modes.We see that the long range modes (a) and (c) are symmetric modes.The short range modes (b) and (d) are anti-symmetric modes.The mode profiles are similar to the profiles of the dielectric waveguide modes, but have significant deviations from the dielectric waveguide modes near the metal layer due to the surface plasmon effect.Using the same technique, we searched the guided plasmon wave modes for varying intermediate dielectric layer thickness starting from zero to one micron.The gold metal film was fixed at 20 nm.Fig. 5 (a) shows the real part of the mode index of the guided plasmon wave modes for different intermediate layer thickness.We see that the structure supports two guided modes for the intermediate layer thickness less than 465 nm.The third guided mode (dashed green line) starts when the intermediate layer thickness reaches 465 nm.The structure starts to support the fourth guided mode (red dot line) when the intermediate layer thickness reaches 543 nm.As the intermediate dielectric layer thickness continuously increases, we expect that more guided modes will be supported.Fig. 5 (b) shows the imaginary part of the mode index of the four guided plasmon modes for varying intermediate layer thickness.
Figure 6 shows the propagation distances versus the intermediate dielectric layer thickness for these guided plasmon modes.We see that when the intermediate layer thickness is approaching to a critical thickness, the propagation distance of the longest range plasmon mode goes to infinite.For the intermediate layer thickness of 543 nm, the complex mode index of the longest range mode (red dot line) is 1.4500012-7.91467x10-7 .The propagation distance of this mode is 155.844mm.This propagation distance is approximately 42 times of the propagation distance of the long range plasmon wave mode guided by the gold film of same thickness without intermediate dielectric layers.It is well known that without the metal film, this waveguide structure supports dielectric waveguide modes.Dielectric waveguide modes are due to the total internal reflections from the boundaries of the high and the low index materials.Depending on the thickness of high index material layer and the contrast of the index of refraction, the dielectric waveguide can support only the fundamental symmetric mode with the Cosine function mode profile in the high index region or support additional higher order modes with the Sine and Cosine function mode profiles.With the presence of the metal film, the modes we have found here are due to the two effects: the total internal reflection effect and the surface plasmon effect.The mode profiles shown in Fig. 4 (b), (c), and (d

Lower dielectric constant intermediate layer
Then we chose the lower dielectric constant silicon dioxide (ε 1 = ε 3 =(1.45) 2 at 1.55 micron) as the intermediate layer dielectric material and the high dielectric constant silicon nitride (ε 0 = ε 4 =4.0 at 1.55 micron) as the top and bottom cladding material . The thickness of the gold metal layer is kept the same at d 2 = 20 nm.We searched the guided plasmon wave modes for different intermediate silicon dioxide dielectric layer thickness.Fig. 7 (a) shows the real part of the mode index for different intermediate dielectric layer thickness.The solid line curve is the real part of the mode index of the long range mode.The dashed line curve is the real part of the mode index of the short range mode.From Fig. 7 (a), we see that two guided modes are supported when the intermediate dielectric layer thickness is smaller than a critical thickness.Fig. 7 (b) shows the imaginary part of the mode index versus the low index intermediate layer thickness for the long range mode (solid line curve) and the short range mode (dashed line curve).We plotted the magnetic field mode profiles (H y ) of two guided plasmon modes for the intermediate layer thickness of 10.5 nm.The long range symmetric mode magnetic field profile is shown in Fig. 9 (a).Fig. 9 (b) is the close look of the mode profile near the metal layer.We see that the symmetric mode has a small notch inside the metal layer.This small notch is due to the electromagnetic field absorption inside the metal material.The short range magnetic field profile is shown in Fig. 10

Summary
We studied the guided plasmon wave modes of a thin gold metal film in inhomogeneous layered dielectric materials with the reflection pole method.The inhomogeneous layered plasmon waveguide structure consists of a finite thickness gold metal film and two identical dielectric intermediate layers on both sides of the metal film in a homogeneous dielectric material.The intermediate dielectric layers have a different dielectric constant from the surrounding uniform dielectric medium.First, we studied the guided plasmon wave modes by choosing the high dielectric constant material as the intermediate layer material and the low dielectric constant material as the top and bottom homogeneous cladding dielectric material.We found that multiple guided long range modes can be supported by this kind of structure and one of them can propagate sever orders of magnitude longer than the long range mode of the same thickness thin metal film in the homogeneous dielectric medium.As the intermediate layer thickness approaches to a critical thickness, the propagation distance of this long range mode goes to infinite.Without the metal film, this structure supports the dielectric waveguide modes, which are based on the principle of the total internal reflection.With the metal film, the modes we have found here are due to the two effects: the total internal reflection from the dielectric-dielectric boundaries and the surface plasmon oscillation.Then, we studied the guided modes by choosing the low dielectric constant material as the intermediate layer material and the high dielectric constant material as the top and bottom cladding dielectric medium.As we increase the low index intermediate layer thickness from zero, the propagation distance of the symmetric long range plasmon wave mode increases.The travel range of the long range mode goes to infinite when the intermediate layer approaches to a critical thickness.In this later case, the structure without metal film cannot support any guided mode since the intermediate dielectric layer material has the lower index of refraction than the surrounding homogeneous dielectric material.Both the long range mode and the short range mode are purely due to the surface plasmon effect.We believe that this ultra-long range surface plasmon waveguide structure with low index intermediate layers will have many applications since the symmetric mode profile is close to the Gaussian shape, which enables reasonably large couple efficiency from Gaussian beam laser to the plasmon wave mode.In summary, our studies have shown that the propagation distance of one guided plasmon wave mode of a thin metal film in inhomogeneous layered dielectric materials can be extended several orders of magnitude longer than the long range mode of the metal film of the same thickness in homogeneous dielectric medium.The propagation distance of this very long range mode goes to infinite as the intermediate layer thickness approaches to the critical thickness.
By applying these boundary conditions, we have

Fig. 1 .
Fig. 1.The plasmon waveguide structure of a thin metal film in layered inhomogeneous dielectric materials.

(
$15.00 USD Received 7 September 2006; revised 3 November 2006; accepted 17 November 2006 . The guided modes correspond to the resonances of the layered waveguide structure, which happen at the poles of the reflection coefficient R. The poles of the reflection coefficient are the solutions of of the real part index )(n of the propagation constant.The plasmon wave guided modes correspond to the rapid changes of the phase of 22 t .According to Bode plot theory[20], a peak of the derivative of the phase of 22 t corresponds to the real part index of a guided mode, and the full width at the half maximum of the phase derivative curve is the imaginary part of mode index.Using the reflection pole method, we searched the plasmon wave modes for a gold (Au) metal film of 20 nm thickness embedded in the silicon dioxide dielectric material ( 45) 2 ) at the 1.55 micron wavelength (d 1 = d 3 =0, d 2 =20 nm).The dielectric constant of gold film at 1.55 micron[21] is ε 2 =-114.925-11.0918j.Fig.2(a) is the plot of the phase of 22 t versus the real part of the propagation constant index.Fig.2 (b) is the derivative of the phase of22

Fig. 2 .
Fig. 2. (a) The phase of t 22 versus the real part of the mode index; (b) the derivative of the phase of t 22 versus the real part of the mode index.

First
we chose a thin gold (Au) film as the metal layer and silicon nitride (ε 1 = ε 3 =4.0)as the intermediate layer material and silicon dioxide (ε 0 = ε 4 = (1.45) 2 ) as the top and bottom cladding of the gold metal layer is d 2 = 20 nm.For the intermediate dielectric layer thickness )

Fig. 3 .
Fig. 3. (a) The phase of t 22 versus the real part of the mode index; (b) the derivative of the phase of t 22 versus the real part of the mode index.

Fig. 4 .
Fig. 4. The magnetic field Hy (in the unit of A/m) profiles of the guided plasmon modes.(a) is the longest range mode; (b) is one short range mode; (c) is the second long range mode; and (d) is the second short range mode.

Fig. 5 .
Fig. 5. (a) The real part of the mode index for varying intermediate dielectric layer thickness; (b) the imaginary part of the mode index for varying intermediate dielectric layer thickness.

Fig. 6 .
Fig. 6.Propagation distances of the guided modes versus intermediate layer thickness for the gold film thickness of 20 nm.

Fig. 7 .
Fig. 7. (a) The real part of the guided mode index versus the intermediate layer thickness; (b) the imaginary part of the guided mode index versus the intermediate layer thickness.

Figure 8 10 .Fig. 8 .
Figure8shows the propagation distance versus different intermediate layer thickness for the two guided modes.One mode is the long range mode (solid line curve).Another is the short range mode (dashed line curve).When the intermediate layer approaches to a critical thickness, the propagation distance of the long range mode goes to infinite.The critical intermediate layer thickness is approximately at nm 7 .10.At the intermediate layer thickness (a).Fig.10 (b) is the close look of short range magnetic profile near the metal region.The anti-symmetric mode magnetic field reaches zero in the center of the metal layer.Without the metal film, the waveguide structure cannot support any dielectric waveguide mode since the intermediate dielectric layer has the lower index of refraction than the cladding dielectric material.Both the long range mode and the short range mode here are purely due to the surface plasmon effect.The extended ultra-long travel range of the symmetric mode is due to the redistribution of the electromagnetic energy caused by the presence of the low index intermediate layers.More electromagnetic energy propagates in the dielectric media.

Fig. 9 .Fig. 10 .
Fig. 9. Magnetic field mode profile (in the unit of A/m) of the long range mode at the intermediate layer thickness d=10.5 nm.

Table 1 .
Guided plasmon wave modes for the intermediate layer thickness of 1000 nm