## Abstract

We study the surface plasmon (SP) resonance energy of isolated spherical Ag nanoparticles dispersed on a silicon nitride substrate in the diameter range 3.5–26 nm with monochromated electron energy-loss spectroscopy. A significant blueshift of the SP resonance energy of 0.5 eV is measured when the particle size decreases from 26 down to 3.5 nm. We interpret the observed blueshift using three models for a metallic sphere embedded in homogeneous background material: a classical Drude model with a homogeneous electron density profile in the metal, a semiclassical model corrected for an inhomogeneous electron density associated with quantum confinement, and a semiclassical nonlocal hydrodynamic description of the electron density. We find that the latter two models provide a qualitative explanation for the observed blueshift, but the theoretical predictions show smaller blueshifts than observed experimentally.

## Reviewed Publication:

de Abajo Javier Garcia

## 1 Introduction

Surface plasmons are collective excitations of the electron gas in metallic structures at the metal/dielectric interface [1]. The ability to concentrate light with SPs [2] and to enhance light-matter interaction on a subwavelength scale enables few and even single-molecule spectroscopy when the size of the metallic structures is decreased to a few nanometers [3]. These collective excitations are usually well-described by the classical Drude model for nanoparticles with dimensions of tens of nanometer and larger [1]. In the quasistatic limit, i.e., when the wavelength of the exciting electromagnetic wave considerably exceeds the dimensions of the structure, the local-response Drude theory predicts that the resonance energy of localized SPs is independent of the size of the nanostructure [4], and that the field enhancement created in the gap between two metallic nanostructures diverges for vanishing gap size [5]. These predictions are however in conflict both with earlier [6–9] and with more recent experimental results, which have shown a size dependency of the localized SP resonance in noble metal nanoparticles in the size range of 1–10 nm [10] and pronounced deviations for dimer geometries [11, 12].

This dependence of the SP resonance on the size of noble metal nanostructures is believed to be a signature of quantum properties of the free-electron gas. With decreasing sizes of the nanoparticles, the quantum wave nature of the electrons is theoretically expected to manifest itself in the optical response due to the effects of quantum confinement [13–17], quantum tunneling [17–20], as well as nonlocal response [21–27]. Nonlocal effects are a direct consequence of the inhomogeneity of the electron gas, which arises due to the quantum wave nature and the many-body properties of the electron gas.

The recent developments in analytical scanning transmission electron microscopes (STEM) equipped with a monochromator and electron energy-loss spectroscopy (EELS) [28] give the possibility of accessing the near-field energy distribution of the plasmon resonance of individual nanoparticles on a subnanometer scale with an energy resolution better than 0.2 eV. This method has been used for the imaging of surface plasmons in many different metallic nanostructures [10, 29–32]. With STEM EELS it is possible to correlate the structural and chemical information on the nanometer scale, such as the shape and the presence of organic ligands, with the spectral information of the SP resonance of single isolated nanoparticles. STEM EELS is thus perfectly suited to probe and access plasmonic nanostructures and SP resonances at length scales where quantum mechanics is anticipated to become important.

In this paper we report the experimental study of the SP resonance of chemically grown single Ag nanoparticles dispersed on 10 nm thick Si_{3}N_{4} membranes with STEM EELS. Our measurements present a significant blueshift of the SP resonance energy from 3.2 to 3.7 eV for particle diameters ranging from 26 down to 3.5 nm. Our results also confirm very recent experiments made with Ag nanoparticles on different substrates using different STEM operating conditions [10], thereby strengthening the interpretation that the blueshift is predominantly associated with the tight confinement of the plasma and the intrinsic quantum properties of the electron gas itself rather than having an extrinsic cause.

We compare our experimental data to three different models: a purely classical local-response Drude model which assumes a constant electron density profile in the metal nanoparticle, a semiclassical local-response Drude model where the electron density is determined from the quantum mechanical problem of electrons moving in an infinite spherical potential well [16], and finally, a semiclassical model based on the hydrodynamic description of the motion of the electron gas which takes into account nonlocal response through the internal quantum kinetics of the electron gas in the Thomas-Fermi (TF) approximation [33, 34]. We find good qualitative agreement between our experimental data and the two semiclassical models, thus supporting the anticipated nonlocal nature of SPs of Ag nanoparticles in the 1–10 nm size regime. The experimentally observed blueshift is however significantly larger than the predictions by the two semiclassical models.

