## Abstract

We studied the far-field optical response of supported gold-silica-gold nanosandwiches using spectroscopic ellipsometry, reflectance and transmittance measurements. Although transmittance data clearly shows that the gold nanodisks in the sandwich structure interact very weakly, oblique reflectance spectra of s- and p-polarized light show clearly asymmetric line-shapes of the Fano type. However, all experimental results are very well described by modeling the gold nanodisks as oblate spheroids and by employing a 2 × 2 scattering matrix formulation of the Fresnel coefficients provided by an island film theory. In particular, the Fano asymmetry can be explained in terms of interference between the scattered waves from the decoupled nanodisks in the spectral range limited by their respective plasmon resonances. We also show that the reflectance and ellipsometry spectra can be described by a three-layer system with uniaxial effective dielectric functions.

©2012 Optical Society of America

## 1. Introduction

Metal-dielectric composite nanostructures support confined electron oscillations driven by electromagnetic waves, i.e. localized surface plasmon resonances. Bringing together two or more metal nanoparticles, the plasmon resonances interact and coupled resonances emerge [1], as has been shown for several stacked configurations comprised of strips [2,3], cut-wires [4–6], and nanodisks [7–11] supported on a substrate. Several of these metal-dielectric-metal nanostructures show symmetric and asymmetric coupled resonances, where the latter has been interpreted as producing optical magnetism. In particular, the optical response of nanostructures comprised by metallic nanodisks in-tandem pairs (nanosandwiches) has been subject of several experimental and theoretical reports [12–15]. Besides of optical magnetism in nanosandwiches, Fano resonances have been reported in their normal incidence reflectance spectra and has been explained by using numerical simulations [7,8]. The Fano resonance is characterized by an asymmetric line-shape that results from the interference between a broad continuum spectrally overlapped by a narrower resonance. The Fano resonance phenomenon is commonly found in diverse metallic nanostructures [5–8,16–18] and the subject has been reviewed recently [19,20].

Generally, the optical response of plasmonic nanostructures is analyzed by numerical methods, often at high computational costs, and analytical descriptions involving the influence of electromagnetic interactions on the resonance line shape are scarce [21–24]. Furthermore, most of the previous reports are based on measurements of transmittance and/or reflectance at normal incidence. Owing to the intrinsic anisotropy of the nanostructures, such analysis is limited because only the in-plane component of the polarizability is probed. Therefore, appropriate modeling of optical measurements at oblique incidence is required to explain the observed phenomena identifying the individual contributions in composite nanostructures. Understanding the optical properties of weakly interacting nanodisks in a nanosandwich structure can serve as a first step in analysis of more complex anisotropic structures.

In this work, we investigate the optical response of gold-silica-gold nanosandwiches fabricated by hole-mask colloidal lithography, which produces short-range ordered structures covering large areas [25]. Spectroscopic ellipsometry and irradiance measurements at oblique and normal incidence are used to probe the in-plane as well as the out-of-plane components. Experimental evidence of weak coupling and interference effects are discussed in Sec. 3.1. The gold nanodisks in the sandwich structure are modeled as oblate spheroids in Sec. 3.2 and the polarizabilities in a modified long wavelength approximation [26] including interaction with the substrate are calculated. In Sec. 3.3 generalities of the island film theory [27] and the modified Fresnel coefficients in the Abelès 2 × 2 scattering matrix formalism are described. The latter formalism is applied in Sec. 3.4 to investigate interference effects and a comparison is made between model-calculated data and experimental spectra. Finally, in Sec. 3.5 the homogenization of the individual uniaxial layers is explored using effective dielectric functions of the Yamaguchi type [28].

## 2. Experimental

Arrays of gold-silica-gold nanosandwiches were fabricated on glass substrates by hole-mask colloidal lithography [25]. Details of thickness and diameter control during fabrication of nanosandwiches can be found elsewhere [10,11]. The fabrication process produces truncated cone-shaped structures schematically represented in the inset of Fig. 1 . In the present work the bottom and top gold nanodisks have average diameters of 170 and 120 nm, respectively. The two gold nanodisks are 20 nm thick and the thickness of the silica spacer is 40 nm.

