Tunable plasmon resonance of a touching gold cylinder arrays

Similar documents
Tuning of plasmon resonance through gold slit arrays with Y-shaped channels

Electromagnetic Absorption by Metamaterial Grating System

Surface Plasmon Resonance in Metallic Nanocylinder Array

SCATTERING OF ELECTROMAGNETIC WAVES ON METAL NANOPARTICLES. Tomáš Váry, Juraj Chlpík, Peter Markoš

Strong plasmon coupling between two gold nanospheres on a gold slab

Metamaterials & Plasmonics

Extinction properties of a sphere with negative permittivity and permeability

Surface Plasmon Polariton Assisted Metal-Dielectric Multilayers as Passband Filters for Ultraviolet Range

U-Shaped Nano-Apertures for Enhanced Optical Transmission and Resolution

Superconductivity Induced Transparency

Propagation of Surface Plasmon Polariton in the Single Interface of Gallium Lanthanum Sulfide and Silver

A Study on the Suitability of Indium Nitride for Terahertz Plasmonics

FDTD Simulations of Surface Plasmons Using the Effective Permittivity Applied to the Dispersive Media

Surface Plasmon Polaritons on Metallic Surfaces

Lecture 10 Light-Matter Interaction Part 4 Surface Polaritons 2. EECS Winter 2006 Nanophotonics and Nano-scale Fabrication P.C.

A Novel Design of Photonic Crystal Lens Based on Negative Refractive Index

MODAL ANALYSIS OF EXTRAORDINARY TRANSMISSION THROUGH AN ARRAY OF SUBWAVELENGTH SLITS

Backscattering enhancement of light by nanoparticles positioned in localized optical intensity peaks

Effect of Paired Apertures in a Periodic Hole Array on Higher Order Plasmon Modes

SURFACE PLASMONS AND THEIR APPLICATIONS IN ELECTRO-OPTICAL DEVICES

Density of modes maps for design of photonic crystal devices

Comparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses

Taking cascaded plasmonic field enhancement to the ultimate limit in silver nanoparticle dimers S. Toroghi* a, P. G. Kik a,b

Optimizing the performance of metal-semiconductor-metal photodetectors by embedding nanoparticles in the absorption layer

Spring 2009 EE 710: Nanoscience and Engineering

Photonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique

Single Emitter Detection with Fluorescence and Extinction Spectroscopy

SUPPLEMENTARY INFORMATION

Plasmon-suppressed vertically-standing nanometal structures

Nanophysics: Main trends

ABSTRACT 1. INTRODUCTION

Negative Index of Refraction in Optical Metamaterials

Chapter 2 Surface Plasmon Resonance

Chapter 5. Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice

Biosensing based on slow plasmon nanocavities

Quantum Information Processing with Electrons?

Tooth-shaped plasmonic waveguide filters with nanometeric. sizes

Multiple Fano Resonances Structure for Terahertz Applications

Demonstration of Near-Infrared Negative-Index Materials

Routing of Deep-Subwavelength Optical Beams and Images without Reflection and Diffraction Using Infinitely Anisotropic Metamaterials

Surface plasmon waveguides

Publication II Wiley Periodicals. Reprinted by permission of John Wiley & Sons.

Design of Uniform Fiber Bragg grating using Transfer matrix method

Title. Author(s)Nagasaki, Akira; Saitoh, Kunimasa; Koshiba, Masanori. CitationOptics Express, 19(4): Issue Date Doc URL.

Polarization anisotropic transmission through metallic Sierpinski-Carpet aperture array

Ankara, Turkey Published online: 20 Sep 2013.

Grating-coupled transmission-type surface plasmon resonance sensors based on dielectric and metallic gratings

Light trapping in thin-film solar cells: the role of guided modes

Research on the Wide-angle and Broadband 2D Photonic Crystal Polarization Splitter

Multiple extraordinary optical transmission peaks from evanescent coupling in perforated metal plates surrounded by dielectrics

Optimization of enhanced absorption in 3D-woodpile metallic photonic crystals

The Dielectric Function of a Metal ( Jellium )

Surface plasmon resonance based refractive index sensor for liquids

SENSITIVITY ENHANCEMENT OF A D-SHAPE SPR-POF LOW-COST SENSOR USING GRAPHENE

