# Surface Plasmon States in Inhomogeneous Media at Critical and Subcritical Metal Concentrations

Authors: Katyayani Seal; Dentcho A. Genov
Journal: International Journal of Optics
 (2012)

---

## Abstract

Semicontinuous metal-dielectric films are composed of a wide range of metal clusters of various geometries—sizes as well as structures. This ensures that at any given wavelength of incident radiation, clusters exist in the film that will respond resonantly, akin to resonating nanoantennas, resulting in the broad optical response (absorption) that is a characteristic of semicontinuous films. The physics of the surface plasmon states that are supported by such systems is complex and can involve both localized and propagating plasmons. This chapter describes near-field experimental and numerical studies of the surface plasmon states in semicontinuous films at critical and subcritical metal concentrations and evaluates the local field intensity statistics to discuss the interplay between various eigenmodes.


---

## Body

## 1. Introduction

Random metal-dielectric composites, such as nanoscaled semicontinuous metallic structures, exhibit unique optical and electrical transport properties [1–26]. One of the most remarkable properties of such composites is the enhancement of electric and magnetic fields in the visible and infrared spectral ranges [1, 8–22]. These materials show promise in wide ranging photonics technologies including applications in negative index materials [27–29], biosensing [30, 31], and nonlinear optics [10] spectroscopy [32, 33] broadband optical amplifiers and optical absorbers [10, 34, 35].The synthesis of metal-dielectric composites typically involves the deposition of metal particles onto a dielectric substrate [14, 21–23]. At low metal concentrations, mutually separated nanometer-sized grains are formed at random locations. Self-similar clusters form as the metal surface concentration increases [1, 3, 10]. At a certain metal concentration (pc), that is, the percolation threshold, the metal clusters interconnect to form an infinite “backbone cluster”, setting off an insulator-to-metal phase transition. This is marked by a sharp drop in dc resistivity [14, 23] and anomalous absorption at visible and near-infrared wavelengths [1, 10, 11, 14–19, 21–23]. At even higher metal concentrations, the sample becomes mostly metallic with dielectric voids, ultimately resulting in a uniform metal film.From both theoretical predictions and experimental corroboration, it has been observed that under illumination with light at optical and near infrared frequencies, strong local field enhancement of up to 104 can be achieved on the surface of two-dimensional metal-dielectric composites [1, 10–22, 36]. This light trapping effect is facilitated by the excitation of collective electron oscillations, that is, surface plasmons (SP). Surface plasmon modes are morphological resonances governed by the structural inhomogeneities inherent in the fractal clusters formed within the random metal-dielectric films. The observed local fields appear to be strongly localized [1, 10, 11, 17–19] with relatively high Q-factors dependent on the intrinsic dissipation in the metal. This phenomenon has been observed in rough and semicontinuous metal films [8, 9, 15, 16, 23]. Notably, due to the geometry of the clusters, the spatial distribution of the field maxima is sensitive to the polarization, wavelength and angle of incidence of the applied field [1, 2, 10]. Thus the metal clusters in a semicontinuous film that support SP modes may be thought of as resonating nanoantennas with a range of resonant frequencies and Q-factors. The accumulation and subsequent dissipation of electromagnetic energy in these localized SP modes result in the anomalous absorption characteristic of semicontinuous metal films [1, 2, 5–10, 14–22].The near-field intensity distribution across a semicontinuous metal film near the “percolation threshold” has been observed to be extremely inhomogeneous with giant local field maxima [15, 16, 20–22], that vary not only in size by orders of magnitude from subwavelength to multiple wavelengths, but also vary strongly in intensity. In the early theoretical studies on percolating metal-dielectric composites, an observed similarity between the Kirchhoff Hamiltonian used to model random metal dielectric systems and the Anderson Hamiltonian [37] led to the expectation that all surface plasmon eigenmodes are localized [10]. Follow-up numerical simulations subsequently suggested that a planar disordered nanosystem can simultaneously have localized and delocalized states [18]. This concept of “inhomogeneous localization” was originally introduced for metal fractal aggregates, where self-similarity in the structural geometry can result in a mode consisting of multiple hot spots spread over the entire system [10]. Later calculations have shown that in an infinitely large semicontinuous metal film, delocalized SP modes exist only at one point in the eigenvalue spectrum and their measure is zero [19]. The conceptual origin for these delocalized states lies in the short-range correlations in the governing Kirchhoff Hamiltonian, which are absent in the Anderson Hamiltonian [37]. Remarkably, even though they have a zero measure, delocalized modes can affect the local field statistics by modifying the critical indices for the high-order field moments [19, 22] and thus influence the optical properties. The interplay between localized and delocalized states in semicontinuous films and their role in the transport properties therein remain a topic of interest and speculation.For disordered systems with a range of possible eigenmodes and transport paths, studies of field and intensity statistics are a valuable avenue to uncover the system’s dominant response and the nature of wave transport. In disordered metallic nanostructures in particular, the intensity statistics are more complex because of the excitation of surface plasmons and resulting inhomogeneous localization [18]. In the study of random metal-dielectric composites, both numerically [2, 14, 17–19] and experimentally [5–7, 15, 16, 19–23] obtained near-field intensity statistics are often used to characterize the local field fluctuations in various metal aggregates, fractal structures, and rough and semicontinuous metal films.This chapter provides a comprehensive view of the interplay between localized and delocalized plasmon excitations in metal semicontinuous films through the analysis of experimentally and numerically obtained near-field intensity statistics.