Spectroscopic ellipsometry (SE) measurements were performed in the 300-1300 nm wavelength range at angles of incidence between 40 and 70° 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 at angles of incidence ranging from 20 to 70°. For ellipsometry and reflectance measurements the back-side of the substrate was roughened. Also, transmittance (s-polarized) spectra at normal incidence and 25° angle of incidence were acquired. The refractive index of the glass substrate was determined from ellipsometric data and the extinction coefficient was subsequently 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 [29].

## 3. Results and discussion

#### 3.1 Weakly coupled plasmon resonances

Figure 1(a) shows a comparison between transmittance (*T*) spectra measured at normal incidence, for the nanosandwich sample and a layer of 20 nm thick silica-capped gold nanodisks. We reported on the optical response of the latter sample in a previous work [30]. The structure and relevant parameters of the two samples are schematically represented in the insets of Fig. 1. It can be observed that the spectrum of the sandwiched structure consists of two plasmon resonances with nearly symmetric line-shapes; one weak resonance at 620 nm and a stronger mode at 810 nm. This result is in contrast to previously published results for nanosandwiches constructed from strongly coupled nanodisks. Such structures typically exhibit two hybridized resonances, that is, a strong high frequency mode and a weak low-frequency mode that correspond to the in-phase and anti-phase motion of the individual dipolar particle plasmons [7–11]. Furthermore, it is clear from Fig. 1(a) that the *T* spectra of the nanosandwiches and of the 170 nm diameter nanodisks are essentially identical for wavelengths longer than ~700 nm. As shown in Fig. 1(b), the s-polarized (*T _{s}*) transmittance spectra at an angle of incidence of 25° are nearly identical to the ones taken at normal incidence, in particular, the symmetric line-shape of the nanosandwich spectrum is retained. The results demonstrate that the plasmon resonances of the composite nanosandwich structure are well approximated by the dipole modes of the individual 120 nm and 170 nm nanodisks and, consequently, that plasmon hybridization is weak. This result is expected when considering that the 40 nm vertical separation between the disks is considerably larger than the 20 nm thickness of the individual nanodisks in the sandwich, as well as by considering the large difference in nanodisk aspect ratio, leading to nondegenerate modes.

Figure 1(c) compares reflectance spectra for the two samples discussed above. Since reflection measurements at near normal incidence are not possible in the experimental setup used, we instead show s-polarized reflectance (*R*_{s}) spectra at an angle of incidence of 20°. In the case of nanosandwich spectrum, one clearly observes asymmetric Fano-like line-shapes, in particular for the short-wavelength mode. A similar asymmetric line-shape was found for p-polarized reflectance (not shown). Notice that the maximum of *R _{s}* in Fig. 1(c) at 810 nm is almost identical to that of the silica-capped gold nanodisks, similar to the case of transmittance in Figs. 1(a) and 1(b).

#### 3.2 Nanodisk polarizabilities

According to the interpretation of the transmission spectra in Fig. 1, it is instructive to analyze the individual optical response of the decoupled nanodisks in the sandwiched structure. For that, and as in our previous work [30], the nanodisks are modeled as oblate spheroids with the major (minor) semi-axis *a _{j}* (

*c*) parallel (perpendicular) to the substrate. The sub-index

_{j}*j*= 1(2) identifies the nanodisk at the bottom (top) as shown in the insets of Figs. 2(b) and 2(c), in the latter

*s*is the thickness of the silica spacer. Thus, the spheroids optical response is expressed in terms of two polarizabilities: one in-plane

*α*

_{j},_{||}and one out-of-plane

*α*

_{j}_{,⊥}which are given by [27,30],

*ε*and

_{a}*ε*are the ambient and nanodisks dielectric functions, respectively. For the in-plane component, the effective depolarization factor in Eq. (1) is given by [30],

_{m}*l*

_{j}_{,||}is the depolarization factor in the electrostatic limit and

*L*

_{j}_{,}

*accounts for the interaction with the substrate using the spheroidal field [31], viz.,*

_{sub}*η*=

_{j}*a*/

_{j}*c*is the spheroid aspect ratio and

_{j}*d*is the distance between the spheroid and its image in the substrate with dielectric function

