Universal dispersion of surface plasmons in flat nanostructures.

Dimensionality has a significant impact on the optical properties of solid-state nanostructures. For example, dimensionality-dependent carrier confinement in semiconductors leads to the formation of quantum wells, quantum wires and quantum dots. While semiconductor properties are governed by excitonic effects, the optical response of metal nanostructures is dominated by surface plasmons. Here we find that, in contrast to excitonic systems, the mode dispersions in plasmonic structures of different dimensionality are related by simple scaling rules. Employing electron energy loss spectroscopy, we show that the modes of silver nanodisks can be scaled to the surface and edge modes of extended silver thin films. We thereby introduce a general and intuitive ordering scheme for plasmonic excitations with edge and surface modes as the elementary building blocks.

Light can be confined to subdiffraction volumes by binding light fields to nanostructures. Besides being of fundamental interest, this topic holds promise for a variety of photonics applications, such as optical communication and storage1 or quantum optics2. Surface plasmons (<span class="Disease">SPs) at <span class="Chemical">metal surfaces, that is, resonant collective electron oscillations coupled to an evanescent light field, play a central role in this quest34. SPs display highly confined modes with accordingly high mode densities that thus interact particularly strong with light. In addition, the spatial and spectral properties of SPs can be easily tuned by controlling the structure geometry. On one hand, this enables to tailor the SP field pattern. On the other hand, the SP resonances can be set within a wide range of visible and infrared frequencies5. Improved experimental and theoretical characterization methods have recently considerably widened our insight in the fundamental SP properties. By using scanning optical near-field microscopy, both the intensity and phase patterns of nanoparticle SP modes were imaged in the visible6 and infrared spectral ranges7. Unprecedented spatial resolution is provided by scanning transmission electron microscopes (STEMs) that, in combination with electron energy loss spectroscopy (EELS)8, can probe the (partial) local density of states of plasmonic modes910. This approach was used to map SP modes in nanoparticles1112, nanowires13 and coupled nanoparticles14, enabling even full three-dimensional (3D) reconstruction of field patterns1516. In particular, STEM–EELS can probe optically dark modes1718, that is, modes with a vanishing net dipole moment, thereby complementing optical methods to map the full modal spectrum of a plasmonic system.

Here we use STEM–EELS on flat <span class="Chemical">metallic nanostructures, as typically produced with electron-beam lithography, to show that plasmonic excitations can be reduced to surface and edge modes as the two elementary building blocks. We probe the specific SP modes sustained by the flat surface and the edge of a <span class="Chemical">silver thin film, comparing an extended system with a straight edge and a nanoparticle (nanodisk) with a radius as small as 50 nm. We find that the nanosdisk modes correspond perfectly to those of the extended film, within measurement error, on applying simple scaling rules, and support our finding with a Drude-like model. In analogy to the hybridization model of plasmon modes19, this allows for a simple and intuitive ordering scheme of SPs. In what follows, we systematically work through the SP modes of varying dimensionality (quasi-2D film-like modes, quasi-1D edge-like modes and quasi-0D particle modes), and finally submit our results to the unifying scaling laws.

Results

Film plasmons

We start with a 50-nm-thick extended <span class="Chemical">silver film on a 15-nm <span class="Chemical">Si3N4 membrane (Fig. 1a, left). The EEL spectrum measured with the focused electron beam of the STEM shows three peaks in the energy range from 1.5 to 5 eV (Fig. 1b, red spectrum). Besides the well-known bulk plasmon peak (B) at 3.75 eV20, we find two peaks at lower energies that, as we show in the following, are due to film SP modes. Assuming the general case of a coupling of the SP modes at both film interfaces, we label the lower and higher energy modes FA and FS, respectively, with the subscripts standing for the antisymmetric (A) and symmetric (S) nature of the transversal magnetic field component. As evident from the calculated field profiles in Fig. 1c, the field strength maxima of FS and FA are in the silver/vacuum and Si3N4/silver interfaces, respectively, with only weak mutual coupling for the chosen film thickness. We corroborate this interpretation with the EEL spectrum acquired from a 50-nm-thick silver film covered with a 22-nm-thick SiO2 layer (Fig. 1a, right). As illustrated by the pink spectrum in Fig. 1b, the SP mode FA (as well as the bulk mode B) is not affected by the cover layer. In contrast, the FS mode at the upper silver interface exhibits a distinct energy shift from 3.6 to 3.3 eV due to the presence of the SiO2 layer.

Figure 1

Thin-film SP modes measured by EELS.