## 2. Synthesis of Semicontinuous Films and Near-Field Statistics

Semicontinuous silver films on glass substrates are typically synthesized by thermal evaporation [5–9] or by pulsed laser deposition [16]. Transmission electron microscope (TEM) images from a silver-glass semicontinuous film grown by pulsed laser deposition are shown in Figure 1(a). From the images, it is evident that these samples are composed of individual silver grains of average size about 20–30 nm. An increase in deposition time (surface concentration of silver) induces a structural transition from isolated metal grains (p<0.4) to interconnected metal clusters (p=0.6) and finally to a nearly continuous metal film with dielectric voids (p>0.8) [16]. For these samples, the percolation threshold was found to be at pc=0.65 from previous work [16]. The 2D structural correlation functions were found to be isotropic. The correlation radius, defined as the half width at half maximum (HWHM) of the structural correlation function, is close to the grain size ~20 nm.Figure 1
TEM micrographs of silver semicontinuous films on glass substrates showing the metal nanostructure geometry for various metal concentrations in the first row. The second row shows the inhomogeneous distribution of the near-field intensity for corresponding samples.Near-field optical microscopy (NSOM) is a powerful tool to image the near-field intensity distribution over the surface of a sample [15, 16]. Figure  1(b) shows NSOM images for silver on glass semicontinuous films with metal concentration p=0.36, 0.45, 0.65, 0.75, and 0.83 at a probe wavelength of 543 nm. The samples were illuminated by the evanescent field (in the total internal reflection geometry) of “p” polarized laser light. The projection of the incident beam’s k-vector on the film plane (x-yplane), k∥, is along the vertical y axis. The local optical signal was collected by a tapered, uncoated optical fiber of tip radius of about 50 nm. The tip-to-sample distance, controlled by shear-force feedback, was approximately 10 nm. The tip resolution is estimated to be of the order of 150 nm from the smallest features in the near-field images [16, 22]. The images exhibit inhomogeneous intensity distributions with intensity maxima (or optical “hot spots”) of various sizes. At p=0.36, there is an indication of interference fringes parallel to the horizontal x axis, which are weakened at p=0.45, completely lost at p=0.65, and reappear at p=0.83.To understand the role of the eigenmodes supported by various clusters in the semicontinuous film, it is instructive to study the near-field intensity statistics for different metal concentrationsp. Specifically, the variance var(I)=〈I(x,y)2〉/〈I(x,y)〉2-1 and higher-order moments M2n=〈I(x,y)n〉/I0n can be obtained from the NSOM images, where I(x,y) is the intensity over the sample surface as a function of x and y coordinates. Figure 2 shows the experimentally obtained values of var(I) as a function of Δp=p-pc at a probe wavelength of 633 nm calculated from the NSOM images similar to the ones in Figure  1(b). Each data point was obtained from approximately 300,000 near-field intensity values. The error bar was evaluated as the standard error of data collected from several NSOM images of different parts of the samples. Counterintuitively, the intensity fluctuations do not maximize at pc(Δp=0), where scattering is the strongest. Instead, a local minimum appears at pc. Away from pc,var(I) has two maxima at Δp=-0.2 and Δp=+0.1. A similar dependence of var(I) on p is also observed at lower probe wavelengths (543 nm). However, the absolute value of var(I) at λ=543nm is lower than that at λ=633 nm, indicating that the intensity fluctuations are stronger at longer wavelengths.Figure 2
Near-field intensity variance var(I) as a function of silver concentration Δp=p-pc, obtained experimentally (solid line with triangles) and numerically with (dotted line) and without (dashed line) spatial averaging [22, Figure 1].The variance of intensity fluctuations is proportional to the inverse participation ratio (IPR) of the local field, often used to describe the spatial extent of eigenstates [2]. In the experimental set up described, where the incident field excites a number of SP modes, IPR gives a qualitative measure of the surface area occupied by hot spots. Similar to var(I), IPR has a minimum at pc. Larger values of IPR at |p-pc|=0.1 and 0.2 indicate that a smaller percentage of the film surface is occupied by hot spots. This is consistent with the direct observation that away from pc, there are fewer hot spots but each is brighter.The experimentally observed dip in IPR suggests a weakening of the SP localization strength atpc. For the calculation of the high-order moments of intensity enhancement M2n, the incident light intensity I0 obtained from the near-field measurement of a clean glass substrate can be used as a standard normalizing factor. Figure 3(a) is a plot of M2n versus Δp for n=2,3,4, and 5 obtained from the same NSOM images (at a probe wavelength of 633 nm) as were used in Figure 2. The error bars were obtained by the same method as for var(I) in Figure 2. The M2n exhibit a prominent dip near pc, which becomes more pronounced for higher n. This observation once more confirms that the local field fluctuations are weakened at the percolation threshold.High-order moments of near-field intensity enhancementM2n as a function of silver concentration Δp=p-pc, obtained experimentally (a) and numerically with (b) and without (c) spatial averaging  [22, Figure 2].
(a)(b)(c)