## 2 Materials and methods

The nanoparticles are grown chemically following the method described in Ref. [35] and subsequently stabilized in an aqueous solution with borohydride ions. The mean size of the nanoparticles is 12 nm with a very broad size distribution ranging from 3 to 30 nm. The nanoparticle solution is dispersed on a 10 nm thick commercially available Si_{3}N_{4} membrane (TEMwindows.com), which has a refractive index of approximately *n*≈2.1 [36]. To characterize our nanoparticles we have used an aberration-corrected STEM FEI Titan (www.FEI.com) operated at 120 kV with a probe diameter of approximately 0.5 nm, and convergence and collection angles of 15 mrads and 17 mrads, respectively. The Titan is equipped with a monochromator allowing us to perform EELS with an energy resolution of 0.15±0.05 eV. We systematically performed EELS measurements at the surface and in the middle of each nanoparticle. The EELS spectra were taken with an exposure time of 90 ms to avoid beam damage as much as possible. To improve the signal-to-noise ratio we accumulated 10–15 spectra for each measurement point. We observed no evidence of damage after each measurement.

The experimental data were analyzed with the aid of commercially available software (Digital Micrograph) and three different methods were used to reconstruct and remove the zero-loss peak (ZLP): the first method is the reflected tail (RT) method, where the negative-energy half part of the ZLP is reflected about the zero-energy axis to approximate the ZLP at positive energies, while the second method is based on fitting the ZLP to the sum of a Gaussian and a Lorentzian functions. The third method is to pre-record the ZLP prior to each set of EELS measurements. All three methods yielded consistent results.

The energies of the SP resonance peaks were determined by using a nonlinear least-squares fit of our data to Gaussian functions. The error in the resonance energy is given by the 95 % confidence interval for the estimate of the position of the center of the Gaussian peak. Nanoparticle diameters were determined by calculating the area of the imaged particle and assigning to the area an effective diameter by assuming a perfect circular shape. The error bars in the size therefore correspond to the deviation from the assumption of a circular shape, which is estimated as the difference between the largest and smallest diameter of the particle.

## 3 Theory

In the following theoretical analysis our hypothesis is that the blueshift of the SP resonance energy is related to the properties of the electron density profile in the metal nanoparticle. Therefore, we use three different approaches to model the electron density of the Ag nanoparticle. In all three approaches, we calculate the optical response and thereby also the resonance energies of the nanoparticle through the quasistatic polarizability *α* of a sphere embedded in a homogeneous background dielectric with permittivity *ε*_{B}. With this approach, we make two implicit assumptions: the first is that we can neglect retardation effects and the second is that we can neglect the symmetry-breaking effect of the substrate. We have validated the quasistatic approach by comparing to fully retarded calculations [37], which shows excellent agreement in the particle size range we consider. The effect of the substrate will be taken into account indirectly by determining an effective homogeneous background permittivity *ε*_{B} using the average resonance frequency of the largest particles (2*R*>20 nm) as the classical limit.

The first, and simplest, approach is to assume a constant free-electron density *n*_{0} in the metal particle, which drops abruptly to zero outside the particle. This assumption is the starting point of the classical local-response Drude model for the response of the Ag nanoparticle, where the polarizability is given by the Clausius-Mossotti relation, which is well-known to be size independent for subwavelength particles. The classical local-response polarizability *α*_{L} is [1]

where *R* is the radius of the particle and

*ε*(

_{∞}*ω*), not included in the plasma response of the free-electron gas itself.

The second approach is to correct the standard approximation in local-response theory of a homogeneous electron density profile by using insight from the quantum wave nature of electrons to model the electron density profile and take into account the quantum confinement of the electrons. For nanometer-sized spheres, the classical polarizability given by the Clausius-Mossotti relation must be altered to take into account an inhomogeneous electron density. In Ref. [16], it is shown that in general the local-response polarizability for a sphere embedded in a homogeneous material is given as

now with a spatially varying Drude permittivity [16, 17]

