Nonlinear pulsed excitation of high-<i>Q</i> optical modes of plasmonic nanocavities

doi:10.1364/OE.18.017165

« Show journal navigation

Nonlinear pulsed excitation of high-Q optical modes of plasmonic nanocavities

Claudiu G. Biris and Nicolae C. Panoiu »View Author Affiliations

Optics Express, Vol. 18, Issue 16, pp. 17165-17179 (2010)
http://dx.doi.org/10.1364/OE.18.017165

View Full Text Article

Acrobat PDF (1522 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Description of the numerical method
3. Results and discussions
4. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Multimedia
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (49)
Cited By
Figures (8)
Multimedia (2)
Metrics

Abstract

We present a comprehensive theoretical and numerical analysis of the physical mechanisms pertaining to the nonlinear pulsed excitation of optical modes in plasmonic cavities made of metallic nanowires. Our analysis is based on extensive numerical simulations carried out both in the frequency and time domains. The numerical algorithm used in our study is based on the multiple scattering method and allows us to include in our analysis the effects of both the surface and bulk nonlinear polarizations generated at the second harmonic (SH). In particular, we investigate the physical roperties of plasmonic modes excited at the SH as the result of the interaction of femtosecond optical pulses with plasmonic nanocavities. We show that such cavities have two distinct types of modes, namely, plasmonic cavity modes and multipole plasmon modes generated via the hybridization of modes of single nanowires. Our analysis reveals that the properties of the latter modes depend only weakly on the cavity geometry, whereas the lifetime and quality factor of plasmonic cavity modes vary considerably with the system parameters.

1. Introduction

In recent years, we have witnessed a paradigmatic shift in the research methodologies used in studying the physical properties of electromagnetic media. This development can be primarily traced to the advent of the research field of metamaterials and its applications to optics and materials science. Metamaterials are a new class of electromagnetic media, whose basic unit cell is artificially engineered so as the optical medium possesses pre-designed macroscopic physical properties tailored for specific applications. In particular, by engineering the material parameters and geometry of the unit cell at a scale that is comparable to or smaller than the operating wavelength, a suite of new materials with remarkable properties have recently emerged. Examples of such optical media include materials that are magnetically active at terahertz and optical frequencies [1

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004). [CrossRef] [PubMed]

], materials with negative index of refraction [2

R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental Verification of a Negative Index of Refraction,” Science 292, 77–79 (2001). [CrossRef] [PubMed]

, 3

S. Zhang, W. Fan, N. C. Panoiu, K. J. Malloy, R. M. Osgood, and S. R. J. Brueck, “Experimental Demonstration of Near-Infrared Negative-Index Metamaterials,” Phys. Rev. Lett. 95, 137404 (2005). [CrossRef] [PubMed]

, 4

V. M. Shalaev, W. Cai, U. K. Chettiar, H-K. Yuan, A. K. Sarychev, V. P. Drachev, and A. V. Kildishev, “Negative index of refraction in optical metamaterials,” Opt. Lett. 30, 3356–3358 (2005). [CrossRef]

], and homogeneous or periodic media with low-index of refraction [5

R. W. Ziolkowski, “Propagation in and scattering from a matched metamaterial having a zero index of refraction,” Phys. Rev. E 70, 046608 (2004). [CrossRef]

, 6

N. C. Panoiu, R. M. Osgood, S. Zhang, and S. R. J. Brueck, “Zero-n bandgap in photonic crystal superlattices,” J. Opt. Soc. Am. B 23, 506–513 (2006). [CrossRef]

, 7

S. Kocaman, R. Chatterjee, N. C. Panoiu, J. F. McMillan, M. B. Yu, R. M. Osgood, D. L. Kwong, and C. W. Wong, “Observation of Zeroth-Order Band Gaps in Negative-Refraction Photonic Crystal Superlattices at Near-Infrared Frequencies,” Phys. Rev. Lett. 102, 203905 (2009). [CrossRef] [PubMed]

It is well known that metallic structures support localized and propagating electromagnetic excitations, also known as surface plasmon polariton (SPP) modes [8

U. Kreibig and M. Vollmer, Optical Properties of Metal Clusters (Springer-Verlag, Berlin, 1995).

, 9

A. V. Zayats, I. I. Smolyaninov, and A. A. Maradudin, “Nano-optics of surface plasmon polaritons,” Phys. Rep. 408, 131–314 (2005). [CrossRef]

]. As a consequence of the resonant excitation of these modes the electromagnetic field in the vicinity of metal surfaces is strongly enhanced and therefore the macroscopic properties of metal-based metamaterials are significantly affected [10

W. Fan, S. Zhang, N. C. Panoiu, A. Abdenour, S. Krishna, R. M. Osgood, K. J. Malloy, and S. R. J. Brueck, “Second Harmonic Generation from a Nanopatterned Isotropic Nonlinear Material,” Nano Lett. 6, 1027–1030 (2006). [CrossRef]

, 11

J. A. H. van Nieuwstadt, M. Sandtke, R. H. Harmsen, F. B. Segerink, J. C. Prangsma, S. Enoch, and L. Kuipers, “Strong Modification of the Nonlinear Optical Response of Metallic Subwavelength Hole Arrays,” Phys. Rev. Lett. 97, 146102 (2006). [CrossRef] [PubMed]

, 12

L. Cao, N. C. Panoiu, and R. M. Osgood, “Surface second-harmonic generation from surface plasmon waves scattered by metallic nanostructures,” Phys. Rev. B 75, 205401 (2007). [CrossRef]

, 13

Y. Zeng, W. Hoyer, J. J. Liu, S. W. Koch, and J. V. Moloney, “Classical theory for second-harmonic generation from metallic nanoparticles,” Phys. Rev. B 79, 235109 (2009). [CrossRef]

, 14

I. I. Smolyaninov, A. V. Zayats, A. Gungor, and C. C. Davis, “Single-Photon Tunneling via Localized Surface Plasmons,” Phys. Rev. Lett. 88, 187402 (2002). [CrossRef] [PubMed]

, 15

X. W. Wang, G. C. Schatz, and S. K. Gray, “Ultrafast pulse excitation of a metallic nanosystem containing a Kerr nonlinear material,” Phys. Rev. B 74, 195439 (2006). [CrossRef]

]. This effect has been successfully employed in several applications, i.e., surface-enhanced Ramman scattering [16

K. Kneipp, Y. Wang, H. Kneipp, L. T. Perelman, I. Itzkan, R. R. Dasari, and M. S. Feld, “Single Molecule Detection Using Surface-Enhanced Raman Scattering (SERS),” Phys. Rev. Lett. 78, 1667 (1997). [CrossRef]

], nanoscale optical antennae [17

T. H. Taminiau, R. J. Moerland, F. B. Segerink, L. Kuipers, and N. F. van Hulst, “ λ/4 Resonance of an Optical Monopole Antenna Probed by Single Molecule Fluorescence,” Nano Lett. 7, 28 (2007). [CrossRef] [PubMed]

], and metallic nanotips for near-field optical microscopy [18

B. Knoll and F. Keilmann, “Near-field probing of vibrational absorption for chemical microscopy,” Nature (London) 399, 134–137 (1999). [CrossRef]

, 19

R. M. Roth, N. C. Panoiu, M. M. Adams, R. M. Osgood, C. C. Neacsu, and M. B. Raschke, “Resonant-plasmon field enhancement from asymmetrically illuminated conical metallic-probe tips,” Opt. Express 14, 2921–2913 (2006). [CrossRef] [PubMed]

]. In addition, enhancement of the second harmonic (SH) field intensity as large as 10⁸ has been observed when employing metallic nanoparticles [20

S. I. Bozhevolnyi, J. Beermann, and V. Coello, “Direct Observation of Localized Second-Harmonic Enhancement in Random Metal Nanostructures,” Phys. Rev. Lett. 90, 197403 (2003). [CrossRef] [PubMed]

]. It is thus clear that achieving a complete haracterization of nonlinear effects occurring at the nanoscale represents an important step towards designing future subwavelength active nanodevices, which can be operated at low optical power.

One important property of localized SPPs is that their physical properties, such as, resonance frequency, optical loss, field enhancement, are strongly dependent on the geometry and material parameters of the corresponding plasmonic nanostructures. Therefore, similar to the effective, macroscopic optical properties of metamaterials, the linear and nonlinear optical response of plasmonic nanostructures can be effectively tailored to specific applications by, e.g., molding the spatial distribution of the optical near-field so as to provide pre-defined functionalities. For example, previous work has shown that plasmonic nanostructures can be used to significantly enhance the optical absorption (efficiency) of solar cells, detectors, and other photovoltaic devices [21

T. Ishi, J. Fujikata, K. Makita, T. Baba, and K. Ohashi, “Si nano-photodiode with a surface plasmon antenna,” Jpn. J. Appl. Phys. 44, L364–L366 (2005). [CrossRef]

, 22

Z. Yu, G. Veronis, S. Fan, and M. L. Brongersma, “Design of midinfrared photodetectors enhanced by surface plasmons on grating structures,” Appl. Phys. Lett. 89, 151116 (2006). [CrossRef]

, 23

R. D. R. Bhat, N. C. Panoiu, S. R. J. Brueck, and R. M. Osgood, “Enhancing the signal-to-noise ratio of an infrared photodetector with a circular metal grating,” Opt. Express 16, 4588–4596 (2008). [CrossRef] [PubMed]

, 24

D. M. Schaadt, B. Feng, and E. T. Yu, “Enhanced semiconductor optical absorption via surface plasmon excitation in metal nanoparticles,” Appl. Phys. Lett. 86, 063106 (2005). [CrossRef]

, 25

N. C. Panoiu and R. M. Osgood, “Enhanced optical absorption for photovoltaics via excitation of waveguide and plasmon-polariton modes,” Opt. Lett. 32, 2825–2827 (2007). [CrossRef] [PubMed]

], as well as subwavelength, all-optical control of the optical power flow in active nanodevices [26

N. C. Panoiu and R. M. Osgood, “Subwavelength Nonlinear Plasmonic Nanowire,” Nano Lett. 4, 2427–2430 (2004). [CrossRef]

, 27

F. Ye, D. Mihalache, B. Hu, and N. C. Panoiu, “Subwavelength Plasmonic Lattice Solitons in Arrays of Metallic Nanowires,” Phys. Rev. Lett. 104, 106802 (2010). [CrossRef] [PubMed]

]. Moreover, because of the evanescent character of localized SPPs, they are ideal tools for achieving subwavelength confinement of the optical field, and as such, SPPs can be instrumental in designing ultra-compact devices, such as nanolasers or laser arrays [28

B. Min, E. Ostby, V. Sorger, E. Ulin-Avila, L. Yang, X. Zhang, and K. Vahala, “High-Q surface-plasmon-polariton whispering-gallery microcavity,” Nature (London) 457, 455–459 (2009). [CrossRef] [PubMed]

, 29

M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature (London) 460, 1110–1112 (2009). [CrossRef] [PubMed]

, 30

M. P. Nezhand, A. Simic, O. Bondarenko, B. Slutsky, A. Mizrahi, L. Feng, V. Lomakin, and Y. Fainman, “Room-temperature subwavelength metallo-dielectric lasers,” Nature Phot. 4, 395–399 (2010). [CrossRef]

, 31

K. Yu, A. Lakhani, and M. C. Wu, “Subwavelength metal-optic semiconductor nanopatch lasers,” Opt. Express 18, 8790–8799 (2010). [CrossRef] [PubMed]

] and optical microcavities [32

V. Berger, “Second-harmonic generation in monolithic cavities”, J. Opt. Soc. Am. B 14, 1351–1360 (1997). [CrossRef]

, 33

D. G. Gusev, I. V. Soboleva, M. G. Martemyanov, T. V. Dolgova, A. A. Fedyanin, and O. A. Aktsipetrov, “Enhanced second-harmonic generation in coupled microcavities based on all-silicon photonic crystals,” Phys. Rev. B 68, 233303 (2003). [CrossRef]

, 34

G. Kozyreff, J. L. D. Juarez, and J. Martorell, “Whispering-gallery-mode phase matching for surface second-order nonlinear optical processes in spherical microresonators,” Phys. Rev. A 77, 043817 (2008). [CrossRef]

, 35

Y. Xu, M. Han, A. B. Wang, Z. Liu, and J. R. Heflin, “Second order parametric processes in nonlinear silica microspheres,” Phys. Rev. Lett. 100, 163905 (2008). [CrossRef] [PubMed]

]. In this connection, a central issue is to devise plasmonic structures which support localized SPPs with low optical losses (modes with large Q factor). A promising approach to address this challenge, which we explore in this work, is to reduce the radiative losses by employing cavity-shaped nonlinear plasmonic structures rather than plasmonic nanoparticles.

Fig. 1. Schematic of the system geometry.

In this paper we investigate the characteristics of linear and nonlinear localized SPPs excited upon the interaction of ultra-short optical pulses with plasmonic structures made of twodimensional (2D) distributions of metallic nanowires. In our study, we employ the multiple scattering matrix (MSM) algorithm for calculating the field dynamics and the spectral characteristics of the optical field, at both the fundamental frequency (FF) and the SH. Unlike our previous study [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

], in which we have investigated the spectral optical response of plasmonic nanostructures, in this paper we analyze via time-domain techniques the optical properties of plasmonic cavities. Our study reveals the existence at the SH of two markedly different types of plasmon resonances, namely, geometry independent multipole plasmon modes, which correspond to the excitation of weakly interacting modes of single cylinders, and geometry dependent plasmonic cavity modes. In addition, we show that the main parameter characterizing the cavity modes, namely, the Q factor, strongly depends on the structure and shape of the cavity. Our results suggest that by carefully designing the system geometry, specifically, the separation distance between the scatterers, the radiative losses can be greatly reduced, thus being possible to design plasmonic cavities with extremely large Q factor. The paper is organized as follows. Section 2 presents the plasmonic structures we investigate and provides a brief overview of the numerical method employed in our analysis. Section 3 contains a detailed discussion of the results obtained for a series of plasmonic cavities with different geometries. The final section summarizes our findings and outlines the main conclusions of our work.

2. Description of the numerical method

We begin our brief outline of the numerical approach used in our study with a description of the geometry of the system we investigate. It consists of an ensemble of N parallel, infinitely long cylinders, C_j , j = 1, 2, …, N, embedded in a background medium with electric permittivity ε_b and magnetic permeability µ_b (in our calculations we assume µ_b = µ ₀). The location and physical characteristics of the cylinders are defined by their polar coordinates (r_j,φ_j ), radius R_j , electric permittivity ε_j (w), magnetic permeability µ_j , and surface second-order susceptibility, $\hat{χ}$ _s,j ⁽²⁾, their axis being perpendicular to the xy-plane (see Fig. 1). Furthermore, the position of the center of the k-th cylinder, specified in a coordinate system with the origin in O_j , is defined by the polar coordinates (r^j _k ,φ^j _k ), whereas the position of an arbitrary point P, defined with respect to the coordinate systems with the origin in O and O_j , is specified by the polar coordinates (r_P ,φ_P ) and (r^j _P,φ^j _P ), respectively. In what follows, we assume that the wave vector k ₀ of the incoming pulse lies in the xy-plane, the angle between k ₀ and the x-axis being ϕ ₀. Moreover, we consider only the case of TE polarized light since in the case of the TM polarization no SH is generated. Under these conditions, the only nonzero field components are H_z , E_x , and E_y .

In order to determine the temporal evolution of the field generated upon the interaction between an incoming optical pulse and an assembly of metallic nanowires we use a three-step approach [37

J. A. Sanchez-Gil and A. A. Maradudin, “Dynamic near-field calculations of surface-plasmon polariton pulses resonantly scattered at sub-micron metal defects,” Opt. Express 12, 883–894 (2004). [CrossRef] [PubMed]

]. First, the incoming pulse is Fourier transformed in a series of monochromatic plane waves, which have the same wave vector but different amplitude and frequency. Then, following the formalism introduced in Ref. [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

], the linear and nonlinear scattering problems are solved for each monochromatic component of the input pulse. Finally, the optical field in the time domain is reconstructed by inverse Fourier transforming the results obtained at the previous step. Thus, following this algorithm, the magnetic field of the input pulse is given by

H_{z}^{inc} (r, t) = e^{i (k_{0} \cdot r - ω_{0} t)} \frac{1}{2 π} \int_{- \infty}^{\infty} {\tilde{H}}_{z}^{inc} (\bar{ω}) e^{- i \bar{ω} t} d \bar{ω},

(1)

where ω ₀ is the pulse carrier frequency and H̃^inc _z is the Fourier transform of the pulse envelope.In our calculations we have considered input pulses with Gaussian envelopes,

H_{z}^{inc} (t) = H_{0} e^{- \frac{t^{2}}{{2 T}_{0}^{2}}} e^{- i ω_{0} t},

(2)

where H ₀ is the pulse amplitude and the pulse duration T ₀ is related to the pulse full-width at half-maximum (FWHM) by

T_{FWHM} = 2 \sqrt{\ln 2} T_{0}

. With these definitions, the amplitude of the plane waves contained in the input pulse can be written as

{\tilde{H}}_{z}^{inc} (ω - ω_{0}) = H_{0} \sqrt{π} T_{0} e^{- \frac{T_{0}^{2} {(ω - ω_{0})}^{2}}{4}},

(3)

where the physical angular frequency ω = $\bar{ω}$ + ω ₀. For the sake of simplicity, we have assumed that the center of the pulse would reach the origin of the coordinate system at t = 0.

The second step of our method consists in solving the scattering problem for a monochromatic plane wave, at both the FF and SH. To this end, we used a numerical method based on the MSM formalism, which has been successfully employed in the study of linear scattering from arrays of cylinders [38

E. Centeno and D. Felbacq, “Rigorous vector diffraction of electromagnetic waves by bidimensional photonic crystals,” J. Opt. Soc. Am. A 17, 320–327 (2000). [CrossRef]

] and nonlinear scattering from a single cylinder [39

J. I. Dadap, “Optical second-harmonic scattering from cylindrical particles,” Phys. Rev. B 78, 205322 (2008). [CrossRef]

, 40

C. I. Valencia, E. R. Mendez, and B. S. Mendoza, “Second-harmonic generation in the scattering of light by an infinite cylinder,” J. Opt. Soc. Am. B 21, 36–44 (2004). [CrossRef]

], a 2D object of arbitrary shape [41

C. I. Valencia, E. R. Mendez, and B. S. Mendoza, “Second-harmonic generation in the scattering of light by two-dimensional particles,” J. Opt. Soc. Am. B 20, 2150–2161 (2003). [CrossRef]

], and arrays of cylinders embedded in a nonlinear medium [42

E. Centeno and D. Felbacq, “Second-harmonic emission in two-dimensional photonic crystals,” J. Opt. Soc. Am. B 23, 2257–2264 (2006). [CrossRef]

]. Thus, we use a modified version of the nonlinear MSM algorithm [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

], which allows one to include in the theoretical formalism the effects induced by both the surface and bulk nonlinear polarizations. This feature is particularly important in the case of centrosymmetric media, such as metals, since in this case the local, bulk nonlinear polarization vanishes and therefore surface nonlinear effects become dominant [43

J. Rudnick and E. A. Stern, “Second-harmonic Radiation from Metal Surfaces,” Phys. Rev. B 4, 4274 (1971). [CrossRef]

, 44

J. I. Dadap, J. Shan, K. B. Eisenthal, and T. F. Heinz, “Second-Harmonic Rayleigh Scattering from a Sphere of Centrosymmetric Material,” Phys. Rev. Lett. 83, 4045 (1999). [CrossRef]

]. Thus, the numerical algorithm consists of two steps. First, the scattered field at the FF is calculated and it is then used to compute the total nonlinear polarization at the SH. This source polarization is then used to calculate the SH field.

Without being a constraint imposed by the MSM formalism, we assume in our study that all cylinders are made of the same material, namely, Ag. Moreover, we assume that the permittivity of the metal is described by the Drude model,

ε (ω) = ε_{0} [1 - \frac{ω_{p}^{2}}{ω (ω + i v)}]

, where ω_p = 1.35× 10¹⁶ rad/s is the plasma frequency and ν = 2.73×10¹³ s⁻¹ is the damping frequency [45

M. A. Ordal, R. J. Bell, R. W. Alexander, Jr, L. L. Long, and M. R. Querry “Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W.,” Appl. Opt. 24, 4493–4499 (1985). [CrossRef] [PubMed]

In what follows we present a brief description of the MSM formalism, a comprehensive description of this numerical method being provided in Ref. [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

]. Thus, the magnetic field of an incoming plane wave can be expanded in a Fourier-Bessel series as,

H_{z}^{inc} (r, φ) = Σ_{m = - \infty}^{\infty} a_{m} J_{m} (k_{0} r) e^{im φ},

(4)

where J_m denotes the Bessel functions of the first kind and

a_{m} = H_{ω} e^{- im (\frac{π}{2} + ϕ_{0})}

are the Fourier-Bessel coefficients, with H_ω the field amplitude at the frequency ω. Similarly, the field scattered from the cylinder C_j can be expressed as:

H_{z, j}^{sc} (r, φ) = Σ_{m = - \infty}^{\infty} b_{mj} H_{m}^{(2)} (k_{0} r_{P}^{j}) e^{im φ_{P}^{j}},

(5)

where H_m ⁽²⁾ is the Hankel function of the second kind and b_mj are the expansion coefficients of the scattered field; they represent the main unknowns of the linear scattering problem. Once these coefficients are calculated, the total field at a point P can be determined by using:

H_{z}^{tot} (P) = Σ_{m = - \infty}^{\infty} a_{m} J_{m} (k_{0} r_{P}) e^{im φ_{P}} + Σ_{j = 1}^{N} Σ_{m = - \infty}^{\infty} b_{mj} H_{m}^{(2)} (k_{0} r_{P}^{j}) e^{im φ_{P}^{j}} .

(6)

By introducing the column vectors

a_{j} = {e^{- {ik}_{0} r_{j} \cos (φ_{j} - ϕ_{0})} a_{m}}

and b _j = {b_mj } we can relate the two sets of expansion coefficients, for the incoming and scattered fields, via the scattering matrix, S _j, of the cylinder C_j (for more details, see Ref. [38

E. Centeno and D. Felbacq, “Rigorous vector diffraction of electromagnetic waves by bidimensional photonic crystals,” J. Opt. Soc. Am. A 17, 320–327 (2000). [CrossRef]

]):

Σ_{k = 1}^{N} [δ_{kj} I - (1 - δ_{kj}) S_{j} T_{jk}] b_{k} = S_{j} a_{j}, j = 1, 2, . . ., N,

(7)

where I is the identity matrix and T _jk is the coupling matrix between the cylinders C_j and C_k . The matrices S _j can be easily found by solving the scattering problem for an isolated cylinder, with the appropriate boundary conditions. Then, by introducing the single column vectors B = {b _j} and G = {S _j a _j}, we can reduce the Eq. (7) to a single matrix equation:

S_{ω} B = G,

(8)

where S _ω is the scattering matrix of the whole system; from the Eq. (7) it is found to be:

S_{ω} = (\begin{matrix} I & - S_{1} T_{12} & - S_{1} T_{13} & \dots \\ - S_{2} T_{21} & I & - S_{2} T_{23} & \dots \\ - S_{3} T_{31} & - S_{3} T_{32} & I & \dots \\ ⋮ & ⋮ & ⋮ & ⋱ \end{matrix}) .

(9)

Since G and T _jk are known if the system characteristics, such as the material and geometrical parameters, are specified, the linear scattering problem has been reduced to the calculation of S _ω and finding the solution of the Eq. (8). These calculations yield the coefficients b_mj , which in turn can be used to determine the total field at the FF via the Eq. (6).

The second step in the nonlinear MSM algorithm consists of finding the field at the SH. In order to do this, one first has to determine the nonlinear polarization at the SH, which represents the source of the SH field. In the case of centrosymmetric media, such as noble metals considered in our study, the nonlinear polarization is the sum between a dipole-allowed surface contribution, P _s(r;2ω), which is strongly localized within a thin interfacial layer at the surface of the metal, and a nonlocal, bulk nonlinear polarization, P _b(r;2ω), generated inside the metal by electric quadrupoles and magnetic dipoles [46

T. F. Heinz, “Second-Order Nonlinear Optical Effects at Surfaces and Interfaces,” in Nonlinear Surface Electromagnetic Phenomena, edited by H. E. Ponath and G. I. Stegeman (Elsevier,Amsterdam, 1991); pp. 353–416.

]. The surface polarization is defined as

P_{s} (r; 2 ω) = ε_{0} {\hat{χ}}_{s}^{(2)} (r; 2 ω) : E (r; ω) E (r; ω) δ (r - r_{s}),

(10)

where r _s defines the surface of the metal, χˆ_s ⁽²⁾ is the surface second-order susceptibility tensor, and the Dirac function models the surface characteristic of this polarization. Due to the symmetries of a homogeneous and isotropic interface the surface nonlinear susceptibility tensor has only three independent components. In the case of Ag their values are χˆ ⁽²⁾ _s,⊥⊥⊥ = 2.79×10⁻¹⁸ m²/V, χˆ ⁽²⁾ _s,∥∥⊥ = χˆ ⁽²⁾ _s,∥⊥∥ = 3.98×10⁻²⁰ m²/V, and χˆ ⁽²⁾ _s,⊥∥∥ = 0 [47

D. Krause, C. W. Teplin, and C. T. Rogers, “Optical surface second harmonic measurements of isotropic thin-film metals: Gold, silver, copper, aluminum, and tantalum,” J. App. Phys. 96, 3626–3635 (2004). [CrossRef]

For an isotropic centrosymmetric medium the bulk nonlinear polarization is given by [48

N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, “Optical Second-Harmonic Generation in Reflection from Media with Inversion Symmetry,” Phys. Rev. 174, 813–822 (1968). [CrossRef]

P_{b} (2 ω) = α [E (ω) \cdot \nabla] E (ω) + β E (ω) [\nabla \cdot E (ω)] + γ \nabla [E (ω) \cdot E (ω)],

(11)

where α, β, and γ are material dependent physical constants. Assuming that the electrons in the metal are described by the free-electron gas model, these bulk nonlinear coefficients are

α = 0; β = ε_{0} \frac{e}{{2 m}_{0} ω^{2}}; γ = \frac{β}{4} [1 - ε_{r} (ω)],

(12)

with e and m ₀ being the electron charge and mass, respectively, and ε_r (ω) = ε(ω)/ε ₀.

Using an approach similar to that used to determine the field at the FF, the total field at the SH, at a point P outside the cylinder C_j , can be written as [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

H_{z, j}^{tot} (P, Ω) = Σ_{m = - \infty}^{\infty} (a_{Ω, mj} + b_{Ω, mj}) H_{m}^{(2)} ({kr}_{P}^{j}) e^{im φ_{P}^{j}}

+ Σ_{\underset{k \neq j}{k = 1}}^{N} Σ_{m, q = - \infty}^{\infty} T_{jk, mq} (a_{Ω, qk} + b_{Ω, qk}) J_{m} ({kr}_{P}^{j}) e^{im φ_{P}^{j}} .

(13)

where

k = \sqrt{μ_{0} ε_{b}} Ω

is the wave vector at the SH. In the Eq. (13) the index j signifies that H^tot _z,j is calculated in the coordinate system with the origin in O_j , whereas the coefficients a _Ω,mj and b _Ω,mj are the Fourier-Bessel coefficients of the series expansion of the field generated by the nonlinear polarization and the scattered field, respectively. Note that the coefficients a _Ω,j = {a _Ω,mj} are fully determined by the nonlinear polarization and therefore the main unknowns of the nonlinear scattering problem are the coefficients b _Ω,j = {b _Ω,mj}. The two sets of expansion coefficients satisfy the following system of matrix equations:

Σ_{\underset{k \neq j}{k = 1}}^{N} S_{Ω, j} T_{jk} (a_{Ω, k} + b_{Ω, k}) = (a_{Ω, j} + b_{Ω, j}) - g_{Ω, j}^{self}, j = 1, 2, . . ., N,

(14)

where g ^self _Ω,j are the source coefficients of the nonlinear field generated by a single cylinder; their rather cumbersome expressions can be found in Ref. [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

]. By denoting the vectors A _Ω ={a _Ω,j},B _Ω = {b _Ω,j}, and G ^self _Ω = {g ^self _Ω,j} we can reduce the system (14) to a single matrix equation:

S_{Ω} B_{Ω} = - S_{Ω} A_{Ω} + G_{Ω}^{self} .

(15)

Finally, by solving this linear system one finds the coefficients b _Ω,j and, subsequently, via the Eq. (13), the total field at the SH at any arbitrary position.

The final step of our algorithm consists in calculating the optical field in the time domain. To achieve this, one simply inverse Fourier transforms the field spectra obtained in the preceding step. Specifically, accurate description of the optical field dynamics has been achieved by dividing the total temporal domain, T = 70√2T ₀, in N_FT = 8192 equal time intervals. Typical execution time when employing 64 Intel® Xeon cores is about 1 hour, provided that the field variables are computed at a single spatial point. This procedure allows one to obtain a complete description of the dynamics of the optical near- and far-field, at both the FF and the SH. As a final remark on the numerical algorithm used here, we want to stress that it allows one to easily incorporate nonlinear surface effects, which would prove to be a challenging task for alternative methods based on the finite-difference time-domain discretization of the field variables.

3. Results and discussions

In this section we will illustrate how the numerical method just described can be applied to investigate the physical characteristics of linear and nonlinear localized SPP modes excited by sub-picosecond pulses upon their interaction with plasmonic nanocavities. In particular, we explore the relation between the geometry and material parameters of the plasmonic cavities and the main optical properties of the localized SPP modes. In the presentation of our main results we will focus on the optical properties of the localized modes excited at the SH as the main conclusions derived in this case also apply to the modes observed at the FF. In addition, since there is no incoming pulse at the SH, some physical quantities characterizing the localized SPP modes, such as their Q factor, can be much easier calculated if the optical field at the SH is analyzed. This approach can also be relevant for a series of potential technological applications, such as sensing or optical detection, as the optical signal generated at the SH is spectrally well separated from the incoming and scattered fields at the FF.

To begin with, we illustrate in Fig. 2 the generic characteristics of the linear and nonlinear interaction between an incoming optical pulse and a SPP cavity mode. One convenient approach for identifying these modes relies on the spectra of the absorption cross-section at the FF, σ_a (ω), and the corresponding absorption spectra at the SH, ∑_a(2ω). To be more specific, it is expected that at the resonance frequency of the localized SPP modes the spectra of the optical absorption presents resonances, as at these specific frequencies the interaction of the optical near-field with the metallic nanowires is enhanced. At the FF the absorption spectra are calculated by using the scattering coefficients b _j, in conjunction with the optical theorem, whereas at the SH we simply compute the total Joule dissipation power by integrating over the metal regions the density of the dissipated power,

P_{abs} = \frac{1}{2} ℜ 𝔢 (J \cdot E^{*})

, with J = σ_j E being the conduction current density and σ_j the conductivity of the cylinder C_j (see Ref. [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

] for more details pertaining to these calculations). In the case of the Drude model, σ = ε ₀ ω ² _p/(v −iω).

The top panels in Fig. 2, which corresponds to a plasmonic cavity made of 6 metallic nanowires, show that the spectra of the absorption cross sections at the FF and SH contain a series of spectral peaks. In addition, the spectral location of these peaks changes with the separation distance between adjacent nanowires, d, a variation of d of less than 100 nm leading to a spectral shift of the resonances at the FF of almost 300 nm. Additional information pertaining to the optical properties of plasmonic cavity modes is revealed by the temporal dynamics of the optical near-field. Thus, the bottom panels in Fig. 2 (and the accompanying movies) reveal that the field at the SH remains trapped in the cavity long after the initial pulsed excitation at the FF has passed through the cavity. Since the nonlinear response of the metal is assumed to be instantaneous, it can be inferred that a nonlinear cavity mode with a significant lifetime and, implicitly, large Q factor, is formed in the cavity. In addition, the plasmonic character of this mode is evident from the spatial distribution of the near-field, namely, the field has large values at the metal surfaces and decays steeply towards the center of the plasmonic cavity. In what follows, we present a more detailed analysis of these plasmonic cavity modes, as well as modes of a different physical nature, namely, multipole plasmon modes.

Fig. 2. Top panels show logarithmic plots of the absorption cross sections, calculated for a plasmonic cavity containing 6 nanowires. The legend indicates the separation distance, in nanometers. Bottom panels show a snapshot of the temporal evolution of the intensity of the electric field. The left and right panels correspond to the FF (Media 1) and SH (Media 2), respectively. The plasmonic cavity consists of Ag cylinders with R = 200 nm and d = 60 nm. The wavelength at the FF is λ_FF = 858 nm and the angle of incidence is ϕ ₀ = 90°.

3.1. Localized plasmon modes in coupled cylindrical nanowires

We begin by considering a series of different cavity geometries, namely, cylinder distributions containing 2, 3, and 4 cylinders. Figure 3 summarizes the main results pertaining to these three geometries. As in the previous case, the spectra of the absorption cross section at the SH show a series of sharp peaks, which suggests the existence of SPP modes. Importantly, the resonance frequency of these modes does not depend on the number of cylinders in the distribution of scatterers, which means that these modes are formed primarily due to the excitation of optical modes in each of the metallic nanowires. This conclusion is supported by the field distributions presented in the bottom panels of the Fig. 3. Thus, these plots clearly show that the modes at λ_SH = 578 nm correspond to dipole (cut-off) modes of the nanowires. Similar field distributions, shown in Fig. 4, demonstrate that the resonances at smaller wavelength (λ_SH = 336 nm) correspond to quadrupole modes of the nanowires. Figure 4 also shows that, as expected, multipole resonances are not only excited in plasmonic structures containing a small number of scatterers but that in fact their existance is a generic phenomenon. However, when the number of scatterers increases the amplitude of the optical modes excited in each nanowire varies with its location in the 2D nanowire assembly, especially when the wavelength becomes comparable to the size of the plasmonic structure. It should be noted that similar resonant modes are exited at the FF, too, but they do not appear in the absorption spectra since they are “buried” in the background generated due to the absorption of the input pulse. These modes, however, can be identified as resonance peaks in the spectra of the scattering cross section [36

C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]

]. Moreover, note that the SH field in the region in-between adjacent cylinders is small, which means that for this separation distance the “hybridization” effects due to the interaction between the optical modes excited in adjacent cylinders are weak.

Fig. 3. Top panels show logarithmic plots of the absorption cross sections. The legends indicate the number of cylinders. Bottom panels present the amplitude of the electric field at the SH, for Ag cylinders with R = 200 nm and d = 20 nm, at λ_SH = 578 nm and ϕ ₀ = 0 (for better visualization, we plot the fourth-order square root of the field amplitude).

To characterize the influence of the system geometry on the resonance frequencies of the multipole plasmon modes we have determined the absorption spectra for different values of the angle of incidence ϕ ₀ and separation distance d. The dispersion plots corresponding to the absorption spectra at the SH are presented in the Fig. 5. These absorption spectra clearly indicate that the resonance frequencies of the multipole plasmon modes are almost independent on the system parameters, supporting therefore the conclusion that they correspond to optical modes excited in each of the metallic nanowires. More specifically, the resonance frequencies of these modes remain unchanged even if the distance between the adjacent metallic nanowires is decreased to a value as small as 1 nm. On the other hand, a significant increase of the optical absorption occurs when the separation distance becomes smaller than a few nanometers, an effect explained by the field enhancement observed in the region separating adjacent nanowires. Figure 5 also shows that the spectral width of the resonances varies with the angle of incidence ϕ ₀. This dependence suggests that, as expected, the strength of the interaction between the input wave and the plasmonic structure and, consequently, the magnitude of the field generated at the SH, changes with the angle of incidence. In addition, it can be seen from Fig. 5 that as the separation distance between the nanowires increases the width of the spectral resonances decreases, and effect that is explained by the fact that the strength of the coupling between the modes excited in adjacent nanowires decreases with the separation distance.

Fig. 4. Distribution of the amplitude of the electric field (top panels) and the real part of the magnetic field (bottom panels) at the SH, calculated for three different plasmonic cavities made of Ag cylinders with R = 200 nm and separation distance d = 20 nm. The wavelength at the SH is λ_SH = 336 nm and ϕ ₀ = 0 (for better visualization, in the case of the electric field, we plot the fourth-order square root of the field amplitude).

Fig. 5. Logarithmic plots of absorption cross section spectra at the SH vs. the angle of incidence ϕ ₀ (top) and separation distance d (bottom), determined for a three-cylinder geometry. The top and bottom panels correspond to d = 60 nm and ϕ ₀ = 0, respectively.

3.2. Plasmonic cavity modes

While the analysis of localized multipole plasmon modes can provide a valuable insight into the contribution of each individual scatterer to the overall optical response of the plasmonic structure, it does not reveal the complete picture of the interaction between optical pulses and plasmonic cavities. To be more specific, our analysis reveals that plasmonic structures containing a larger number of metallic nanowires support additional plasmonic modes, which have a different physical origin as compared to that of the multipole plasmon modes. In order to illustrate this conclusion, we present in Fig. 6 the absorption spectra, at both the FF and the SH, of cavity-shaped plasmonic structures containing 4, 6, and 8 metallic nanowires. As in the previous cases these spectra present a series of spectral peaks, which correspond to resonances of the plasmonic system. By inspecting the field profiles corresponding to these spectral peaks we found that besides the multipole plasmon modes similar to those supported by plasmonic structures with a smaller number of metallic nanowires there are additional, markedly different type of modes, which we call plasmonic cavity modes. These localized plasmon modes are formed due to the coherent response of the whole cavity. For example, as expected, in the case of the four-cylinder cavity the absorption spectrum at the SH presents two resonance peaks at λ_SH = 578 nm and λ_SH = 336 nm, spectral peaks that correspond to the dipole and quadrupole plasmon modes, respectively. The absorption spectra at the SH have, however, additional resonance peaks at λ_SH = 321 nm, λ_SH = 429 nm, and λ_SH = 333 nm. Since the resonance wavelength changes significantly with the number of nanowires forming the plasmonic cavity, it can be concluded that these optical modes are determined by the coherent response of the whole structure. The field profiles presented in Fig. 6 further support this conclusion, by showing that at the corresponding resonance wavelengths the optical field is not confined only to the region surrounding each nanowire but spreads inside the plasmonic cavity.

Fig. 6. Top panels show logarithmic plots of the absorption cross section at the FF and the SH. The legend indicates the number of cylinders forming the plasmonic cavity. Bottom panels present the distribution of the amplitude of the electric field at the SH, for Ag cylinders with R=200 nm and separation distance d=60 nm. From left to right, the wavelength at the SH is λ_SH = 321 nm, λ_SH = 429 nm, and λ_SH = 333 nm.

We want to point out that the plasmonic cavity modes investigated here are similar to whispering-gallery modes recently observed in plasmonic structures with a different geometry [28

], the main difference being that in our case the angular momentum of the plasmonic cavity modes is equal to zero. Indeed, the angular momentum of the incident beam is zero and therefore the angular momentum of the excited modes must be zero, too. One additional important feature of the plasmonic cavity modes presented in Fig. 6 is that they do not appear as resonances in the scattering cross section spectra. Therefore, they are dark-plasmon modes that do not couple with the radiation continuum [49

M. Liu, T. W. Lee, S. K. Gray, P. Guyot-Sionnest, and M. Pelton, “Excitation of Dark Plasmons in Metal Nanoparticles by a Localized Emitter,” Phys. Rev. Lett. 102, 107401 (2009). [CrossRef] [PubMed]

], and as a result the corresponding radiative losses are suppressed. Nevertheless, in our case these modes are excited via the nonlinear polarization generated at the SH, which acts as localized dipole sources. These results suggest that the nonlinear polarization at the SH can be used to excite subradiant (low loss) propagating modes formed in chains of metallic nanoparticles, via the near-field resonant coupling of single-particle dark-plasmon modes. As we will prove later, the suppression of the radiative losses also leads to a considerable increase of the lifetime of the plasmonic cavity modes and, consequently, makes it possible to design plasmonic cavities with very large Q factor.

As it can be seen in Fig. 3, for a separation distance of d = 20 nm the four-cylinder cavity does not have a plasmonic cavity mode but such a mode exists for d = 60 nm. This observation provides further evidence that the characteristics of the plasmonic cavity modes are strongly influenced by the geometry of the cavity. In order to explore this dependence in more detail, we will focuss in what follows on the optical properties of the plasmonic cavity formed by placing metallic nanowires at the corners of a hexagon. In making this choice we were primarily guided by the fact that the cavity modes of this structure are very well defined and, as we will see in the next sub-section, they have a large Q factor.

Fig. 7. The same as in Fig. 5 but for a hexagonal plasmonic cavity.

Because the dispersion spectra of the absorption cross section represent a powerful tool for investigating the properties of localized SPP modes, we have calculated these spectra for the hexagonal plasmonic cavity. The results, plotted in Fig. 7, provide further insight into the specific properties of plasmonic cavity modes. As expected, because these modes do not depend on the optical coupling between the incoming wave and the plasmonic cavity, they are independent on the angle of incidence ϕ ₀. On the other hand, the separation distance between adjacent cylinders does have a notable effect on the spectral location of the resonance wavelength of the plasmonic cavity modes. Thus, the resonance wavelength of the cavity mode increases almost linearly with the separation distance, a wavelength shift of almost 150 nm being observed when the separation distance changes by about 100 nm. In addition, it can be seen that, as in the case of multipole plasmon modes, the spectral width of the resonance decreases as d increases; however, as we will demonstrate later, in the case of plasmonic cavity modes this behavior is determined by the interplay between the radiative and absorption losses. As we will show in the next sub-section, other parameters characterizing the plasmonic cavity modes, such as the Q factor, have a more intricate dependence on the geometry of the cavity. A general feature, however, of these modes is that their optical properties can be easily tailored by modifying the shape of the cavity. Moreover, similar to the case of plasmonic cavities containing a smaller number of metallic nanowires, high optical absorption is observed for a separation distance approaching 1 nm. Again, this effect is attributable to the strong electromagnetic field generated at the surface of the nanowires for small inter-cylinder separation distance.

The strong dependence of the resonance frequency of plasmonic cavity modes on the distance between cylinders or other geometrical parameters can have important applications to sensing or photovoltaic devices. Specifically, the plasmonic cavity can be viewed as playing the role of an optical antenna that collects and concentrates the signal carried by the input pulse into a reduced volume, making it possible to increase the signal-to-noise ratio and/or speed of a detector. These plasmonic cavities can also be employed in the design of lasers with subwavelength size, as has in fact been recently demonstrated [30

, 31

K. Yu, A. Lakhani, and M. C. Wu, “Subwavelength metal-optic semiconductor nanopatch lasers,” Opt. Express 18, 8790–8799 (2010). [CrossRef] [PubMed]

]. In particular, the Q factor of plasmonic cavities employed in laser applications plays a crucial role in determining the performance of such nanolasers. Consequently, in what follows we study in more detail the dependence of the Q factor of the plasmonic cavity modes on the parameters defining the plasmonic structure.

3.3. Time domain analysis

The main physical quantity that describes the temporal response of an optical mode is the Q factor (or, equivalently, its lifetime). In order to calculate this important parameter that characterizes a plasmonic cavity mode we proceeded as follows. First, we illuminated the cavity with an optical pulse of sub-picosecond duration and subsequently we recorded the optical field, at both the FF and SH, at an arbitrary location inside the cavity. If the carrier frequency of the input optical pulse is close to a resonance frequency of a plasmonic cavity mode the asymptotic temporal evolution of the optical field can be represented by an exponential dependence,

E (t) = E_{0} e^{- \frac{ω_{r} t}{Q}},

(16)

where ω_r is the resonance frequency of the plasmonic cavity mode. It should be noted that the relation (16) is independent of the location of the point in the cavity where the field is measured, the value of the carrier frequency of the input optical pulse (as long as it is close to ω_r ), and the duration of the incident optical pulse, a conclusion that is fully verified by our numerical simulations. Therefore, the Q factor of the optical mode can be easily determined by calculating the slope of the line representing the linear fit of the semi-logarithmic temporal dependence of the computed field inside the cavity.

We have calculated the Q factor of the plasmonic cavity mode of the the hexagonal plasmonic cavity and the main results are summarized in Fig. 8. As expected, when the carrier frequency of the input optical pulse is close to the resonance frequency of the cavity mode the temporal evolution of the normalized field inside the cavity follows the exponential decay described by the Eq. (16). From this asymptotic dependence we have derived the value of the Q factor of the cavity mode and the corresponding lifetime, τ = Q/ω_r .

Our calculations show that the lifetime of the plasmonic cavity mode increases from τ = 380.23 fs at d = 20 nm to τ = 582.2 fs at d = 100 nm. If the separation distance is further increased, the lifetime begins to decrease. Moreover, the Q factor of the optical mode follows the same dependence on the separation distance, its maximum value, Q_max = 2294, being reached for d = 113 nm. Note that this extremely large value is more than twice as large as the Q factor of recently observed plasmonic whispering-gallery modes [28

], although it should be mentioned that these modes were observed in three-dimensional plasmonic cavities. This dependence of the Q factor on the separation distance between the metallic nanowires is somewhat surprising because one would expect that the radiative losses increase with the separation distance and therefore the Q factor should monotonously decrease as d increases. In order to explain this apparent contradiction, we decompose the Q factor as

Fig. 8. a) Normalized electric field inside the cavity vs. time. The numbers in the legend refer to the separation distance in nanometers. b) The dependence of the Q factor on the separation distance d: circles represent simulation results whereas the dotted line is provided as a guide to the eye. The cavity consists of Ag cylinders with R = 200 nm and the incoming pulse has T ₀ = 283 fs.

\frac{1}{Q} = \frac{1}{Q_{abs}} + \frac{1}{Q_{rad}},

(17)

where Q_abs is determined by the absorption in the metallic nanowires and Q_rad is due to radiative losses. The Q_abs factor can be calculated from the absorption spectra by fitting with a Lorentzian the region the spectrum corresponding to the resonance peak. Thus, Q_abs = ω_r /∆ω, where ∆ω is the spectral width of the Lorentzian. By using this procedure we found that for small values of the separation distance Q_abs ≈ Q, which means that the losses of the plasmonic cavity mode are primarily due to the absorption loss in the metal. This result also supports the conclusion that the plasmonic cavity mode is a dark-plasmon mode, whose radiative losses are suppressed. As d increases, the field confinement decreases, and therefore the absorption losses decrease. Consequently, the Q factor of the optical mode increases. If the separation distance is further increased, the optical field begins to leak more easily out of the cavity, the cavity effects become weaker, and consequently the radiative losses start to dominate. As a result, the Q factor of the mode begins to decrease. This scenario predicts that there is a separation distance for which the Q factor reaches a maximum value, a prediction which is fully verified by the results presented in Fig. 8b. This analysis also suggests that by minimizing the optical losses associated to plasmonic cavity modes it is possible to optimize considerably the efficiency of these plasmonic cavities, a property that can have important implications to the development of efficient subwavelength nanolasers.

4. Conclusions

In conclusion, we have presented and discussed the main optical properties of linear and nonlinear localized SPP modes excited upon the interaction between ultra-short optical pulses and nanocavities made of metallic nanowires. Our numerical analysis, based on the MSM formalism, has revealed that plasmonic cavities support two distinct types of localized SPP modes, namely, multipole plasmon modes that are the result of the hybridization of coupled plasmon modes supported by each metallic nanowire of the plasmonic cavity, and plasmonic cavity modes, which can be viewed as the coherent optical response of the entire assembly of metallic nanowires. We have also demonstrated that this dichotomy in the physical origin of these optical modes is responsible for their markedly different optical properties. For example, whereas the properties of the multipole plasmon modes depend almost exclusively on the size of the individual nanowires, the geometrical and material parameters of the plasmonic cavity strongly influence the characteristics of the plasmonic cavity modes. In particular, we have demonstrated that this feature can be used effectively to design plasmonic cavities with extremely large Q factors. As a consequence, we believe that the results presented in this work could lead to a better understanding of linear and nonlinear optical effects at the nanoscale and, equally important, will have significant technological applications to nano-sensors, ultra-small detectors, subwavelength lasers, and other active photonic nanodevices.

Acknowledgements

The authors acknowledge the use of the UCL Legion High Performance Computing Facility, and associated support services, in the completion of this work. This work was supported by the UK Engineering and Physical Sciences Research Council (ESPRC), under Grant No. EP/G030502/1.

References and links

1.	S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004). [CrossRef] [PubMed]
2.	R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental Verification of a Negative Index of Refraction,” Science 292, 77–79 (2001). [CrossRef] [PubMed]
3.	S. Zhang, W. Fan, N. C. Panoiu, K. J. Malloy, R. M. Osgood, and S. R. J. Brueck, “Experimental Demonstration of Near-Infrared Negative-Index Metamaterials,” Phys. Rev. Lett. 95, 137404 (2005). [CrossRef] [PubMed]
4.	V. M. Shalaev, W. Cai, U. K. Chettiar, H-K. Yuan, A. K. Sarychev, V. P. Drachev, and A. V. Kildishev, “Negative index of refraction in optical metamaterials,” Opt. Lett. 30, 3356–3358 (2005). [CrossRef]
5.	R. W. Ziolkowski, “Propagation in and scattering from a matched metamaterial having a zero index of refraction,” Phys. Rev. E 70, 046608 (2004). [CrossRef]
6.	N. C. Panoiu, R. M. Osgood, S. Zhang, and S. R. J. Brueck, “Zero-n bandgap in photonic crystal superlattices,” J. Opt. Soc. Am. B 23, 506–513 (2006). [CrossRef]
7.	S. Kocaman, R. Chatterjee, N. C. Panoiu, J. F. McMillan, M. B. Yu, R. M. Osgood, D. L. Kwong, and C. W. Wong, “Observation of Zeroth-Order Band Gaps in Negative-Refraction Photonic Crystal Superlattices at Near-Infrared Frequencies,” Phys. Rev. Lett. 102, 203905 (2009). [CrossRef] [PubMed]
8.	U. Kreibig and M. Vollmer, Optical Properties of Metal Clusters (Springer-Verlag, Berlin, 1995).
9.	A. V. Zayats, I. I. Smolyaninov, and A. A. Maradudin, “Nano-optics of surface plasmon polaritons,” Phys. Rep. 408, 131–314 (2005). [CrossRef]
10.	W. Fan, S. Zhang, N. C. Panoiu, A. Abdenour, S. Krishna, R. M. Osgood, K. J. Malloy, and S. R. J. Brueck, “Second Harmonic Generation from a Nanopatterned Isotropic Nonlinear Material,” Nano Lett. 6, 1027–1030 (2006). [CrossRef]
11.	J. A. H. van Nieuwstadt, M. Sandtke, R. H. Harmsen, F. B. Segerink, J. C. Prangsma, S. Enoch, and L. Kuipers, “Strong Modification of the Nonlinear Optical Response of Metallic Subwavelength Hole Arrays,” Phys. Rev. Lett. 97, 146102 (2006). [CrossRef] [PubMed]
12.	L. Cao, N. C. Panoiu, and R. M. Osgood, “Surface second-harmonic generation from surface plasmon waves scattered by metallic nanostructures,” Phys. Rev. B 75, 205401 (2007). [CrossRef]
13.	Y. Zeng, W. Hoyer, J. J. Liu, S. W. Koch, and J. V. Moloney, “Classical theory for second-harmonic generation from metallic nanoparticles,” Phys. Rev. B 79, 235109 (2009). [CrossRef]
14.	I. I. Smolyaninov, A. V. Zayats, A. Gungor, and C. C. Davis, “Single-Photon Tunneling via Localized Surface Plasmons,” Phys. Rev. Lett. 88, 187402 (2002). [CrossRef] [PubMed]
15.	X. W. Wang, G. C. Schatz, and S. K. Gray, “Ultrafast pulse excitation of a metallic nanosystem containing a Kerr nonlinear material,” Phys. Rev. B 74, 195439 (2006). [CrossRef]
16.	K. Kneipp, Y. Wang, H. Kneipp, L. T. Perelman, I. Itzkan, R. R. Dasari, and M. S. Feld, “Single Molecule Detection Using Surface-Enhanced Raman Scattering (SERS),” Phys. Rev. Lett. 78, 1667 (1997). [CrossRef]
17.	T. H. Taminiau, R. J. Moerland, F. B. Segerink, L. Kuipers, and N. F. van Hulst, “ λ/4 Resonance of an Optical Monopole Antenna Probed by Single Molecule Fluorescence,” Nano Lett. 7, 28 (2007). [CrossRef] [PubMed]
18.	B. Knoll and F. Keilmann, “Near-field probing of vibrational absorption for chemical microscopy,” Nature (London) 399, 134–137 (1999). [CrossRef]
19.	R. M. Roth, N. C. Panoiu, M. M. Adams, R. M. Osgood, C. C. Neacsu, and M. B. Raschke, “Resonant-plasmon field enhancement from asymmetrically illuminated conical metallic-probe tips,” Opt. Express 14, 2921–2913 (2006). [CrossRef] [PubMed]
20.	S. I. Bozhevolnyi, J. Beermann, and V. Coello, “Direct Observation of Localized Second-Harmonic Enhancement in Random Metal Nanostructures,” Phys. Rev. Lett. 90, 197403 (2003). [CrossRef] [PubMed]
21.	T. Ishi, J. Fujikata, K. Makita, T. Baba, and K. Ohashi, “Si nano-photodiode with a surface plasmon antenna,” Jpn. J. Appl. Phys. 44, L364–L366 (2005). [CrossRef]
22.	Z. Yu, G. Veronis, S. Fan, and M. L. Brongersma, “Design of midinfrared photodetectors enhanced by surface plasmons on grating structures,” Appl. Phys. Lett. 89, 151116 (2006). [CrossRef]
23.	R. D. R. Bhat, N. C. Panoiu, S. R. J. Brueck, and R. M. Osgood, “Enhancing the signal-to-noise ratio of an infrared photodetector with a circular metal grating,” Opt. Express 16, 4588–4596 (2008). [CrossRef] [PubMed]
24.	D. M. Schaadt, B. Feng, and E. T. Yu, “Enhanced semiconductor optical absorption via surface plasmon excitation in metal nanoparticles,” Appl. Phys. Lett. 86, 063106 (2005). [CrossRef]
25.	N. C. Panoiu and R. M. Osgood, “Enhanced optical absorption for photovoltaics via excitation of waveguide and plasmon-polariton modes,” Opt. Lett. 32, 2825–2827 (2007). [CrossRef] [PubMed]
26.	N. C. Panoiu and R. M. Osgood, “Subwavelength Nonlinear Plasmonic Nanowire,” Nano Lett. 4, 2427–2430 (2004). [CrossRef]
27.	F. Ye, D. Mihalache, B. Hu, and N. C. Panoiu, “Subwavelength Plasmonic Lattice Solitons in Arrays of Metallic Nanowires,” Phys. Rev. Lett. 104, 106802 (2010). [CrossRef] [PubMed]
28.	B. Min, E. Ostby, V. Sorger, E. Ulin-Avila, L. Yang, X. Zhang, and K. Vahala, “High-Q surface-plasmon-polariton whispering-gallery microcavity,” Nature (London) 457, 455–459 (2009). [CrossRef] [PubMed]
29.	M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature (London) 460, 1110–1112 (2009). [CrossRef] [PubMed]
30.	M. P. Nezhand, A. Simic, O. Bondarenko, B. Slutsky, A. Mizrahi, L. Feng, V. Lomakin, and Y. Fainman, “Room-temperature subwavelength metallo-dielectric lasers,” Nature Phot. 4, 395–399 (2010). [CrossRef]
31.	K. Yu, A. Lakhani, and M. C. Wu, “Subwavelength metal-optic semiconductor nanopatch lasers,” Opt. Express 18, 8790–8799 (2010). [CrossRef] [PubMed]
32.	V. Berger, “Second-harmonic generation in monolithic cavities”, J. Opt. Soc. Am. B 14, 1351–1360 (1997). [CrossRef]
33.	D. G. Gusev, I. V. Soboleva, M. G. Martemyanov, T. V. Dolgova, A. A. Fedyanin, and O. A. Aktsipetrov, “Enhanced second-harmonic generation in coupled microcavities based on all-silicon photonic crystals,” Phys. Rev. B 68, 233303 (2003). [CrossRef]
34.	G. Kozyreff, J. L. D. Juarez, and J. Martorell, “Whispering-gallery-mode phase matching for surface second-order nonlinear optical processes in spherical microresonators,” Phys. Rev. A 77, 043817 (2008). [CrossRef]
35.	Y. Xu, M. Han, A. B. Wang, Z. Liu, and J. R. Heflin, “Second order parametric processes in nonlinear silica microspheres,” Phys. Rev. Lett. 100, 163905 (2008). [CrossRef] [PubMed]
36.	C. G. Biris and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]
37.	J. A. Sanchez-Gil and A. A. Maradudin, “Dynamic near-field calculations of surface-plasmon polariton pulses resonantly scattered at sub-micron metal defects,” Opt. Express 12, 883–894 (2004). [CrossRef] [PubMed]
38.	E. Centeno and D. Felbacq, “Rigorous vector diffraction of electromagnetic waves by bidimensional photonic crystals,” J. Opt. Soc. Am. A 17, 320–327 (2000). [CrossRef]
39.	J. I. Dadap, “Optical second-harmonic scattering from cylindrical particles,” Phys. Rev. B 78, 205322 (2008). [CrossRef]
40.	C. I. Valencia, E. R. Mendez, and B. S. Mendoza, “Second-harmonic generation in the scattering of light by an infinite cylinder,” J. Opt. Soc. Am. B 21, 36–44 (2004). [CrossRef]
41.	C. I. Valencia, E. R. Mendez, and B. S. Mendoza, “Second-harmonic generation in the scattering of light by two-dimensional particles,” J. Opt. Soc. Am. B 20, 2150–2161 (2003). [CrossRef]
42.	E. Centeno and D. Felbacq, “Second-harmonic emission in two-dimensional photonic crystals,” J. Opt. Soc. Am. B 23, 2257–2264 (2006). [CrossRef]
43.	J. Rudnick and E. A. Stern, “Second-harmonic Radiation from Metal Surfaces,” Phys. Rev. B 4, 4274 (1971). [CrossRef]
44.	J. I. Dadap, J. Shan, K. B. Eisenthal, and T. F. Heinz, “Second-Harmonic Rayleigh Scattering from a Sphere of Centrosymmetric Material,” Phys. Rev. Lett. 83, 4045 (1999). [CrossRef]
45.	M. A. Ordal, R. J. Bell, R. W. Alexander, Jr, L. L. Long, and M. R. Querry “Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W.,” Appl. Opt. 24, 4493–4499 (1985). [CrossRef] [PubMed]
46.	T. F. Heinz, “Second-Order Nonlinear Optical Effects at Surfaces and Interfaces,” in Nonlinear Surface Electromagnetic Phenomena, edited by H. E. Ponath and G. I. Stegeman (Elsevier,Amsterdam, 1991); pp. 353–416.
47.	D. Krause, C. W. Teplin, and C. T. Rogers, “Optical surface second harmonic measurements of isotropic thin-film metals: Gold, silver, copper, aluminum, and tantalum,” J. App. Phys. 96, 3626–3635 (2004). [CrossRef]
48.	N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, “Optical Second-Harmonic Generation in Reflection from Media with Inversion Symmetry,” Phys. Rev. 174, 813–822 (1968). [CrossRef]
49.	M. Liu, T. W. Lee, S. K. Gray, P. Guyot-Sionnest, and M. Pelton, “Excitation of Dark Plasmons in Metal Nanoparticles by a Localized Emitter,” Phys. Rev. Lett. 102, 107401 (2009). [CrossRef] [PubMed]