## 3. Numerical Model of a Semicontinuous Film

In order to understand the experimental findings, it is useful to perform numerical simulations to simulate the semicontinuous film. In the visible frequency regime, metal-dielectric composites can be viewed as a system of lumpedRLCcircuits [2, 10], and the bond percolation model can be employed to describe their electromagnetic response. The bond model has a percolation threshold pc=0.5, which differs slightly from the experimental observation. This deviation, however, should not preclude a direct mapping between experiment and theory because the scaling properties of random metal films do not depend on the absolute value of pc [3].The inhomogeneous metal dielectric composite is a random system of metal nanoparticles incorporated within a dielectric host and can be modeled as a strongly coupled lumpedR-L-C network. Provided the feature sizes along the direction of illumination are small compared to the incident wavelength λ=2πc/ω, the impinging electric field E⃗0 couples with the metal dielectric composite quasistatically, with a local electric potential φ(ω) given as
(1)∇⃗⋅[σ(ω,r⃗)∇⃗φ(r⃗)]=E⃗0⋅∇⃗σ(ω,r⃗),
where the local conductivity σ(ω,r⃗) is a discontinuous random function across the composite taking discrete values σm=-iωεm for the metal particles and σd=-iωεd for the dielectric voids with probabilities p and 1-p, respectively [14]. Equation (1) describes a generalized Kirchhoff problem and can be solved numerically by discretization of the electric potential on a cubic lattice, with the metal and dielectric constituents represented by randomly distributed metal/dielectric bonds connecting each lattice site. The resulting linear system of equations is written in short form Ĥ·φ⃗=Ŵ·Ê0, where Ĥ is a complex value Kirchhoff Hamiltonian (KH), φ⃗={φi,j}is a vector containing the values of the potential at each lattice side (i,j), and Ŵ is the incident field projection operator. The Kirchhoff Hamiltonian is a random matrix with diagonal elements Hnn=∑σ{i,j}n, where the summation is over all bond conductivities connecting the nth site with its nearest neighbors, and correlated off-diagonal elements proportional to the respective bond conductivities σ{i,j}n. Traditional numerical methods such as Gaussian eliminations require O(L6) floating point operations and O(L4) operational memory to achieve solution (here L is the lateral size of the composite). For the system sizes of interest (a few μm2 or up to a million metal particles), these algorithms are not feasible even with high performance computing. To overcome this computational problem, a block elimination method has been developed that splits the Kirchhoff Hamiltonian on L×L block matrixes and performs triangulation on these separated block of matrixes [2]. This effectively allows the reduction of this two-dimensional problem to an equivalent 4/3-dimensional problem, resulting in dramatic improvement in performance with solution obtained in only O(L4) floating point operations and operational memory requirement of the order of O(L3).The numerical results based on the Block Elimination procedure are shown in Figures2, 3(b), and 3(c), where the local field distribution is calculated for 2μm×2μm-sized semicontinuous films of different silver concentrations, followed by calculations of var(I) and M2n. Detailed information on the method of calculation used is provided in [2, 10]. At the percolation threshold a highly inhomogeneous local field distribution is observed, with peak intensities surpassing the incident laser light intensity by factors of 103–104.The experimentally obtained local field enhancements are considerably lower in comparison because of the finite spatial resolution of NSOM. For optimal comparison with the experimental data, the local fields at a distance equal to the tip-sample separation (~10 nm for the NSOM configuration used) from the sample surface can be estimated numerically and then averaged using the NSOM aperture function (in this case, corresponding to a tip resolution of ~150 nm). Numerical data for var(I) and M2n before and after the spatial averaging are provided in Figures 2, 3(b), and 3(c). These curves were obtained from 100 realizations of the semicontinuous film surface for each p. As shown in Figure 2, the var(I) calculated without spatial averaging exhibits a wider and deeper dip near pc. The spatial averaging not only reduces the value of var(I) but also partially smooths out the dip at pc. After the resolution limit of the NSOM is taken into account, the numerical data closely mimic the experimental var(I). However, some of the numerical values fall outside the experimental error bars because of inevitable approximations used in the calculations such as uncertainty in the value of dielectric constant of silver, the tip-sample distance, and the spatial averaging parameter. In Figure 3, the calculated local field moments M2n exhibit a p-dependence similar to the measured one, that is,M2n reaches a local minimum at pc and has two distinct maxima at |Δp|≤0.3. The above numerical uncertainties also contribute to the discrepancies between experimental and numerical (spatially averaged) data shown in Figure 2.