FDTD Analysis on Optical Confinement Structure with Electromagnetic Metamaterial

Lecture 10: Surface Plasmon Excitation. 5 nm

TUNABLE MULTI-CHANNEL FILTERING USING 1-D PHOTONIC QUANTUM WELL STRUCTURES

Directional emitter and beam splitter based on self-collimation effect

The observation of super-long range surface plasmon polaritons modes and its application as sensory devices

Negative refraction of photonic and polaritonic waves in periodic structures

Long-Wavelength Optical Properties of a Plasmonic Crystal

Light transmission through a single cylindrical hole in a metallic film

Nanomaterials and their Optical Applications

Light Localization in Left-Handed Media

GENERALIZED SURFACE PLASMON RESONANCE SENSORS USING METAMATERIALS AND NEGATIVE INDEX MATERIALS

Large omnidirectional band gaps in metallodielectric photonic crystals

l* = 109 nm Glycerol Clean Water Glycerol l = 108 nm Wavelength (nm)

Modeling of Kerr non-linear photonic components with mode expansion

Plasmonics: elementary excitation of a plasma (gas of free charges) nano-scale optics done with plasmons at metal interfaces

Extraordinary optical transmission through multi-layered systems of corrugated metallic thin films

Enhancing and suppressing radiation with some permeability-near-zero structures

Coherent thermal emission by excitation of magnetic polaritons between periodic strips and a metallic film

Polarization control and sensing with two-dimensional coupled photonic crystal microcavity arrays. Hatice Altug * and Jelena Vučković

Attenuation of Surface Plasmon Intensity by Transverse and Longitudinal. Slits

Left-handed materials: Transfer matrix method studies

Design of a Multi-Mode Interference Crossing Structure for Three Periodic Dielectric Waveguides

Monolayer Black Phosphorus

Super-reflection and Cloaking Based on Zero Index Metamaterial

ORE Open Research Exeter

APPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION AND REFLECTION BY A CONDUCTING SHEET

Photonic/Plasmonic Structures from Metallic Nanoparticles in a Glass Matrix

Electromagnetic scattering from multiple sub-wavelength apertures in metallic screens using the surface integral equation method

Collective effects in second-harmonic generation from plasmonic oligomers

Redirection of flexural waves in platonic crystal slabs

Simulated Study of Plasmonic Coupling in Noble Bimetallic Alloy Nanosphere Arrays

Transmission resonances on metallic gratings with very narrow slits

Electromagnetic Metamaterials

Author(s) Tamayama, Y; Nakanishi, T; Sugiyama. Citation PHYSICAL REVIEW B (2006), 73(19)

Microcavity enhanced directional transmission through a subwavelength plasmonic slit

Nanocomposite photonic crystal devices

transmission reflection absorption

Energy transport in metal nanoparticle plasmon waveguides

Transmission of Light Through Small Elliptical Apertures (Part 2)

FINITE DIFFERENCE TIME DOMAIN SIMULATION OF LIGHT TRAPPING IN A GaAs COMPLEX STRUCTURE

Enhanced Transmission by Periodic Hole. Arrays in Metal Films

Low Losses Left Handed Materials Using Metallic Magnetic Cylinders.