_{j}*ε*. The effective depolarization factor for the out-of-plane polarizability in Eq. (1) is

_{s}*L*

_{j}_{,⊥}=

*l*

_{j}_{,⊥}-2

*L*

_{j}_{,}

*where*

_{sub}*l*

_{j}_{,⊥}= 1-2

*l*

_{j}_{,||}. The expressions of

*l*

_{j}_{,||}and

*L*

_{j}_{,}

*in Eqs. (3) and (4) explicitly show their dependence on the spheroid aspect ratio*

_{sub}*η*=

_{j}*a*/

_{j}*c*. The third and fourth terms of Eq. (2) correspond to dynamic depolarization and radiation damping [30], respectively, where

_{j}*k*= 2

*π*/

*λ*is the wave number of the incident light.

Figure 2 shows the nanodisks polarizabilities, *α _{j},*

_{||}and

*α*

_{j}_{,⊥}, calculated with Eqs. (1)-(4), using the nominal parameters of the nanosandwich structure:

*a*

_{1}= 85 nm,

*a*

_{2}= 60 nm,

*c*

_{1}=

*c*

_{2}= 10 nm, and

*s*= 40 nm. It is clear from Fig. 2 that $\left|{\alpha}_{1,\left|\right|}\right|>\left|{\alpha}_{2,\left|\right|}\right|$ and $\left|{\alpha}_{1,\perp}\right|>\left|{\alpha}_{2,\perp}\right|$ which is mainly due to the larger volume of the bottom nanodisk. As was previously reported, the interaction with the substrate enhances (and shifts to longer wavelengths) the plasmon resonance of supported gold nanodisks [30]. In the present case, the distances

*d*between the nanodisks and their respective images are

_{j}*d*

_{1}= 2

*c*

_{1}and

*d*

_{2}= 2(

*c*

_{2}+

*s*+ 2

*c*

_{1}) for the bottom and top nanodisks, respectively. Since the interaction with the substrate decreases with

*d*according to Eq. (4),

_{j}*L*

_{1,}

*>*

_{sub}*L*

_{2,}

*, implying a larger decrease of the effective depolarization factor*

_{sub}*L*

_{1,||}in Eq. (2) for the nanodisk at the bottom. The latter result and the dependence on

*L*

_{1,||}in Eq. (1) show how the interaction with the substrate contributes to the result $\left|{\alpha}_{1,\left|\right|}\right|>\left|{\alpha}_{2,\left|\right|}\right|$.

In the case of the out-of-plane polarizabilities, an important contribution to the result $\left|{\alpha}_{1,\perp}\right|>\left|{\alpha}_{2,\perp}\right|$ resides on the fact that the nanodisk at the bottom is sandwiched by the silica spacer and the substrate (inset of Fig. 1). Thus, for p-polarized incident light, screening charges are induced in these materials due to the normal component of the incident electric field. Therefore, for this component the surrounding ambient differs from air and in Eq. (1), *ε _{a}* was taken as that of fused silica in the calculation of the out-of-plane polarizability

*α*

_{1,⊥}. It should be noted that the spectral position of the maxima in Im{

*α*

_{1},

_{||}} and Im{

*α*

_{2},

_{||}} at 810 and 620 nm, respectively, are coincident with the plasmon resonances in the

*T*spectrum in Fig. 1(a). In Fig. 2(c), the steep rise of Im{

*α*

_{j}_{,⊥}} at wavelengths shorter than 500 nm can be attributed to the onset of interband transitions in gold. Furthermore, in Fig. 2(b) it is important to note that

*α*

_{1},

_{||}and

*α*

_{2},

_{||}differ in phase in the wavelength range between 620 and 810 nm.

#### 3.3 Island film theory and 2 × 2 matrix Abelès formalism

The optical response of systems comprised by particles (islands) attached to a flat substrate has been successfully explained by the island film theory [27]. This theory is based on the concept of excess currents, charge densities and fields that modify the boundary conditions at the surface. For non-magnetic systems and in absence of spatial dispersion, the effect of those excess quantities is accounted for through constitutive relations in terms of two surface susceptibilities, one parallel (*γ*) and the other perpendicular (*β*) to the substrate. Thus, modified Fresnel coefficients are obtained, which for light polarized perpendicular *r*_{s} and parallel *r*_{p} to the plane of incidence are given by [27],