## 4. A Qualitative Understanding of the Variations in Localization Strength

The observed decrease of var(I) and M2n at |Δp|≥0.3 can be qualitatively explained by considering the resonant nature of SP excitation. According to percolation theory [3], the largest metal cluster in the random composite has a size lm(p)=apcν/|p-pc|ν, where a is the average metal particle size and ν is a critical exponent. Resonant SP modes can be excited by incident light of a wavelength λ if lm(p)≥lr=a(|εm′|/εd)ν/(s+t), where εm′+iεm′′ and εd are the dielectric constants of the metal particles and dielectric host, respectively [10]. The critical exponents s and t correspond to the scaling of the composite’s static conductivity and dielectric constant, respectively, with the sample size. For the planar systems under consideration,s=ν=t=4/3 [3]. Thus the SP resonance wavelength in a cluster of size l≥lr scales as λr≈λp(l/a)ν/(s+t)≈λpl/a, where λp is the plasma wavelength [10]. This unique property of localized SPs is distinct from resonances in the retarded (nonquasistatic) regime where λr≈l.When the metal concentrationp differs considerably from pc, there is a decrease in the maximum cluster size lm. If lm(p)<lr, SP modes in the semicontinuous metal film can no longer be excited by the external field, leading to a dramatic decrease in the local electromagnetic response. The cutoff metal concentrations are estimated from the condition lr=lm(p±) as p±=pc[1±(εd/|εm′|)1/(s+t)]. At λ=633nm, Δp±=|p±-pc|≈0.3. For concentrations in the range |Δp|>Δp±, both var(I) and M2n are expected to decrease sharply. Estimates for the cutoff concentrations Δp agree well with the experimental and numerical results shown in Figures 2 and 3. The presence of cutoff metal concentrations also manifests itself as an anomalous absorption band in semicontinuous metal films [23]. Specifically, when κ=|εm′′/εm′|≪1, the buildup of giant local fields and storage of electromagnetic energy in the films leads to strong light absorption, which occurs for metal concentrations inside the band pc-Δp-<p<pc+Δp+ [26]. The quantity Δp±, calculated as |εm′|/εd/(|εm′|+εd) from effective medium theory, gives values close to the cutoff concentration estimates based on the resonant cluster size.