OCIS Codes
(190.4350) Nonlinear optics : Nonlinear optics at surfaces
(240.6680) Optics at surfaces : Surface plasmons
(140.3945) Lasers and laser optics : Microcavities

ToC Category:
Optics at Surfaces

History
Original Manuscript: June 21, 2010
Revised Manuscript: July 16, 2010
Manuscript Accepted: July 21, 2010
Published: July 28, 2010

Citation
Claudiu G. Biris and Nicolae C. Panoiu, "Nonlinear pulsed excitation of high-Q optical modes of plasmonic nanocavities," Opt. Express 18, 17165-17179 (2010)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-18-16-17165

Sort: Author | Year | Journal | Reset

References

S. Linden, C. Enkrich, M. Wegener, J. Zhou, T. Koschny, and C. M. Soukoulis, “Magnetic Response of Metamaterials at 100 Terahertz,” Science 306, 1351–1353 (2004). [CrossRef] [PubMed]
R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental Verification of a Negative Index of Refraction,” Science 292, 77–79 (2001). [CrossRef] [PubMed]
S. Zhang, W. Fan, N. C. Panoiu, K. J. Malloy, R. M. Osgood, and S. R. J. Brueck, “Experimental Demonstration of Near-Infrared Negative-Index Metamaterials,” Phys. Rev. Lett. 95, 137404 (2005). [CrossRef] [PubMed]
V. M. Shalaev, W. Cai, U. K. Chettiar, H.-K. Yuan, A. K. Sarychev, V. P. Drachev, and A. V. Kildishev, “Negative index of refraction in optical metamaterials,” Opt. Lett. 30, 3356–3358 (2005). [CrossRef]
R. W. Ziolkowski, “Propagation in and scattering from a matched metamaterial having a zero index of refraction,” Phys. Rev. E 70, 046608 (2004). [CrossRef]
N. C. Panoiu, R. M. Osgood, S. Zhang, and S. R. J. Brueck, “Zero-n bandgap in photonic crystal superlattices,” J. Opt. Soc. Am. B 23, 506–513 (2006). [CrossRef]
S. Kocaman, R. Chatterjee, N. C. Panoiu, J. F. McMillan, M. B. Yu, R. M. Osgood, D. L. Kwong, and C. W. Wong, “Observation of Zeroth-Order Band Gaps in Negative-Refraction Photonic Crystal Superlattices at Near-Infrared Frequencies,” Phys. Rev. Lett. 102, 203905 (2009). [CrossRef] [PubMed]
U. Kreibig, and M. Vollmer, Optical Properties of Metal Clusters (Springer-Verlag, Berlin, 1995).
A. V. Zayats, I. I. Smolyaninov, and A. A. Maradudin, “Nano-optics of surface plasmon polaritons,” Phys. Rep. 408, 131–314 (2005). [CrossRef]
W. Fan, S. Zhang, N. C. Panoiu, A. Abdenour, S. Krishna, R. M. Osgood, K. J. Malloy, and S. R. J. Brueck, “Second Harmonic Generation from a Nano-patterned Isotropic Nonlinear Material,” Nano Lett. 6, 1027–1030 (2006). [CrossRef]
J. A. H. van Nieuwstadt, M. Sandtke, R. H. Harmsen, F. B. Segerink, J. C. Prangsma, S. Enoch, and L. Kuipers, “Strong Modification of the Nonlinear Optical Response of Metallic Subwavelength Hole Arrays,” Phys. Rev. Lett. 97, 146102 (2006). [CrossRef] [PubMed]
L. Cao, N. C. Panoiu, and R. M. Osgood, “Surface second-harmonic generation from surface plasmon waves scattered by metallic nanostructures,” Phys. Rev. B 75, 205401 (2007). [CrossRef]
Y. Zeng, W. Hoyer, J. J. Liu, S. W. Koch, and J. V. Moloney, “Classical theory for second-harmonic generation from metallic nanoparticles,” Phys. Rev. B 79, 235109 (2009). [CrossRef]
I. I. Smolyaninov, A. V. Zayats, A. Gungor, and C. C. Davis, “Single-Photon Tunneling via Localized Surface Plasmons,” Phys. Rev. Lett. 88, 187402 (2002). [CrossRef] [PubMed]
X. W. Wang, G. C. Schatz, and S. K. Gray, “Ultrafast pulse excitation of a metallic nanosystem containing a Kerr nonlinear material,” Phys. Rev. B 74, 195439 (2006). [CrossRef]
K. Kneipp, Y. Wang, H. Kneipp, L. T. Perelman, I. Itzkan, R. R. Dasari, and M. S. Feld, “Single Molecule Detection Using Surface-Enhanced Raman Scattering (SERS),” Phys. Rev. Lett. 78, 1667 (1997). [CrossRef]
T. H. Taminiau, R. J. Moerland, F. B. Segerink, L. Kuipers, and N. F. van Hulst, “λ /4 Resonance of an Optical Monopole Antenna Probed by Single Molecule Fluorescence,” Nano Lett. 7, 28 (2007). [CrossRef] [PubMed]
B. Knoll, and F. Keilmann, “Near-field probing of vibrational absorption for chemical microscopy,” Nature 399, 134–137 (1999). [CrossRef]
R. M. Roth, N. C. Panoiu, M. M. Adams, R. M. Osgood, C. C. Neacsu, and M. B. Raschke, “Resonant-plasmon field enhancement from asymmetrically illuminated conical metallic-probe tips,” Opt. Express 14, 2921–2913 (2006). [CrossRef] [PubMed]
S. I. Bozhevolnyi, J. Beermann, and V. Coello, “Direct Observation of Localized Second-Harmonic Enhancement in Random Metal Nanostructures,” Phys. Rev. Lett. 90, 197403 (2003). [CrossRef] [PubMed]
T. Ishi, J. Fujikata, K. Makita, T. Baba, and K. Ohashi, “Si nano-photodiode with a surface plasmon antenna,” Jpn. J. Appl. Phys. 44, L364–L366 (2005). [CrossRef]
Z. Yu, G. Veronis, S. Fan, and M. L. Brongersma, “Design of midinfrared photodetectors enhanced by surface plasmons on grating structures,” Appl. Phys. Lett. 89, 151116 (2006). [CrossRef]
R. D. R. Bhat, N. C. Panoiu, S. R. J. Brueck, and R. M. Osgood, “Enhancing the signal-to-noise ratio of an infrared photodetector with a circular metal grating,” Opt. Express 16, 4588–4596 (2008). [CrossRef] [PubMed]
D. M. Schaadt, B. Feng, and E. T. Yu, “Enhanced semiconductor optical absorption via surface plasmon excitation in metal nanoparticles,” Appl. Phys. Lett. 86, 063106 (2005). [CrossRef]
N. C. Panoiu, and R. M. Osgood, “Enhanced optical absorption for photovoltaics via excitation of waveguide and plasmon-polariton modes,” Opt. Lett. 32, 2825–2827 (2007). [CrossRef] [PubMed]
N. C. Panoiu, and R. M. Osgood, “Subwavelength Nonlinear Plasmonic Nanowire,” Nano Lett. 4, 2427–2430 (2004). [CrossRef]
F. Ye, D. Mihalache, B. Hu, and N. C. Panoiu, “Subwavelength Plasmonic Lattice Solitons in Arrays of Metallic Nanowires,” Phys. Rev. Lett. 104, 106802 (2010). [CrossRef] [PubMed]
B. Min, E. Ostby, V. Sorger, E. Ulin-Avila, L. Yang, X. Zhang, and K. Vahala, “High-Q surface-plasmon-polariton whispering-gallery microcavity,” Nature 457, 455–459 (2009). [CrossRef] [PubMed]
M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature 460, 1110–1112 (2009). [CrossRef] [PubMed]
M. P. Nezhand, A. Simic, O. Bondarenko, B. Slutsky, A. Mizrahi, L. Feng, V. Lomakin, and Y. Fainman, “Room temperature subwavelength metallo-dielectric lasers,” Nat. Photonics 4, 395–399 (2010). [CrossRef]
K. Yu, A. Lakhani, and M. C. Wu, “Subwavelength metal-optic semiconductor nanopatch lasers,” Opt. Express 18, 8790–8799 (2010). [CrossRef] [PubMed]
V. Berger, “Second-harmonic generation in monolithic cavities,” J. Opt. Soc. Am. B 14, 1351–1360 (1997). [CrossRef]
D. G. Gusev, I. V. Soboleva, M. G. Martemyanov, T. V. Dolgova, A. A. Fedyanin, and O. A. Aktsipetrov, “Enhanced second-harmonic generation in coupled microcavities based on all-silicon photonic crystals,” Phys. Rev. B 68, 233303 (2003). [CrossRef]
G. Kozyreff, J. L. D. Juarez, and J. Martorell, “Whispering-gallery-mode phase matching for surface second-order nonlinear optical processes in spherical microresonators,” Phys. Rev. A 77, 043817 (2008). [CrossRef]
Y. Xu, M. Han, A. B. Wang, Z. Liu, and J. R. Heflin, “Second order parametric processes in nonlinear silica microspheres,” Phys. Rev. Lett. 100, 163905 (2008). [CrossRef] [PubMed]
C. G. Biris, and N. C. Panoiu, “Second harmonic generation in metamaterials based on homogenous centrosymmetric nanowires,” Phys. Rev. B 81, 195102 (2010). [CrossRef]
J. A. Sanchez-Gil, and A. A. Maradudin, “Dynamic near-field calculations of surface-plasmon polariton pulses resonantly scattered at sub-micron metal defects,” Opt. Express 12, 883–894 (2004). [CrossRef] [PubMed]
E. Centeno, and D. Felbacq, “Rigorous vector diffraction of electromagnetic waves by bidimensional photonic crystals,” J. Opt. Soc. Am. A 17, 320–327 (2000). [CrossRef]
J. I. Dadap, “Optical second-harmonic scattering from cylindrical particles,” Phys. Rev. B 78, 205322 (2008). [CrossRef]
C. I. Valencia, E. R. Mendez, and B. S. Mendoza, “Second-harmonic generation in the scattering of light by an infinite cylinder,” J. Opt. Soc. Am. B 21, 36–44 (2004). [CrossRef]
C. I. Valencia, E. R. Mendez, and B. S. Mendoza, “Second-harmonic generation in the scattering of light by two-dimensional particles,” J. Opt. Soc. Am. B 20, 2150–2161 (2003). [CrossRef]
E. Centeno, and D. Felbacq, “Second-harmonic emission in two-dimensional photonic crystals,” J. Opt. Soc. Am. B 23, 2257–2264 (2006). [CrossRef]
J. Rudnick, and E. A. Stern, “Second-harmonic Radiation from Metal Surfaces,” Phys. Rev. B 4, 4274 (1971). [CrossRef]
J. I. Dadap, J. Shan, K. B. Eisenthal, and T. F. Heinz, “Second-Harmonic Rayleigh Scattering from a Sphere of Centrosymmetric Material,” Phys. Rev. Lett. 83, 4045 (1999). [CrossRef]
M. A. Ordal, R. J. Bell, R. W. Alexander, Jr., L. L. Long, and M. R. Querry, “Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W,” Appl. Opt. 24, 4493–4499 (1985). [CrossRef] [PubMed]
T. F. Heinz, “Second-Order Nonlinear Optical Effects at Surfaces and Interfaces,” in Nonlinear Surface Electromagnetic Phenomena, edited by H. E. Ponath and G. I. Stegeman (Elsevier, Amsterdam, 1991); pp. 353–416.
D. Krause, C. W. Teplin, and C. T. Rogers, “Optical surface second harmonic measurements of isotropic thin-film metals: Gold, silver, copper, aluminum, and tantalum,” J. Appl. Phys. 96, 3626–3635 (2004). [CrossRef]
N. Bloembergen, R. K. Chang, S. S. Jha, and C. H. Lee, “Optical Second-Harmonic Generation in Reflection from Media with Inversion Symmetry,” Phys. Rev. 174, 813–822 (1968). [CrossRef]
M. Liu, T. W. Lee, S. K. Gray, P. Guyot-Sionnest, and M. Pelton, “Excitation of Dark Plasmons in Metal Nanoparticles by a Localized Emitter,” Phys. Rev. Lett. 102, 107401 (2009). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper