## Abstract

The linear and nonlinear optical properties of a composite containing interacting spherical silver nanoparticles embedded in a dielectric host are studied as a function of interparticle separation using three dimensional frequency domain simulations. It is shown that for a fixed amount of metal, the effective third-order nonlinear susceptibility of the composite *χ*
^{(3)}(*ω*) can be significantly enhanced with respect to the linear optical properties, due to a combination of resonant surface plasmon excitation and local field redistribution. It is shown that this geometry-dependent susceptibility enhancement can lead to an improved figure of merit for nonlinear absorption. Enhancement factors for the nonlinear susceptibility of the composite are calculated, and the complex nature of the enhancement factors is discussed.

© 2009 OSA

## 1. Introduction

Recently a surge of interest has occurred in the theory, fabrication and characterization of metamaterials: micro- and nano-structured materials with properties significantly different from their constituents, for example negative index materials [1–3] and cloaking devices [4,5]. One specific thrust focuses on the design of highly nonlinear optical materials [6] based on nanocomposites. A possible approach to obtaining a large nonlinear optical response involves taking advantage of the strong local electric fields that can be achieved in and around plasmon resonant metal nanoparticles. Since the enhancement of the third-order nonlinear susceptibility scales with the fourth power of the electric field [7–9], these resonantly enhanced local fields can dramatically increase the nonlinear response of a composite compared to that of its constituent materials. Such plasmon enhanced metal-dielectric composites have many potential optical applications; for example for enhanced two-photon fluorescence [10–12] or as optical switches [13–15], and nonlinear optical absorption [16,17]. The latter application in particular requires materials that provide significant linear transmission, while exhibiting large optical absorption under high incident irradiance. Several studies have considered the nonlinear effective medium properties of such metal-dielectric composites either experimentally or theoretically. Reports indicate that while the addition of metal leads to significant linear absorption, the composite nonlinear response can be increased with respect to the linear absorption, for example, by increasing the volume filling fraction of the nanoparticles [18,19], by increasing host refractive index [20,21], and by increasing the aspect ratio of spheroidal metal particles [22]. However, no systematic studies appear to exist on the effect of interparticle interactions in metal dielectric composites.

In the present study we investigate the nonlinear optical properties of a nanocomposite consisting of spherical silver nanoparticles arranged in a regular lattice, with a lattice spacing that is sufficiently small to prevent diffractive effects. The linear and nonlinear optical properties of these composites are studied as a function of interparticle interaction, which is varied by changing the lattice spacing while maintaining a fixed metal volume fraction. The effective optical properties are determined using numerical simulations of the three dimensional electric field distribution in the composite under plane wave excitation. We show that changing the nanoparticle arrangement significantly increases the complex nonlinear refractive index with respect to the linear absorption. The feasibility of experiments demonstrating the presence of a surface plasmon enhanced host nonlinear optical response is discussed.

## 2. Theory

The composite linear dielectric function *ε _{c}* of an isotropic composite with a position dependent isotropic dielectric function $\epsilon \left(\omega ,\overrightarrow{r}\right)$ can be obtained based on a known linear electric field distribution of the form $\overrightarrow{E}\left(\omega ,\overrightarrow{r}\right){e}^{-i\omega t}$ within a volume V using the following relation [23]:

Here $\overrightarrow{E}\left(\omega ,\overrightarrow{r}\right)$ is the position dependent complex electric field in the composite, the notation ${\u3008\mathrm{...}\u3009}_{V}$ represents spatial averaging over volume *V*, and the notation ${\overrightarrow{E}}^{2}$ represents taking the dot product $\overrightarrow{E}\cdot \overrightarrow{E}$. Note that this dot product yields a complex number with a phase angle that is twice that of the original electric field vector, which will become important in this study. For simplicity of notation, frequency and position arguments are omitted in the remainder of this manuscript. The definition of the effective dielectric function given in Eq. (1) follows from the requirement of equal energy density in the effective medium and the composite. Note that identical results can be obtained by using the more common assumption of equal electric displacement in the effective medium and in the composite, as was confirmed in numerical evaluation of simulated electric field data using both approaches (data not shown). In the case of a dilute random distribution of isolated isotropic spherical nanoparticles in an isotropic host, Eq. (1) leads to the well-known Maxwell Garnett result [24]. Here we consider binary composites containing spherical inclusions of a material with an isotropic dielectric function *ε _{in}* embedded in a host material with an isotropic dielectric function