(a) Sketch of the sample geometry, (1) a 50-nm-thick silver film on a 15-nm-thick Si3N4 substrate, (2) the same with an additional 22-nm-thick SiO2 layer. (b, top) Calculated dispersion relations of antisymmetric (solid lines) and symmetric (dashed lines) SP modes SPA and SPS, respectively, for both sample geometries (1), red and (2), pink. (b, bottom) EEL spectra acquired on both systems (1) and (2). (c) Calculated magnetic field profile for FA and FS for sample geometry (1).

The spectral position of the film mode peaks can be assigned to the asymptotic branches of the corresponding dispersion relations. This is illustrated in the upper panel of Fig. 1b showing calculated dispersion relations in a k(E) plot (k wavenumber, E energy) for the layer systems with and without <span class="Chemical">SiO2 cover layer for both symmetric and antisymmetric coupling. The asymptotic SP energies at large k-values correspond well to the EEL peak energies (Fig. 1b, lower panel) that reflect the SP density of states. We observe, however, a small systematic shift of the peak maximum towards lower energy values, which we attribute to the strong SP damping at large k-values21 and the accordingly reduced interaction strength with the electron beam.

Edge plasmons

Next, we analyse a 30-nm-thick <span class="Chemical">silver film with lithographically defined straight edges and an overall dimension of 100 × 100 μm2 on a 15-nm-thick <span class="Chemical">Si3N4 membrane. Great care was taken to produce defect-free edges for which the reduced film thickness turned out advantageous, while leading to no significant increase in coupling between both film interfaces as compared with the case in Fig. 1. We note that a thickness of 30 nm is maintained for all considered structures in the following. While a thickness dependence is not included in the scaling laws that we discuss, we show in Supplementary Fig. 1 that such a dependence is only weak for thickness range relevant for flat lithographed structures (≥30 nm). A sketch and a TEM image of the sample are shown in Fig. 2a. The EEL spectrum acquired on the film (integrated along the red line in Fig. 2a, distance to edge 45 nm) as plotted in Fig. 2b fits well with the corresponding spectrum in Fig. 1b. The higher relative intensity of the FA mode that is mainly localized at the substrate/silver interface is due to the reduced silver thickness and further illustrates our correct mode assignment. In the EEL spectrum taken outside the silver film at a distance of 15 nm from the edge (integrated along the blue line in Fig. 2a), neither the bulk plasmon peak nor the two F peaks are present. However, two new peaks emerge at loss energies of 2.3 and 3.4 eV. In the following we argue that these peaks are due to antisymmetric and symmetric SP edge modes, which we label EA and ES. These are indeed distinct peaks that do not evolve due to spectral shifts from the film mode peaks when moving across the film edge, as evidenced by a series of spectra shown in Supplementary Fig. 2 and discussed in Supplementary Note 1.

Figure 2

Film versus edge plasmons.

(a) Sketch (left) and TEM image (right, scale bar in upper corner is 500 nm) of the edge of a 30-nm-thick silver film on a 15-nm-thick Si3N4 substrate. The dashed arrows in the sketch indicate the electron trajectories transversing (red) and bypassing (blue) the film, the red and blue solid lines in the image mark the sites where the spectra in b were measured. (b) EEL spectra acquired on the film 45 nm away from the edge (red) and 15 nm outside the edge (blue), as marked in a.

Edge plasmon dispersion

As we now discuss for the EA mode, the edge plasmon dispersion relation can be deduced by probing laterally confined film edges. The length of the edges was 0.3, 0.95, 1.25 and 1.55 μm while the film extension in the orthogonal direction was kept constant at 100 μm. Figure 3a exemplarily shows a TEM image of the 0.95-μm-long edge (for the other edge lengths see Supplementary Fig. 3) superimposed with EEL maps, revealing standing wave patterns due to edge plasmon reflection at the film corners. The energy ranges as indicated for each map were chosen from the spectrum in Fig. 3b that shows the EEL signal integrated along the edge. The individual peaks EA (m being the linear plasmon mode order) are fit with Gaussian functions and the energy windows for the EEL maps in Fig. 3a are indicated. Mode orders higher than 5 are not taken into account, because the corresponding peaks are energetically too close to be resolved experimentally. To extract the mode wavenumbers k=2π/λ from the data, cross-cuts taken from the EEL maps are plotted in Fig. 3c and fit by sinusoidal functions describing standing waves with wavelength λ. The vicinity of the lateral film boundaries (marked by the dashed lines) was excluded from the fitting procedure due to obvious edge- and reflection phase-induced effects distorting the signal. We now combine the data to a dispersion relation by plotting, with blue symbols in Fig. 3d, the peak energies retrieved from the <span class="Disease">fits in Fig. 3b versus the according wavenumber values (Fig. 3c). We note that the estimated measurement inaccuracies are smaller than the symbol size, as discussed in detail in Supplementary Figs 3–5 and Supplementary Note 2. The symbols trace a function that is significantly separated from the film dispersion relation (FA) plotted in red, which is calculated for an extended 30-nm-thick <span class="Chemical">silver film on a 15-nm-thick Si3N4 membrane. The grey line is plotted as a guide to the eye only, as the implementation of the calculated edge mode dispersion22 goes beyond the scope of this paper. We interpret our results as evidencing a specific plasmon edge mode that in a confined geometry can be well described in a standing wave picture, along the lines of former work by Gu et al.23, Nelayah et al.24 and Novotny25.

Figure 3

Plasmon modes of straight edges and their dispersion relation.

(a) TEM image (scale bar, 500 nm) of a laterally structured 30-nm-thick silver film with a 0.95-μm-long lower edge superimposed with EEL maps acquired in the indicated energy ranges. The maps show standing wave patterns corresponding to antisymmetric plasmons of linear mode order m, denoted by . (b) EEL spectrum (red line) integrated along the lower film edge shown in a. The energy windows of the EEL maps in a are marked by the dashed boxes, the Gaussian curves (black lines) are fits to the experimental EEL peaks. (c) Intensity cross-cuts taken from a, the dashed lines mark the lateral film extension. The parts plotted in black (experimental data) are fit by a sinusoidal describing a standing wave (red lines). (d) Dispersion relation of edge plasmons (blue symbols), the grey line is a guide to the eye. The error bar size is discussed in detail in Supplementary Figs 3–5 and Supplementary Note 2. The red curve plots the calculated FA dispersion relation for an extended 30-nm-thick silver film on a 15-nm-thick Si3N4 membrane.

The case of nanodisks

Having identified the film and edge modes of a semi-infinite thin film, the question of the applicability of this concept to geometries of different dimensionality, in particular to nanoparticles smaller or comparable with the light wavelength, arises. Indeed, specific plasmonic excitations with radial symmetry (‘breathing’ modes) on <span class="Chemical">silver nanodisks have recently been identified as antisymmetric film <span class="Disease">SPs confined to the disk geometry18. In this case the SP wavenumber is scaled by the radial confinement, with k=2π/λ=2πn/d, where λ is the SP wavelength, d is the disk diameter and n is the radial mode order. Deducing the wavenumbers from disks with diameters in the 100–800 nm range (disk height 30 nm, substrate 15-nm-thick Si3N4) and combining them with the according EEL peak energies yield the data set plotted by the red symbols in Fig. 4. The data perfectly fit within measurement accuracy to the calculated dispersion relation of the antisymmetric SP (FA) of a 30-nm-thick film on a 15-nm-thick Si3N4 substrate18.

Figure 4

Nanodisk modes interpreted as edge and film plasmons.

(a) Data points retrieved from the breathing modes and from the dipolar/multipolar modes of nanodisks are plotted by the red and blue symbols, respectively. The error bar size is discussed in detail in Supplementary Figs 3–5 and Supplementary Note 2. The nanodisks are 30 nm high with diameters of 100–800 nm, the substrate is a 15-nm-thick Si3N4 membrane. The red curve is the calculated dispersion relation of the FA mode of a 30-nm-thick silver film on a 15-nm-thick Si3N4 substrate. The grey line is the guide to the eye copied from Fig. 3d. The orange symbols show the two in-plane dipolar excitations of an elliptical nanoparticle with axes lengths of 100 and 200 nm. (b) TEM bright-field images of both particle types are shown (scale bar, 100 nm).

Besides the breathing modes, the nanodisks sustain dipolar and multipolar plasmon excitations. These modes are confined to the particle boundaries with a specific number of field antinodes along the particle circumference18. It is thus tempting to ask for their correspondence to the modes on the straight film edge we have analysed in Fig. 3. Accordingly, we assume the disk circumference dπ to determine the possible SP wavenumbers k=2π/λ=2ℓ/d for an angular mode order ℓ. The mode energies extracted from the peaks in the EEL spectra versus the so-defined wavenumbers are plotted in Fig. 4 as blue symbols. We observe that all data points closely follow the dispersion relation experimentally found for the EA mode in Fig. 3c, from where the grey line in Fig. 4 is copied. We thus find that the dipolar and multipolar particle modes can be simply scaled to the plasmon modes of a straight edge. This result might appear unexpected at the first view. Generally, it could be expected that bent edges with a curvature radius as small as 50 nm would yield a dispersion relation different from that of a straight edge. In particular for small particles an influence of the geometry-dependent coupling of charges across the particle should play a role.

To shed light on this point, we apply a simple quasi 1D Drude-like model accounting for electrons confined to a straight or circular wire, as discussed in more detail in Supplementary Fig. 6 and Supplementary Note 3. For the straight wire and a given wavenumber k, the plasmon frequencies ω are given by the expression