Here, *n*(*r*) is the electron density in the metal nanoparticle. Clearly, if *n*(*r*)=*n*_{0} we arrive at the classical Clausius-Mossotti relation Eq. (1) as expected. To determine the density profile in this local-response model, we follow the approach of Ref. [16] and assume that the free electrons move in an infinite spherical potential well. The approach just outlined of a local-response theory with an inhomogeneous electron density is very similar to the theoretical model used in Ref. [10] for explaining their experimental results. It should be noted that any effects due to electron spill-out and quantum tunneling are neglected in all of the approaches that we consider.

The third and final approach is to compare our experimental data with a linearized nonlocal hydrodynamic model in which the electron density is allowed to deviate slightly from the constant electron density used in classical local-response theories [22, 38–40]. The dynamics of the electron gas is governed by the semiclassical hydrodynamic equation of motion [25, 26, 34], which results in an inhomogeneous electron density profile. The nonlocal hydrodynamic polarizability *α*_{NL}(*ω*) is exactly given as

and these results constitute our nonlocal-response generalization of the Clausius-Mossotti relation of classical optics. Here,

is the wave vector of the additional longitudinal wave allowed to be excited in the hydrodynamic nonlocal theory [25, 34], and*j*

_{1}is the spherical Bessel function of first order. Finally, within TF theory where

*ν*

_{F}is the Fermi velocity [34]. We emphasize that for

*β*→0, the local-response Drude result is retrieved, since

*δ*

_{NL}→0 and Eq. (4) simplifies to the classical Clausius-Mossotti relation Eq. (1).

The SP resonance energy follows theoretically from the Fröhlich condition, i.e., we must consider the poles of Eq. (4). For sufficiently small blueshifts and neglecting damping, the resonance frequency can be approximated by

where the first term is the common size-independent local-response Drude result for the SP resonance that also follows from Eq. (1), and the second term gives the size-dependent blueshift due to nonlocal corrections. At this stage, we note that a 1/(2*R*) dependence was experimentally observed in Refs. [6, 7] using optical spectroscopy. However, Eq. (6) reveals, besides a 1/(2*R*) dependence, that there is a delicate interplay in the blueshift between the material parameters of the metal, through *ε*_{∞}(*ω*) and *β*, and the background medium *ε*_{B}. Furthermore, Eq. (6) shows that the blueshift can be enhanced with a large-permittivity background medium.

## 4 Results

Figures 1(A–C) display STEM images of Ag nanoparticles with diameters of 15.5, 10.0, and 5.5 nm, respectively. The images show that no chemical residue is left from the synthesis and that the particles are faceted. We find that approximately 70% of the studied nanoparticles have a relative size error (i.e., the ratio of the size error bar to the particle diameter) below 20% (determined from the 2D STEM images), verifying that the shape of the nanoparticles is to a first approximation overall spherical (see Supplementary Figure 1). On a subset of the particles, thickness measurements using image recordings at different tilt angles were performed, revealing information about the shape of the nanoparticle in the third dimension. Such 3D investigations confirmed that the shape is overall spherical, but however could not be completed for all particles due to stability issues: the positions of tiny nanoparticles fluctuate under too long exposure of the electron beam, thus preventing accurate determination of the shape of the nanoparticle in the third dimension perpendicular to the substrate.

Figures 1(D–F) display raw normalized EELS data, acquired on Ag nanoparticles with diameters 14.1, 9.8, and 6.6 nm, respectively. The peaks correspond to the excitation of the SP. When the diameter of the nanoparticle decreases, the SP resonance clearly shifts progressively to higher energies. Figures 1(D–F) also display that the amplitude and linewidth of the SP resonances can vary from particle to particle (with the same size) and at times show narrowing instead of the expected broadening of the resonance for decreasing nanoparticle sizes [6, 13, 14]. This is for example seen in the linewidths in Figures 1(D–F) which seem to decrease with size. However, as will be explained in more detail in the next paragraph, we did not find a systematic trend of the linewidths in our EELS measurements probably due to the shape variations in our ensemble of nanoparticles.