*ε*. We limit our study to composites in which the inclusions are arranged in a rectangular lattice, and consider optical excitation with electric fields aligned with one of the principal axes of the lattice. In this specific case, Eq. (1) provides only a single diagonal element of the dielectric tensor of the anisotropic effective medium.

_{h}Analogous to the approach shown for the composite dielectric function, the effective complex third order nonlinear susceptibility of an isotropic composite *χ _{c}*

^{(3)}(

*ω*) can be determined based on a known linear electric field distribution for the composite using the following relation [9,23,25]:

*χ*

_{c}^{(3)}(

*ω*) must be represented by a complex fourth rank tensor. In this study we derive the effective

*χ*

_{c}^{(3)}only for excitation along one of the principal axes of the simulated structure, here the x-direction. In this specific case Eq. (2) provides only the tensor component

*χ*

_{c,xxxx}^{(3)}, which in the following will be denoted

*χ*

_{c}^{(3)}. Given this simplified notation, it is important to keep in mind that

*the calculated susceptibilities do not represent an effective isotropic nonlinear susceptibility*, but instead represent the calculated

*χ*

_{c}^{(3)}for excitation with a specific electric field direction with respect to a highly symmetric structure.

For a binary composite in which the inclusion and host materials have an isotropic Kerr-type third order nonlinear susceptibility *χ _{in}*

^{(3)}and

*χ*

_{h}^{(3)}respectively, it can easily be shown that Eq. (2) can be written in the following form:

*f*is the volume fraction of the inclusion in the composite,

_{in}*f*is the volume fraction of the host material given by

_{h}*f*= 1 -

_{h}*f*, and the factors

_{in}*g*represent susceptibility enhancement factors that satisfy the following relation:

_{j}^{(3)}*j*represents either the inclusion (‘in’) or the host (‘h’). It is important to note that the numerator is averaged over a limited volume

*V*indicating the volume of the inclusion or the host respectively. From Eq. (3) the factors

_{j}*g*can be seen to represent the enhancement of the third order susceptibility contribution from host and inclusion, relative to the value expected based on a homogeneous electric field distribution throughout the volume. Note that another common definition of the enhancement factor considers the volume fraction

_{j}^{(3)}*f*as part of the enhancement factor, see e.g. Ref [8]. Although the latter approach certainly provides valid results, we favor the definition shown in Eq. (3) as it more clearly highlights the contribution of the electric field redistribution to the

*χ*enhancement, even for components with a small fill-fraction.

^{(3)}Based on the calculated composite nonlinear susceptibility *χ _{c}*

^{(3)}one can obtain the composite nonlinear refractive index

*η*=

_{2,c}*n*+

_{2,c}*iκ*where

_{2,c}*n*and

_{2,c}*κ*represent the real and imaginary parts of the nonlinear refractive index. In highly absorptive composites the conversion from

_{2,c}*χ*

_{c}^{(3)}to the nonlinear refractive index

*η*must take into account linear loss in the composite [26] which has been shown to lead to the relation ${\eta}_{2,c}=(3/(4{\epsilon}_{0}c|{\eta}_{c}{|}^{2}))(1-i({\kappa}_{c}/{n}_{c})){\chi}_{c}^{(3)}$ where

_{2}*ε*is the permittivity of vacuum,

_{0}*c*is the speed of light in vacuum,

*η*is the complex linear refractive index of the composite, and

_{c}*n*and

_{c}*κ*are the real and imaginary parts of

