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

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.

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 (𝑝_𝑐), 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 10⁴ 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 (𝑝<0.4) to interconnected metal clusters (𝑝=0.6) and finally to a nearly continuous metal film with dielectric voids (𝑝>0.8) [16]. For these samples, the percolation threshold was found to be at 𝑝_𝑐=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 𝑝=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-y plane), 𝐤_‖, is along the vertical 𝑦 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 𝑝=0.36, there is an indication of interference fringes parallel to the horizontal 𝑥 axis, which are weakened at 𝑝=0.45, completely lost at 𝑝=0.65, and reappear at 𝑝=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 concentrations 𝑝. Specifically, the variance var(𝐼)=⟨𝐼(𝑥,𝑦)²⟩/⟨𝐼(𝑥,𝑦)⟩²−1 and higher-order moments 𝑀_2𝑛=⟨𝐼(𝑥,𝑦)^𝑛⟩/𝐼₀^𝑛 can be obtained from the NSOM images, where 𝐼(𝑥,𝑦) is the intensity over the sample surface as a function of x and y coordinates. Figure 2 shows the experimentally obtained values of var(𝐼) as a function of Δ𝑝=𝑝−𝑝_𝑐 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 𝑝_𝑐(Δ𝑝=0), where scattering is the strongest. Instead, a local minimum appears at 𝑝_𝑐. Away from 𝑝_𝑐,var(𝐼) has two maxima at Δ𝑝=−0.2 and Δ𝑝=+0.1. A similar dependence of var(𝐼) on 𝑝 is also observed at lower probe wavelengths (543 nm). However, the absolute value of var(𝐼) 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 Δ𝑝=𝑝−𝑝𝑐, 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(𝐼), IPR has a minimum at 𝑝_𝑐. Larger values of IPR at |𝑝−𝑝_𝑐|=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 𝑝_𝑐, there are fewer hot spots but each is brighter.

The experimentally observed dip in IPR suggests a weakening of the SP localization strength at 𝑝_𝑐. For the calculation of the high-order moments of intensity enhancement 𝑀_2𝑛, the incident light intensity 𝐼₀ 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 𝑀_2𝑛 versus Δ𝑝 for 𝑛=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(𝐼) in Figure 2. The 𝑀_2𝑛 exhibit a prominent dip near 𝑝_𝑐, which becomes more pronounced for higher 𝑛. This observation once more confirms that the local field fluctuations are weakened at the percolation threshold.

(a)

(b)

(c)

Figure 3 