Figure 2 displays the resonance energy of the SP as a function of the diameter of the nanoparticles. A significant blueshift of the SP resonance of 0.5 eV is observed when the nanoparticle diameter decreases from 26 to 3.5 nm. This trend is in good agreement with the results shown in Ref. [10], despite the difference in the substrate and the STEM operating conditions, a strong indication that the blueshift of Ag nanoparticles is robust to extrinsic variations. Another prominent feature in Figure 2 is the scatter of resonance energies at a fixed particle diameter. We mainly attribute the spread in resonance energies at a given particle size to shape variations in our ensemble of nanoparticles (see Supplementary Material). Slight deviations from perfect circular shape in the STEM images will result in a delicate dependency on the location of the electron probe and give rise to splitting of SP resonance energies due to degeneracy lifting. In this regard, we also note that even a perfectly circular particle on a 2D STEM image may still possess some weak prolate or oblate deformation in the third dimension, resulting in a departure from spherical shape. Calculations using the local response model show that a 20% deformation of a sphere into an oblate or prolate spheroid results in a 0.4 eV spread in resonance energy (see Supplementary Figure 2), which is approximately the spread in resonance energy we observe for particles larger than 10 nm. Furthermore, shape deviations may also impact the linewidth of the SP resonance, since the electron probe can excite the closely-spaced non-degenerate resonance energies simultaneously, which may appear as a single broadened peak. This broadening mechanism could explain the apparent linewidth narrowing for decreasing particle size seen in Figures 1(D–F). However, we cannot rule out that other effects beyond shape deviations contribute to the spread of resonance energies and impact the SP resonance linewidth. These could for example be the facets or the particle-to-substrate interface [42].

Along with the EELS measurements in Figure 2, we show Eq. (1) for the local-response Drude model (red line) and the semiclassical local-response model Eq. (3) (blue line). Furthermore, the nonlocal relation of Eq. (3) (green solid line) and the approximate nonlocal relation of Eq. (6) (green dashed line) are also depicted, and we see that Eq. (6) is accurate for particle sizes

nm.Due to the narrow energy range in consideration (~3.0–3.9 eV), we approximate *ε*_{∞}(*ω*) as a second-order Taylor polynomial based on the frequency-dependent values given for Ag in Ref. [41]. We find *ε*_{∞}(*ω*)=(59.8+*i*55.1)(*ω/ω*_{P})* ^{2}*-(40.3+

*i*42.4)(

*ω/ω*

_{P})+(10.5+i8.6). Since the refractive index of the Si

_{3}N

_{4}substrate varies hardly (

*n*≈2.1) in the narrow energy range we consider [36], we assume that the background permittivity

*ε*

_{B}is constant and determine it by approximating the average resonance energy of the largest particles (2

*R*>20 nm) as the classical limit, i.e., the first term of Eq. (6).

It is known that local Drude theory produces size-independent resonance frequencies of subwavelength particles, but this theory is clearly inadequate to describe the measurements of Figure 2. The nonlocal quasistatic hydrodynamic model predicts a blueshift in agreement with the experimental EELS measurements. Interestingly, the measured blueshift is even larger than predicted. We also see that the local-response model with an inhomogeneous electron density profile shows a similar trend as the nonlocal hydrodynamic model, indicating that these two different models describe very similar physical effects. The oscillations in the resonance energy in the inhomogeneous local-response model seen for small particle diameter are due to small variations in the density profile with decreasing size (i.e., discrete changes in the number of electrons), as also stated in Ref. [10].

The inhomogeneous local-response model and the nonlocal hydrodynamic model, when applied to a sphere in a homogeneous background medium, agree qualitatively with the EELS measurements. However, they do not provide the full picture. One of the probable issues arising is that the substrate is taken into account indirectly through a homogeneous background medium, a state-of-the-art procedure [10] which however may not be adequate to describe the effects of the presence of a dielectric substrate. It has been shown that the dielectric substrate modifies the absorption spectrum of an isolated sphere [43] and also the waveguiding properties of nanowires [31, 44, 45]. In an attempt to include the symmetry breaking effect of the substrate in our theoretical analysis, we apply a simple image charge model. The main effect of the substrate in this picture stems from the interaction of the dipole mode of the nanoparticle with the induced dipole mode in the substrate [46–48]. However, we find that such a dipole-dipole model for the substrate is inadequate to describe the large blueshift observed experimentally (see Supplementary Material). Indeed, it has been shown that the induced image charges in the substrate can make the contributions of higher order multipoles in the nanoparticle important [49], and it has also been observed theoretically that higher order multipoles produce larger blueshifts in the nonlocal hydrodynamic model (Figure 2 in Ref. [50]). The impact of the substrate on the electron density inhomogeneity and thereby the SP resonance energy depends on the thickness and refractive index of the substrate, which may explain the quantitative agreement between theory and experiment reported in Ref. [10], since thinner substrates with smaller refractive indexes were used in their experiments. In order to completely address this issue, one would need to go beyond the dipole-dipole model for the substrate, thus future 3D EELS simulations taking nonlocal effects and/or inhomogeneous electron densities into account would be needed.

Another complementary explanation in the context of the inhomogeneity of the free-electron density could be the combined contribution of both the inhomogeneous static equilibrium electron density and nonlocality. It is well-known that the static equilibrium electron density is inhomogeneous, even in a semi-infinite metal [51], due to Friedel oscillations and the electron spill-out effect at the metal surface. The Friedel oscillations are modeled in the local quantum-confined model given by Eq. (3) while nonlocality is neglected, and *vice versa* in the nonlocal hydrodynamic model given by Eq. (3). As seen in Figure 2, the two effects separately give rise to similar-sized blueshifts, suggesting that the contribution of both effects simultaneously could add up to the significantly larger experimentally observed blueshift. Simply put, an extension of the nonlocal hydrodynamic model to include an inhomogeneous equilibrium free-electron density could produce a larger blueshift, which may be in accordance with the experimental observations. Furthermore, such a model could also take into account the electron spill-out effect, which in free-electron models has been argued to produce a redshift of the SP resonance [21, 50, 52–54], describing adequately simple metals. In contrast, it has also been shown that the spill-out effect in combination with the screening from the *d* electrons gives rise to the blueshift seen in Ag nanoparticles [55].

Additional size effects such as changes of the electronic band structure of the smallest nanoparticles, which are considerably more difficult to take into account, also impact the shift in SP resonance energy [6].

## 5 Conclusion

We have investigated the surface plasmon resonance of spherical silver nanoparticles ranging from 26 down to 3.5 nm in size with STEM EELS and observed a significant blueshift of 0.5 eV of the resonance energy. We have compared our experimental data with three different models based on the quasistatic optical polarizability of a sphere embedded in a homogeneous material. Two of the models, a nonlocal hydrodynamic model and a generalized local model, incorporate an inhomogeneity of the electron density induced by the quantum wave nature of the electrons. These two different models produce similar results in the SP resonance energy and describe qualitatively the blueshift observed in our measurements. Although our exact hydrodynamic generalization of the Clausius-Mossotti relation predicts a nonlocal blueshift that grows fast [as 1/(2*R*)] when decreasing the diameter and increases even faster for the smallest particles (2*R*<10 nm), the observed blueshifts are nevertheless larger than predicted.

The quantitative agreement between the two different theoretical models and the discrepancy with the larger observed blueshift suggest that a more detailed theoretical description of the system is needed to fully understand the influence of the substrate and the effect of the confinement of free electrons on the SP resonance shift in silver nanoparticles. On the experimental side, further EELS studies of other metallic materials and on different substrates could unveil the mechanism behind the size dependency of the SP resonance of nanometer scale particles.

We thank S. I. Bozhevolnyi for directing our attention to the theoretical model in Ref. [16] and G. Toscano for fruitful discussions. The Center for Nanostructured Graphene is sponsored by the Danish National Research Foundation, Project DNRF58. The A. P. Møller and Chastine Mc-Kinney Møller Foundation is gratefully acknowledged for the contribution toward the establishment of the Center for Electron Nanoscopy. N. Stenger acknowledges financial support by a Lundbeck Foundation Grant Nbr. R95-A10663.

## References

[1] Maier SA. Plasmonics: fundamentals and applications. New York: Springer; 2007.10.1007/0-387-37825-1Search in Google Scholar

[2] Schuller JA, Barnard ES, Cai W, Jun YC, White JS, Brongersma ML. Plasmonics for extreme light concentration and manipulation. Nat Mater 2010;9:193–204.10.1038/nmat2630Search in Google Scholar