_{c}*η*. The thus obtained nonlinear refractive index can be used to determine the nonlinear absorption coefficient of the composite

_{c}*β*according to

_{c}*β*= (4π/

_{c}*λ*)

*κ*with λ the free space wavelength. Finally, using the calculated linear and nonlinear optical properties one can obtain the frequency dependent figure of merit for nonlinear absorption given by

_{2,c}*β*/

*α*. This figure of merit reflects the fact that high nonlinear absorption requires both a high

*β*value as well as a low linear absorption to enable a long interaction length. This figure of merit is especially relevant in metal dielectric composites, where plasmon enhanced linear absorption will significantly affect the figure of merit.

## 3. Simulation geometry

For all simulations, we consider arrays of spherical silver nanoparticles embedded in a host with a frequency independent real refractive index of 1.5. This index is chosen as a typical value of commonly used host materials at visible frequencies, similar to that of for example SiO_{2} (n = 1.45), many organic polymers (n ~1.4-1.6), and soda lime glass (n = 1.52). For the linear silver dielectric properties we use a surface scattering corrected Drude model fit of available literature values [27] given by *ε _{in} = ε_{∞} - ω_{p}^{2}*/(

*ω*(

^{2}+ iω*Γ*)) where

_{0}+ Γ_{s}*ε*= 5.451,

_{∞}*ω*= 1.474 × 10

_{p}^{16}rad/s, and the bulk electron scattering rate

*Γ*= 8.354 × 10

_{0}^{13}s

^{−1}. The surface scattering rate

*Γ*is given by

_{s}*Γ*/

_{s}= Aυ_{f}*r*= 2.8 × 10

^{14}s

^{−1}, where

*A*is set to 1,

*υ*= 1.39 × 10

_{f}^{6}m/s is the Fermi velocity in silver, and

*r =*5 nm is the radius of the particle. This particle radius allows the use of sufficiently small interparticle separations to consider the composite as an effective medium in the wavelength range of interest. Note that this choice makes surface scattering the dominant contribution to the total electron damping.

Three-dimensional frequency domain electromagnetic simulations were performed using Microwave Studio [28]. For all simulations, the incident electric field is polarized along the *x* direction and the light propagates along the *z* direction. Periodic arrays of nanoparticles are simulated by considering a rectangular unit cell that contains a single 10 nm diameter Ag sphere centered in the simulation volume. The length (*L _{z}*) of the unit cell is held fixed at 30 nm, while the height (

*L*) and width (

_{x}*L*) are both varied while maintaining a constant unit cell volume of 30 nm × 24 nm × 24 nm. These dimensions correspond to a fixed inclusion volume fill fraction

_{y}*f*= 0.03. An infinite array of nanoparticles is simulated by taking advantage of symmetry planes, setting the tangential electric field components to zero in the

_{in}*y*-

*z*plane at x = 0 and at x =

*L*/2, and setting tangential magnetic field components to zero in the

_{x}*x*-

*z*plane at y = 0 and at y =

*L*/2. In this way only a quarter of the unit cell needs to be considered in the simulations, reducing the required simulation time. Open boundary conditions are used along the z-direction. At the chosen length

_{y}*L*of the simulation volume all particle-related near fields were found to be contained within the unit cell. Tetrahedral meshing is employed to simulate physically realistic electric field distributions that are virtually free of grid artifacts while using a reasonable number of mesh cells. The minimum cell size used was smaller than 1 nm on a side. The simulated linear electric field distributions were analyzed using Eqs. (1) and 2 to compute the effective linear and nonlinear composite properties for x-polarized excitation. To study the effect of interparticle interactions on the optical properties, simulations were performed for five different geometries with dimensions

_{z}*L*×

_{x}*L*of 12.5 × 46 nm

_{y}^{2}, 15 × 38.4 nm

^{2}, 24 × 24 nm

^{2}, 38.4 × 15 nm

^{2}, and 46 × 12.5 nm

