Double resonance surface enhanced Raman scattering substrates: an intuitive coupled oscillator model

Double resonance surface enhanced Raman scattering substrates: an intuitive coupled oscillator model Yizhuo Chu, Dongxing Wang, Wenqi Zhu, and Kenneth B. Crozier* School of Engineering and Applied Sciences, Harvard University, 29 Oxford St. Cambridge, Massachusetts 02138, USA *kcrozier@seas.harvard.edu Abstract: The strong coupling between localized surface plasmons and surface plasmon polaritons in a double resonance surface enhanced Raman scattering (SERS) substrate is described by a classical coupled oscillator model. The effects of the particle density, the particle size and the SiO 2 spacer thickness on the coupling strength are experimentally investigated. We demonstrate that by tuning the geometrical parameters of the double resonance substrate, we can readily control the resonance frequencies and tailor the SERS enhancement spectrum. 2011 Optical Society of America OCIS codes: (240.6680) Surface plasmons; (240.6695) surface-enhanced Raman scattering. References and links 1. S. A. Maier, Plasmonics: Fundamentals and Applications (Springer, 2007). 2. W. L. Barnes, A. Dereux, and T. W. Ebbesen, Surface plasmon subwavelength optics, Nature 424(6950), 824 830 (2003). 3. K. A. Willets and R. P. Van Duyne, Localized surface plasmon resonance spectroscopy and sensing, Annu. Rev. Phys. Chem. 58(1), 267 297 (2007). 4. J. Zuloaga and P. Nordlander, On the energy shift between near-field and far-field peak intensities in localized plasmon systems, Nano Lett. 11(3), 1280 1283 (2011). 5. K. L. Kelly, E. Coronado, L. L. Zhao, and G. C. Schatz, The optical properties of metal nanoparticles: the influence of size, shape, and dielectric environment, J. Phys. Chem. B 107(3), 668 677 (2003). 6. 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(9), 1667 1670 (1997). 7. A. D. McFarland, M. A. Young, J. A. Dieringer, and R. P. Van Duyne, Wavelength-scanned surface-enhanced Raman excitation spectroscopy, J. Phys. Chem. B 109(22), 11279 11285 (2005). 8. T. H. Reilly, S. Chang, J. D. Corbman, G. C. Schatz, and K. L. Rowlen, Quantitative evaluation of plasmon enhanced Raman scattering from nanoaperture arrays, J. Phys. Chem. C 111(4), 1689 1694 (2007). 9. P. L. Stiles, J. A. Dieringer, N. C. Shah, and R. P. Van Duyne, Surface-enhanced Raman spectroscopy, Ann. Rev. Anal. Chem. 1(1), 601 626 (2008). 10. D. L. Jeanmaire and R. P. Van Duyne, Surface Raman spectroelectrochemistry. Heterocyclic, aromatic, and aliphatic-amines adsorbed on anodized silver electrode, J. Electroanal. Chem. 84(1), 1 20 (1977). 11. Y. Chu, M. G. Banaee, and K. B. Crozier, Double-resonance plasmon substrates for surface-enhanced Raman scattering with enhancement at excitation and stokes frequencies, ACS Nano 4(5), 2804 2810 (2010). 12. J. B. Jackson, S. L. Westcott, L. R. Hirsch, J. L. West, and N. J. Halas, Controlling the surface enhanced Raman effect via the nanoshell geometry, Appl. Phys. Lett. 82(2), 257 259 (2003). 13. C. E. Talley, J. B. Jackson, C. Oubre, N. K. Grady, C. W. Hollars, S. M. Lane, T. R. Huser, P. Nordlander, and N. J. Halas, Surface-enhanced Raman scattering from individual au nanoparticles and nanoparticle dimer substrates, Nano Lett. 5(8), 1569 1574 (2005). 14. M. G. Banaee and K. B. Crozier, Mixed dimer double-resonance substrates for surface-enhanced Raman spectroscopy, ACS Nano 5(1), 307 314 (2011). 15. M. G. Banaee and K. B. Crozier, Gold nanorings as substrates for surface-enhanced Raman scattering, Opt. Lett. 35(5), 760 762 (2010). 16. J. R. Lakowicz, Radiative decay engineering 5: metal-enhanced fluorescence and plasmon emission, Anal. Biochem. 337(2), 171 194 (2005). 17. K. Kneipp, H. Kneipp, I. Itzkan, R. R. Dasari, and M. S. Feld, Surface-enhanced Raman scattering and biophysics, J. Phys. Condens. Matter 14(18), R597 R624 (2002). 18. Y. Chu and K. B. Crozier, Experimental study of the interaction between localized and propagating surface plasmons, Opt. Lett. 34(3), 244 246 (2009). 19. Y. Chu, E. Schonbrun, T. Yang, and K. B. Crozier, Experimental observation of narrow surface plasmon resonances in gold nanoparticle arrays, Appl. Phys. Lett. 93(18), 181108 (2008). (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14919

20. B. Luk yanchuk, N. I. Zheludev, S. A. Maier, N. J. Halas, P. Nordlander, H. Giessen, and C. T. Chong, The Fano resonance in plasmonic nanostructures and metamaterials, Nat. Mater. 9(9), 707 715 (2010). 21. N. Liu, L. Langguth, T. Weiss, J. Kästel, M. Fleischhauer, T. Pfau, and H. Giessen, Plasmonic analogue of electromagnetically induced transparency at the Drude damping limit, Nat. Mater. 8(9), 758 762 (2009). 22. S. Zhang, D. A. Genov, Y. Wang, M. Liu, and X. Zhang, Plasmon-induced transparency in metamaterials, Phys. Rev. Lett. 101(4), 047401 (2008). 23. A. Ghoshal and P. G. Kik, Theory and simulation of surface plasmon excitation using resonant metal nanoparticle arrays, J. Appl. Phys. 103(11), 113111 (2008). 24. D. Wang, T. Yang, and K. B. Crozier, Optical antennas integrated with concentric ring gratings: electric field enhancement and directional radiation, Opt. Express 19(3), 2148 2157 (2011). 25. S. Mukherjee, H. Sobhani, J. B. Lassiter, R. Bardhan, P. Nordlander, and N. J. Halas, Fanoshells: nanoparticles with built-in Fano resonances, Nano Lett. 10(7), 2694 2701 (2010). 26. Q. F. Xu, S. Sandhu, M. L. Povinelli, J. Shakya, S. H. Fan, and M. Lipson, Experimental realization of an onchip all-optical analogue to electromagnetically induced transparency, Phys. Rev. Lett. 96(12), 123901 (2006). 27. S. M. Ma, H. Xu, and B. S. Ham, Electromagnetically-induced transparency and slow light in GaAs/AlGaAs multiple quantum wells in a transient regime, Opt. Express 17(17), 14902 14908 (2009). 28. C. L. Garrido Alzar, M. A. G. Martinez, and P. Nussenzveig, Classical analog of electromagnetically induced transparency, Am. J. Phys. 70(1), 37 (2002). 29. Y. S. Joe, A. M. Satanin, and C. S. Kim, Classical analogy of Fano resonances, Phys. Scr. 74(2), 259 266 (2006). 30. S. E. Harris, Electromagnetically induced transparency, Phys. Today 50(7), 36 42 (1997). 31. W. R. Holland and D. G. Hall, Frequency shift of an electric-dipole resonance near a conducting surface, Phys. Rev. Lett. 52(12), 1041 1044 (1984). 32. I. P. Radko, S. I. Bozhevolnyi, G. Brucoli, L. Martin-Moreno, F. J. Garcia-Vidal, and A. Boltasseva, Efficiency of local surface plasmon polariton excitation on ridges, Phys. Rev. B 78(11), 115115 (2008). 1. Introduction Noble metal nanostructures exhibit unique optical properties due to the excitation of localized surface plasmons (LSPs) and/or surface plasmon polaritons (SPPs) [1 5]. These phenomena can lead to strong electromagnetic field enhancement near the nanostructure surfaces, and this has been employed for enhanced fluorescence, surface-enhanced Raman scattering (SERS), and biosensing [6 17]. The coupling between plasmonic phenomena provides additional interesting properties such as resonance shift, hybridized modes and Fano line shape [18 25]. The effects of coupling are attracting considerable interest across a variety of systems, ranging from coupled microresonators to coupled quantum wells [26,27]. It has been shown that the strong coupling between a broad radiative mode and a narrow dark mode results in a pronounced Fano line shape [20 22, 26 29]. Analogies can be drawn between these phenomena and electromagnetically induced transparency (EIT) in an atomic system [30]. This phenomenon can be modeled by a system with two classical oscillators coupled by a weak spring [28, 29]. EIT-like properties have been demonstrated using plasmonic nanoparticles in which the broad radiative mode is represented by the dipole plasmon mode and the narrow dark mode is represented by quadrupole plasmon mode [20 22]. We recently investigated strong coupling between LSPs and SPPs in a double resonance plasmonic structure consisting of a gold nanoparticle array, a SiO 2 spacer and a continuous gold film [11, 18]. Two hybridized modes were observed and the structure was employed for SERS [11]. We demonstrated that the double resonance structure can provide SERS enhancement factor of nearly 10 9, which was more than two orders of magnitude larger than the largest enhancement we measured for gold nanoparticles on glass substrates [11]. In that work [11], however, we did not develop a model. The plasmon resonance frequencies of the structure are critical for its performance as a SERS substrate, but an intuitive model for how they follow from the structure parameters was lacking. They therefore had to be predicted by finite difference time domain (FDTD) simulations. Similarly, an unusual extinction cross section line shape was observed, but could not be physically interpreted. Here, we address this by studying the use of a classical coupled oscillator model to describe the coupling between LSPs and SPPs. We investigate the influence of the geometrical parameters on the coupling strength. We find that changes on the nanoparticle density, the nanoparticle size and the SiO 2 spacer thickness lead to the change on the coupling strength, therefore lead to the shift of the hybridized resonance frequencies. For the application of SERS, we demonstrate that by (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14920

choosing the geometrical parameters of the device appropriately, different portions of the Raman spectrum can be selectively enhanced. 2. Method A schematic diagram of the double resonance SERS substrate is shown in Fig. 1(a). It consists of a gold nanoparticle array, an SiO 2 spacer and a continuous gold film. The fabrication method is described in Ref. 18. Fig. 1(b) shows a SEM image of a fabricated device. The gold nanoparticles are cylinders with thickness of 40 nm. The nanoparticle array is a rectangular lattice with overall extent of about 100 μm 100 μm. We fabricate devices with different periods in the range from 400 nm to 900 nm. The gold film is 100 nm thick. The thickness of the SiO 2 spacer is varied from 20 nm to 60 nm. The extinction cross sections are determined using the method described in Ref.18. The fabricated devices are illuminated with collimated and polarized white light from a Xe lamp at normal incidence. The polarization is parallel to one of the axes of the rectangular lattice, i.e. the x axis in Fig. 1(a). The reflection is collected by a long working distance objective. The reflection spectra are normalized by the spectra of the area without the gold nanoparticle pattern, i.e. just the SiO 2 spacer and underlying gold film. The gold film is sufficiently thick to make the transmission negligible. The extinction cross section associated with each gold nanoparticle can be calculated by the relation C ext = (1-R) Λ x Λ y, where R is reflectance, and Λ x and Λ y are the periods along x and y axes, respectively. Fig. 1. (a) Schematic diagram of the double resonance SERS substrate. Actual device has hundreds of periods along x and y axes; only a small portion of it depicted here. (b) Scanning electron micrographs (SEM) of a fabricated device. (c) Extinction cross section spectrum of a double resonance SERS substrate. SiO 2 thickness is 24 nm, gold nanoparticle diameter is 120 nm and array periods along x and y axes are 780 nm. Black curve: experimental data. Red circles: fitted coupled oscillator model. As shown in Fig. 1(a), the double resonance substrate is illuminated at normal incidence. The LSPs on the gold nanoparticles can be directly excited by the illumination. The reciprocal lattice vector provided by the gold nanoparticle array enables coupling to the SPP, which exists at the interface between the SiO 2 spacer and the gold film. Our previous work showed that the LSP and SPP resonances exhibit an anti-crossing behavior, with a splitting between resonances due to strong coupling [11,18]. The black curve in Fig. 1(c) shows the measured extinction spectrum of a typical device. The gold nanoparticle diameter is 120 nm. The periods Λ x and Λ y are both 780 nm. The thickness of the SiO 2 is 24 nm, as measured by ellipsometry. Two hybridized resonances with similar line width can be seen. Here, we utilize the classic coupled oscillator model to describe the coupling between the LSP and SPP. The LSP in our structure is represented by oscillator 1, which is driven by an external source E(t) = E 0 e -iωt, where E 0 and ω are the amplitude and the frequency of the electric field, respectively. The SPP mode is represented by oscillator 2, and can be excited only through the coupling between the two oscillators. The coupling between the two oscillators is provided by scattering by the gold nanoparticle array. The motion of the charges of oscillators 1 and 2 are given as follows: (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14921

dr() t dr() t ee e + + = (1) 2 iωt 1 1 2 0 γ 2 1 ω1r1() t gr2() t dt dt me 2 dr2() t dr2() t 2 + γ 2 2 + ω2r2() t gr1() t = 0 (2) dt dt where ω 1 and ω 2 are the resonance frequencies of the uncoupled LSP mode and the uncoupled SPP mode, respectively. r 1 and r 2 are the displacements of an electron from its equilibrium position. m e is the effective mass of the electrons. The anti-crossing occurs when ω 1 ω 2. γ 1 and γ 2 are the damping coefficients, and related to the loss of the oscillators. Since the LSP is a radiative mode and the SPP is a dark mode, we have γ 2 << γ 1. Both γ 1 and γ 2 are << ω 1. g is defined as the coupling coefficient, and we have g<< ω 1 2. As we discuss below in detail, the strength of the coupling between the two oscillators determines the splitting between the hybridized resonances. The energy dissipation can be written as: 2 2 2 2 2 πie E0 ωω ( 2 ω iγω 2 ) P( ω) = (3) 2 2 2 2 2 me [( ω1 ω iγω 1 )( ω2 ω iγω 2 ) g ] To find the parameters {ω 1, ω 2, γ 1, γ 2, g}, we fit the experimental data (black curve in Fig. 1(c)) to Eq. (3). The result is shown as the curve comprising red circles in Fig. 1(c). Good agreement between the experimental spectrum and the fitting curve can be seen. 3. Experimental results 3.1 Effects of geometrical parameters on the coupling strength We first investigate the influence of the nanoparticle density on the coupling strength. A series of structures is fabricated. The SiO 2 thickness is 24 nm. For all the structures, the diameter of the gold nanoparticles is 120 nm. The period Λ x along the x axis, i.e. parallel to the direction of the polarization of the incidence is fixed at 780 nm. To change the particle density, we vary the period Λ y along the y axis of the array from 400 nm to 900 nm. For normal incidence, for a particular frequency, the first order SPP mode occurs when the wave vector satisfies the equation k spp = 2π/Λ x where k spp is the wave vector of the SPP for that frequency. Thus, for a given SiO 2 thickness, the first order SPP resonance frequency is only determined by the period along the x axis. Therefore the frequency of the uncoupled SPP mode is fixed for all the structures, though the period along the y axis changes. The measured extinction spectra of three structures are shown in Fig. 2(a). It can be seen that the splitting between the two resonances is larger when the period along the y axis is smaller. The parameters of the coupled oscillator model are fitted to the experimental data, with the results also shown Fig. 2(a) as black, red and green circles. The coupling coefficients g are calculated and plotted in Fig. 2(b) as a function of Λ y. It can be seen that the coupling coefficient decreases with decrease in particle density, i.e. with increase in the y period Λ y. This is due to the fact that, for the structure with larger particle density, there are more scattering centers. For larger particle densities, more energy from the incident illumination is coupled into the SPP and more of the SPP s energy is radiated by scattering by the nanoparticles. The resonance frequencies of uncoupled LSP and SPP, the damping coefficients and the coupling coefficient are listed in Table 1. It can be seen that, for the fabricated structures, ω 1 is slightly larger than ω 2, indicating that the uncoupled resonance frequency of the LSP is slightly larger than that of the SPP mode. The dampling coefficient γ 1 is about one order of magnitude larger than γ 2, indicating that oscillator 1 is bright radiative mode. (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14922

Fig. 2. (a) Extinction cross section spectra of three devices with different periods along y axis (Λ y). Experimental measurements: solid curves. Coupled oscillator model: circles. For all the devices, the SiO 2 spacer is 24 nm thick, the diameter of the gold nanoparticles is 120 nm and the period Λ x of the array along the x axis is 780 nm. Black curve: Λ y = 400 nm. Red curve: Λ y = 600 nm. Green curve: Λ y = 900 nm. (b) Coupling coefficient g as a function of array period Λ y. Table 1. Fitting Parameters for the Devices with Different Periods Along y Axis Λ y (nm) ω 1( 10 15 ) γ 1( 10 14 ) ω 2( 10 15 ) γ 2( 10 13 ) g( 10 29 ) 400 2.373 2.215 2.301 1.616 4.843 500 2.385 1.910 2.306 1.805 4.362 600 2.375 1.620 2.308 1.278 4.006 780 2.345 1.883 2.310 1.246 3.484 900 2.373 1.794 2.311 1.604 3.171 We next investigate the influence of the particle size. A series of structures is fabricated. The periods along x and y axes are fixed at 780 nm. The LSP resonance position mainly depends on the particle s dimension along the same direction of the polarization, i.e. the x axis. As shown in Fig. 3(a), the nanoparticles are chosen to have a rod-shaped cross section in the xy plane, comprising a rectangle with a semicircle at each of its two ends along the x axis. The diameter of each semicircle matches the side length of the rectangle along the y axis. The total length of the nanorod along the x axis, including the semicircular ends, is fixed at 120 nm. The width dy of nanorod, i.e. its dimension along the y axis, is varied from 40 nm to 100 nm. Figure 3(b) shows the extinction spectra of three structures for which the nanorod widths dy are 50, 80, and 100 nm. A SEM image of a fabricated device with nanorod array is shown in Fig. 3(c).One can see that the splitting of the two resonance increases as the nanoparticle width increases, indicating that larger particles provide greater coupling. The coupling coefficients calculated by fitting the experimental spectra according to the coupled oscillator model are plotted in Fig. 3(d). The resonance frequencies of uncoupled LSP and SPP, the damping coefficients and the coupling coefficient are listed in Table 2. (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14923

Fig. 3. (a) Schematic illustration of nanoparticle cross section. (b) Extinction cross section spectra of three devices with different nanoparticle cross sections. Experimental measurements: solid curves. Coupled oscillator model: circles. For all the devices, the SiO 2 is 24 nm thick, the length of the gold nanoparticles along the x axis is 120 nm and the periods of the array along both the x and y axes are 780 nm. The dimension of the nanoparticle along the y axis is denoted as dy. Black curve: dy = 50 nm. Red curve: dy = 80 nm. Green curve: dy = 100 nm. (c) Scanning electron micrographs (SEM) of fabricated device consisting of gold nanorod array. (d) Coupling coefficient of double resonance structures as a function of nanoparticle width (dy). Table 2. Fitting Parameters for the Devices with Different Nanoparticle Width dy dy (nm) ω 1( 10 15 ) γ 1( 10 14 ) ω 2( 10 15 ) γ 2( 10 13 ) g( 10 29 ) 40 2.300 1.529 2.313 1.882 2.606 50 2.332 1.683 2.311 1.609 2.706 60 2.356 1.344 2.311 1.398 2.895 70 2.363 1.511 2.312 1.752 3.038 80 2.368 1.620 2.312 1.915 3.143 90 2.362 1.359 2.311 1.555 3.270 100 2.356 1.399 2.311 1.634 3.389 Finally, we investigate the influence of the SiO 2 spacer thickness. A series of structures with different SiO 2 thicknesses is fabricated. The periods along the x and y axes are fixed at 780 nm. With increasing SiO 2 spacer thickness, the frequencies of both the uncoupled LSP mode and the uncoupled SPP mode change. For the LSP mode, the interaction of the excited dipole on the nanoparticle and its image dipole becomes weaker for a thicker SiO 2 spacer, resulting in a blueshift of the LSP resonance [31]. On the other hand, the resonance of the uncoupled SPP mode has a slight redshift due to the increase of the effective index of the space above the gold film. To match the frequencies of the uncoupled LSP and SPP modes, (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14924

with the period of the array fixed, the nanoparticle length along the x axis has to be increased as the SiO 2 spacer thickness is increased. The extinction spectra of two structures with different SiO 2 thicknesses are shown in Fig. 4(a). For the structure with SiO 2 thickness of 21 nm, the nanoparticle length dx along the x-axis is 120 nm. For the structure with SiO 2 thickness of 56 nm, the nanoparticle length dx is 155 nm. It is obvious that the splitting of the two resonances for the device with the 56 nm SiO 2 spacer is larger than that for the device with the 21 nm SiO 2 spacer, indicating that the coupling strength for the structure with the thicker SiO 2 spacer is larger. However, as discussed above, the increasing coupling strength could come from the nanoparticles being larger. To understand the influence of the SiO 2 thickness on the coupling strength, we make another series of structures. For these structures, the periods along x and y axes are fixed at 780 nm, and the nanoparticle length dx along the x axis is fixed at 150 nm, so that the effect of the nanoparticle size can be eliminated. Here, we assume that the coupling between the two oscillators is independent of the resonance frequencies of the uncoupled oscillators. The extinction spectra of the devices with SiO 2 thicknesses of 21 nm and 56 nm are shown in Fig. 4(b) and 4(c). By fitting the extinction spectra to the coupled oscillator model, the coupling coefficient is found and plotted in Fig. 4(d). The fitting curves are shown as circles in Fig. 4(b) and 4(c). The features around 600 nm in Fig. 4(b) and around 650 nm in Fig. 4(c) are the second order SPP modes. They are not considered in the coupled oscillator model, and therefore contribute to differences between the coupled oscillator model fit and the experimental data. The fitting parameters are listed in Table 3. It can be seen that with an increasing thickness of SiO 2, the coupling coefficient decreases when the SiO 2 thickness is smaller than 30 nm. On the other hand, when the SiO 2 thickness is larger than 30 nm, the coupling coefficient increases with thickness. To explain this phenomenon, we investigate the coupling efficiency of a dielectric grating coupler on top of a continuous gold film numerically, using the FDTD method [32]. The dielectric grating is a 2D square array. The diameters of the dielectric nanoparticles are 150 nm and the x- and y- periods of the array are 780 nm. The spacer between the dielectric particle layer and the gold film is either air or SiO 2. The structures are illuminated by a polarized plane wave at normal incidence. The electric field component in the direction perpendicular to the surface of the structure, E z, is recoded by a time monitor (indicated by black dot in Fig. 5) placed on the surface of the gold film at the 1/4 of period away from the center of the dielectric particle to measure amplitude of the first SPP mode on the top surface of the gold film. The first order SPP mode occurs around 800 nm. When the thickness of the spacer changes, the resonance wavelength of the first order SPP mode undergoes a slight shift due to the change of the local refractive index. For each structure, we find E z over a range of illumination wavelengths. Its peak value is denoted as E zspp and plotted in Fig. 5 as a function of the spacer thickness. The amplitude of the SPP decreases with thickness when the thickness is smaller than 20 nm. When the thickness is between 20 nm and 120 nm, the amplitude increases with thickness. When the thickness is greater than 120 nm, the amplitude of the SPP decreases with thickness. It can be seen that the coupling coefficient of the double resonance structure (Fig. 4(d)) also follows this pattern. We conclude that the variation in coupling coefficient vs spacer thickness for the double resonance structure arises from the dependence of the efficiency with which SPPs on a gold film are excited by an array of scatterers on the distance from the scatterers to the film. (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14925

Fig. 4. (a) Measured extinction cross section spectra of double resonance devices. h denotes the thickness of the SiO2 spacer. dx denotes the length of the gold nanoparticles along x axis. For all devices, the periods of the array along the x and y axes are 780 nm. Black curve: h = 21nm, dx = 120nm. Red curve: h = 56nm, dx = 155nm. (b) Extinction cross section spectra for device with 21 nm thick SiO2 spacer. Black line: experimental results. Red circles: coupled oscillator model fit. Nanoparticles have dx = 150nm. (c) Extinction cross section spectra for device with 56 nm thick SiO2 spacer. Black line: experimental results. Red circles: coupled oscillator model fit. Nanoparticles have dx = 150nm. (d) Coupling coefficient of double resonance structures as a function of SiO2 thickness. Table 3. Fitting Parameters for the Devices with Different Thicknesses of SiO2 thickness (nm) 21 31 39 48 56 #146801 - $15.00 USD (C) 2011 OSA ω1( 1015) γ1( 1014) ω2( 1015) γ2( 1013) g( 1029) 2.089 2.210 2.262 2.284 2.312 1.359 1.129 1.222 1.465 1.457 2.308 2.286 2.258 2.229 2.204 0.812 1.190 1.875 1.673 1.549 5.083 4.858 4.994 5.046 5.403 Received 29 Apr 2011; revised 6 Jun 2011; accepted 6 Jun 2011; published 19 Jul 2011 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14926

Fig. 5. Perpendicular component of electric field E zspp on the surface of gold film as a function of spacer thickness. Position at which electric field E zspp is monitored is indicated by black dot and arrow. Dielectric grating has period of 780 nm along x and y axes. Structure is illuminated at normal incidence with a plane wave. Each structure consists of a 2D SiO 2 grating on top of a spacer and a continuous gold film. Black dots: the spacer is air. Red dots: the spacer is SiO 2. 3.2 SERS enhancement spectra We have shown above that the double resonance substrate can be modeled by coupled oscillators, and that the splitting between the resonances can be modified by appropriate choice of the substrate s geometric parameters. We now show that this enables portions of the SERS spectrum to be selectively enhanced. SERS measurements are performed on benzenthiol monolayers formed on two double resonance substrates. The first has an SiO 2 spacer thickness of 23 nm, gold nanoparticle diameter of 133 nm and x- and y-axis periods of 780 nm. The second has an SiO 2 spacer thickness of 60 nm, gold nanoparticle diameter of 150 nm and x- and y-axis periods of 760 nm. For both structures, the geometrical parameters are chosen so ensure that the LSP and SPP are strongly coupled. The extinction spectra are found by measuring the reflection from the samples using an optical microscope equipped with a spectrometer [11]. The objective lens has a numerical aperture of 0.12. The illumination, therefore, does not comprise the normally incidence plane waves depicted in Fig. 1(a), but the range of angular spread is relatively small due to the modest numerical aperture. The extinction spectra measured for the two substrates are shown in Fig. 6(a) and 6(b) as red curves. For both substrates, two hybridized resonances with approximate equal intensity and line width can be seen. The small features located between the two primary resonances are due to the angular spread of the illumination [11]. Benzenthiol molecules are coated on the two structures and form monolayers on the surfaces of the gold nanoparticles. This is performed by immersing the devices in a 3mM solution of benzenethiol in ethanol for 1 hour, rinsing with neat ethanol and then blow drying with nitrogen. SERS measurements are carried out using a confocal Raman microscope equipped with a 5 objective (NA = 0.12), with excitation at λ = 783nm. The incident laser power at sample is ~20 mw and the acquisition time for each spectrum is 5 seconds. The black curves in Fig. 6(a) and 6(b) show the SERS spectra obtained from the two structures. The excitation laser wavelength is indicated with a dashed line in each spectrum, and is located near the high energy resonance peak for both structures. It can be seen, however, that the low energy resonance of the substrate with the 60 nm thick SiO 2 spacer is red shifted from that of the substrate with the 23nm SiO 2 spacer. This phenomenon is due to the fact that the coupling strength is larger for the 60 nm SiO 2 structure than that for the 23 nm SiO 2 structure, as discussed above. The SERS spectrum of the substrate with 23 nm thick SiO 2 has a distinctly strong signal around 1100cm 1 corresponding to the low energy resonance at 855 nm, while the SERS spectrum of the substrate with 60 nm thick SiO 2 has a strong signal around 1500 cm 1 corresponding to the low energy resonance at 890 nm. The background around 1350cm 1 (at 875nm) is from the SiO 2 layer. One can see (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14927

that with different designs of the double resonance SERS substrate, the different portions of SERS spectrum can be selectively enhanced. 4. Conclusions Fig. 6. Measured extinction (left hand axes) and SERS (right hand axes) spectra of double resonance structures on which monolayers of benzenthiol are formed. (a) Double resonance structure with 23 nm thick SiO 2 spacer, gold nanoparticles with diameters of 133 nm and x- and y-axis periods of 780 nm. (b) Double resonance structure with 60 nm thick SiO 2 spacer, gold nanoparticles with diameters of 150 nm and x- and y-axis periods of 760 nm. Dashed line: wavelength of the excitation laser (λ = 783 nm). Strong coupling between LSPs and SPPs in double resonance SERS substrates has been explained by a classical coupled oscillator model. We investigated the influence of the geometrical parameters of the substrate on the coupling strength between LSPs and SPPs. We found that by increasing the particle density and the particle size, the coupling between the LSPs and SPPs becomes larger, leading to enhanced splitting between the two hybridized resonances. The SiO 2 spacer thickness was found to play an important role in determining the coupling strength. We showed that proper choice of the substrate geometrical parameters enables selective enhancement of different portions of the SERS spectra to be achieved. Acknowledgments The authors thank Eric Diebold, Dr Paul Peng and Prof. Eric Mazur for use of their Raman confocal microscope and thank Dr. Tal Ellenbogen, Dr. Mohamad G. Banaee and Kai Wang for helpful discussions. This work was supported by the Defense Advanced Research Projects Agency (DARPA) SERS S&T Fundamentals program under grant no. FA9550-08-1-0285, issued by the Air Force Office of Scientific Research, by the DARPA N/MEMS S&T Fundamentals program under grant no. N66001-10-1-4008 issued by the Space and Naval Warfare Systems Center Pacific (SPAWAR), and by the National Science Foundation under grant no. ECCS-0747560 (CAREER award). Fabrication work was carried out at the Harvard Center for Nanoscale Systems, which is supported by the National Science Foundation. (C) 2011 OSA 1 August 2011 / Vol. 19, No. 16 / OPTICS EXPRESS 14928