[3] Kneipp K. Surface-enhanced Raman scattering. Phys Today 2007;60:40–6.10.1063/1.2812122Search in Google Scholar

[4] Wang F, Shen YR. General properties of local plasmons in metal nanostructures. Phys Rev Lett 2006;97:206806.10.1103/PhysRevLett.97.206806Search in Google Scholar

[5] Romero I, Aizpurua J, Bryant GW, García de Abajo FJ. Plasmons in nearly touching metallic nanoparticles: singular response in the limit of touching dimers. Opt Express 2006;14:9988–99.10.1364/OE.14.009988Search in Google Scholar

[6] Kreibig U, Genzel L. Optical absorption of small metallic particles. Surf Sci 1985;156:678–700.10.1016/0039-6028(85)90239-0Search in Google Scholar

[7] Charlé K-P, Schulze W, Winter B. The size dependent shift of the surface-plasmon absorption-band of small spherical metal particles. Z Phys D 1989;12:471–5.10.1007/BF01427000Search in Google Scholar

[8] Ouyang F, Batson P, Isaacson M. Quantum size effects in the surface-plasmon excitation of small metallic particles by electron-energy-loss spectroscopy. Phys Rev B 1992;46:15421–5.10.1103/PhysRevB.46.15421Search in Google Scholar PubMed

[9] Berciaud S, Cognet L, Tamarat P, Lounis B. Observation of intrinsic size effects in the optical response of individual gold nanoparticles. Nano Lett 2005;5:515–8.10.1021/nl050062tSearch in Google Scholar PubMed

[10] Scholl JA, Koh AL, Dionne JA. Quantum plasmon resonances of individual metallic nanoparticles. Nature 2012;483:421–7.10.1038/nature10904Search in Google Scholar PubMed

[11] Ciracì C, Hill RT, Mock JJ, Urzhumov Y, Fernández-Domínguez AI, Maier SA, Pendry JB, Chilkoti A, Smith DR. Probing the ultimate limits of plasmonic enhancement. Science 2012;337:1072–4.10.1126/science.1224823Search in Google Scholar

[12] Kern J, Grossmann S, Tarakina NV, Häckel T, Emmerling M, Kamp M, Huang J-S, Biagioni P, Prangsma J, Hecht B. Atomic-scale confinement of resonant optical fields. Nano Lett 2012;12:5504–9.10.1021/nl302315gSearch in Google Scholar

[13] Genzel L, Martin TP, Kreibig U, Dielectric function and plasma resonance of small metal particles. Z Phys B 1975;21:339–46.10.1007/BF01325393Search in Google Scholar

[14] Kraus WA, Schatz GC. Plasmon resonance broadening in small metal particles. J Chem Phys 1983;79:6130.10.1063/1.445794Search in Google Scholar

[15] Halperin WP, Quantum size effects in metal particles. Rev Mod Phys 1986;58:533–606.10.1103/RevModPhys.58.533Search in Google Scholar

[16] Keller O, Xiao M, Bozhevolnyi S. Optical diamagnetic polarizability of a mesoscopic metallic sphere: transverse self-field approach. Opt Comm 1993;102:238–44.10.1016/0030-4018(93)90389-MSearch in Google Scholar

[17] Öztürk ZF, Xiao S, Yan M, Wubs M, Jauho A-P, Mortensen NA. Field enhancement at metallic interfaces due to quantum confinement. J Nanophot 2011;5:051602.10.1117/1.3574159Search in Google Scholar

[18] Zuloaga J, Prodan E, Nordlander P. Quantum description of the plasmon resonances of a nanoparticle dimer. Nano Lett 2009;9:887–91.10.1021/nl803811gSearch in Google Scholar PubMed

[19] Mao L, Li Z, Wu B, Xu H. Effects of quantum tunneling in metal nanogap on surface-enhanced Raman scattering. Appl Phys Lett 2009;94:243102.10.1063/1.3155157Search in Google Scholar

[20] Esteban R, Borisov AG, Nordlander P, Aizpurua J. Bridging quantum and classical plasmonics with a quantum-corrected model. Nat Commun 2012;3:825.10.1038/ncomms1806Search in Google Scholar PubMed

[21] Ljungbert A, Lundqvist S. Non-local effects in the optical absorption of small metallic particles. Phys Rev Lett 1985;156:839–44.10.1016/0039-6028(85)90256-0Search in Google Scholar

[22] García de Abajo FJ. Nonlocal effects in the plasmons of strongly interacting nanoparticles, dimers, and waveguides. J Phys Chem C 2008;112:17983–7.10.1021/jp807345hSearch in Google Scholar

[23] David C, García de Abajo FJ. Spatial nonlocality in the optical response of metal nanoparticles. J Phys Chem C 2012;115:19470–5.10.1021/jp204261uSearch in Google Scholar

[24] Aizpurua J, Rivacoba A. Nonlocal effects in the plasmons of nanowires and nanocavities excited by fast electron beams. Phys Rev B 2008;78:035404.10.1103/PhysRevB.78.035404Search in Google Scholar

[25] Raza S, Toscano G, Jauho A-P, Wubs M, Mortensen NA. Unusual resonances in nanoplasmonic structures due to nonlocal response. Phys Rev B 2011;84:121412(R).10.1103/PhysRevB.84.121412Search in Google Scholar

[26] Toscano G, Raza S, Jauho A-P, Mortensen NA, Wubs M. Modified field enhancement and extinction in plasmonic nanowire dimers due to nonlocal response. Opt Express 2012;20:4176.10.1364/OE.20.004176Search in Google Scholar PubMed

[27] Fernández-Domínguez AI, Wiener A, García-Vidal FJ, Maier SA, Pendry JB. Transformation-optics description of nonlocal effects in plasmonic nanostructures. Phys Rev Lett 2012;108:106802.10.1103/PhysRevLett.108.106802Search in Google Scholar PubMed

[28] García de Abajo FJ. Optical excitations in electron microscopy. Rev Mod Phys 2010;82:209–75.10.1103/RevModPhys.82.209Search in Google Scholar

[29] Nelayah J, Kociak M, Stephan O, García de Abajo FJ, Tence M, Henrard L, Taverna D, Pastoriza-Santos I, Liz-Marzan LM, Colliex C. Mapping surface plasmons on a single metallic nanoparticle. Nat Phys 2007;3:348–53.10.1038/nphys575Search in Google Scholar

[30] Koh AL, Bao K, Khan I, Smith WE, Kothleitner G, Nordlander P, Maier SA, McComb DW. Electron energy-loss spectroscopy (EELS) of surface plasmons in single silver nanoparticles and dimers: influence of beam damage and mapping of dark modes. ACS Nano 2009;3:3015–22.10.1021/nn900922zSearch in Google Scholar PubMed

[31] Nicoletti O, Wubs M, Mortensen NA, Sigle W, van Aken PA, Midgley PA. Surface plasmon modes of a single silver nanorod: an electron energy loss study. Opt Express 2011;19:15371.10.1364/OE.19.015371Search in Google Scholar PubMed

[32] Koh AL, Fernández-Domínguez AI, McComb DW, Maier SA, Yang JKW. High-resolution mapping of electron-beam-excited plasmon modes in lithographically defined gold nanostructures. Nano Lett 2011;11:1323–30.10.1021/nl104410tSearch in Google Scholar PubMed

[33] Bloch F. Bremsvermögen von Atomen mit mehreren Elektronen. Z Phys A 1933;81:363–76.10.1007/BF01344553Search in Google Scholar

[34] Boardman A. Electromagnetic surface modes. Hydrodynamic theory of plasmon-polaritons on plane surfaces. Chichester: John Wiley and Sons; 1982.Search in Google Scholar

[35] Mulfinger L, Solomon SD, Bahadory M, Jeyarajasingam A, Rutkowsky SA, Boritz C. Synthesis and study of silver nanoparticles. J Chem Educ 2007;84:322–5.10.1021/ed084p322Search in Google Scholar

[36] Bååk T. Silicon oxynitride; a material for GRIN optics. Appl Opt 1982;21:1069–72.10.1364/AO.21.001069Search in Google Scholar PubMed