^{2}, corresponding to an increasing edge-to-edge interparticle separation along the x-direction of 2.5 nm, 5 nm, 14 nm, 28.4 nm, and 36 nm. These choices of longitudinal spacing

*L*(along the incident electric field) and transverse spacing

_{x}*L*(normal to the incident electric field) lead to significant changes of the interparticle interaction while maintaining a constant metal fill fraction. The selected interparticle separations are all significantly smaller than half the shortest optical wavelength considered, and consequently diffractive effects do not affect the optical properties of these composites.

_{y}## 4. Results and discussion

Figure 1(a)
shows the computed electric field distribution *E _{x}* in the x-y plane for

*L*= 24 nm (square arrangement of nanoparticles) under excitation with x-polarized light with a field strength of 1 V/m at a frequency close to the nanoparticle surface plasmon resonance (

_{x}*ω*= 4.64 × 10

^{15}rad/s). The vertical scale indicates the enhancement of the lateral electric field strength (

*E*) relative to the incident field, and the field is shown at the phase angle for which maximum field enhancement is observed. Contour lines are spaced 2.5 V/m apart. The field distribution around the 10 nm diameter particle appears dipolar, and shows significant field enhancement (~15 × ) near the nanoparticle surface. Figure 1(b) shows the corresponding

_{x}*E*distribution for a longitudinal spacing

_{x}*L*= 12.5 nm (2.5 nm edge-to-edge spacing) excited near the collective plasmon resonance of these interacting particles (

_{x}*ω*= 4.0 × 10

^{15}rad/s). The maximum field enhancement is seen to be significantly increased to ~45 × at this reduced longitudinal interparticle spacing due to near-field interaction between adjacent particles. Consequently, an enhancement of the nonlinear optical response is expected in the regions of high field for such closely spaced geometries. Note that the geometry-related increase in peak electric field amplitude occurs predominantly in the small volume between the particles.

Figure 2
shows the linear absorption coefficient as calculated using Eq. (1) based on the simulated three-dimensional electric field data, evaluated for several frequencies of the incident plane wave. Each of the curves in Fig. 2 corresponds to one of the five interparticle spacings considered. The locations of the symbols indicate the frequency samples at which the field distribution and the corresponding optical properties were evaluated. For comparison, the dashed line labeled ‘MG limit’ shows the analytical result obtained using Maxwell Garnett (MG) theory for this fill fraction. The observed absorption peaks are due to the resonant excitation of approximately dipolar plasmon modes on the metal nanoparticles, leading to resonantly enhanced energy dissipation. As the longitudinal interparticle separation is reduced, a red-shift is observed in the location of the plasmon resonance compared to the Maxwell Garnett result. Conversely, at large longitudinal interparticle spacing (and small transverse interparticle spacing) a blue-shift is observed. These are well known effects that can be understood in terms of near-field interactions between neighboring metal nanoparticles [29]. Consequently these resonance must be interpreted as *collective* plasmon resonances of the lattice, rather than as individual nanoparticle plasmon resonances. The maximum linear absorption coefficient is seen to increase slightly as the longitudinal interparticle spacing is reduced. Finally, the analytical Maxwell Garnett result is seen to lie close to the numerically computed absorption coefficient for *L _{x}* = 24 nm (square arrangement), indicating that interparticle interactions are minimal in this specific geometry.

Figure 3
shows the complex geometry-dependent nonlinear susceptibility enhancement factor *g _{h}*

^{(3)}calculated based on the simulated three-dimensional electric field data using Eq. (4), for different frequencies of the incident plane wave. Figure 3(a) shows the phase angle of the complex enhancement factor, and Fig. 3(b) shows the magnitude of the enhancement on a logarithmic scale. The dashed curve shows the corresponding results calculated using an analytic expression for the third order nonlinear susceptibility of composites consisting of non-interacting particles [8]. The inset in Fig. 3(b) shows the real and imaginary part of the enhancement factor on a linear scale for the non-interacting particles. Note that the analytical curve closely resembles our numerically obtained results for

