Optical response of supported gold nanodisks

It is shown that the ellipsometric spectra of short range ordered planar arrays of gold nanodisks supported on glass substrates can be described by modeling the nanostructured arrays as uniaxial homogeneous layers with dielectric functions of the Lorentz type. However, appreciable deviations from experimental data are observed in calculated spectra of irradiance measurements. A qualitative and quantitative description of all measured spectra is obtained with a uniaxial effective medium dielectric function in which the nanodisks are modeled as oblate spheroids. Dynamic depolarization factors in the long-wavelength approximation and interaction with the substrate are considered. Similar results are obtained calculating the optical spectra using the island-film theory. Nevertheless, a small in-plane anisotropy and quadrupolar coupling effects reveal a very complex optical response of the nanostructured arrays. ©2011 Optical Society of America OCIS codes: (240.2130) Ellipsometry and polarimetry; (240.6680) Surface plasmons; (160.1190) Anisotropic optical materials; (160.1245) Artificially engineered materials; (160.4760) Optical properties. References and links 1. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley & Sons, 1983). 2. N. V. Voshchinnikov and V. G. Farafonov, “Optical properties of spheroidal particles,” Astrophys. Space Sci. 204(1), 19–86 (1993). 3. M. E. Stewart, C. R. Anderton, L. B. Thompson, J. Maria, S. K. Gray, J. A. Rogers, and R. G. Nuzzo, “Nanostructured plasmonic sensors,” Chem. Rev. 108(2), 494–521 (2008). 4. N. G. Khlebtsov and L. A. Dykman, “Optical properties and biomedical applications of plasmonic nanoparticles,” J. Quant. Spectrosc. Radiat. Transf. 111(1), 1–35 (2010). 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. V. K. S. Hsiao, Y. B. Zheng, B. K. Juluri, and T. J. Huang, “Light-driven plasmonic switches based on Au nanodisk arrays and photoresponsive liquid crystals,” Adv. Mater. (Deerfield Beach Fla.) 20(18), 3528–3532 (2008). 7. M. Svedendahl, S. Chen, A. Dmitriev, and M. Käll, “Refractometric sensing using propagating versus localized surface plasmons: a direct comparison,” Nano Lett. 9(12), 4428–4433 (2009). 8. H. Fredriksson, Y. Alaverdyan, A. Dmitriev, C. Langhammer, D. S. Sutherland, M. Zäch, and B. Kasemo, “Hole–mask colloidal lithography,” Adv. Mater. (Deerfield Beach Fla.) 19(23), 4297–4302 (2007). 9. T. Pakizeh, M. S. Abrishamian, N. Granpayeh, A. Dmitriev, and M. Käll, “Magnetic-field enhancement in gold nanosandwiches,” Opt. Express 14(18), 8240–8246 (2006). 10. A. Dmitriev, T. Pakizeh, M. Käll, and D. S. Sutherland, “Gold-silica-gold nanosandwiches: tunable bimodal plasmonic resonators,” Small 3(2), 294–299 (2007). 11. T. Pakizeh, A. Dmitriev, M. S. Abrishamian, N. Granpayeh, and M. Käll, “Structural asymmetry and induced optical magnetism in plasmonic nanosandwiches,” J. Opt. Soc. Am. B 25(4), 659–667 (2008). 12. T. Pakizeh, C. Langhammer, I. Zorić, P. Apell, and M. Käll, “Intrinsic Fano interference of localized plasmons in Pd nanoparticles,” Nano Lett. 9(2), 882–886 (2009). 13. B. Sepúlveda, J. B. González-Díaz, A. García-Martín, L. M. Lechuga, and G. Armelles, “Plasmon-induced magneto-optical activity in nanosized gold disks,” Phys. Rev. Lett. 104(14), 147401 (2010). 14. P. Hanarp, M. Käll, and D. S. Sutherland, “Optical properties of short range ordered arrays of nanometer gold disks prepared by colloidal lithography,” J. Phys. Chem. B 107(24), 5768–5772 (2003). #146347 $15.00 USD Received 21 Apr 2011; revised 20 May 2011; accepted 28 May 2011; published 7 Jun 2011 (C) 2011 OSA 20 June 2011 / Vol. 19, No. 13 / OPTICS EXPRESS 12093 15. C. Langhammer, Z. Yuan, I. Zorić, and B. Kasemo, “Plasmonic properties of supported Pt and Pd nanostructures,” Nano Lett. 6(4), 833–838 (2006). 16. C. Langhammer, B. Kasemo, and I. Zorić, “Absorption and scattering of light by Pt, Pd, Ag, and Au nanodisks: absolute cross sections and branching ratios,” J. Chem. Phys. 126(19), 194702 (2007). 17. Y. B. Zheng, B. Kiraly, S. Cheunkar, T. J. Huang, and P. S. Weiss, “Incident-angle-modulated molecular plasmonic switches: a case of weak exciton-plasmon coupling,” Nano Lett. 11(5), 2061–2065 (2011). 18. Y. B. Zheng, B. K. Juluri, L. Lin Jensen, D. Ahmed, M. Lu, L. Jensen, and T. J. Huang, “Dynamical tuning of plasmon-exciton coupling in arrays of nanodisk-J-aggregate complexes,” Adv. Mater. (Deerfield Beach Fla.) 22(32), 3603–3607 (2010). 19. T. Yang and K. B. Crozier, “Surface plasmon coupling in periodic metallic nanoparticle structures: a semianalytical model,” Opt. Express 16(17), 13070–13079 (2008). 20. T. Yang and K. B. Crozier, “Dispersion and extinction of surface plasmons in an array of gold nanoparticle chains: influence of the air/glass interface,” Opt. Express 16(12), 8570–8580 (2008). 21. K. B. Crozier, E. Togan, E. Simsek, and T. Yang, “Experimental measurement of the dispersion relations of the surface plasmon modes of metal nanoparticle chains,” Opt. Express 15(26), 17482–17493 (2007). 22. E. Ringe, J. M. McMahon, K. Sohn, C. Cobley, Y. Xia, J. Huang, G. C. Schatz, L. D. Marks, and R. P. Van Duyne, “Unraveling the effects of size, composition, and substrate on the localized surface plasmon resonance frequencies of gold and silver nanocubes: a systematic single-particle approach,” J. Phys. Chem. C 114(29), 12511–12516 (2010). 23. C. P. Burrows and W. L. Barnes, “Large spectral extinction due to overlap of dipolar and quadrupolar plasmonic modes of metallic nanoparticles in arrays,” Opt. Express 18(3), 3187–3198 (2010). 24. S. Berthier, “Anisotropic effective medium theories,” J. Phys. I 4(2), 303–318 (1994). 25. T. Yamaguchi, S. Yoshida, and A. Kinbara, “Optical effect of the substrate on the anomalous absorption of aggregated silver films,” Thin Solid Films 21(1), 173–187 (1974). 26. V. A. Fedotov, V. I. Emel’yanov, K. F. MacDonald, and N. I. Zheludev, “Optical properties of closely packed nanoparticle films: spheroids and nanoshells,” J. Opt. A, Pure Appl. Opt. 6(2), 155–160 (2004). 27. D. Bedeaux and J. Vlieger, Optical Properties of Surfaces (Imperial College Press, 2001). 28. R. Lazzari and I. Simonsen, “GRANFILM: a software for calculating thin-layer dielectric properties and Fresnel coefficients,” Thin Solid Films 419(1-2), 124–136 (2002). 29. H. Wormeester, E. S. Kooij, A. Mewe, S. Rekveld, and B. Poelsema, “Ellipsometric characterisation of heterogeneous 2D layers,” Thin Solid Films 455–456, 323–334 (2004). 30. H. Wormeester, A.-I. Henry, E. S. Kooij, B. Poelsema, and M. P. Pileni, “Ellipsometric identification of collective optical properties of silver nanocrystal arrays,” J. Chem. Phys. 124(20), 204713 (2006). 31. J. M. Flores-Camacho, L. D. Sun, N. Saucedo-Zeni, G. Weidlinger, M. Hohage, and P. Zeppenfeld, “Optical anisotropies of metal clusters supported on a birefringent substrate,” Phys. Rev. B 78(7), 075416 (2008). 32. M. Losurdo, M. Bergmair, G. Bruno, D. Cattelan, C. Cobet, A. Martino, K. Fleischer, Z. Dohcevic-Mitrovic, N. Esser, M. Galliet, R. Gajic, D. Hemzal, K. Hingerl, J. Humlicek, R. Ossikovski, Z. V. Popovic, and O. Saxl, “Spectroscopic ellipsometry and polarimetry for materials and systems analysis at the nanometer scale: state-ofthe-art, potential, and perspectives,” J. Nanopart. Res. 11(7), 1521–1554 (2009). 33. K. Johansen, H. Arwin, I. Lundström, and B. Liedberg, “Imaging surface plasmon resonance sensor based on multiple wavelengths: sensitivity considerations,” Rev. Sci. Instrum. 71(9), 3530–3538 (2000). 34. R. M. A. Azzam and N. M. Bashara, Ellipsometry and Polarized Light (North Holland, 1977). 35. G. Xu, Y. Chen, M. Tazawa, and P. Jin, “Influence of dielectric properties of a substrate upon plasmon resonance spectrum of supported Ag nanoparticles,” Appl. Phys. Lett. 88(4), 043114 (2006). 36. C. A. Foss, G. L. Hornyak, J. A. Stockert, and C. R. Martin, “Template-synthesized nanoscopic gold particles: optical spectra and the effects of particle size and shape,” J. Phys. Chem. 98(11), 2963–2971 (1994). 37. M. Meier and A. Wokaun, “Enhanced fields on large metal particles: dynamic depolarization,” Opt. Lett. 8(11), 581–583 (1983). 38. A. Moroz, “Depolarization field of spheroidal particles,” J. Opt. Soc. Am. B 26(3), 517–527 (2009). 39. R. Esteban, R. Vogelgesang, J. Dorfmüller, A. Dmitriev, C. Rockstuhl, C. Etrich, and K. Kern, “Direct near-field optical imaging of higher order plasmonic resonances,” Nano Lett. 8(10), 3155–3159 (2008). 40. R. Boyack and E. C. Le Ru, “Investigation of particle shape and size effects in SERS using T-matrix calculations,” Phys. Chem. Chem. Phys. 11(34), 7398–7405 (2009). 41. I. O. Sosa, C. Noguez, and R. G. Barrera, “Optical properties of metal nanoparticles with arbitrary shapes,” J. Phys. Chem. B 107(26), 6269–6275 (2003). 42. C. E. Román-Velázquez, C. Noguez, and R. G. Barrera, “Substrate effects on the optical properties of spheroidal nanoparticles,” Phys. Rev. B 61(15), 10427–10436 (2000). 43. R. Lazzari, S. Roux, I. Simonsen, J. Jupille, D. Bedeaux, and J. Vlieger, “Multipolar plasmon resonances in supported silver particles: the case of Ag/α-Al2O3 (0001),” Phys. Rev. B 65(23), 235424 (2002). 44. Y. Ekinci, A. Christ, M. Agio, O. J. F. Martin, H. H. Solak, and J. F. Löffler, “Electric and magnetic resonances in arrays of coupled gold nanoparticle in-tandem pairs,” Opt. Express 16(17), 13287–13295 (2008).