## 5. Density of States and Localization Length from the SP Eigenproblem

To further understand the behavior of var(I) and M2n at pc, it is instructive to solve the SP eigenproblem. In order to model the local electric response of a random system comprised of a large number of metal particles, one can work in the quasistatic limit. The general solution of the current conservation (1) can be described completely by examining the underlying eigenproblem [10, 19](2)∇⇀⋅[σ(r⃗)∇⇀ΨΛ(r⃗)]=ΛΨΛ(r⃗).Reduction of (2) on a square lattice leads to a system of linear equations ĤKH|n〉=Λn|n〉, where again the Kirchhoff Hamiltonian ĤKH=∑nhnn|n〉〈n|+∑n≠khnk|n〉〈k| is a random matrix with correlated diagonal and off-diagonal elements and |n〉 are the SP eigenstates [5–7]. The existence of short range correlations in the SP eigenproblem makes it substantially different than the corresponding quantum mechanical case where the random Hamiltonians are noncorrelated [19]. The correlations in the KH are due to the local current conservation and exist for all spatial dimensions.The statistical properties of the solution to the SP eigenproblem can be investigated in terms of the density of statesρ(λ) and SP localization lengths ξ(λ). These are plotted for different metal concentrations in Figures 4(a) and 4(b), respectively. At the band center, and concentrations close to percolation, the density of states exhibits a singularity as shown in Figure 4(a). In the first approximation the density of states seems to diverge as a power law ρ(Λ)~|Λ|-γ at p=pc (red curve), where γ=0.14±0.01 [17] is a critical exponent. The localization length for each eigenmode (presented in Figure 4(b)) is also calculated using the gyration radius  ξ2(p,Λ)=∫(r⃗-<r⃗>Λ)2|ΨΛ(r⃗,p)|2dr⃗, where <r⃗>Λ=∫r⃗|ΨΛ(r⃗,p)|2dr⃗ is the “mass center” of the nth mode and the integration is performed over the film surface. Similar to the density of states, there is a singularity at the band center. At pc, the localization length diverges logarithmically for Λ→0, but also can be fitted with a power law ξ(Λ)~|Λ|-α, where α=0.15±0.02.The surface plasmons density of states (DOS) (a) and localization length (b) are calculated for an ensemble of 100 metal-dielectric films each with size1μm. Four different metal filling fractions are investigated: p=0.1 (yellow color), p=0.2 (green color), p=0.3 (blue color), and p=0.5 (red color). The singularity at the band center and its dependence on the metal concentration is closely inspected in the insets [24, Figure 2].
(a)(b)Relying on the SP eigenproblem (2), it is possible to estimate the optical response of the random metal films for an incident frequency ωr corresponding to the single particle resonance Reεm(ωr)=-εd. Specifically, the local potential is expanded over the SP eigenspace [10] to obtain an integral relationship for the local field moments(3)M2n=1S|E0|2n〈∫|E(r⃗,p)|2ndr⃗〉≅∫-∞∞ρ(p,Λ)[a/ξ(p,Λ)]2(2n-1)(Λ2+κ2)ndΛ,
where κ=|εm′′/εm′|≪1 is the intrinsic loss factor, S is the sample area, and the brackets correspond to ensemble averaging over large numbers of system realizations.The integral in (3) can be analyzed as follows. Close to the band center (|Λ|≪1), the properties of the generalized SP eigenstates can be modeled (see inserts in Figure 4) by simplified expressions for the SP DOS fp(ρ,λ)=ρ(p,λ)/ρ(pc,λ)≈1 and the SP localization length fξ(p,Λ)=ξ(p,Λ)/ξ(pc,Λ)≈g(p)+θ(Λ)(1-g(p)), where θ(Λ) is a step function and g(p)<1. Since only the SP eigenmodes in the vicinity of the band center have a dominant contribution to the moments, the p dependence from the integral (3) can then be taken out to give M2n(p)≈M2n(pc)[1+g4(1+n)(p)],where M2n(pc) are the local field moments at percolation [19]. Using g(p)<1, the limiting condition M2n(p)>M2n(pc)g4(1-n)(p)>M2n(pc) is obtained, which implies that at pc and for n>1, the local field moments have a minimum. This simple reasoning based on the general properties of the SP eigenstates provides a qualitative explanation of the observed dip in the field moments as per Figure 3.From Figure4(a), it is also clear that for p<pc the localization lengths of the excited SP eigenstates, predominantly those with Λ<0, is lower as compared to the critical case (p=pc). This suggests that at noncritical metal concentrations the role of the localized states in the local field statistics becomes stronger. Moreover, since current conservation holds in all cases, one may expect maxima in ρ(λ) and ξ(λ) at noncritical metal concentrations as well, as shown in Figures 4(a) and 4(b). With the decrease of p, a shift in the DOS is observed with an increase of the number of eigenstates with positive eigenvalues. This phenomenon corresponds to the fact that the eigenproblem maps the underlying film geometry with positive eigenvalues corresponding to eigenstates situated mostly in the dielectric host. It should be pointed out that although Figure 4 presents only the results for p<pc, the properties of the SP eigenstates at high and low metal concentrations are equivalent due to symmetry. Thus, the relationships ρ(Λ,pc+Δp)=ρ(Λ,pc-Δp) and ξ(Λ,pc+Δp)=ξ(Λ,pc+Δp) hold for any given Δp≤pc.