*L*= 24 nm in terms of magnitude, shape, and resonance frequency, showing that inter-particle interactions only weakly affect the nonlinear response in this particular geometry. For all geometries a significant nonlinear susceptibility enhancement (

_{x}*g*

^{(3)}> 1) is observed across a large frequency bandwidth near the plasmon resonance, with maximum enhancement occurring near the surface plasmon resonance of the structure. As the longitudinal interparticle spacing is reduced, the magnitude of the enhancement is seen to increase by more than an order of magnitude. Apparently the field enhancement obtained at small interparticle spacing results in a net increase of

*g*

^{(3)}, despite the fact that the enhancement occurs only within a small volume. While interparticle interaction does contribute to the enhancement, the field enhancement obtained due to the plasmon resonance provides the main contribution to the observed enhancement factors. The resonant nature of the enhancement has an important consequence: as can be seen in Fig. 3(a), the phase angle of the enhancement factor increases from approximately 0° to 360° as the frequency is increased from below the plasmon resonance frequency to well above the resonance frequency. This behavior can be understood by realizing that near resonance the local electric fields inside and just outside the nanoparticles occur with respectively a −90° and + 90° phase delay with respect to the incident field. As can be seen in Eq. (4), the phase of the numerator of the enhancement factor scales with E

^{2}, and consequently the 90° phase difference of the linear fields on resonance lead to a 180° phase difference of the numerator in Eq. (4) compared to the phase of the average field. The complex nature of the enhancement indicates that in the presence of metal nanoparticles, a Kerr-type positive nonlinear refractive host can act as a nonlinear absorber, a negative nonlinear refractive medium, or a saturable absorber, depending on the frequency used [16,17]. Note that in experiments on metal-dielectric composites these effects may be overshadowed by nonlinearities introduced by the metal nanoparticles.

Figure 4
shows the phase (Fig. 4(a)) and magnitude (Fig. 4(b)) of the geometry-dependent nonlinear susceptibility enhancement factor *g _{in}*

^{(3)}of the inclusion, based on the same simulated three-dimensional electric field data used to generate Fig. 3. Note that while the phase dependence on frequency appears virtually identical to that observed in Fig. 3(a), the curves shown in Fig. 4(a) are in fact independently calculated values based on the

*internal*electric fields, as opposed to the external electric fields that were used to generate Fig. 3(a). The dashed curves represent the corresponding analytical result for this fill fraction. This analytical result relies on a different equation than the one used in Fig. 3, since it also considers the internal fields [8]. In the limit of weak interparticle interaction (

*L*= 24 nm), the magnitude of the enhancement

_{x}*g*

_{in}^{(3)}is seen to be significantly larger than that of

*g*

_{h}^{(3)}observed in Fig. 3. This is due to fact that the internal field enhancement that contributes to

*g*

_{in}^{(3)}occurs throughout the entire volume of the particle, while the external field enhancement that contributes to

*g*

_{h}^{(3)}occurs only in a small fraction of the host volume. In stark contrast with the observations made in Fig. 3, the magnitude of the nonlinear susceptibility enhancement of the inclusion is found to be nearly independent of interparticle separation. This very different behavior is due to the fact that the internal field distribution is relatively unaffected by changes in the interparticle separation.

Using the calculated linear effective dielectric function and nonlinear susceptibility enhancement factors *g _{j}*

^{(3)}it is now possible to evaluate the influence of the geometry-dependent field enhancement on the figure of merit for nonlinear absorption. The figure of merit is given by

*β*/

_{c}*α*, where

_{c}*β*is the nonlinear absorption coefficient of the composite. The figure of merit indicates the inverse of the irradiance required to achieve a nonlinear absorption coefficient equal to the linear absorption coefficient. For the following analysis the silver nonlinear response is approximated by a Kerr-type response with

_{c}*χ*

_{in}^{(3)}=

*i**10

^{−10}esu (1.75

*i*× 10

^{−17}V

^{2}/m

^{2}), while the host is assumed to be nonlinearly refractive with