In Eqs. (5)-(7), *n _{a}* ( =

*ε*

_{a}^{1/2}) and

*n*( =

_{s}*ε*

_{s}^{1/2}) are the refractive indices of the ambient and substrate, respectively;

*θ*and

_{i}*θ*are the incidence and refraction angles, respectively, related by Snell’s law

_{t}*n*sin

_{a}*θ*=

_{i}*n*sin

_{s}*θ*. The transmission coefficients for s- and p-polarization are given by [27],

_{t}*D*and

_{s}*D*are the denominators in Eq. (5) and (6), respectively.

_{p}In the case of non-interacting particles the constitutive coefficients are *γ = ρα*_{||} and *β = ρα*_{⊥}/*ε _{a}*

^{2}, where

*ρ*is the number of particles per unit area. Using this formalism, we previously studied the optical response of supported gold nanodisks approximated as oblate spheroids, Fig. 3(a) [30].

Because the composite structure of the nanosandwiches, it is necessary to consider the Abelès 2 × 2 matrix formalism of the island film theory which leads to the interface matrix [32],

*r*and

_{f}*t*are the reflection and transmission coefficients of the hypothetical free-standing island film schematically represented in Fig. 3(b). Because the same medium is above and below the island film, the coefficients

_{f}*r*and

_{f}*t*for p- and s-polarizations in Eq. (10) are obtained from Eqs. (5)-(9) making

_{f}*n*=

_{a}*n*and

_{s}*θ*=

_{i}*θ*, that is,

_{t}The 2 × 2 matrix formalism allows straightforward calculations of the optical response of stratified systems. In fact, Eqs. (5)-(9) follow from the ratios $r={S}_{21}/{S}_{11}$ and $t=1/{S}_{11}$ where S* _{ij}* are elements of the scattering matrix $S={I}_{f}{I}_{as}$ (for the corresponding polarization) being

**I**

*the standard Abelès ambient-substrate interface matrix [33].*

_{as}In the present work Eqs. (11)-(13) are of particular interest because they allow for an analysis of the individual optical response of the free-standing island layers, that is, with polarizabilities (*α _{j,}*

_{||},

*α*

_{j,}_{⊥}) corresponding to the nanodisks at the bottom (

*j*= 1) and at the top (

*j*= 2) of the sandwich structure. For illustration purposes, we show in Fig. 4 only the normal incidence transmittance ${T}_{j}={\left|{t}_{f,j}\right|}^{2}$ and oblique (

*θ*= 50°) reflectance ${R}_{j}^{p,s}={\left|{r}_{f,j}^{p,s}\right|}^{2}$ spectra of the surface susceptibilities ${\gamma}_{j}=\rho {\alpha}_{j,\left|\right|}$ and ${\beta}_{j}=\rho {\alpha}_{j,\perp}/{\epsilon}_{a}^{2}$ with

_{i}*ρ*= 4.5 μm

^{−2}. Similar results were obtained for other angles of incidence. As expected, the

*T*,

_{j}*R*

_{j}^{p}and

*R*

_{j}^{s}spectra peak at the corresponding spectral position of the maximum in Im{

*α*

_{1},

_{||}} or Im{

*α*

_{2},

_{||}} shown in Fig. 2(a). Nevertheless, the individual spectra of

*T*,

_{j}*R*, and

_{j}^{p}*R*overlap to different extent in the range of 600 to 800 nm where the complex in-plane polarizabilities ${\alpha}_{j,\left|\right|}=\left|{\alpha}_{j,\left|\right|}\right|\mathrm{exp}\left(i{\phi}_{j}\right)$ in Fig. 2 present the largest phase difference Δ

_{j}^{s}*φ*=

*φ*

_{1}-

*φ*

_{2}as is shown in Fig. 4(d). Also, in Fig. 4(d) it is shown the spectral dependence of the phase difference of reflection coefficients ${r}_{f,j}^{p,s}=\left|{r}_{f,j}^{p,s}\right|\mathrm{exp}\left(i{\phi}_{j}^{p,s}\right)$ in Eq. (10), for p- and s-polarization ($\Delta {\phi}^{p,s}={\phi}_{1}^{p,s}-{\phi}_{2}^{p,s}$) calculated at an angle of incidence of 50°. Clearly, the phase differences show similar wavelength dependences. Results similar to those in Fig. 4(d) were obtained also for reflection spectra calculated for other angles of incidence but such large phase differences were not found for the transmission coefficients.

#### 3.4 Interference effects in free-standing and supported two-island layers

In this section, we first investigate interference effects in the hypothetical system shown in Fig. 3(c) which is comprised by two free-standing gold nanodisks layers, separated by a distance *h* and immersed in an ambient of refractive index *n _{a}*. The scattering matrix of this system is $S={I}_{f2}{L}_{h}{I}_{f1}$ where the top and bottom nanodisks layers are represented with interface matrices

**I**

_{f}_{2}and

**I**

_{f}_{1}of the type given in Eq. (10) and

**L**

*is the standard layer matrix ${L}_{h}=diag\left[\mathrm{exp}\left(-i\varphi \right),\mathrm{exp}\left(+i\varphi \right)\right]$ with $\varphi =kh{n}_{a}\mathrm{cos}{\theta}_{i}$ and*

_{h}*h*=

*s*+

*c*

_{1}+

*c*

_{2}. This phase shift is introduced by the electromagnetic waves traveling back and forth in the region bounded by the two interfaces of the nanosandwich structure as represented in Fig. 3(d). Thus, the reflectance of this system for p- and s-polarization ${R}^{p,s}={\left|{S}_{21}^{p,s}/{S}_{11}^{p,s}\right|}^{2}$ is given by,

*δ*is the phase of ${t}_{2}^{(p,s)2}-{r}_{2}^{(p,s)2}$. The cosine term in the numerator of Eq. (14) produces interference effects. The transmittance is obtained from ${T}_{d}^{p,s}={\left|1/{S}_{11}^{p,s}\right|}^{2}={T}_{1}^{p,s}{T}_{2}^{p,s}/{D}^{p,s}$ where

^{p,s}*D*is the denominator in Eq. (14) and no interference effects are expected in transmission spectra.

^{p,s}Figure 5(a)
shows the calculated spectrum at normal incidence of the two interfaces system and clearly *T _{d}*≈

*T*

_{1}

*T*

_{2}. However, in Figs. 5(b) and 5(c) the reflectance spectra

*R*and

^{p}*R*calculated with Eq. (14) for p- and s-polarization, respectively, exhibit strong destructive interference at λ = 660 nm coincident with what is experimentally observed in Fig. 1(c). In the case of p-polarization the reflectance decreases with increasing angle of incidence from 20 to 70°, as indicated by the arrow in Fig. 5(b), because the contribution of

^{s}*α*

_{j},_{||}decreases continuously with increasing angle of incidence. On the contrary,

*R*increases and broadens at grazing angles of incidence. The fact that different interference effects arose depending on the observation conditions (reflection or transmission) has previously been reported in other plasmonic systems [16]. The Fano-like interference in Figs. 5(b) and 5(c) results from the coherent addition of electromagnetic waves reflected from the nanodisks layers and determined by the phase Φ

^{s}*in Eq. (14) whose spectral dependence for p-polarization is shown in Fig. 5(d). Actually,*

^{p,s}*δ*is very small and the interference is dominated by the other terms in Φ

^{p,s}*. Similar results were obtained for s-polarization (not shown).*

^{p,s}The formalism in the last paragraph is readily extended for supported gold nanosandwiches. For this system the scattering matrix is $S={I}_{f2}{L}_{h}{I}_{f1}{L}_{a1}{I}_{as}$ where **I*** _{as}* is the normal Abelès ambient-substrate interface matrix [33]. Figure 6
shows a comparison between simulated and experimental spectra. As can be observed for

*R*

_{s}in Figs. 6(a) and 6(b) at

*θ*= 50° and 70°, respectively, all the features in the experimental spectra are reproduced by the simulation, including the destructive Fano interference at λ = 660 nm. A close inspection of the

_{i}*R*

_{s}spectra in Figs. 6(a) and 6(b) reveals that they are less broad than the calculated spectra shown in Fig. 5(c). This observation indicates that interference effects associated with reflections from the substrate are also present. On the other hand, it can be noticed that the experimental

*R*

_{p}spectra are well described by the model-calculated data, including the inverted line-shape at an angle of incidence of 70°, which is not seen in Fig. 5(c) because it is caused by the presence of the substrate. Furthermore, the model-calculated

*R*

_{p}spectrum in Fig. 6(b) shows weaker constructive interference and the minimum at about λ = 800 nm is slightly blue-shifted. It is clear that modeling the nanodisks as oblate spheroids represents some limitations. In particular, the sharpness of edges and interfaces in the real structure (which obviously have been neglected) strongly modify the screening effects mainly on the out-of-plane components (

*α*

_{j}_{,⊥}) controlling the p-reflectance and an increasing effect with the angle of incidence is expected.

Figure 6(c) shows the normal incidence overall transmittance of light exiting from the back side of the substrate calculated according to,

*t*is the transmission coefficient of the substrate-ambient back interface (the reflection from this interface is neglected);

_{sa}*α*and

_{s}*d*are the absorption coefficient and thickness of the substrate, respectively. Finally, Figs. 6(d) and 6(e) show the experimental and model-calculated ellipsometric spectra for Ψ and Δ which are defined by the complex-valued ratio between the reflection coefficients for light polarized parallel

_{s}*r*

_{p}and perpendicular

*r*

_{s}to the plane of incidence [33],

The importance of reporting ellipsometric data resides on the fact that they represent different illumination and measuring conditions, which are incidence of simultaneous p- and s-polarized waves and measurement of their intensity ratio and phase difference. As shown in Figs. 6(d) and 6(e) there is a very good agreement between experimental and model-calculated data sets at angles of incidence of 50 and 70°, thus supporting to the reliability of the model employed as well as the use of the nominal structural parameters.

#### 3.5 Homogenization of nanosandwiches layers

The results of the previous sections clearly show that homogenization of the nanosandwich layer is not possible, that is, the internal structure of the sample cannot be neglected. Nevertheless, previously we reported homogenization of supported nanodisks using the effective medium Yamaguchi expressions [30]. Therefore, this approach represents a homogenization alternative to be explored for the layers comprising the nanosandwich. In this case, the sample is represented as a three-layer system on glass which is depicted in the inset of Fig. 7(c)
. In this system, each layer has uniaxial anisotropy with effective in-plane (*ε _{j}*

_{,}

*) and out-of-plane (*

_{||}*ε*

_{j}_{,⊥}) components of the complex dielectric tensor of the Yamaguchi type [28,30],

*q*is the volume fraction occupied by the nanodisks with volume

_{j}*V*in layer

_{j}*j*( = 1,2) and polarizabilities

*α*

_{j}_{,(||,⊥)}of Eq. (1) and Fig. 2. Figure 7 shows the in-plane (

*ε*

_{j}_{,}

*) and out-of -plane (*

_{||}*ε*

_{j}_{,⊥}) components of the effective complex dielectric tensor for the bottom (

*j*= 1) and top (

*j*= 2) nanodisks layers. As expected, $\left|{\epsilon}_{1,\left|\right|}\right|>\left|{\epsilon}_{2,\left|\right|}\right|$ and $\left|{\epsilon}_{1,\perp}\right|>\left|{\epsilon}_{2,\perp}\right|$ because of the same reasons as discussed above in the case of the nanodisks polarizabilities in Fig. 2. For the spacer layer the polarizabilities in Eqs. (17) and (18) were calculated with

*ε*=

_{m}*ε*in Eq. (1) and geometric depolarization factors

_{silica}*l*

_{2}_{,||}and ${l}_{2,\perp}=1-2{l}_{2,\left|\right|}$using the nominal values of the structural parameters. The resulting tensor components are close to unity.

Figure 8
shows the experimental and calculated spectra of *R _{s}*,

*R*, Ψ, and Δ using the air-three uniaxial layers-glass system and the 4 × 4 matrix Berreman method [33]. The film thicknesses for the layers were the nominal nanodisks thicknesses namely, 20, 40, and 20 nm and the volume fractions

_{p}*q*were 0.075, 0.064, and 0.037 for layer 1, spacer layer and layer 2, respectively. As can be observed in Figs. 8(a) and 8(b), the calculated data reproduce the Fano line-shape of the

_{j}*R*and

_{s}*R*spectra quite well at the two angles of incidence

_{p}*θ*= 40 and 70°, respectively. The experimental and simulated spectra of the ellipsometric angles Ψ and Δ at the two angles of incidence are shown in Figs. 8(c) and 8(d), respectively. The homogenized layers model description of the experimental spectra in Fig. 8 is clearly similar to the one based on the island film theory in Fig. 6. In particular, both models give an excellent representation of both plasmon resonances and the experimentally observed Fano interference.

_{i}## 4. Conclusions

We have presented a theoretical analysis of experimental spectra of a layer of gold nanosandwiches comprised of two weakly coupled 20 nm thick gold nanodisks separated by a 40 nm silica spacer. It was found that modified Fresnel coefficients based on an island film theory and the Abelès 2 × 2 scattering matrix formalism accurately describe all experimental features observed, including ellipsometry data and spectral interference effects of the Fano type. The surface susceptibilities required by the island film theory were modeled by approximating the nanodisks as oblate spheroids while taking into account the interaction with the substrate. The asymmetric Fano line-shape observed in the reflectance spectra results from coherent superposition and interference of electromagnetic waves being reflected by the individual nanodisks layers and the substrate. Recently, pronounced Fano interference between localized plasmons and interface reflection was reported for layers of nanodisks made by hole-mask colloidal lithography [18]. However, in that case, the measurements were performed with light incident from the glass side. In particular, the strongest Fano effect was observed for angles close to the critical angle, that is, close to the transition from a spectral reflectance peak, similar to the spectral features discussed here, and a spectral dip observed under total internal reflection. In the present case, the Fano effect is instead dominated by interference between the light emerging from the two nanodisk layers comprising the sandwich structure. The strongest interference is found for wavelengths in between the plasmon resonances of the individual nanodisks, corresponding to out-of-phase oscillations of the decoupled nanodisks plasmons. As an alternative approach, we show that the measured spectra can also be well described through homogenization of the individual layers in the nanosandwiches array using dielectric functions of the Yamaguchi type.

## Acknowledgments

AMG acknowledges the support of Conacyt-Mexico and of the “Fondo Sectorial de Investigación para la Educación” grant No. 103385. The Knut and Alice Wallenberg foundation is acknowledged for support to instrumentation. MK and AD acknowledge support from the Swedish Foundation for Strategic Research through the grant RMA08 Functional Electromagnetic Metamaterials.

## References and links

**1. **N. J. Halas, S. Lal, W. S. Chang, S. Link, and P. Nordlander, “Plasmons in strongly coupled metallic nanostructures,” Chem. Rev. **111**(6), 3913–3961 (2011). [CrossRef] [PubMed]

**2. **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**(24), 3356–3358 (2005). [CrossRef] [PubMed]

**3. **H.-K. Yuan, U. K. Chettiar, W. Cai, A. V. Kildishev, A. Boltasseva, V. P. Drachev, and V. M. Shalaev, “A negative permeability material at red light,” Opt. Express **15**(3), 1076–1083 (2007). [CrossRef] [PubMed]

**4. **E. Pshenay-Severin, A. Chipouline, J. Petschulat, U. Hübner, A. Thünnerman, and T. Pertsch, “Optical properties of metamaterials based on asymmetric double-wire structures,” Opt. Express **19**(7), 6269–6283 (2011). [CrossRef] [PubMed]

**5. **A. Christ, Y. Ekinci, H. H. Solak, N. A. Gippius, S. G. Tikhodeev, and O. J. F. Martin, “Controlling the Fano interference in a plasmonic lattice,” Phys. Rev. B **76**(20), 201405 (2007). [CrossRef]

**6. **A. Christ, O. J. F. Martin, Y. Ekinci, N. A. Gippius, and S. G. Tikhodeev, “Symmetry breaking in a plasmonic metamaterial at optical wavelength,” Nano Lett. **8**(8), 2171–2175 (2008). [CrossRef] [PubMed]

**7. **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). [CrossRef] [PubMed]

**8. **Y. Jeyaram, S. K. Jha, M. Agio, J. F. Löffler, and Y. Ekinci, “Magnetic metamaterials in the blue range using aluminum nanostructures,” Opt. Lett. **35**(10), 1656–1658 (2010). [CrossRef] [PubMed]

**9. **Q. Li and G. P. Wang, “Tunable photonic metamaterials in the near infrared frequencies,” Opt. Express **18**(13), 14123–14128 (2010). [CrossRef] [PubMed]

**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). [CrossRef] [PubMed]

**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). [CrossRef]

**12. **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). [CrossRef] [PubMed]

**13. **C. Tserkezis, N. Papanikolaou, G. Gantzounis, and N. Stefanou, “Understanding artificial optical magnetism of periodic metal-dielectric-metal layered structures,” Phys. Rev. B **78**(16), 165114 (2008). [CrossRef]

**14. **L. Qiu, S. Wang, H. Liu, T. Li, S. Zhu, and X. Zhang, “Magnetic plasmon in coupled nanosandwich structure,” J. Opt. Soc. Am. B **28**(7), 1655–1660 (2011). [CrossRef]

**15. **J. Parsons, E. Hendry, J. R. Sambles, and W. L. Barnes, “Localized surface-plasmon resonances and negative refractive index in nanostructured electromagnetic metamaterials,” Phys. Rev. B **80**(24), 245117 (2009). [CrossRef]

**16. **B. Gallinet and O. J. F. Martin, “Influence of electromagnetic interactions on the line shape of plasmonic Fano resonances,” ACS Nano **5**(11), 8999–9008 (2011). [CrossRef] [PubMed]

**17. **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). [CrossRef] [PubMed]

**18. **M. Svedendahl and M. Käll, “Fano interference between localized plasmons and interface reflections,” ACS Nano **6**(8), 7533–7539 (2012). [CrossRef] [PubMed]

**19. **A. E. Miroshnichenko, S. Flach, and Y. S. Kivshar, “Fano resonances in nanoscale structures,” Rev. Mod. Phys. **82**(3), 2257–2298 (2010). [CrossRef]

**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). [CrossRef] [PubMed]

**21. **B. Gallinet and O. J. F. Martin, “Ab initio theory of Fano resonances in plasmonic nanostructures and metamaterials,” Phys. Rev. B **83**(23), 235427 (2011). [CrossRef]

**22. **N. A. Mirin, K. Bao, and P. Nordlander, “Fano resonances in plasmonic nanoparticle aggregates,” J. Phys. Chem. A **113**(16), 4028–4034 (2009). [CrossRef] [PubMed]

**23. **A. Pannipitiya, I. D. Rukhlenko, M. Premaratne, H. T. Hattori, and G. P. Agrawal, “Improved transmission model for metal-dielectric-metal plasmonic waveguides with stub structure,” Opt. Express **18**(6), 6191–6204 (2010). [CrossRef] [PubMed]

**24. **X. Piao, S. Yu, S. Koo, K. Lee, and N. Park, “Fano-type spectral asymmetry and its control for plasmonic metal-insulator-metal stub structures,” Opt. Express **19**(11), 10907–10912 (2011). [CrossRef] [PubMed]

**25. **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). [CrossRef]

**26. **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). [CrossRef]

**27. **D. Bedeaux and J. Vlieger, *Optical Properties of Surfaces* (Imperial College Press, 2001).

**28. **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). [CrossRef]

**29. **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). [CrossRef]

**30. **A. Mendoza-Galván, K. Järrendahl, A. Dmitriev, T. Pakizeh, M. Käll, and H. Arwin, “Optical response of supported gold nanodisks,” Opt. Express **19**(13), 12093–12107 (2011). [CrossRef] [PubMed]

**31. **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). [CrossRef]

**32. **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). [CrossRef]

**33. **R. M. A. Azzam and N. M. Bashara, *Ellipsometry and Polarized Light* (North Holland, , 1977).