Supplementary Figure S1 Anticrossing and mode exchange between D1 (Wood's anomaly)

Theoretical studies on the transmission and reflection properties of metallic planar fractals

arxiv: v2 [cond-mat.other] 20 Nov 2008

The Broadband Fixed-Angle Source Technique (BFAST) LUMERICAL SOLUTIONS INC

Transcription:

J. At. Mol. Sci. doi: 10.4208/jams.091511.101811a Vol. 3, No. 3, pp. 252-261 August 2012 Tunable plasmon resonance of a touching gold cylinder arrays Geng-Hua Yan a, Yan-Ying Xiao a, Su-Xia Xie b, and Hong-Jian Li a,b, a College of Physics Science and Technology, Central South University, Changsha 410083, China b College of Materials Science and Engineering, Central South University, Changsha 410083, China Received 15 September 2011; Accepted (in revised version) 18 October 2011 Published Online 18 February 2012 Abstract. We investigate the plasmon resonance of a touching gold cylinder arrays surrounded with air or dielectric. We show that the plasmon resonance is tunable by modulating the cylinder radius, the cylinder layer and refractive index of the dielectric. It is found resonance peak blue-shifts and splits as the gold cylinder radius reduces, width of the resonance peak gets much smaller as layer increases, and resonant peak red-shifts noticeably as dielectric refractive index increases. Based on electric field distributions at the resonance wavelengths we reveal the mechanism of the transmission enhancement. These phenomena are helpful for the design of potential optics devices. PACS: 73.20.Mf, 78.20.-e, 78.67.-n Key words: touching gold cylinders, surface plasmon, electric field distributions 1 Introduction Metallic nano-structures have attracted a great attention as they can provide many intriguing properties for application in fields like near-field microscopy and spectroscopy, nanoscale photonic devices [1], biological applications, and biosensors [2]. Controllable and tunable surface plasmon resonance of metallic nanopaticales [3], nanoshells [4], nanorods [5], and nanocylinders [6] are investigated recently. It is found that the electromagnetic field can be enhanced in the vicinity of the nanostructures [7,8]. Parameters such as size [9], shape [10], and the surrounding dielectric [11,12] of the structure are investigated in order to clarify the optical properties. It is reported there is an extraordinary optical transmission for a metallic slab perforated with hole arrays [13], even for a continuous metallic Corresponding author. Email address: lihj398@yahoo.com.cn (H. J. Li) http://www.global-sci.org/jams 252 c 2012 Global-Science Press

G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 253 structure with periodically nanostructured on the surface [14]. If a series of nanocylinders are arranged parallelly touching each other periodically, forming a continuous structure, high enhancement transmission is expected. In this paper, we investigate plasmon resonance of a touching gold cylinder arrays surrounded with air or dielectric. We show that the plasmon resonance is tunable by modulating the cylinder radius, the cylinder layer and refractive index of the dielectric. It is found that resonance peak blue-shifts and splits as the gold cylinder radius reduces, width of the resonance peak gets much smaller as layer increases, and resonant peak red shifts noticeably as dielectric refractive index increases. By analyzing electric field distributions at the resonance wavelengths, we reveal the mechanism of the transmission enhancement. These properties are helpful for the design of potential optics devices. 2 Model and theory We apply a periodic array of gold cylinders as shown in Fig. 1, in order to demonstrate the optical property. Transmission spectra and electric fields are simulated by using the finite-difference time-domain (FDTD) method [15, 16]. In the entire work, we investigate a lattice of the periodic cylinder structure with 600 nm length in x-direction and 800 nm length in y-direction; cylinder radius is noted as r, and the thickness of the dielectric is noted as h, and the spatial and temporal steps are set at x= y=lnm and t= x/2c (c is the velocity of light in vacuum), and we send a Gaussian single pulse of light with a wide frequency profile. Periodic boundary conditions are imposed on the left and right surfaces, and the perfectly matched layers (PML) are used on the top and bottom of the lattice [17]. Since there are no gaps between cylinders in the structure, all effects related to mechanism of light transmission directly through the slits are absent, and the resonant tunneling through a metal film is the only mechanism responsible for the enhanced transmission. In all calculations below, light illuminates on top of the structure at a normal incidence polarizing along x-direction. When a plane wave is incident normally on the interface between a metal and a dielectric, in order to excite surface plasmon waves, the phase-matching condition as follows should be obeyed by the incident wave vector [13, 18] ksp = k0 sinθ+i Gx, (1) where k 0 is the wave vector of incident light and θ is the incident angle. Gx is the Bragg vector and Gx =2π/p, i is an integer, which is the mode indices. ksp is the surface plasmon wave vector whose length is εd ε ksp = k0 m (2) ε d +ε m ε d and ε m are permittivities of the dielectric and the metal, respectively.

254 G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 Figure 1: A schematic of the x y cross section of the investigated lattice of touching gold cylinder array. Light illuminates on top of the structure in a normal incidence and polarizes along x-direction. The metallic structure is made of gold, and the dielectric constant of gold is described by the Drude model [19] ε m (ω)=1 ω2 p ω 2 +iγω, (3) where ω p is the plasma frequency and γ is the collision frequency, related to energy loss. The parameters used in the Drude model for gold are ω p = 1.374 10 16 s 1 and γ=4.08 10 13 s 1, taken from Ref. [20]. 3 Results and discussion Fig. 2 shows the zeroth order transmission, reflection, and absorption spectra as a function of wavelength λ for different cylinder radii of the bare touching gold cylinder structure. The absorption spectra are obtained via A =1 T R, where A, T, and R are the absorbance, transmittance, and reflectance, respectively. It is worth noting that peak value, peak shape and center peak wavelength of the resonance spectra varied noticeably as the radius alters. From Fig. 2 it can be seen that the center wavelengths of the transmission resonance peaks for the structure with cylinder radius r=150 nm, r=75 nm, and r=50 nm are λ=642 nm, λ=450 nm, and λ=432 nm, respectively, here the length of the investigated structure is set as 600 nm invariably. With radius decreasing, transmission is intensified, and the transmission peak splits into two peaks (when the cylinder radius decreases to r=60 nm in the lattice, the transmission peak starts splitting) and the split becomes greater as the cylinder number increases. It is reported that transmission can be enhanced in corrugated metallic films without holes by Bonod et al. [21], which is recently demonstrated experimentally [14]. Since

G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 255 Figure 2: Transmission, reflection, and absorption spectra of the bare touching gold cylinder array for different cylinder radii: (a) r= 150 nm, (b) r= 75 nm, and (c) r= 50 nm. the cylinders touch each other in a line, the system can be treated as a continuous film with periodic corrugations on both surfaces. The corrugations play no significant role in the light transmission through the film directly, but only serve as a periodical perturbation to ensure the coupling between the incident wave and the surface plasmons at the two interfaces of the cylinder and air. When compared with grating hole array, high transmission through film requires smaller thickness. This is natural since propagation in gold is attenuated stronger than in the holes. For a structure with larger cylinder radius, the interaction between the surface plasmon polariton modes on the up and bottom interfaces is weaker, most electromagnetic energy is absorbed and reflected at the transmission resonant wavelength, while for a structure with smaller cylinder radius, the interaction between the surface plasmon polariton modes on the opposite interfaces is stronger, so the transmission peak is magnified and splitting, and the surface plasmon resonance would happen in a high energy state, so the transmission peak shifts to shorter

256 G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 wavelength range as the cylinder radius gets smaller, it is obtained from Fig. 2 that continuous the transmission wavelength shifts λ=210 nm to the shorter wavelength region as the cylinder radius is reduced from 150 nm to 50 nm. In addition, there is a low transmission resonance peak at a shorter wavelength region, which exhibits a similar trend of the high transmission peak. We treat it as the high order coupling transmission resonance. It is found that a peak in a transmission curve always corresponds to a minimum in the reflection and absorption curves at the same frequency as shown in Fig. 2, and the reflection and the absorption at the high transmission wavelength are reduced as the transmission strengthened, which means electromagnetic radiation at this wavelength range is mostly transmitted but reflected and absorbed little. As the cylinder radius becomes smaller, the dips of the reflection and absorption get even lower. Furthermore, we can see a clear blue-shift in the resonance dip of the reflectance and absorbance spectra as the cylinder radius is reduced. The reflectance and absorbance spectra clearly indicate the resonant features. In Fig. 3 we show the field distributions associated with the main transmission resonances of the periodic structure with different cylinder radii reported in Fig. 2. The x/y components of the electric field have completely different distributions at the main resonance. For the E x component (polarized perpendicular to the wave propagating), the electric field is subject to an enhancement near the contact point of cylinders. Wave nodes of E x confined near the cylinder junctions with the same charge sign on the top and bottom spacing. When a normal incidence matches the frequency of a standing-wave surface plasmon mode, a surface plasmon mode on the top surface is excited and strong field is built up inside the groove schemed by two adjacent cylinders, which in turn resonantly excites the standing-wave surface plasmon mode on the bottom surface of the cylinder structure, then a strong electric field is built up in the grooves of the bottom surface, which finally emits the radiation downwards into air. As a result, resonant transmission happens at the frequencies of standing-wave surface plasmon polariton modes. This Figure 3: The distribution of the electric field intensity E x and E y at transmission wavelengths λ=642 nm, λ=450 nm and λ=432 nm, for the bare touching gold cylinder structure with cylinder radii r=150 nm, r=75 nm, and r= 50 nm, respectively, corresponding to Fig. 2(a), 2(b) and 2(c).

G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 257 process is effectively photon tunneling through touching gold cylinders via resonantly exiting surface plasmon modes. Resonant excitation of a surface plasmon mode in the groove between cylinders results in strong field enhancement and reflection and absorption minimum at the frequency of the surface plasmon mode. For the E y component (polarized parallel to the wave propagating vector), the electric field exhibits an opposite symmetry property of the field modes in x and y directions with +/- charges as a form of quadruple on each cylinder surface, and each cylinder remains neutral. Figure 4: Comparison of transmission spectra of touching and un-touching structure composed of gold cylinders with radius r= 75 nm. Cylinder separation is marked as s=0 and s=50 nm. Touching gold cylinders periodically arranged form a special kind of corrugated gold film. It has the enhanced transmission property, which is the same as a corrugated metallic film. It is found the transmission, reflection and absorption of this touching structure depend on the cylinder radius sensitively. Additionally, it is convenient and practicable to tune the transmission by adjusting cylinder radius in the realization of experiment. In our study, we focus our attention on the optical transmission of the structure with gold touching nanocylinder array. It is found that the structure with touching cylinders has a good filter property, in other words, there is only one enhanced resonance mode in the transmission spectra. However, for the un-touching one, there are many resonance modes in the transmission spectrum as shown in Fig. 4 and Fig. 5. In Fig. 4, we simulate the comparison of transmission spectra of touching and un-touching structures by varying the separation between cylinders, while keeping cylinder radius unchanged. In Fig. 5, we depict the difference of transmission characteristics for the touching and un-touching cylinder structure by modifying the cylinder radius with central position of each cylinder being immobilized. It is known that the surface plasmon is very sensitive to the refractive index in the

258 G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 Figure 5: Transmission spectra for the structure changes from a touching one to a un-touching one by decreasing cylinder radius from r=75 nm to r=55 nm in step of 10 nm with central position of each cylinder being immobilized. vicinity of the metal surface of the periodic hole array. Here, we investigate the transmission property of the periodic touching gold cylinders with radius r = 100 nm embedded in a dielectric with thickness of h=600 nm in y direction. Transmission spectra for the structure with single layer touching gold cylinders embedded in the center of the dielectric with refractive index increasing from n=1 to n=2 in steps of 0.5 are shown in Fig. 6(a). We can see clearly the transmission peak red shifts noticeably with the dielectric Figure 6: Transmission through the periodic touching gold cylinder array embedded in a material with refractive index varying from n=1 to 2 in steps of 0.5: (a) for the single layer structure, and (b) for the double-layer structure.

G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 259 Figure 7: Transmission through 1 to 3 layers of periodic touching gold cylinder arrays: (a) for the structure with bare gold cylinders and (b) for the structure with cylinders embedded in a dielectric with refractive index n=1.5. Figure 8: The distribution of the electric field intensity E x and E y for the bare structure with one, two, and three layers of touching gold cylinders at transmission wavelengths λ=518 nm, λ=518 nm and λ=510 nm, respectively, corresponding to spectra in Fig. 7(a). refractive index increasing. In other words, the surface plasmon resonant modes depend on the refractive index of the dielectric covered on the cylinder surface sensitively, so that the surface plasmon resonance can be used to design a potential device such as a sensor. Maybe the detected peak is quite broad as shown in Fig. 6(a), to optimize the sensitivity of the potential sensor, it is necessary to narrow the peak width. On the purpose, we investigate the transmission spectra for the structure composed of two layers of touching gold cylinders close to each other embedded in a dielectric with refractive index from n=1 to n=2 as shown in Fig. 6(b). With the increase of the refractive index, the resonance peak redshifts gradually with a decrease of the peak value. It is found that full

260 G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 width at half maximum of the resonance peak is narrowed significantly for the structure with two layers compared to that with a single layer, which meets our expectation well. So structure with double-layer cylinders embedded in a dielectric has a better sensitivity. We have calculated the transmission spectra of the bare structure with different layers of touching gold cylinders as shown in Fig. 7(a) and the structure with different layers of touching gold cylinders embedded in a dielectric with refractive index n=1.5 as shown in Fig. 7(b). It can be seen clearly from both figures that transmission resonance peak shifts little for both structures as cylinder layer increases, however, the resonance occurs at a much longer wavelength region for the structure embedded in a dielectric than that for the bare structure. In addition, with the layer increasing, the full width at half maximum of the resonance peak decreases for both the bare and embedded structures. Particularly, for the bare structure, the resonant maximum of the structure with double-layer cylinders is almost the same as that with a single layer as shown in Fig. 7(a). Such a structure is able to provide a sharp filtering. Peak value of the embedded structure is reduced as the layer of cylinders increases as shown in Fig. 7(b). In Fig. 8 we show the cross sections of the field distribution as a function of the cylinder layer at the transmission resonance wavelength λ=518 nm, λ=518 nm and λ=510 nm for a normal illumination corresponding to the spectra in Fig. 7(a). The field distribution of E x is homogeneous for structures with different layers, and the field amplitude changes insignificantly, which according well to the resonance transmission property as shown in Fig. 7(a). E y distribution has both plus and minus charges on each cylinder, which is associated with strongly confined charges of opposite signs as a form of quadruple on each cylinder surface. 4 Conclusions In conclusion, light tunneling via surface plasmon modes through a periodic touching gold cylinder array has been investigated by using FDTD method. A significant enhancement of the light transmission is provided by the cylinder structure, and the enhancement can be modulated by the cylinder radius, cylinder layer, and refractive index of the dielectric around cylinders. Resonance peak blue shifts and splits as the gold cylinder radius reduces, width of the resonance peak gets much smaller as layer increases, and resonant peak red shifts noticeably as dielectric refractive index increases. The electric field distributions are calculated to explicate the transmission mechanism. These results are helpful to design a good filter device. Acknowledgments. This work was funded by the Research Fund for the Doctoral Program of Higher Education of China under Grant No. 20100162110068 and the National Natural Science Foundation of China under Grant No. 11164007.

G. H. Yan, Y. Y. Xiao, S. X. Xie, and H. J. Li / J. At. Mol. Sci. 3 (2012) 252-261 261 References [1] P. N. Prasad, Nanophotonics (Wiley, Hoboken, New York, 2004). [2] D. A. Schultz, Curr. Opin. Biotechnol. 14 (2003)13. [3] J. Aizpurua, G. W. Bryant, L. J. Richter, et al., Phys. Rev. B 71 (2005) 235420. [4] K. Höflich, U. Gösele, and S. Christiansen, J. Chem. Phys. 131 (2009) 164704. [5] K. G. Thomas, S. Barazzouk, B. I. Ipe, et al., J. Phys. Chem. B 108 (2004) 13066. [6] S. Seguia, J. L. Gervasoni, and N. R. Arista, Surf. Sci. 601 (2007) 4169. [7] C. Bohren and D. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, New York, 1983). [8] U. Kreibig and M. Vollmer, Optical Properties of Metal Clusters (Springer, Berlin, 1995). [9] O. Stranik, R. Nooney, C. McDonagh, et al., Plasmonics 2 (2007)15. [10] M. H. Rashid, R. R. Bhattacharjee, and T. K. Mandal, J. Phys. Chem. C 111 (2007) 9684. [11] H. Takele, H. Greve, C. Pochstein, et al., Nanotechnology 17 (2006) 3499. [12] D. P. Lyvers, J. M. Moon, A. V. Kildishev, et al., ACS Nano. 2 (2008) 2569. [13] T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, et al., Nature 391 (1998) 667. [14] B. F. Bai, L. F. Li, and L. J. Zeng, Opt. Lett. 30 (2005) 2360. [15] K. S. Yee, IEEE Trans. Ant. Prop. 14 (1996) 302. [16] A. Taflove, Computational Electrodynamics: the Finite Difference Time Domain Method, 2nd ed. (Artech House, Norwood, 2002). [17] J. P. Berenger, J. Comput. Phys. 114 (1994) 185. [18] H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings (Springer- Verlag, Berlin, 1988). [19] P. Drude, Ann. Phys. 306 (2009) 566. [20] M. A. Ordal, L. L. Long, R. J. Bell, et al., Appl. Opt. 22(1983)1099. [21] N. Bonod, S. Enoch, L. Li, et al., Opt. Express 11 (2003) 482.