*χ*

_{h}^{(3)}= 10

^{−14}esu (1.75 × 10

^{−21}V

^{2}/m

^{2}). For clarity of presentation these values are assumed to be frequency independent near the nanoparticle resonance. Note that while metal nonlinear optical properties are often described in terms of an effective Kerr nonlinearity, the underlying physical mechanisms include non-Kerr type effects such as Fermi smearing and thermal nonlinearities. Consequently, the metal contribution to the composite nonlinearity in real-world experiments may not reproduce the exact functional form derived here, however the predicted trends are expected to be observable.

The total figure of merit for nonlinear absorption can be shown to be separable into a host contribution and an inclusion contribution. Figure 5(a)
shows the host contribution to the figure of merit of the composite as a function of geometry. Despite the fact that a nonlinear *refractive* host is assumed, a positive figure of merit for nonlinear *absorption* is observed across a significant frequency range due to the complex nature of the susceptibility enhancement factor near the plasmon resonance of the structure. As the longitudinal interparticle spacing is decreased from 46 nm to 12.5 nm, the peak figure of merit is seen to increase by more than an order of magnitude. This increase follows from the observed large increase in the effective nonlinear susceptibility (Fig. 3) due to near-field coupling between adjacent particles accompanied by additional field enhancement and confinement (Fig. 1), combined with the relatively small change in peak absorption coefficient as the interparticle spacing is modified (Fig. 2). These results demonstrate that while the presence of the metal nanoparticles does introduce significant absorption, the nonlinear optical absorption performance of a thin metal-dielectric composite can be improved significantly by modifying the spatial distribution of the metal. Figure 5(b) shows the contribution of the inclusion to the overall figure of merit of the composite as for the same geometries. The entirely imaginary *χ _{in}*

^{(3)}assumed for the metal is seen to lead to a

*negative*figure of merit for nonlinear absorption, indicative of a composite that exhibits saturable absorption. Experimental studies of metal-dielectric composites indeed show saturable absorption near the plasmon resonance [17,30]. As the longitudinal interparticle spacing is decreased, the magnitude of the figure of merit is seen to decrease slightly due to the observed weak increase in linear absorption as the longitudinal interparticle spacing is decreased (Fig. 2), and partly due to an additional frequency dependence introduced by the factor 1/

*λ*in the conversion from

*χ*

^{(3)}to

*β*.

It is important to note that the specific example shown in Fig. 5 assumes a weak nonlinear contribution, using equations derived to first order in nonlinear polarization. This assumption will cease to be valid in cases where a significant change in total absorption is achieved relative to the already present metal-induced absorption. Consequently the complete evaluation of plasmon-enhanced nonlinear absorbers will require calculations that go beyond the first order perturbation approach. Secondly, the high irradiances required to achieve measurable changes in the absorption will induce significant heating and possibly melting and sintering of metal nanoparticles, adding an experimental challenge to studies of these effects, particularly for relatively slow measurements involving nanosecond pulses. Note also that for the values of *χ _{in}*

^{(3)}and

*χ*

_{h}^{(3)}chosen in this example the metal nonlinear response provides by far the dominant contribution to the figure of merit. While the values used here were chosen mainly for illustrative purposes, their magnitudes are not unrealistic, and consequently the experimental observation of plasmon enhancement of nonlinear refractive or absorptive host properties are expected to be challenging. A potentially mitigating factor is that our choice of nanoparticle diameter as discussed previously leads to a significant predicted contribution of surface scattering to the total electronic damping. Consequently experimental studies using larger particles may produce a significantly stronger nonlinear optical response due to a reduction in surface scattering. Further optimization of particle arrangement and particle shape may lead to an experimentally observable enhancement of host-related nonlinear optical refractive and absorptive effects in real metal-dielectric composite materials.

## 5. Summary and conclusions

