Abstract
The fact that surface-induced damping rate of surface plasmon polaritons (SPPs) in metal nanoparticles increases with the decrease of particle size is well known. We show that this rate also increases with the degree of the mode confinement, hence damping of the higher order nonradiative SPP modes in spherical particles is greatly enhanced relative to damping of the fundamental (dipole) SPP mode. Since higher order modes are the ones responsible for quenching of luminescence in the vicinity of metal surfaces, the degree of quenching increases resulting in a substantial decrease in the amount of attainable enhancement of the luminescence.
© 2015 Optical Society of America
1. Introduction
It is well known that emission properties of nanometer-scale objects, such as atoms, molecules, or quantum dots placed near metal surfaces get modified. This modification stems from interaction between the emitting object and collective oscillations of free-electrons in the metal called surface plasmon polaritons (SPP’s) [1, 2]. SPPs are characterized by strong electric fields localized near metal surfaces, and if the dimensions of the metal nanoparticles are sub-wavelength, so is the spatial extent of the localized field [3]. Enhancements of radiative decay rates [4], fluorescence [5–9], Raman scattering [10–12] and others near metal nanoparticles have been demonstrated by numerous groups. Spectacular enhancement of Raman scattering by metal nanoparticles with far more orders of magnitude [10–13] has since made surface enhanced Raman scattering (SERS) a practical sensing technology, but when it comes to the luminescence enhancement the progress has been less impressive [14]. Not only the emission efficiency suffers from non-radiative decay (damping) of SPP modes [15,16], but the situation is exacerbated by the process of life-time quenching [17] in which the energy of the emitting object (“emitter”) gets coupled into the higher-order so-called “dark” SPP modes that do not radiate and hence dissipates inside the metal [18].
The modes are “dark” because the electric field in them undulates (changes sign) on much shorter scales than the wavelength of radiation, hence they are poorly coupled to the radiative modes. This is the case for the large wavevector propagating SPP modes on the metal-dielectric interface, or higher order localized SPP modes in metal nanoparticles, which are the subject of this work. The dark modes can of course be observed using electron loss spectroscopy [19] or near field optical imaging [20]. Previously in [18,21,22] we have developed the theory of luminescence quenching due to excitation of higher order modes. We have shown that as emitting object approaches the metal surface, the coupling into higher modes gets enhanced relative to the coupling into the radiative lowest order (dipole) mode. As larger fraction of the energy gets coupled from the emitter into the dark modes, the luminescence gets quenched. Since the resonant frequencies of higher order modes are blue-shifted relative to the dipole mode, the coupling (and hence quenching) critically depends on the detuning and linewidth of SPP resonances.
In the first-order approximation SPP mode can be treated as a harmonic oscillator whose Lorentzian linewidth is equal to the damping rate that has two constituents –radiative decay rate that depends on geometry (and is vanishingly small for all non-radiative higher-order modes), and the non-radiative rate specific to the metal, the same rate that enters the Drude expression for the dielectric constant of the metal . Using the language of classical physics, this “bulk” damping rate represents ohmic losses in the metal. Using the “quantum” language, one can describe it as a combination of processes of an SPP (or a photon) annihilation and creation of one and more electron-hole pairs in the metal. Since both energy and momentum must be conserved in this process, a mechanism that compensates for the momentum mismatch must be involved. This can be a scattering by phonons, defects, and impurities, or electron-electron scattering [23] and all the SPP modes experience the same “bulk” non-radiative damping . Yet as early as 1980’s [24] it was demonstrated that as metallic objects become smaller there appears a new size-dependent contribution to the non-radiative damping which, classically speaking, can be associated with the electron scattering on the surface of nanoparticle. Using this phenomenological description Kreibig and Vollmer [24] introduced the “surface collision scattering rate” , where is a characteristic dimension of nanoparticle and is a Fermi velocity. On the quantum level this “surface collision damping” [25] is nothing but the Landau damping [26–28], a process in which large wave-vector components of the electro-magnetic field can directly create electron-hole pairs in the metal without need for any additional momentum-matching mechanisms. The same phenomenon can also be described by the “nonlocality”, i.e., spatial dispersion of the metal dielectric constant [29,30]. While simple phenomenological picture of surface collision damping relates it to the size of the metal nanoparticles, i.e., to the electron confinement, quantum description [31] shows that the damping is determined by the photon confinement, i.e. the spatial extent of the electric field inside the metal nanoparticle. For the fundamental (dipole) SPP mode, the spatial extent of the field inside the nanoparticle is equal to its size, so two approaches give exactly the same result. In the higher order modes, however, the field is confined near the surface, and the spatial extent of the mode of the -th order decreases, roughly as shown in Section 2 below. It is then reasonably to expect the surface damping rate of the -th order mode to increase as where is the radius of the metal nanosphere, causing the higher order modes to broaden in the spectral domain. That in turn strongly affects the rate of coupling of energy from the emitters placed near the surface and has impact on quenching. In this work we evaluate this impact and come to conclusion that surface collisions considerably increase the quenching and reduce the amount of attainable enhancement of the luminescence.
2. Theory of confinement-dependent damping
Consider the spherical metal nanoparticle of radius embedded in a medium with dielectric constant shown in Fig. 1(a). Specifically, we shall treat the nano-structure in which either Au or Ag nanospheres are embedded in the wide bandgap GaN. Each nanoparticle supports a large number of SPP modes with resonant frequencies, ranging from the for the lowest order radiating dipole mode to for the highest order modes that behave as surface plasmons propagating along the nanoparticle surface.
If the polar axis z is associated with the polarization of the emitter as in Fig. 1(a), only modes with have non-vanishing field along this axis,
Where and are unit vectors in polar coordinates, andis a Legendre polynomial. The radial dependence of the electric field for different modes is shown in Fig. 1(b) along z axis () where the energy of each mode has been normalized to . Since the higher order mode occupies progressively smaller volume [18] the amplitudes of higher order modes increase with and if an emitter is placed close enough to the metal surface to interact with these dark high -modes the energy will be very efficiently transferred into these modes and subsequently dissipated. But in order for that to happen, just having spatial overlap is not enough as the spectral overlap is also required, as shown in Fig. 1(c) for Ag nanosphere where each mode is represented by Lorentzian with full width at half maximum (FWHM) equal to damping which for mode is equal to where is the radiative decay rate of the nanoparticle proportional to the nanoparticle volume [18], and for all other modes . All the modes have the same bulk damping rate , but the surface collisions are far more prevalent in the higher order modes. One can make a very rough estimate of these rates by evaluating the mean penetration depth of the radial field into the metal aswhich should result in surface collision damping rate of where according to [24] is constant of the order of unity.To obtain a more precise value of the dependence of on the order of the mode we consider the Lindhard formula for the longitudinal dielectric constant of the metal [31]
where is a plasma frequency and is the Fermi velocity of the metal. For the small wavevectors the dielectric constant is real, but once , the logarithm in (3) becomes imaginary number, and the dielectric constant becomes complex, acquiring imaginary part . Thus in the spatial Fourier transform of the -th mode taken inside the metal the longitudinal spectral components with will get absorbed, as shown in Fig. 2(a), and the surface collision contribution to the imaginary part of the effective dielectric constant of this mode will becomewherewith similar expressions for and , and the power spectrum has been normalized as . Using Drude formula, the surface contribution to imaginary part of dielectric constant can be related to the effective surface scattering rate as . Normalizing the wavevector to the diameter of the nanoparticle, and introducing we obtain from Eq. (4) , where ,And the dimensionless is normalized as Fermi velocity in gold or silver is . For the frequencies in the visible range (say 500 nm) hence for we have . If longitudinal field is confined on the scale then its power spectrum is expected to decay at large values of as and integration will yield that does not depend on and. To check this proposition we performed numerical calculations of for a wide range of (corresponding to the range of from 2 to 100nm) and found that indeed these coefficients are constants and hence their sum is just a number, plotted in Fig. 2(b). It is clear that the coefficient can be very well approximated as and the collision scattering damping as which is only slightly different from the phenomenological estimate made from Eq. (2). The result of this surface damping adding on to the bulk damping rate leads to increase of Lorentzian FWHM of various SPP modes as shown in Fig. 1(d) with higher order ones being progressively broadened.It should be noted that results similar to one shown in Fig. 2(a) have been recently obtained [32] using the hydrodynamic model [30,33] in which Landau damping contribution is represented by a phenomenological diffusion term. In the hydrodynamic model electron density is co-ordinate dependent while in our model the damping occurs in the reciprocal space. Yet both models predict suppression of the large wavevector components of electron wavefunction and therefore lead to similar results. Applying the theory derived here to more complicated, non-spherical particles will allow a more comprehensive comparison between two theories.
3. Impact on the quenching coefficient
In comparison with Fig. 1(c), the excitation spectra of different SPP modes of Ag nanosphere with surface collision broadening taken into account as shown in Fig. 1(d) shows that higher order modes get broadened considerably. For example, spectrum of mode broadens by a factor of 5 and its overlap with the emission line of the emitter, centered at frequency increases accordingly. The coupling of the emitter’s energy into the l-th order SPP mode is characterized by the Purcell factor [22],
where Lorentzian line shape of the mode is . Since the emission line is typically resonant with the dipole mode, , it is clear that for as long as , i.e. practically for all the dark modes with , the increase in damping causes increase in Purcell factor . For the l = 1 bright dipole mode on the other hand, increase in damping causes reduction of . The changes caused by increased surface collision damping is well described by the “quenching” ratio [18], i.e., the ratio of summation of all Purcell factors associated with all the dark modes () to the Purcell factor for the “bright” dipole mode (). The quenching factors for gold and silver nanospheres with radius nm are plotted versus the separation distance of the emitter in Figs. 3(a) and 3(b) respectively. The “Bulk” curves in each figure are describing the quenching under assumption that only the “bulk” non-radiative damping is contributing. The curves “Surface” include only the surface contribution to damping and the curves “Bulk + Surface” correspond to the realistic description of damping that includes both bulk and surface contributions.Other than the obvious quenching effect being stronger for smaller separation between the metal nanoparticle and the emitter, quenching increases when surface collision damping is taken into account. While the quenching ratio for Ag is generally smaller than that of Au because of its smaller ohmic losses, the impact of surface damping is more significant for Ag given the fact that both nanoparticles of equal size have similar surface collision damping rate , but this rate represents a high fraction of the total damping rate for Ag. Another interesting observation is that even if we had ways to magically eliminate the “bulk” ohmic losses in the metal, we would still have been stuck with surface damping and quenching would still take place because this surface effect is inherently associated with the field localization of the various orders of the SPP modes.
4. Optimization of luminescence enhancement
Once the quenching factor had been determined one can evaluate the luminescence enhancement factor for an emitter with an original radiative efficiency placed in the vicinity of meal nanoparticle [18]
where is the out-coupling efficiency of dipole mode. All the factors (, , and ) depend on the nanosphere radius . Hence one can perform an optimization to determine the ideal size of the nanosphere that is small enough to have a small effective mode volume assuring large Purcell factor of the dipole SPP mode in (7), yet is still sufficiently large to assure a decent radiative out-coupling efficiency of this mode. Furthermore, according to (7) if the quenching through high order modes could be neglected, it would have been always better to have emitters placed as close as possible to the metal nanoparticle. But as shown in previous section the quenching also increases as the separation distance d gets shorter. Hence the distance can also be optimized to achieve a favorable compromise of good coupling into the dipole mode and adequate suppression of the luminescence quenching by higher order modes. The dependence of the enhancement by Au nanosphere on and is shown in inset of Fig. 4(a) for the emitter with . The results of optimization of luminescence enhancement for Au and Ag nanospheres are shown in Fig. 4(a) and (b) respectively. The optimization is performed for a wide range of emitter radiative efficiencies. The optimal sphere radius and separation at which the maximum attainable enhancement is obtained for each are shown in Figs. 4(c)-4(f). The three graphs in each plot correspond to the same three cases of damping as in Fig. 3.As expected [16], the maximum attainable enhancement deteriorates quite rapidly with the increase of . With the inclusion of surface collision damping, the attainable enhancement is reduced for both Au and Ag nanoparticles, but the impact is much stronger for Ag for the same reason as explained in section 3. Also, the optimized values of radius and separation distances also increase once surface collision damping is taken into account. The increase in is easy to understand – the surface collision damping is inversely proportional to the radius. The modest increase in obviously reduces coupling into the broadened higher order modes, while keeping coupling into the dipole mode reasonably high.
As one can see that for noble metal spheres the changes imposed by the surface collision damping are significant. For Ag the quenching ratio [Fig. 3(b)] is increased by an order of magnitude and the maximum enhancement [Fig. 4(b)] is reduced by a factor of 3.5. The reason for this is that even when Ag sphere is embedded in high index dielectric the dipole SPP resonance is near 2.3eV where the bulk damping rate in Ag [34] since that energy is close to the onset of interband transitions. Since the surface collision rate is roughly , in the 10 nm radius nanosphere this rate is comparable to the bulk rate for the quadrupole mode . If, however, one uses nanoshells [35] or dimers [36] the SPP resonances can be moved to the red part of the spectrum, say 700nm where and obviously the influence of surface collision damping will grow four-fold leading to much larger changes in quenching ratios.
If one considers other SPP structures, whether it is a simple interface SPP, an SPP in the metal grooves [37], or other design, there are always some highly confined modes whose Fourier spectrum contains larger fraction of large-k spectral components - these modes are not coupled well with the radiating modes and the surface collisions will also lead to additional fluorescence quenching in them. One should expect a similar increase in quenching for the strong coupling regime [38] as surface collision damping will make the effective coupling weaker.
Once again, it should be noted that even in the perfect world where one could eliminate all “bulk” ohmic loss of any metal, the enhancement would never be “perfect”, as will be limited by the surface damping as shown in Figs. 4(a) and 4(b) for both Au and Ag nanoparticles, respectively.
5. Conclusion
In conclusion, we have shown that high order SPP modes in metal nanoparticles are subject to much stronger damping than the lowest order dipole mode. This additional damping, roughly proportional to the order of the mode is caused by the confinement of the higher order mode in the vicinity of the surface. Increased damping causes spectral broadening of the higher order SPP resonances and increases the probability of their excitation by the luminescing object placed near the nanoparticle. Since the higher order modes are dark the energy gets trapped and eventually dissipated in them and the luminescence experiences additional quenching. As result, maximum enhancement of luminescence attainable in the vicinity of small metal nanoparticles gets reduced by a factor of 3 or 4 in the blue green range and even larger factors in the red part of spectrum.
Acknowledgments
JK acknowledges fruitful discussions with Dr. P. Noir. GS acknowledges support from Air Force Office of Scientific Research (AFOSR) (FA9550-14-1-0196 Dr. Gernot Pomrenke, Program Manager).
References and links
1. M. I. Stockman, “Nanoplasmonics: past, present, and glimpse into future,” Opt. Express 19(22), 22029–22106 (2011). [CrossRef] [PubMed]
2. S. A. Maier, Plasmonics: Fundamentals and Applications (Springer, 2007).
3. W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature 424(6950), 824–830 (2003). [CrossRef] [PubMed]
4. K. Okamoto, I. Niki, A. Scherer, Y. Narukawa, T. Mukai, and Y. Kawakami, “Surface plasmon enhanced spontaneous emission rate of InGaN∕GaN quantum wells probed by time-resolved photoluminescence spectroscopy,” Appl. Phys. Lett. 87(7), 071102 (2005). [CrossRef]
5. S. Kühn, U. Håkanson, L. Rogobete, and V. Sandoghdar, “Enhancement of single-molecule fluorescence using a gold nanoparticle as an optical nanoantenna,” Phys. Rev. Lett. 97(1), 017402 (2006). [CrossRef] [PubMed]
6. L. Novotny, “Single molecule fluorescence in inhomogeneous environments,” Appl. Phys. Lett. 69(25), 3806 (1996). [CrossRef]
7. S. Kühn, G. Mori, M. Agio, and V. Sandoghdar, “Modification of single molecule fluorescence close to a nanostructure: radiation pattern, spontaneous emission and quenching,” Mol. Phys. 106(7), 893–908 (2008). [CrossRef]
8. P. Bharadwaj and L. Novotny, “Spectral dependence of single molecule fluorescence enhancement,” Opt. Express 15(21), 14266–14274 (2007). [CrossRef] [PubMed]
9. R. Bardhan, N. K. Grady, J. R. Cole, A. Joshi, and N. J. Halas, “Fluorescence enhancement by Au nanostructures: nanoshells and nanorods,” ACS Nano 3(3), 744–752 (2009). [CrossRef] [PubMed]
10. S. Nie and S. R. Emory, “Probing single molecules and single nanoparticles by surface-enhanced Raman scattering,” Science 275(5303), 1102–1106 (1997). [CrossRef] [PubMed]
11. A. M. Michaels, M. Nirmal, and L. E. Brus, “Surface enhanced Raman spectroscopy of individual Rhodamine 6G molecules on large Ag nanocrystals,” J. Am. Chem. Soc. 121(43), 9932–9939 (1999). [CrossRef]
12. Z. Wang, S. Pan, T. D. Krauss, H. Du, and L. J. Rothberg, “The structural basis for giant enhancement enabling single-molecule Raman scattering,” Proc. Natl. Acad. Sci. U.S.A. 100(15), 8638–8643 (2003). [CrossRef] [PubMed]
13. E. C. Le Ru and P. G. Etchegoin, Principles of Surface Enhanced Raman Spectroscopy and Related Plasmonic Effects (Elsevier, Amsterdam, 2009).
14. M. Bakker, V. P. Drachev, Z. Liu, H.-K. Yuan, R. H. Pedersen, A. Boltasseva, J. Chen, J. Irudayaraj, A. V. Kildishev, and V. M. Shalaev, “Nanoantenna array-induced fluorescence enhancement and reduced lifetimes,” New J. Phys. 10(12), 125022 (2008). [CrossRef]
15. R. Carminati, J.-J. Greffet, C. Henkel, and J. M. Vigoureux, “Radiative and non-radiative decay of a single molecule close to a metallic nanoparticle,” Opt. Commun. 261(2), 368–375 (2006). [CrossRef]
16. G. Sun, J. B. Khurgin, and R. A. Soref, “Practical enhancement of photoluminescence by metal nanoparticles,” Appl. Phys. Lett. 94(10), 101103 (2009). [CrossRef]
17. E. Dulkeith, A. C. Morteani, T. Niedereichholz, T. A. Klar, J. Feldmann, S. A. Levi, F. C. van Veggel, D. N. Reinhoudt, M. Möller, and D. I. Gittins, “Fluorescence quenching of dye molecules near gold nanoparticles: radiative and nonradiative effects,” Phys. Rev. Lett. 89(20), 203002 (2002). [CrossRef] [PubMed]
18. G. Sun, J. B. Khurgin, and C. C. Yang, “Impact of high-order surface Plasmon modes of metal nanoparticles on enhancement of optical emission,” Appl. Phys. Lett. 95(17), 171103 (2009). [CrossRef]
19. F.-P. Schmidt, H. Ditlbacher, U. Hohenester, A. Hohenau, F. Hofer, and J. R. Krenn, “Dark plasmonic breathing modes in silver nanodisks,” Nano Lett. 12(11), 5780–5783 (2012). [CrossRef] [PubMed]
20. R. Esteban, R. Vogelgesang, J. Dorfmüller, A. Dmitriev, C. Rockstuhl, C. Etrich, and K. Kern, “Direct near-field optical imaging of higher order plasmonic resonances,” Nano Lett. 8(10), 3155–3159 (2008). [CrossRef] [PubMed]
21. G. Sun, J. B. Khurgin, and D. P. Tsai, “Comparative analysis of photoluminescence and Raman enhancement by metal nanoparticles,” Opt. Lett. 37(9), 1583–1585 (2012). [CrossRef] [PubMed]
22. G. Sun and J. B. Khurgin, “Origin of giant difference between fluorescence, resonance, and nonresonance Raman scattering enhancement by surface plasmons,” Phys. Rev. A 85(6), 063410 (2012). [CrossRef]
23. J. B. Khurgin, “How to deal with the loss in plasmonics and metamaterials,” Nat. Nanotechnol. 10(1), 2–6 (2015). [CrossRef] [PubMed]
24. U. Kreibig and M. Vollmer, Optical Properties of Metal Clusters (Springer-Verlag, Berlin, 1995).
25. A. V. Uskov, I. E. Protsenko, N. A. Mortensen, and E. P. O’Reilly, “Broadening of plasmonic resonance due to electron collisions with nanoparticle boundary: a quantum mechanical consideration,” Plasmonics 9(1), 185–192 (2014). [CrossRef]
26. C. Yannouleas and R. A. Broglia, “Landau damping and wall dissipation in large metal clusters,” Ann. Phys. 217(1), 105–141 (1992). [CrossRef]
27. R. A. Molina, D. Weinmann, and R. A. Jalabert, “Oscillatory size dependence of the surface plasmon linewidth in metallic nanoparticles,” Phys. Rev. B 65(15), 155427 (2002). [CrossRef]
28. Z. Yuan and S. Gao, “Landau damping and lifetime oscillation of surface plasmons in metallic thin films studied in a jellium slab model,” Surf. Sci. 602(2), 460–464 (2008). [CrossRef]
29. F. J. García de Abajo, “Nonlocal effects in the plasmons of strongly interacting nanoparticles, dimers, and waveguides,” J. Phys. Chem. C 112(46), 17983–17987 (2008). [CrossRef]
30. N. A. Mortensen, S. Raza, M. Wubs, T. Søndergaard, and S. I. Bozhevolnyi, “A generalized non-local optical response theory for plasmonic nanostructures,” Nat. Commun. 5, 3809 (2014). [CrossRef] [PubMed]
31. J. Lindhard, “On the properties of a gas of charged particles,” Mat. Fys. Medd. Dan. Vid. 28, 1–57 (1954).
32. S. Raza, S. I. Bozhevolnyi, M. Wubs, and N. Asger Mortensen, “Nonlocal optical response in metallic nanostructures,” J. Phys. Condens. Matter 27(18), 183204 (2015). [CrossRef] [PubMed]
33. T. Christensen, W. Yan, S. Raza, A. P. Jauho, N. A. Mortensen, and M. Wubs, “Nonlocal response of metallic nanospheres probed by light, electrons, and atoms,” ACS Nano 8(2), 1745–1758 (2014). [CrossRef] [PubMed]
34. P. B. Johnson and R. W. Christy, “Optical constants of noble metals,” Phys. Rev. B 6(12), 4370–4379 (1972). [CrossRef]
35. F. Tam, G. P. Goodrich, B. R. Johnson, and N. J. Halas, “Plasmonic enhancement of molecular fluorescence,” Nano Lett. 7(2), 496–501 (2007). [CrossRef] [PubMed]
36. O. L. Muskens, V. Giannini, J. A. Sánchez-Gil, and J. Gómez Rivas, “Optical scattering resonances of single and coupled dimer plasmonic nanoantennas,” Opt. Express 15(26), 17736–17746 (2007). [CrossRef] [PubMed]
37. S. Raza, N. Stenger, A. Pors, T. Holmgaard, S. Kadkhodazadeh, J. B. Wagner, K. Pedersen, M. Wubs, S. I. Bozhevolnyi, and N. A. Mortensen, “Extremely confined gap surface-plasmon modes excited by electrons,” Nat. Commun. 5, 4125 (2014). [CrossRef] [PubMed]
38. A. Delga, J. Feist, J. Bravo-Abad, and F. J. Garcia-Vidal, “Quantum emitters near a metal nanoparticle: strong coupling and quenching,” Phys. Rev. Lett. 112(25), 253601 (2014). [CrossRef] [PubMed]