## 6. Behavior of the Slope of the Higher-Order Moments (M2n)

The increased contribution of the localized states atp≠pc, or alternatively the increase of SP delocalization at pc, is revealed through the dependence of In[M2n(p)] on the order n. According to the scaling theory described above, a linear dependence of the form In[M2n(p)]~2ns(p)] is to be expected, with the slope s(p)=s(pc)+4In[1/g(p)]≥s(pc) increasing for metal concentrations away from the percolation threshold. Figure 5(a) is a log-linear plot of the measured M2n versus order n at pc for two probe wavelengths 543 nm and 633 nm. Clearly, a linear dependence on the field moment n is observed. However, due to the limited NSOM resolution, the slopes obtained experimentally were found to be considerably lower than those obtained numerically. To account for the effect of the finite spatial resolution in the experiment, the calculated local field distribution was first averaged over an area of dimension d, and then M2n was computed to extract the slope of In[M2n(p)] versus n. Figure 5(b) shows that the slope decreases as d increases. When the slope is reduced to the experimental values (marked by horizontal lines), the averaging parameter d is found to be very close to the NSOM resolution limit estimated from the smallest features in the near-field images. The same value of d also leads to good agreement between the calculated var(I) and the experimental data.(a) Log-linear plot of measured high-order moments of near-field intensity enhancementM2n versus n. (b) Slope of log(M2n) over n at the percolation threshold obtained numerically after spatial averaging versus averaging parameter d. The measured slopes are marked by horizontal lines with the standard error given in the legend [22, Figure 3].
(a)(b)To extend this result for metal concentrations away frompc, Figure 6 plots the slope as described above for various metal concentrations. In the calculations, the local field is again spatially averaged to account for the finite NSOM resolution. Figure 6 shows a similar trend with respect to p for both experiment and numerical calculations, with a clear minimum at pc. Although the exact numerical values of the slope at p≠pc cannot be easily interpreted based only on the SP eigenproblem, the observed increase in slope away from pc clearly indicates a transition at pc. The analysis following (3) suggests that stronger SP localization is favored at these metal concentrations, which leads to stronger local field fluctuations. The higher-order moments are more sensitive to those fluctuations, a fact that is manifested by the increase in the slope. It must be pointed out that for p>pc, apart from the quasistatic SP modes caused by the sub-wavelength features in the system, the medium could support propagating surface waves and surface plasmon polaritons [20]. However, the local fields modulations that are associated with those propagating modes are weaker compared to the quasistatic resonances and should have a relatively small effect on the local field distributions. This conclusion is indirectly supported by the good correlation between the numerical calculations and experiment.Figure 6
Slope versus metal concentration for numerical (dashed line) and experimental (solid line and square symbols) data. The wavelengths of the incident light areλ=543nm (green color) and λ=633nm (red color) [24, Figure 3].