The effect of interparticle spacing on the linear and nonlinear optical properties of periodic metal-dielectric nanocomposites was discussed. Under the assumption that the materials exhibit a third-order Kerr type nonlinearity, it is shown that a reduced interparticle spacing along the incident field direction can lead to a significantly increased composite nonlinear optical response at frequencies near the plasmon resonance, while leaving the maximum linear optical absorption largely unaffected. These two findings lead to the observation of an enhanced figure of merit for nonlinear absorption by a nonlinear refractive host.

## Acknowledgements

We would like to thank Prof. David G. Stroud for helpful discussions. This material is based upon work supported by the U. S. Army Research Office under contract/grant number 50372-CH-MUR.

## References and links

**1. **C. M. Soukoulis, S. Linden, and M. Wegener, “Physics. Negative refractive index at optical wavelengths,” Science **315**(5808), 47–49 (2007). [CrossRef]

**2. **R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science **292**(5514), 77–79 (2001). [CrossRef]

**3. **H. J. Lezec, J. A. Dionne, and H. A. Atwater, “Negative refraction at visible frequencies,” Science **316**(5823), 430–432 (2007). [CrossRef]

**4. **J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science **312**(5781), 1780–1782 (2006). [CrossRef]

**5. **D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science **314**(5801), 977–980 (2006). [CrossRef]

**6. **K. Tsuchiya, S. Nagayasu, S. Okamoto, T. Hayakawa, T. Hihara, K. Yamamoto, I. Takumi, S. Hara, H. Hasegawa, S. Akasaka, and N. Kosikawa, “Nonlinear optical properties of gold nanoparticles selectively introduced into the periodic microdomains of block copolymers,” Opt. Express **16**(8), 5362–5371 (2008). [CrossRef]

**7. **J. E. Sipe and R. W. Boyd, “Nonlinear susceptibility of composite optical materials in the Maxwell Garnett model,” Phys. Rev. A **46**(3), 1614–1629 (1992). [CrossRef]

**8. **M. I. Stockman, K. B. Kurlayev, and T. F. George, “Linear and nonlinear optical susceptibilities of Maxwell Garnett composites: Dipolar spectral theory,” Phys. Rev. B **60**(24), 17071–17083 (1999). [CrossRef]

**9. **D. Stroud and P. M. Hui, “Nonlinear susceptibilities of granular matter,” Phys. Rev. B **37**(15), 8719–8724 (1988). [CrossRef]

**10. **H. Ditlbacher, N. Felidj, J. R. Krenn, B. Lamprecht, A. Leitner, and F. R. Aussenegg, “Electromagnetic interaction of fluorophores with designed two-dimensional silver nanoparticle arrays,” Appl. Phys. B **73**(4), 373–377 (2001). [CrossRef]

**11. **A. M. Glass, A. Wokaun, J. P. Heritage, J. G. Bergman, P. F. Liao, and D. H. Olson, “Enhanced two-photon fluorescence of molecules adsorbed on silver particle films,” Phys. Rev. B **24**(8), 4906–4909 (1981). [CrossRef]

**12. **W. Wenseleers, F. Stellacci, T. Meyer-Friedrichsen, T. Mangel, C. A. Bauer, S. J. K. Pond, S. R. Marder, and J. W. Perry, “Five Orders-of-Magnitude Enhancement of Two-Photon Absorption for Dyes on Silver Nanoparticle Fractal Clusters,” J. Phys. Chem. B **106**(27), 6853–6863 (2002). [CrossRef]

**13. **X. Hu, P. Jiang, C. Ding, H. Yang, and Q. Gong, “Picosecond and low-power all-optical switching based on an organic photonic-bandgap microcavity,” Nat. Photonics **2**(3), 185–189 (2008). [CrossRef]

**14. **R. Katouf, T. Komikado, M. Itoh, T. Yatagai, and S. Umegaki, “Ultra-fast optical switches using 1D polymeric photonic crystals,” Photonics Nanostruct. Fundam. Appl. **3**(2-3), 116–119 (2005). [CrossRef]