where e and m are the electron charge and mass, respectively, n0 is the electron density and ν(k)=2K0(ak) is the Fourier transform of the Coulomb potential, with K0 being the modified Bessel function of zeroth order. We have introduced a transversal extension a for the stripe-like wire to avoid divergences of v(k). For the circular wire with radius R, we can again use Equation (1) but have to replace ν(k) with

where we have introduced polar coordi<span class="Gene">nates (x, y)=R(cosϕ, sinϕ) for the integration along the ring circumference. The enumerator accounts for the SP charge density wave with angular order ℓ, the denomi<span class="Gene">nator gives the distance from one point of the ring (here (R,0)) to the other points (a accounts again for effects of the transversal ring extension). Importantly, the dispersion relations for the straight and circular wire practically coincide, irrespective of the softness parameter a (for details see Supplementary Note 3). The ‘magic’ about the straight and circular wire geometries is that in both cases the charge distribution of the plasmonic modes is homogeneous, and the conserved quantities are linear and angular momentum, respectively. In transformation optics, one can map the problem of a thin metallic slab to a metallic cylinder26. Thus, although seemingly different, the two problems are intimately connected through a conformal mapping. The straight wire and ring SP modes of Equations (1) and (2) appear to be equally connected by some hidden symmetry. We note that things are markedly different for geometries with broken symmetry, as we illustrate with an elliptical nanoparticle with in-plane axes 100 and 200 nm long. In this case the two orthogonal dipole modes have different energies, and the corresponding data points (orange symbols in Fig. 4) do not comply with the edge plasmon dispersion relation.

Symmetric modes

Besides the antisymmetric film and edge modes there exists a whole set of symmetric plasmonic modes for our sample geometries. In this case, however, detailed experimental evidence is lacking due to the increased damping at higher mode energies and the spectral overlaps with the bulk mode (Fig. 1). As the observed data indicates that the mode scaling turns out to be analogous to the antisymmetric case as expected, we study this point in a simulation analysis of the complete SP mode spectrum, as detailed in Supplementary Fig. 7 and Supplementary Note 4.

In summary, we have shown that plasmonic modes in thin-film-based (nano)structures can be scaled to film and edge modes. In particular, the mode dispersion of straight edges remains unchanged when large curvatures in the case of nanodisks are considered. This finding is intimately related to the homogeneous charge distribution in symmetric particles. Our results provide a simple ordering scheme of plasmon modes, fostering an intuitive understanding of the full mode spectrum in plasmonic systems and providing guidelines for their design.

Methods

Fabrication

<span class="Chemical">Silver films 30 or 50 nm thick were thermally evaporated on 15-nm-thick <span class="Chemical">Si3N4 membranes. Lateral structuring was done by electron-beam lithography in a RAITH e-line system using a standard lift-off process with a poly(methylmetacrylate) resist18. One sample was covered by an additional 22-nm-thick silicon dioxide layer deposited on top of the silver structures by electron beam evaporation.

EELS measurements

EEL measurements were performed in a STEM (FEI Tecnai F20) with a monochromated 200 keV electron probe of 0.15 eV energy width (full-width-at-half-maximum). EEL spectra were measured with a high resolution <span class="Chemical">Gatan Imaging Filter equipped with a charge-coupled device camera with 2,048 × 2,048 pixel and an energy dispersion of 0.01 eV per channel. The electron beam was scanned over the region of interest and spectra were acquired at each beam position (pixel)27. The scan step size in the spectrum images varied between 2.9 and 27 nm, the exposure times per EEL spectrum were chosen between 0.05 and 0.1 s. STEM images were acquired with a high-angle annular dark-field detector (Fischione Instruments). The EEL spectra were background subtracted by removing the tail of the zero-loss peak that includes directly transmitted and elastically scattered electrons by using a fitted logarithmic function (‘fitted logarithm tail’ model within Digital Micrograph software (<span class="Chemical">Gatan, USA)).

Simulations

Simulations were done with the MNPBEM toolbox28 that relies on the boundary element method approach29. The dispersion relations in Figs 1b, 3d and 4 were calculated with the transfer matrix method30.

Author contributions

F.-P.S. fabricated the samples, performed the experiments and analysed the data. U.H., H.D. and A.H. performed the simulations and contributed to the data analysis. H.D., J.R.K. and F.H. designed the experiments and supervised the project. All authors contributed to writing the manuscript.

Additional information

How to cite this article: Schmidt, F.-P. et al. Universal dispersion of surface plasmons in flat nanostructures. <span class="Gene">Nat. Commun. 5:3604 doi: 10.1038/ncomms4604 (2014).