High-order moments of near-field intensity enhancement 𝑀2𝑛 as a function of silver concentration Δ𝑝=𝑝−𝑝𝑐, obtained experimentally (a) and numerically with (b) and without (c) spatial averaging [22, Figure 2].

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 lumped RLC circuits [2, 10], and the bond percolation model can be employed to describe their electromagnetic response. The bond model has a percolation threshold 𝑝_𝑐=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 𝑝_𝑐 [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 lumped R-L-C network. Provided the feature sizes along the direction of illumination are small compared to the incident wavelength 𝜆=2𝜋𝑐/𝜔, the impinging electric field 𝐸⃗₀ couples with the metal dielectric composite quasistatically, with a local electric potential 𝜑(𝜔) given as ∇⃗⋅𝜎(𝜔,𝑟)⃗∇⃗𝜑(𝑟)⃗=𝐸⃗₀⋅∇⃗𝜎(𝜔,𝑟)⃗, (1) where the local conductivity 𝜎(𝜔,𝑟)⃗ is a discontinuous random function across the composite taking discrete values 𝜎_𝑚=−𝑖𝜔𝜀_𝑚 for the metal particles and 𝜎_𝑑=−𝑖𝜔𝜀_𝑑 for the dielectric voids with probabilities 𝑝 and 1−𝑝, 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 𝐻⋅𝜑⃗=𝑊⋅𝐸₀, where 𝐻 is a complex value Kirchhoff Hamiltonian (KH), 𝜑⃗={𝜑_𝑖,𝑗}is a vector containing the values of the potential at each lattice side (𝑖,𝑗), and 𝑊 is the incident field projection operator. The Kirchhoff Hamiltonian is a random matrix with diagonal elements 𝐻_𝑛𝑛=∑𝜎_{{𝑖,𝑗}𝑛}, where the summation is over all bond conductivities connecting the 𝑛th site with its nearest neighbors, and correlated off-diagonal elements proportional to the respective bond conductivities 𝜎_{{𝑖,𝑗}𝑛}. Traditional numerical methods such as Gaussian eliminations require 𝑂(𝐿⁶) floating point operations and 𝑂(𝐿⁴) operational memory to achieve solution (here 𝐿 is the lateral size of the composite). For the system sizes of interest (a few 𝜇𝑚² 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 𝐿×𝐿 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 𝑂(𝐿⁴) floating point operations and operational memory requirement of the order of 𝑂(𝐿³).

The numerical results based on the Block Elimination procedure are shown in Figures 2, 3(b), and 3(c), where the local field distribution is calculated for 2𝜇𝑚×2𝜇𝑚-sized semicontinuous films of different silver concentrations, followed by calculations of var(I) and 𝑀_2𝑛. 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 10³–10⁴.

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 𝑀_2𝑛 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 𝑝. As shown in Figure 2, the var(I) calculated without spatial averaging exhibits a wider and deeper dip near 𝑝_𝑐. The spatial averaging not only reduces the value of var(I) but also partially smooths out the dip at 𝑝_𝑐. 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 𝑀_2𝑛 exhibit a p-dependence similar to the measured one, that is,𝑀_2𝑛 reaches a local minimum at 𝑝_𝑐 and has two distinct maxima at |Δ𝑝|≤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 𝑀_2𝑛 at |Δ𝑝|≥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 𝑙_𝑚(𝑝)=𝑎𝑝^𝜈_𝑐/|𝑝−𝑝_𝑐|^𝜈, where 𝑎 is the average metal particle size and 𝜈 is a critical exponent. Resonant SP modes can be excited by incident light of a wavelength 𝜆 if 𝑙_𝑚(𝑝)≥𝑙_𝑟=𝑎(|𝜀^′_𝑚|/𝜀_𝑑)^{𝜈/(𝑠+𝑡)}, where 𝜀^′_𝑚+𝑖𝜀_𝑚^′′ and 𝜀_𝑑 are the dielectric constants of the metal particles and dielectric host, respectively [10]. The critical exponents 𝑠 and 𝑡 correspond to the scaling of the composite’s static conductivity and dielectric constant, respectively, with the sample size. For the planar systems under consideration,𝑠=𝜈=𝑡=4/3 [3]. Thus the SP resonance wavelength in a cluster of size 𝑙≥𝑙_𝑟 scales as 𝜆_𝑟≈𝜆_𝑝(𝑙/𝑎)^{𝜈/(𝑠+𝑡)}≈𝜆_𝑝𝑙/𝑎, where 𝜆_𝑝 is the plasma wavelength [10]. This unique property of localized SPs is distinct from resonances in the retarded (nonquasistatic) regime where 𝜆_𝑟≈𝑙.

When the metal concentration 𝑝 differs considerably from 𝑝_𝑐, there is a decrease in the maximum cluster size 𝑙_𝑚. If 𝑙_𝑚(𝑝)<𝑙_𝑟, 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 𝑙_𝑟=𝑙_𝑚(𝑝_±) as 𝑝_±=𝑝_𝑐[1±(𝜀_𝑑/|𝜀^′_𝑚|)^{1/(𝑠+𝑡)}]. At 𝜆=633nm, Δ𝑝_±=|𝑝_±−𝑝_𝑐|≈0.3. For concentrations in the range |Δ𝑝|>Δ𝑝_±, both var(I) and 𝑀_2𝑛 are expected to decrease sharply. Estimates for the cutoff concentrations Δ𝑝 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 𝜅=|𝜀_𝑚^′′/𝜀^′_𝑚|≪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 𝑝_𝑐−Δ𝑝₋<𝑝<𝑝_𝑐+Δ𝑝₊ [26]. The quantity Δ𝑝_±, calculated as |𝜀^′_𝑚|/𝜀_𝑑/(|𝜀^′_𝑚|+𝜀_𝑑) 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 𝑀_2𝑛 at 𝑝_𝑐, 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)

Reduction of (2) on a square lattice leads to a system of linear equations 𝐻_KH|𝑛⟩=Λ_𝑛|𝑛⟩, where again the Kirchhoff Hamiltonian 𝐻_KH=∑_𝑛ℎ_𝑛𝑛|𝑛⟩⟨𝑛|+∑_𝑛≠𝑘ℎ_𝑛𝑘|𝑛⟩⟨𝑘| is a random matrix with correlated diagonal and off-diagonal elements and |𝑛⟩ 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 𝑝=𝑝_𝑐 (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 𝜉²(𝑝,Λ)=∫(𝑟−⃗<𝑟>⃗_Λ)²|Ψ_Λ(𝑟,⃗𝑝)|²𝑑𝑟⃗, where <𝑟⃗>_Λ=∫𝑟|⃗Ψ_Λ(𝑟,⃗𝑝)|²𝑑𝑟⃗ is the “mass center” of the 𝑛th 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 𝑝_c, the localization length diverges logarithmically for Λ→0, but also can be fitted with a power law 𝜉(Λ)∼|Λ|^−𝛼, where 𝛼=0.15±0.02.

(a)

(b)

Figure 4 

The surface plasmons density of states (DOS) (a) and localization length (b) are calculated for an ensemble of 100 metal-dielectric films each with size 1𝜇𝑚. Four different metal filling fractions are investigated: 𝑝=0.1 (yellow color), 𝑝=0.2 (green color), 𝑝=0.3 (blue color), and 𝑝=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].

Relying on the SP eigenproblem (2), it is possible to estimate the optical response of the random metal films for an incident frequency 𝜔_𝑟 corresponding to the single particle resonance Re𝜀_𝑚(𝜔_𝑟)=−𝜀_𝑑. Specifically, the local potential is expanded over the SP eigenspace [10] to obtain an integral relationship for the local field moments𝑀_2𝑛=1𝑆𝐸₀^2𝑛𝐸𝑟,⃗𝑝^2𝑛𝑑𝑟⃗≅^∞_−∞𝜌𝑝,Λ𝑎/𝜉𝑝,Λ^2(2𝑛−1)Λ²+𝜅^2𝑛𝑑Λ, (3) where 𝜅=|𝜀_𝑚^′′/𝜀^′_𝑚|≪1 is the intrinsic loss factor, 𝑆 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 𝑓_𝑝(𝜌,𝜆)=𝜌(𝑝,𝜆)/𝜌(𝑝_𝑐,𝜆)≈1 and the SP localization length 𝑓_𝜉(𝑝,Λ)=𝜉(𝑝,Λ)/𝜉(𝑝_𝑐,Λ)≈𝑔(𝑝)+𝜃(Λ)(1−𝑔(𝑝)), where 𝜃(Λ) is a step function and 𝑔(𝑝)<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 𝑀_2𝑛(𝑝)≈𝑀_2𝑛(𝑝_𝑐)[1+𝑔^4(1+𝑛)(𝑝)],where 𝑀_2𝑛(𝑝_𝑐) are the local field moments at percolation [19]. Using 𝑔(𝑝)<1, the limiting condition 𝑀_2𝑛(𝑝)>𝑀_2𝑛(𝑝_𝑐)𝑔^4(1−𝑛)(𝑝)>𝑀_2𝑛(𝑝_𝑐) is obtained, which implies that at 𝑝_𝑐 and for 𝑛>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 Figure 4(a), it is also clear that for 𝑝<𝑝_𝑐 the localization lengths of the excited SP eigenstates, predominantly those with Λ<0, is lower as compared to the critical case (𝑝=𝑝_𝑐). 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 𝑝, 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 𝑝<𝑝_𝑐, the properties of the SP eigenstates at high and low metal concentrations are equivalent due to symmetry. Thus, the relationships 𝜌(Λ,𝑝_𝑐+Δ𝑝)=𝜌(Λ,𝑝_c−Δ𝑝) and 𝜉(Λ,𝑝_𝑐+Δ𝑝)=𝜉(Λ,𝑝_𝑐+Δ𝑝) hold for any given Δ𝑝≤𝑝_𝑐.

6. Behavior of the Slope of the Higher-Order Moments (𝑀_2𝑛)

The increased contribution of the localized states at 𝑝≠𝑝_𝑐, or alternatively the increase of SP delocalization at 𝑝_𝑐, is revealed through the dependence of In[𝑀_2𝑛(𝑝)] on the order 𝑛. According to the scaling theory described above, a linear dependence of the form In[𝑀_2𝑛(𝑝)]∼2𝑛𝑠(𝑝)] is to be expected, with the slope 𝑠(𝑝)=𝑠(𝑝_𝑐)+4In[1/𝑔(𝑝)]≥𝑠(𝑝_𝑐) increasing for metal concentrations away from the percolation threshold. Figure 5(a) is a log-linear plot of the measured 𝑀_2𝑛 versus order 𝑛 at 𝑝_𝑐 for two probe wavelengths 543 nm and 633 nm. Clearly, a linear dependence on the field moment 𝑛 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 𝑑, and then 𝑀_2𝑛 was computed to extract the slope of In[𝑀_2𝑛(𝑝)] versus 𝑛. Figure 5(b) shows that the slope decreases as 𝑑 increases. When the slope is reduced to the experimental values (marked by horizontal lines), the averaging parameter 𝑑 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 𝑑 also leads to good agreement between the calculated var(I) and the experimental data.

(a)

(b)

Figure 5 

(a) Log-linear plot of measured high-order moments of near-field intensity enhancement 𝑀2𝑛 versus 𝑛. (b) Slope of log(𝑀2𝑛) over 𝑛 at the percolation threshold obtained numerically after spatial averaging versus averaging parameter 𝑑. The measured slopes are marked by horizontal lines with the standard error given in the legend [22, Figure 3].

To extend this result for metal concentrations away from 𝑝_𝑐, 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 𝑝 for both experiment and numerical calculations, with a clear minimum at 𝑝_𝑐. Although the exact numerical values of the slope at 𝑝≠𝑝_𝑐 cannot be easily interpreted based only on the SP eigenproblem, the observed increase in slope away from 𝑝_𝑐 clearly indicates a transition at 𝑝_𝑐. 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 𝑝>𝑝_𝑐, 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 moments 𝑀_2𝑛 correspond to the enhancement of multiphoton optical processes due to the excitation of SP resonances in the substrate. For instance, 𝑀₄ relates to the electromagnetic enhancement of Raman-scattering provided that Raman active molecules are deposited on the metal surface. Since 𝑀_2𝑛(𝑝)≥𝑀_2𝑛(𝑝_𝑐), 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 𝑛th 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.

Acknowledgments

Research partly supported by the U.S. Department of Energy, Basic Energy Sciences, Materials Sciences and Engineering Division, through Oak Ridge National Laboratory, U.S. Department of Energy (KS) and the Louisiana Board of Regents and NSF under Contracts nos. LEQSF(2011-14)-RD-A-18 and NSF(2010)-PFUND-20 (DAG).