**15. **M. Scalora, J. P. Dowling, C. M. Bowden, and M. J. Bloemer, “Optical limiting and switching of ultrashort pulses in nonlinear photonic band gap materials,” Phys. Rev. Lett. **73**(10), 1368–1371 (1994). [CrossRef]

**16. **D. C. Kohlgraf-Owens and P. G. Kik, “Numerical study of surface plasmon enhanced nonlinear absorption and refraction,” Opt. Express **16**(14), 10823–10834 (2008). [CrossRef]

**17. **G. Piredda, D. D. Smith, B. Wendling, and R. W. Boyd, “Nonlinear optical properties of a gold-silica composite with high gold fill fraction and the sign change of its nonlinear absorption coefficient,” J. Opt. Soc. Am. B **25**(6), 945–950 (2008). [CrossRef]

**18. **H. B. Liao, W. Wen, and G. K. L. Wong, “Preparation and characterization of Au/SiO_{2} multilayer composite films with nonspherical Au particles,” Appl. Phys., A Mater. Sci. Process. **80**(4), 861–864 (2005). [CrossRef]

**19. **N. Pinçon, B. Palpant, D. Prot, E. Charron, and S. Debrus, “Third-order nonlinear optical response of Au:SiO_{2} thin films: Influence of gold nanoparticle concentration and morphologic parameters,” Eur. Phys. J. D **19**, 395–402 (2002). [CrossRef]

**20. **O. Maruyama, Y. Senda, and S. Omi, “Non-linear optical properties of titanium dioxide films containing dispersed gold particles,” J. Non-Cryst. Solids **259**(1-3), 100–106 (1999). [CrossRef]

**21. **G. Ma, W. Sun, S.-H. Tang, H. Zhang, Z. Shen, and S. Qian, “Size and dielectric dependence of the third-order nonlinear optical response of Au nanocrystals embedded in matrices,” Opt. Lett. **27**(12), 1043–1045 (2002). [CrossRef]

**22. **J. Jayabalan, A. Singh, R. Chari, and S. M. Oak, “Ultrafast third-order nonlinearity of silver nanospheres and nanodiscs,” Nanotechnology **18**(31), 315704 (2007). [CrossRef]

**23. **O. Levy and D. Stroud, “Maxwell Garnett theory for mixtures of anisotropic inclusions: Application to conducting polymers,” Phys. Rev. B **56**(13), 8035–8046 (1997). [CrossRef]

**24. **J. C. M. Garnett, ““Colours in Metal Glasses and in Metallic Films,” Philos. Trans. R. Soc. London Ser. A **203**(1), 385–420 (1904). [CrossRef]

**25. **H. R. Ma, R. F. Xiao, and P. Sheng, “Third-order optical nonlinearity enhancement through composite microstructures,” J. Opt. Soc. Am. B **15**(3), 1022–1029 (1998). [CrossRef]

**26. **R. del Coso and J. Solis, “Relation between nonlinear refractive index and third-order susceptibility in absorbing media,” J. Opt. Soc. Am. B **21**(3), 640–644 (2004). [CrossRef]

**27. **P. B. Johnson and R. W. Christy, “Optical-Constants Of Noble-Metals,” Phys. Rev. B **6**(12), 4370–4379 (1972). [CrossRef]

**28. **Microwave Studio, Computer Simulation Technology, Darmstadt, Germany.

**29. **S. A. Maier, P. G. Kik, and H. A. Atwater, “Observation of coupled plasmon-polariton modes in Au nanoparticle chain waveguides of different lengths: Estimation of waveguide loss,” Appl. Phys. Lett. **81**(9), 1714–1716 (2002). [CrossRef]

**30. **D. D. Smith, G. Fischer, R. W. Boyd, and D. A. Gregory, “Cancellation of photoinduced absorption in metal nanoparticle composites through a counterintuitive consequence of local field effects,” J. Opt. Soc. Am. B **14**(7), 1625–1631 (1997). [CrossRef]