[37] Ruppin R. Optical properties of a plasma sphere. Phys Rev Lett 1973;31:1434–7.10.1103/PhysRevLett.31.1434Search in Google Scholar

[38] Pendry JB, Aubry A, Smith DR, Maier SA. Transformation optics and subwavelength control of light. Science 2012;337:549–52.10.1126/science.1220600Search in Google Scholar PubMed

[39] Dasgupta BB, Fuchs R. Polarizability of a small sphere including nonlocal effects, Phys Rev B 1981;24:554–61.10.1103/PhysRevB.24.554Search in Google Scholar

[40] Fuchs R, Claro F. Multipolar response of small metallic spheres: nonlocal theory. Phys Rev B 1987;35:3722–7.10.1103/PhysRevB.35.3722Search in Google Scholar

[41] Rakić AD, Djurišić AB, Elazar JM, Majewski ML. Optical properties of metallic films for vertical-cavity optoelectronic devices. Appl Opt 1998;37:5271–83.10.1364/AO.37.005271Search in Google Scholar PubMed

[42] Noguez C. Surface plasmons on metal nanoparticles: the influence of shape and physical environment. J Phys Chem C 2007;111:3806–19.10.1021/jp066539mSearch in Google Scholar

[43] Ruppin R. Optical absorption by a small sphere above a substrate with inclusion of nonlocal effects. Phys Rev B 1992;45:11209–15.10.1103/PhysRevB.45.11209Search in Google Scholar

[44] Li Z, Bao K, Fang Y, Guan Z, Halas NJ, Nordlander P, Xu H. Effect of a proximal substrate on plasmon propagation in silver nanowires. Phys Rev B 2010;82:1.10.1103/PhysRevB.82.241402Search in Google Scholar

[45] Zhang S, Xu H. Optimizing substrate-mediated plasmon coupling toward high-performance plasmonic nanowire waveguides. ACS Nano 2012;6:8128–35.10.1021/nn302755aSearch in Google Scholar

[46] Yamaguchi T, Yoshida S, Kinbara A. Optical effect of the substrate on the anomalous absorption of aggregated silver films. Thin Solid Films 1974;21:173–87.10.1016/0040-6090(74)90099-6Search in Google Scholar

[47] Jain PK, Huang W, El-Sayed MA. On the universal scaling behavior of the distance decay of plasmon coupling in metal nanoparticle pairs: a plasmon ruler equation. Nano Lett 2007;7:2080–8.10.1021/nl071008aSearch in Google Scholar

[48] Novotny L, Hecht B. Principles of nano-optics. New York: Cambridge; 2006.10.1017/CBO9780511813535Search in Google Scholar

[49] Ruppin R. Surface modes and optical absorption of a small sphere above a substrate. Surf Sci 1983;127:108–18.10.1016/0039-6028(83)90402-8Search in Google Scholar

[50] Boardman AD, Paranjape BV. The optical surface modes of metal spheres. J Phys F Met Phys 1977;7:1935.10.1088/0305-4608/7/9/036Search in Google Scholar

[51] Lang ND, Kohn W. Theory of metal surfaces: charge density and surface energy. Phys Rev B 1970;1:4555–68.10.1103/PhysRevB.1.4555Search in Google Scholar

[52] Ascarelli P, Cini M. “Red shift” of the surface plasmon resonance absorption by fine metal particles. Solid State Commun 1976;18:385–8.10.1016/0038-1098(76)90029-6Search in Google Scholar

[53] Ruppin R. Plasmon frequencies of small metal spheres. J Phys Chem Solids 1978;39:233–7.10.1016/0022-3697(78)90048-3Search in Google Scholar

[54] Apell P, Ljungbert Å. Red shift of surface plasmons in small metal particles. Solid State Commun 1982;44:1367–9.10.1016/0038-1098(82)90895-XSearch in Google Scholar

[55] Liebsch A. Surface-plasmon dispersion and size dependence of Mie resonance: silver versus simple metals. Phys Rev B 1993;48:11317–28.10.1103/PhysRevB.48.11317Search in Google Scholar

**Received:**2012-10-09

**Accepted:**2013-02-19

**Published Online:**2013-03-23

**Published in Print:**2013-04-01

©2013 by Science Wise Publishing & De Gruyter Berlin Boston