## 7. Conclusion

The implications of the above statistical analysis in practical applications are wide-ranging. The high-order local field momentsM2n correspond to the enhancement of multiphoton optical processes due to the excitation of SP resonances in the substrate. For instance, M4 relates to the electromagnetic enhancement of Raman-scattering provided that Raman active molecules are deposited on the metal surface. Since M2n(p)≥M2n(pc), one may conclude that contrary to common practice, a highly effective surface enhanced Raman substrate can be achieved at low and high metal concentration. A similar conclusion can be made for the higher-order optical processes such as hyper-Raman scattering, Kerr optical nonlinearity, and nth harmonic generation.To conclude, the near-field intensity statistics in semicontinuous silver films over a range of metal concentrations is effective in providing an overview of the various eigenstates supported by the film and the interplay between localized and delocalized SP eigenmodes. The variance of intensity fluctuations and the high-order moments of intensity enhancement exhibit local minima at the percolation threshold. Experimental observations, through theoretical corroboration, indicate a decrease of the SP localization in metal composites at critical concentrations with maxima in localization at concentrations away from the critical metal concentration. This effect runs opposite to the intuitive notion that stronger interactions lead to increase in localization. A possible explanation for this behavior is the existence of delocalized SP states, due to the local current conservation. Those states are present in the spectra for all metal concentrations and for metal-dielectric systems of any dimension. The possibility for this coexistence of both localized and delocalized SP modes has important ramifications in the understanding and applications of random metal films. Fine tuning the geometrical structure of the metal composites may allow for tuning of the localization properties of the excited SP modes, which may improve the enhancement of various linear and nonlinear optical processes.



---

*Source: 519509-2012-04-02.xml*