Introduction
The optical properties of metallic nanoparticles have been of interest for many years and the scattering of electromagnetic waves by spheres and spheroids has been rigorously calculated before [1,2].One of the most interesting properties resulting from the interaction between electromagnetic waves and metallic nanoparticles is the excitation of collective oscillations of the conduction electrons, i.e. localized surface plasmon resonances (LSPRs).The LSPRs dependence on the size and shape of nanoparticles and the surrounding environment has attracted much attention due to potential applications in different devices [3][4][5][6][7].Of special interest are nanostructured arrays supported on a substrate as those fabricated by hole-mask colloidal lithography [8].This technique allows control of shapes and sizes in the fabrication of an ample variety of nanostructures including disks, ellipses, cones, binary nanoparticles [8] as well composite nanostructures like metal-dielectric-metal structures [9][10][11].Besides the LSPR observed in metallic nanostructures, the scientific literature reports on other interesting phenomena like optical magnetism in nanosandwiches [9][10][11], intrinsic Fano interference of localized plasmons in Pd nanoparticles [12], and plasmon induced magneto-optical activity [13].
Although there are many reports on the optical properties of metallic nanoparticles supported on a substrate, there are still some remaining gaps that deserve attention.First, for large nanoparticles (length scales of several tens or a few hundreds of nm) many of the studies are performed by only analyzing the extinction spectra obtained from transmittance measurements at normal incidence.In particular, the LSPRs of Au, Ag, Pt, and Pd nanodisks short-range ordered arrays have been previously studied by modeling the nanodisks as oblate spheroids in the modified long wavelength approximation [14][15][16].However, under normal incidence the analysis is limited to only the component (in-plane) of the polarizability tensor parallel to the substrate.Furthermore, the effect of the substrate is taken into account by assuming that the spheroids are embedded in an effective medium of refractive index given by the average between the ambient and substrate.On the other hand, both experimental and theoretical optical studies of periodic gold nanoparticle arrays (one-and two-dimensional) have been performed at oblique incidence [17][18][19][20][21]. Therefore, appropriate modeling of optical measurements at oblique incidence is required to retrieve the components of the polarizability tensor.Certainly, electrodynamics calculations are able to extract optical information but often with high computational costs [5,[17][18][19][20][21][22][23].
The simplest model to describe the optical properties of nanostructured arrays supported on a substrate is a homogeneous anisotropic layer.In this case, the components of the dielectric function tensor can be modeled with Lorentz harmonic-type oscillators or anisotropic effective medium theories [24][25][26].For small metallic nanoparticles (length scales < 10 nm) supported on a substrate, in the 2D-efective medium theory of Yamaguchi et al [25] and its improved version [26], the interaction between neighboring particles and with the substrate are taken into account in the dipole approximation.On the other hand, higher-order multipole interactions were extensively investigated by Bedeaux and Vlieger leading to the formulation of the island-film theory [27].This theory is based on the notion of excess quantities and surface susceptibilities providing modified Fresnel coefficients that have successfully explained the optical response of small metallic nanoparticles on a surface [28][29][30][31].At first-order two surface susceptibilities are required: one parallel (γ) and the other perpendicular (β) to the substrate.For non-interacting islands (low coverage), γ and β depend on geometrical details of the islands and the effect of the substrate is taken into account by the charge image technique.As noted above, the 2D-effective medium theory and the island-film theory have not been explored to describe the optical properties of large nanoparticles yet.
Spectroscopic ellipsometry is a well established optical technique for characterization of surfaces, interfaces, and thin films.Recently, its potential use at the nanometer scale has been reviewed [32].In this work, we have investigated the optical response of gold nanodisks arrays by spectroscopic ellipsometry and irradiance measurements at oblique and normal incidence.Thus, both, the in-plane and out-of-plane components are probed.The data are analyzed using two approaches: i) by modeling the nanostructure with a homogenous uniaxial layer with an effective dielectric function of either the Lorentz type (Sec.3.1) or with a modified Yamaguchi effective medium expression that includes dynamic depolarization factors (Sec.3.2).ii) In the second approach, the island-film theory is applied by modeling the in-plane polarizability in the modified long-wavelength approximation (Sec.

Experimental
Gold nanostructured arrays were fabricated on glass substrates by hole-mask colloidal lithography [8].The investigated samples are 40 nm thick Au nanodisks as well as 20 nm thick Au nanodisks covered with a 20 nm thick SiO 2 layer.In both cases, the nanodisks are nearly circular with an average diameter of 170 nm as represented in Figs.1(a) and 1(b).One of the Au nanodisks was capped with a SiO 2 layer in order to investigate the effect of partial embedding on the optical properties of metallic nanoparticles.Also, this structure is part of composite structures like Au-SiO 2 -Au nanosandwiches which have shown very interesting optical properties [9][10][11].Details of thickness and diameter control during fabrication of Au nanodisks can be found elsewhere [7,10,11,14].Spectroscopic ellipsometry (SE) measurements were performed in the 300-1300 nm wavelength range using a variable angle spectroscopic ellipsometer (VASE) of rotating analyzer type equipped with an auto-retarder (J. A. Woollam Co., Inc.).Reflectance spectra were measured for polarizations parallel and perpendicular to the plane of incidence.For measurements in reflection mode the back side of the substrate was roughened and the angles of incidence were between 50° and 70°.Also, transmission ellipsometry and irradiance transmittance (p-polarized) spectra at normal incidence were acquired.The refractive index of the glass substrate was determined from ellipsometric data and then its extinction coefficient was obtained by a point-by-point fitting of transmittance spectra at normal incidence.The complex dielectric function of gold was taken from VASE measurements of a homogeneous gold film [33].
SE measures the change in the polarization state that incident polarized light experiences due to interaction with a sample.In reflection (transmission) mode the change is described by the complex-valued ratio ρ r (ρ t ) between the reflection (transmission) coefficients for light polarized parallel r p (t p ) and perpendicular r s (t s ) to the plane of incidence [34] where ρ r (ρ t ) is expressed in terms of the two ellipsometric angles Ψ and Δ which depend on the microstructure and optical properties of the constituents of the sample.Thus, the construction of an optical model for the reflection (transmission) coefficients r p (t p ) and r s (t s ) allows the determination of parameters like film thicknesses, complex dielectric functions, volume fractions, etc., of a multilayer model representing the sample.Those parameters are obtained in a fitting procedure minimizing the mean square error (MSE) between the model and experimental data, where N is the number of measured Ψ-Δ pairs, M the total number of fit parameters, and σ exp are the standard deviations of the measurements.The superscripts mod and exp indicate model-calculated and experimental data, respectively.The non-linear regression algorithm also provides 90% confidence limits of the fitted parameters.ε a .As a first approach, the optical response of these nanostructured arrays will be described as a continuous layer of thickness d and with a dielectric tensor diag(ε x ,ε y ,ε z ), Fig. 1(c).The origin of this anisotropy is clear from the structure of the arrays because applying an electric field either parallel or perpendicular to the sample surface, the optical response of the nanostructured array will be different.Furthermore, the nearly circular shape of the nanodisks allows the definition of the in-plane (ε x = ε y = ε || ) and out-of -plane (ε z = ε  ) components of a homogeneous uniaxial layer.  .

Homogeneous uniaxial layer approach
where E is the photon energy, ε  is the high frequency dielectric constant; A, B, and E c are the oscillator amplitude, broadening, and central energy, respectively.The fitting was performed with the WVASE32 software (J. A. Woollam Co., Inc.).For the fitting, the layer thickness was fixed to d = 20 nm, corresponding to the thickness of the Au-nanodisks.Otherwise, leaving d as a fitting parameter, its correlation with A || , ε || , and ε  became very large.Thus, the fitting parameters are only those of Eq. ( 3).The obtained values and confidence limits as well as the MSE given by Eq. ( 2) are presented in Table 1.The components of the complex dielectric tensor ε L,|| and ε L, are shown in Fig. 2(c) where clearly a strong anisotropy of the nanostructured array is observed.It should be mentioned that if the total thickness of the nanostructure is considered in the fitting, i.e. d = d silica + d Au = 40 nm, the result is identical to that in Figs.2(a) and 2(b) but with parameters A || = 1.11, ε || = 0.74, A  = 0.03, and ε  = 0.94 whereas the broadening and central energy remained the same as in Table 1.Also, similar results were obtained considering either the weight thicknesses of about 2 nm or any other effective thickness ranging between 2 and 40 nm.Furthermore, it has been previously reported that for 10 nm thick and 110 nm in diameter gold nanodisks array the location of the plasmon resonance is not affected by capping with 20 nm of silica [10].
At this point, it seems natural to conclude that the optical response of the nanodisks array can be given in terms of the homogeneous uniaxial layer 20 nm thick with effective dielectric functions ε L,|| and ε L, of Fig. 2(c).However, in order to investigate whether or not the model truly describes the optical response of the nanostructured array, a comparison should be established with independent measurements as is done in Figs.2(d  that the experimental spectra of oblique p-polarized reflectance (R p ) are well described by the model-calculations.For s-polarization (R s ), presented in Fig. 2(e), the calculated data fairly describe the plasmon resonance but appreciably differences are observed at short wavelengths.
In addition, as the angle of incidence increases the experimental data of R s show a shoulder at about 640 nm and also noticed in the R p spectrum at 70°.The origin of this feature (Q) indicated by the arrows in Figs.2(d) and 2(e) is discussed in Sec.3.5.Another discrepancy is noticed in irradiance transmittance (T) at normal incidence (Fig. 2(f)).In this case, the plasmon resonance obtained from the model-calculations is broader and shifted to longer wavelengths compared to the experimental data.VASE data for the array with 40 nm thick gold nanodisks are shown in Figs.3(a) and 3(b).Once again, an excellent fit was obtained using the nanodisks thickness as the layer thickness.Fitting parameters and confidence limits are presented in Table 1 and the effective dielectric function tensor components obtained are shown in Fig. 3(c).As can be seen in Figs.3(d)-3(f), major discrepancies are obtained between the experimental and model-calculated spectra of R p , R s , and T for this thicker nanodisks array.In Figs.3(d) and 3(e) also is observed a feature at about 600 nm (Q) in the experimental R p spectrum at 70° and in the three R s spectra.In summary, the successful description of the SE data modeling the nanostructured arrays as a homogeneous uniaxial layer does not assure that the dielectric tensor of Lorentz-type obtained represents reliable materials properties.Therefore, other approaches should be used as will be shown below.
Fig. 3. Experimental and best fit ellipsometric spectra of Ψ (a) and Δ (b) for a 40 nm thick gold nanodisks array using the Lorentz oscillator dielectric response tensor components shown in (c).Comparison between experimental and model-calculated data for: (d) oblique reflectance of p-polarized and (e) s-polarized light, as well as (f) normal incidence transmittance.

Dynamic Yamaguchi approach
The homogenization procedure depicted in Fig. 1 within the Yamaguchi effective medium theory lead to effective dielectric function components ε Y,|| and ε Y, given by [25], and, where q is the volume fraction occupied by the nanodisks, F || and F  are effective depolarization factors, including the dipole interaction between particles and with the substrate.Equations ( 4) and ( 5) are applicable for very small particles compared with the wavelength.Since we are dealing with non-interacting nanodisks, F || and F  will depend only on their shape and interaction with the substrate.In their original work, Yamaguchi et al use a point dipole approximation to account for interaction of spherical (or near spherical) particles with the substrate [25].These assumptions impose some limitations to the model which was improved by Fedotov et al by considering the exact local fields of spheroidal particles [26].Furthermore, dynamic depolarization and radiation damping become important for the inplane component due to the relative large diameter of the nanodisks.The nanodisks are thus modeled as oblate spheroids with the major (minor) semi-axis a(c) parallel (perpendicular) to the substrate.We introduce the dynamic depolarization factor F ||,d for the in-plane component as, where l  = (1 + e 2 )(1-e arctan e 1 ) and l || = (1-l  )/2 are the depolarization factors in the electrostatic limit being e = ((a/c) 2 -1) -1/2 [27].The second term in Eqs. ( 6) and ( 7) accounts for the interaction with the substrate according to Eq. ( 28) in [26] and represents a modification of the Yamaguchi expression.The third and fourth terms of Eq. ( 6) correspond to dynamic depolarization and radiation damping, respectively, where k = 2π/λ is the wave number of the incident light with wavelength λ.In other works, these terms were included in a dynamic Maxwell-Garnett expression and were written as k 2 a 2 /3 and 2ik 3 a 3 /9, respectively, leading to only qualitative agreement with measurements [35,36].This is because those terms originally were derived from volume integrals on a sphere [37].In the case of spheroids those integrals lead to the first, third, and fourth terms in Eq. ( 6) up to a dynamic geometrical factor multiplying the dynamic depolarization (~k 2 ) term [38].However, it was shown that this kind of factor, which depends on the aspect ratio, produces a considerable red-shift larger than Tmatrix calculations [38] and herein is taken as the unity.Therefore, the dielectric tensor of non interacting nanodisks is represented by a dynamic modified Yamaguchi expression ε DMY,|| for the in-plane component, Eqs. ( 4) and ( 6), whereas a modified Yamaguchi expression ε MY, , Eqs. ( 5) and (7), is used for the out-of -plane component.
Figure 4 shows a comparison between the experimental and calculated spectra for the silica-capped Au nanodisks, schematically shown in Fig. 1(b).For that, Fresnel coefficients of an air-uniaxial film-substrate system [34] were used with dielectric functions components given by Eqs. ( 4)- (7).The calculations were performed with a FORTRAN code using the parameters: a = 85 nm, c = 10 nm, q = 0.065 and film thickness d = 20 nm.In Fig. 4 it can be observed that this model gives a better description of all measured spectra than by using Lorentz expressions shown in Fig. 2. The quantitative agreement seen in Fig. 4 supports the use of Eq. ( 6) for the description of the LSPR in spheroidal particles.In Fig. 4(f) the transmittance spectrum neglecting the nanodisks-substrate interaction (second term in Eq. ( 6)) is also included which clearly shows that the LSRP then becomes weaker and blue-shifted.This is explained when considering that according to Eq. (4) the resonance is located at wavelengths where Re{ε m } = (1-1/F || )ε a and as the interaction with the substrate decreases F || this condition is fulfilled at longer wavelengths.Capping the Au nanodisks with the SiO 2 layer has not effect on the LSPR because ε a = 1 in Eq. ( 4).Although a very good description of the optical spectra is seen n Fig. 4, some differences can be noticed in the spectral location of the extremes in R p and R s spectra, Figs.4(d Fig. 4. Experimental and calculated spectra for silica-capped gold nanodisks array 20 nm thick using the dynamic Yamaguchi effective dielectric response assuming gold spheroids: (a) Ψ, (b) Δ, (c) effective dielectric function tensor components, (d) reflectance of p-polarized and (e) spolarized light, as well as (f) normal incidence transmittance (subs).In (f) also the spectrum for spheroids non interacting with the substrate (free) is included.
Figure 5 shows a comparison of experimental and calculated spectra using Eqs.( 4)- (7) for the dielectric function of gold nanodisks array with thickness 40 nm.The parameters used in the calculations were a = 88 nm, c = 19 nm, q = 0.065 and film thickness d = 38 nm.Clearly in Figs. 4 and 5 the description of the experimental T, R p , and R s spectra is remarkably improved compared to the Lorentz approach of Figs. 2 and 3.In particular, the width and position of the plasmon resonance in the T spectrum are better described due to the sharper resonance of Im{ε DMY,|| } compared to Im{ε L,|| }.It can be noticed that the two components ε DMY,|| and ε MY, in Figs.4(c) and 5(c) are of the same order as those in Figs.2(c) and 3(c), the former show a different dispersion at short wavelengths which gives the correct spectral dependence of R p , and R s spectra.The latter feature is attributable to the onset of interband transitions in gold.These results indicate that the optical properties of these nanodisks arrays are dependent on their size, shape, properties of the constituent materials, and volume fraction.

Island film theory approach
Another formalism to be explored as a description of the optical response of the nanodisk arrays is the island-film theory.This theory is based on the concept of excess currents, charge densities and fields (E, D, B, and H) developed for the study of the optical properties of a film of islands attached to a flat substrate [27].The island-substrate interface is assumed to be invariant with respect to translation, rotation and reflection in the x-y plane.In this theory, the film is replaced by a dividing surface somewhere within the film separating two homogeneous media, the ambient and the substrate.Then, excess polarization (P s ) and magnetization (M s ) densities are located at that dividing surface in such a way that the electromagnetic fields outside of the original film remain the same.Those excess quantities modify the boundary conditions and for non-magnetic systems in absence of spatial dispersion, the constitutive relations at the dividing surface are found to be: P s || = γ E ||,+ , P s  = β D ,+ , and M s = 0, where the subscript "+" represents the average of the extrapolated values of the corresponding field at the dividing surface.The constitutive relations define the surface susceptibilities γ and β.Thus, for light incident from the ambient, modified Fresnel coefficients r p and r s are found and are given by [27,28], At normal incidence the transmission coefficient of the ambient-island film-substrate system can be calculated from [27,28], In Eqs. ( 8)- (11), ) are the refractive indexes of the ambient and substrate, respectively; θ i and θ t are the incidence and refraction angles, respectively, related by Snell's law n a sinθ i = n s sinθ t .From Eqs. ( 8)-( 11), the well known Fresnel coefficients for the ambient-substrate interface are obtained when γ = β = 0.The overall transmittance of light exiting from the back side of the substrate can be calculated according to, where t sa is to the transmission coefficient of the substrate-ambient back interface (this second reflection is neglected); α s and d s are the absorption coefficient and thickness of the substrate, respectively.
The next step requires the calculation of the two constitutive coefficients γ and β which for identical non-interacting particles (islands) are given by [27,28] where ρ is the number of nanoparticles per unit area; α || and α  are the in-plane and out-ofplane dipolar polarizabilities, respectively.These polarizabilities can be determined by multipole expansion of the potentials in the ambient and substrate and are expressed as a sum of the incident potential, the potential due to the charge distribution induced in the island, and the potential of the image charge distribution in the substrate.By applying the boundary conditions of those potentials at the surface of the substrate, an infinite set of linear equations for the expansion coefficients is obtained.In the practice, that set is numerically solved by neglecting multipoles larger than a suitable order.Also, analytical expressions for α || and α  have been reported for spheres and spheroids (oblate and prolate).Herein, the nanodisks are represented as oblate spheroids and the polarizabilities in Eq. ( 13) can be obtained from the multipole expansion in spheroidal coordinates.In the dipolar approximation α || and α  are given by [27], where a and c are the major and minor semi axes, respectively; ε m is the nanodisks dielectric function; L || and L  are the depolarization factors which for oblate spheroids touching the substrate are given by [27],   where l (||,) are the static depolarization factors and the nanodisks-substrate interaction is accounted through terms multiplied by the factor (ε a -ε s )(ε a + ε s ) 1 (<0 for an air-glass interface).It can be noticed that as a result of this interaction, 2L || + L  1.Also, the interaction with the substrate decreases the depolarization factors resulting in a shifted plasmon resonance to longer wavelengths.
This formalism has been used to analyze the optical response including the ellipsometric spectra of nanocolloidal gold particles on silicon [29], hexagonal arrays of silver nanocrystals on highly oriented pyrolytic graphite [30], and silver clusters supported on α-Al 2 O 3 substrates [31].In these works, the authors corrected the intraband contribution to the metal dielectric function (modifying the relaxation time in the Drude term) because the particle size was smaller than the bulk mean free path.
In the present case, the relative large diameter of the nanodisks makes it necessary to consider dynamic depolarization and radiation damping effects for the in-plane component α || which is renormalized according to [5,15,37], 23 *

11
, 46 and in Eq. ( 13) α || should be substituted by α || *.It can be noticed that in Eqs. ( 8)-( 17) the input parameters are only the lengths of major and minor semi axes a and c, respectively, and the number of nanodisks per unit area ρ.The extinction spectra of supported metallic nanodisks have been analyzed with Eq. ( 17) but an averaged refractive index was used for the surrounding medium [14,15].Figure 6 shows the calculated optical spectra with the island-film theory for the array of SiO 2 capped 20 nm thick gold nanodisks.The calculations were performed with a FORTRAN code with a = 89 nm, c = 9 nm, and ρ = 4.5 μm 2 .It can be noticed that the values of the semi axes a and c are comparable to the nominal radius (85 nm) and half thickness (10 nm), respectively.Also, the value of ρ is in agreement with the volume fraction q introduced in Sec.3.2 because they are related by q = ρ4πa 2 c/3d.As can be seen in Fig. 6, the calculated spectra reproduce qualitatively and quantitatively all main features in the experimental data.The constitutive coefficients γ and β show a strong anisotropy of the in-plane and out-of-plane surface susceptibilities, Fig. 6(c).It can be noticed that both, the dynamic Yamaguchi (Fig. 4) and the island film theory (Fig. 6) approaches, give very similar optical spectra.However, the dynamic Yamaguchi effective medium approach is preferred because the spheroid parameters used, a = 85 nm and c = 10 nm, agree with the nominal radius and half thickness of the nanodisks.In (f) also the spectrum neglecting the spheroids-substrate interaction (free) is included.

In-plane anisotropy
Figure 7 shows the transmission ellipsometry data at normal incidence of the gold-silica nanostructured arrays.The measurements were performed at two different orientations, one from an arbitrary orientation (A) and the other by rotating the sample 90° to an orientation A + 90°.In Fig. 7(a) around the plasmon resonance wavelength of 800 nm, Ψ shows small but significant deviations by less than ± 1° from the isotropic case (45°) whereas in Fig. 7(b) the deviations in Δ from 0° are about ± 2°.These deviations indicate some in-plane anisotropy which exact origin is unknown at present.One possible explanation is by considering that the hole-mask colloidal lithography process produces "nearly" circular shaped nanodisks.In first approximation they can be represented as elliptical nanodisks with in-plane semi axes a and b (see the insert in Fig. 7(a)).Within this approximation, for incident light polarized parallel to a (b) the optical response can be described by an in-plane effective dielectric function ε ||,a (ε ||,b ).
In the general case, when the nanodisks frame (a,b) is rotated by an angle  with respect to the laboratory frame (p,s), the ratio of the transmission coefficients t p and t s in Eq. ( 1) for p and s polarizations depends on the averaged dielectric functions given by ε p = ε ||,a cos for ellipsoidal particles will depend on the four terms considered in Eq. ( 6) as was done for spheroidal particles.In that case, the depolarization factors in the electrostatic limit l a and l b can be obtained from the integral form [1]. Nevertheless, the complete calculation of ε ||,a and ε ||,b is not possible because analytical expressions accounting for interaction with the substrate, dynamic depolarization, and radiation damping for ellipsoidal disks are not available.Another possible origin of the in-plane anisotropy could be due to the exact shape of the nanodisks.That is, depending on the procedure employed for fabrication of the nanodisks, either straight sidewalls or oblique truncation of the top as well as different nanoscopic roughness are obtained [11].Furthermore, given the high sensitivity of  shows an ellipsoidal nanodisk with semi axes (a,b) rotated an angle  from the (p,s) frame.

Quadrupolar coupling
As discussed in Sec.3.2 and 3.3, the main features in the experimental spectra can be quantitatively reproduced by modeling the nanodisks as oblate spheroids in the longwavelength limit.However, both the 2D Yamaguchi model and the island-film theory are unable to reproduce the features of the R s spectra and in the R p (70°) spectrum at about 640 nm in Figs. 2, 4 and 6 and at 600 nm in Figs. 3 and 5.That is because their origin is a quadrupolar effect which cannot be included in the present models.At normal incidence the phase of the incident wave does not change across the diameter of the nanodisks which are homogenously polarized and Eqs. ( 4) and ( 17) apply.On the other hand, increasing the angle of incidence there will be an appreciable phase difference of the incident wave on the two opposite extremes of the nanodisks, see the insert in Fig. 6(e).This phase retardation creates an inhomogeneous polarization field inside the nanodisks and quadrupolar moments can be coupled.It has been reported that for the extinction spectra measured at oblique incidence of 70° for s-polarized light transmittance, a second peak emerge at shorter wavelengths than the plasmon resonance in Au-nanodisks arrays 20 nm thick with diameter of 213 or 355 nm [39].For the larger nanodisks arrays the quadrupolar resonance was imaged using apertureless scanning near-field optical microscopy [39].These retarding effects in large particles coupling quadrupolar resonances have been reported using T-matrix calculations and finite-element modeling for silver spheroids and spheres, respectively [23,40].By using the discrete dipole approximation, the optical spectra of isolated nanoparticles with different shapes have been studied and a richer spectrum was found for particles with sharp edges [41].Also, the excitation of multipolar modes for small spheroids lying on a substrate has been studied using a spectral representation [42].Surface difference reflectivity spectroscopy of small silver nanoparticles have shown that besides particle and substrate interactions multipolar effects arises due to the particle truncation when it is brought in contact with the substrate [43].Therefore, the inhomogeneous incident field at larger angles of incidence and the sharp edges of the nanodisks real shape contribute to the quadrupolar feature.Finally, the results of this work will be valuable for a better understanding of the optical response of more complex nanostructures where nanodisks are the building blocks, for example as in nanosandwiches [9][10][11]44].

Conclusions
The effective optical response of gold-based nanostructured arrays produced by hole-mask colloidal lithography has been studied using two approaches.First, it was found that the ellipsometry spectra in the reflection mode can be described by modeling the nanostructured arrays as homogeneous uniaxial layers with dielectric functions of the Lorentz type.However, model-calculated spectra for oblique reflectance deviate from the experimental data at short wavelengths and the plasmon resonance is red-shifted and broader for normal incidence transmittance.Instead, modeling the nanodisks as oblate spheroids, the modified dynamic Yamaguchi effective medium expression qualitatively and quantitatively describes all measured spectra.The second approach considered was the island-film theory representing the nanodisks as oblate spheroids interacting with the substrate.Also, radiation damping and dynamic depolarization effects for the in-plane polarizability were considered.It was revealed that details of the nanostructured arrays are important to give a more complete description of all measured optical spectra.A small in-plane anisotropy detected by transmission ellipsometry spectra could be due to factors like deviation of circularity of the nanodisks, sidewalls shape, nanoscopic roughness as well as a weak nanodisks interaction.Furthermore, the complexity of the optical response increases because quadrupolar coupling arises at oblique incidence which cannot be included within the models here used.

Figure 1
Figure 1 shows a schematic representation of the investigated samples: (a) 40 nm thick gold nanodisks and (b) SiO 2 -capped 20 nm thick gold nanodisks.In both cases the nanodisks are deposited on a glass substrate with dielectric function ε s in an ambient with dielectric function

Figures 2 (
Figures 2(a) and 2(b) show the experimental VASE data and the best fit using the homogenous uniaxial layer model for the silica-capped Au nanodisks array.The in-plane ε || and out-of -plane ε  components are represented with a single Lorentz oscillator according to the expression, )-2(f).

Fig. 2 .
Fig. 2. Experimental and best fit ellipsometric spectra of Ψ (a) and Δ (b) for a silica capped gold nanodisks array with thickness 20 nm using the Lorentz oscillator dielectric tensor components shown in (c).Comparison between experimental and model-calculated data for: (d) oblique reflectance of p-polarized and (e) s-polarized light, as well as (f) normal incidence transmittance.

#Fig. 5 .
Fig. 5. Experimental and calculated spectra for a gold nanodisks array 40 nm thick using the dynamic Yamaguchi effective dielectric response assuming gold spheroids: (a) Ψ, (b) Δ, (c) effective dielectric function tensor components, (d) reflectance of p-polarized and (e) spolarized light, as well as (f) normal incidence transmittance.

Fig. 7 .
Fig. 7. Transmission ellipsometry spectra of Ψ (a) and Δ (b) at normal incidence for the silica capped Au nanodisks array 20 nm thick at two orientations differing by 90°.The insert in (a)
) and 4(e).The latter were obtained from derivative ; whereas for R s the values are 799 (806), 794 (806), and 792 (806) nm for the angles of incidence 50, 60 and 70°, respectively.This different response to p-or s-polarized light is indicative of some in-plane anisotropy as discussed below. 9) 2  + ε ||,b sin 2  and ε s = ε ||,a sin 2  + ε ||,b cos 2 , respectively.Following the procedure described in Sec.3.2, ε ||,a and ε ||,b transmission ellipsometry (as compared with irradiance transmittance measurements), it could be employed to investigate a likely weak nanodisks interaction by varying the interparticle spacing and analyzing the effect on the in-plane anisotropy.However, such comprehensive study is beyond of the scope of this work.