Long-range surface plasmon polaritons (LRSPPs) are optical surface waves that propagate along a thin symmetric metal slab or stripe over an appreciable length (centimeters). Vigorous interest in LRSPPs has stimulated a large number of studies over three decades spanning a broad topical landscape. Naturally, a good segment of the literature covers fundamentals such as modal characteristics, excitation, and field enhancement. But a large portion also involves the LRSPP in diverse phenomena, including nonlinear interactions, molecular scattering, fluorescence, surface-enhanced Raman spectroscopy, transmission through opaque metal films and emission extraction, amplification and lasing, surface characterization, metal roughness and islandization, optical interconnects and integrated structures, gratings, thermo-, electro- and magneto-optics, and (bio)chemical sensing. Despite the breadth and depth of the research conducted to date, much remains to be uncovered, and the scope for future investigations is broad. We review the properties of the LRSPP, survey the literature involving this wave, and discuss the prospects for applications. Avenues for further work are suggested.
© 2009 Optical Society of America
Surface plasmon polaritons (SPPs) are transverse magnetic (TM) polarized optical surface waves that propagate, typically, along a metal–dielectric interface [1, 2, 3], henceforth termed the single interface. One purely bound (nonradiative) SPP wave (mode) is supported by this structure, with fields that peak at the interface and decay exponentially away into both media. Its fields are associated with a charge density wave in the metal; so the SPP is a coupled excitation (plasmon polariton).
The single-interface SPP exhibits interesting and useful properties such as an energy asymptote in its dispersion curve, very high surface and bulk sensitivities, and subwavelength confinement near its energy asymptote. But it is also characterized by a high attenuation, especially near the energy asymptote, limiting the scope for applications. For a metal bounded by an ideal dielectric, the attenuation is caused primarily by free-electron scattering in the metal and, at short enough operating wavelengths, by interband transitions. Roughness along the interface causes additional attenuation.
A practical way of reducing the SPP attenuation is to use a thin metal film or stripe bounded on all sides by the same dielectric and operating the structure in the so-called long-range SPP (LRSPP) mode. “Long-range” is taken to mean that the LRSPP attenuation is at least a factor of 2 to 3 lower than that of the single-interface SPP, resulting in propagation over a longer distance. Indeed, attenuation reduction factors, or equivalently, range extension factors, greater than 100 have been demonstrated experimentally with the LRSPP. The range extension mitigates an important limitation of the single-interface SPP, but this comes at the expense of reduced confinement. Nevertheless, the extended range may outweigh the reduced confinement potentially enabling better and more competitive applications.
Our purposes in this paper are to highlight the attributes of the LRSPP, review its origins and the body of work that has been conducted on this wave, and to discuss potential applications based on its use. The LRSPP has appeared in a large number of studies, either as the focal point or peripherally within another context. Literary contributions that include consideration of the LRSPP are cited, and, generally, only the salient points involving the LRSPP are described. In many cases, there is more to a cited contribution than this narrow perspective, so the reference list should not be regarded as listing exclusively LRSPP studies. The discussion is at times broadened to include the single-interface SPP, or the short-range (high-attenuation) modes supported by the structures considered, in order to provide context and comparison or to underline important points of differentiation. The references are organized by subsection, then chronologically.
The literature on SPPs in general is of course much broader. The texts edited by Boardman  and by Agranovich and Mills , and the monograph by Raether , review early work and remain important literary landmarks in the field. Papers by Sambles et al. , Welford , and Barnes , serve as good introductions to SPPs, as does a recent textbook by Maier . Good recent reviews of the field include those of Barnes et al. , Zayats et al. , Maier and Atwater, Ozbay , and Atwater . Ebbesen et al.  and Degiron et al.  have recently discussed integrated optical circuitry based on SPPs.
Work conducted on or with the LRSPP can be organized in many ways. The main division adopted here follows the dimensionality of the guiding structures: Section 2 reviews the LRSPP in the metal slab, which provides confinement in one dimension, and Section 3 reviews the LRSPP in the metal stripe, which provides confinement in two dimensions. Subsections 2.1, 3.1 briefly and didactically summarize the properties of the bound modes in these structures in order to provide context and introduce relevant symbols and terminology. Section 4 discusses the prospects for applications, and Section 5 gives concluding remarks.
Throughout this paper, an time dependence is assumed with modes propagating in the direction according to . The complex propagation constant γ in inverse meters expands as , where α and β are the attenuation and phase constants, respectively. The normalized propagation constant is , where is the phase constant of plane waves in free space, the wavelength in free space, and the speed of light in free space. The complex effective index of a mode is then given by , where is the effective index and the normalized attenuation. The mode power attenuation (MPA) in decibels per meter is given by .
The propagation length is the distance from the launch point where the mode power decays by a factor of and is given by . (A less used definition of the propagation length, the power decay length, is also found in the literature and yields a larger length for the same attenuation.)
The group velocity of a mode is given by , and its lifetime is given by .
The relative permittivity is related to the permittivity ε in the usual way, , where is the permittivity of free space. The relative permittivity is related to the optical parameters n, k via . A numerical subscript i identifies these quantities in a specific region. The relative permittivity of a metal at optical wavelengths is decomposed into real and imaginary parts as . For good metals in the infrared and over at least part of the visible (e.g., Au, Ag, Cu, and Al).
2. Metal Slab
2.1. Modes of the Metal Slab
The single interface is sketched in Fig. 1(a), and the structures sketched in Figs. 1(b), 1(c) are important variations. The variation shown as Fig. 1(b) consists of a thin metal film of thickness t and relative permittivity bounded by optically semi-infinite dielectrics (claddings) of relative permittivity and , henceforth referred to as the metal slab. Another variation consists of a thin dielectric film of thickness t and relative permittivity bounded by optically semi-infinite metals of relative permittivity , as sketched in Fig. 1(c), and henceforth referred to as the metal clads. The metal slab is said to be symmetric when and asymmetric otherwise ( or ). The distribution of the main transverse electric field component of the modes is sketched as the red curves over the cross section of each structure. Mode propagation occurs along the axis (upward perpendicular to the page).
In the symmetric metal slab, the bound single-interface SPPs supported by the individual metal–dielectric interfaces at large t, couple as t is reduced, forming two TM-polarized bound supermodes, sometimes termed coupled modes, that exhibit distinct dispersion characteristics and a distinct evolution with structure parameters (ε and t). These supermodes are denoted herein as and for asymmetric bound and symmetric bound, respectively, since their main transverse electric field component varies either asymmetrically or symmetrically across the structure (along the y axis), as sketched in Fig. 1(b). The field component of the modes has the same symmetry as the corresponding field, but the longitudinal electric field has the opposite symmetry. The charge density in the metal linked to the mode has a symmetric distribution over t, as indicated by the pluses in Fig. 1(b), whereas the charge density associated with the mode is asymmetric over t, as indicated by the plus and minus signs.
The LRSPP is the mode of the thin metal slab.
In a symmetric structure with lossless claddings, the attenuation and effective index of the mode decrease smoothly as t is reduced, with its mode fields increasingly expelled from the metal film and penetrating more deeply into the claddings. Indeed, as , the confinement and attenuation of this mode vanish as it evolves smoothly into the vertically polarized TEM (transverse electromagnetic) wave of the background. The LRSPP is the mode of the thin metal slab. The mode exhibits increasing confinement and penetration into the metal with decreasing t and, consequently, increasing attenuation. For large t the and modes are degenerate with the single-interface SPPs supported by the uncoupled top and bottom metal–dielectric interfaces.
The trends are similar in an asymmetric slab except that (i) the mode cuts off below a certain thickness t that depends on the permittivities and operating wavelength, and (ii) with increasing t the mode evolves into the SPP supported by the metal interface with the high-index cladding, while the mode evolves into the SPP at the interface with the low-index cladding. These trends with t are apparent from Fig. 2, which plots the effective index and normalized attenuation of the and modes in an asymmetric slab at assuming Ag for the metal and claddings of relative permittivity and . The symmetric structure may be more convenient for working with the LRSPP because the mode remains nonradiative (purely bound) for , whereas in an asymmetric structure the LRSPP ( mode) remains nonradiative only for t greater than a cutoff thickness.
The and modes exhibit distinct dispersion characteristics if t is small enough. For example, Fig. 3 shows normalized dispersion curves for the and modes propagating along a lossless metal film modeled as a Drude metal and bounded symmetrically by semi-infinite vacuum for three normalized thicknesses. For very thin metal films the dispersion of the and modes are distinct, but for thick films they merge and coincide with the corresponding single-interface SPP. The dispersion curves are asymptotic with the light line as , and they are asymptotic with as , where is the plasma frequency of the Drude model. The limit is termed the SPP energy asymptote. The dispersion curve of the mode is above that of the mode, so a given β occurs at a higher ω for this mode. The confinement increases, tends to zero, and the optical density of states diverges near the energy asymptote. A region of negative occurs for the mode before the energy asymptote. The trends are the same for structures having symmetric dielectric claddings, with the asymptotic limits taking on appropriate values in the dielectric [light line, ; energy, ]. The trends are similar for real metals, except that the attenuation of the modes increases dramatically toward the energy asymptote, and the depth of the asymptote is limited by bendback. The mode is generally long range away from the energy asymptote (and for small t).
The nomenclature used for identifying the and modes varies greatly throughout the literature, and occasionally erroneous assignments are made. The nomenclature that we have adopted follows the integrated optics convention of identifying a mode by features in the spatial distribution of its main transverse electric field ( in the structures of Fig. 1). Other mode nomenclatures are based on identifying features in the longitudinal electric field, in the charge distribution across the metal, or on the mode’s location on a dispersion diagram. Thus, in the literature, one finds the mode sometimes termed the asymmetric mode based on or on charge distribution (Fig. 1(b)), or the mode because it corresponds to the highest curve on a dispersion diagram (Fig. 3). Likewise, the mode has also been termed the symmetric mode or the mode.
The LRSPP can be excited by p-polarized light by using a high-index prism as sketched in longitudinal cross-sectional view in Fig. 4 (case 1, ). In typical experiments, the angle of incidence θ is varied beyond the critical angle of the interface, and the reflected power is monitored by using a detector. The excitation of a mode occurs when the in-plane wavenumber of the incident light equals the mode’s wavenumber β, resulting in a drop of the reflected power. If the reflectance is plotted versus θ, then dips appear for each excited mode as shown in Fig. 5. The mode is excited at a smaller θ than the mode, since the former has a smaller β at the operating ω (Figs. 2, 3). Lower attenuation results in a narrower dip, a signature of the LRSPP. As t increases, the and dips merge into one, because the modes become degenerate and identical to the single-interface SPP. The plot of Fig. 5 is termed an attenuated total reflection (ATR) angular spectrum, since the excitation of modes results in attenuation of the reflected light as measured by the detector.
The LRSPP can also be excited by a TM-polarized optical beam via end-fire coupling as sketched in longitudinal cross-sectional view in Fig. 6. In this arrangement the beam is focused onto the end facet of the structure such that it overlaps well with the fields of the LRSPP. End-fire coupling can be efficient, since the transverse mode fields of the LRSPP are symmetrically distributed over the structure cross section, as are the exciting fields in such arrangements. The technique can be easier to implement than prism coupling, and it eliminates problems associated with prism loading (mode perturbation and unwanted outcoupling); however, it requires access to high-quality end facets, which can be difficult or inconvenient to create in some structures, and it is not β selective; thus all modes that overlap the input beam will be excited to some extent (including radiative modes); so outputs must be interpreted carefully.
The metal clads depicted in Fig. 1(c) also support coupled modes, such as the symmetric mode for which the main transverse electric field distribution is also sketched. Table 1 summarizes modal information for the three structures of Fig. 1, assuming , Ag for the metals, for the dielectrics, for the metal slab, and for the metal clads; is the mode field width. The trade-off between confinement and attenuation across these structures is evident from these data. The LRSPP in the metal slab is at one end of the trade-off, having a low attenuation but also low confinement (smaller and larger mode size ), whereas the symmetric mode of the metal clads is at the other end of the trade-off (higher attenuation and confinement). The single-interface SPP is between them. Neither the metal clads nor the single interface support LRSPPs.
2.2. Origins of the LRSPP [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36]
The LRSPP builds on a number of previous studies involving the single interface and the metal slab, including, for example, the work of Ritchie , Kliewer and Fuchs , Otto [18, 19], Economou , Kretschmann , and Abelès and Lopez-Rios .
The coupled modes in the symmetric metal slab were studied long ago through dispersion computations, for example, such as those produced by Kliewer and Fuchs  for an ideal thin metal film bounded by vacuum (Fig. 3) and by Economou  for a similar structure and a number of variants thereof. These early studies did not include the effects of damping; so the range of mode propagation could not be assessed.
Otto described and demonstrated prism coupling as a means of exciting the single-interface SPP  (the Otto configuration), by introducing a low-index gap between the prism and the metal film (Fig. 4, case 2). Otto also considered, theoretically, prism coupling to the and modes of the metal slab, including damping in the metal, bounded symmetrically by low-index claddings and high-index prisms , and computed the reflection and transmission responses of the system. Interestingly, in Fig. 3 of , he showed a narrow linewidth in the reflection and transmission resonances of the system mediated by the mode ( in the notation of ) at for a Ag slab thick and for a prism-metal spacing of (the largest considered). He also showed this resonance narrowing as the prism-metal spacing increases, a behavior later understood to be characteristic of the LRSPP (e.g., ). He showed theoretically that the transmittance of the system could be very high, 70% in the case of a thick Ag film, and that the transmission was mediated primarily by the mode. (Transmission through a Ag film of this thickness bounded by air is essentially zero.) He suggested using the system in transmission as a polarizer. The arrangement was studied experimentally shortly thereafter , confirming the main findings of Otto , and exploring its possible application as a polarizing spectral filter. (Dragila et al. studied a similar system ).
During the same time, Tien et al.  proposed and demonstrated the prism coupling approach to excite dielectric waveguide modes, and Kretschmann  demonstrated a variant of the Otto configuration  (the Kretschmann–Raether configuration) for exciting the single-interface SPP where the metal film is deposited directly onto the base of the prism (Fig. 4, case 3). (This configuration was also considered by Turbadar years earlier , as pointed out by Welford .) The Otto  and Kretschmann–Raether  configurations are not suitable for exciting LRSPPs, but Otto’s other configuration  is suitable.
Abelès and Lopez-Rios  combined the Otto and Kretschmann–Raether configurations (Fig. 4, case 4) to excite the (essentially) uncoupled single-interface SPPs supported at the opposite interfaces of a highly asymmetric thin metal slab.
Although much had been learned regarding the coupled modes of the metal slab, their attenuation (and range) had remained unexplored until the work of Kovacs [26, 28]. Kovacs computed the propagation constant of the and modes at in In slabs of thickness , bounded symmetrically by semi-infinite (pp. 97–99 of  or Table 1 of ). His computations show that and α of the mode are larger than those of the mode and that they move further apart with decreasing t. Specifically, his results show and α of the mode both increasing with decreasing t, and and α of the mode both decreasing with decreasing t. From his computations, it is also noted that α of the mode for is lower than for . Prism coupling to the coupled modes was also investigated theoretically and experimentally for different In thicknesses .
Kovacs also explored theoretically and experimentally prism coupling to the and modes of thin Ag slabs bounded by identical dielectrics (Fig. 4, case 1; Fig. 5) [26, 27]. The experimental structures comprised an Ag slab thick bounded by cryolite and were excited at by angle scanned prism coupling. These structures allowed the experimenters to unambiguously differentiate and identify the and modes, observing that their angular separation (in θ) decreased with increasing t as expected. Informative plots of the Poynting vector through the structure, and of the current density in the metal, were also given, from which one notes, for instance, that the power flow in the metal has a large component in the direction antiparallel to the direction of modal propagation . Effects caused by the finite thickness of the claddings were also investigated, revealing the importance of the gap s (Fig. 4, case 1) on the coupling efficiency of the modes. The mode was not really long range in these structures because of the large thickness of the metal film, the roughness of the interfaces, and inhomogeneities in the bounding cryolite films.
Thus, many essential features of the LRSPP were uncovered by Kovacs [26, 28]. And by making slight modifications to then-known prism coupling techniques [18, 19, 21, 23, 25], particularly by removing the second prism in Otto’s other geometry , he demonstrated experimentally via optical means the existence of the coupled modes in symmetric metal slabs [26, 27]. (Many of his results  were subsequently summarized in Chap. 4 of .)
At about the same time, Fukui et al.  (in collaboration with G. I. Stegeman) included damping in the metal and computed the lifetime τ of the and modes in thin unsupported Ag films bounded symmetrically by vacuum as a function of the Ag film thickness t. They predicted a propagation length of for the mode on a thick Ag film in vacuum at . The propagation length of the corresponding single-interface SPP is ; so the predicted range extension was a factor of 50. They also showed the mode guided as , with no cutoff thickness being apparent in their results.
Sarid [30, 31] then studied theoretically the and modes in the metal slab bounded by asymmetric dielectrics (, ). From his results, reproduced here as Fig. 2, it is apparent that the LRSPP does not exist below a cutoff thickness, in this case. He modeled the prism coupling arrangement (Fig. 4, case 1), noting that the half-width angle corresponding to the excitation of the LRSPP was very narrow, about 0.004°. He also commented on the criticality of the gap s, pointing out that if it is too large the coupling efficiency suffers, but if it is too small then the prism loads and wipes out the LRSPP. The terminology “long-range” for identifying the low-loss version of the mode originated here.
The structure explored experimentally by Kuwamura et al.  consisted of a thick Ag film deposited on with index-matching oil used as the other cladding. The structure was excited at in a prism-coupled configuration. The measured propagation length of the mode in this structure is approximately longer than that of the SPP on the corresponding single interface.
Craig et al.  reported a propagation length of for the mode supported by a thick Ag film, bounded by glass on one side and index-matching oil on the other, and excited at by using a retroreflecting prism . This measured propagation length is longer than that of the corresponding single-interface SPP.
Quail et al.  characterized a thick Ag film and a thick Al film, both on glass and covered by index-matching oil, at in a prism-coupling configuration. They measured propagation lengths for the mode that are about longer than the corresponding single-interface SPPs.
Dohi et al.  measured a propagation length of for the mode supported by an Ag film about thick on Pyrex and covered with oil having a slightly different index than that of the Pyrex (i.e., in a slightly asymmetric structure ). Their structures were excited at in a prism coupling arrangement. Their range extension is comparable with that achieved by Craig et al. .
2.3. Prism Coupling and Field Enhancement [37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53]
Difficulties with prism coupling to the LRSPP (Fig. 4) were identified early on [27, 30, 31, 32, 33, 34, 35, 36]. Experimental difficulties include controlling the gap spacing s and parallelism between the base of the prism and the metal slab, as well as the angle of incidence of the input light θ, while ensuring that the permittivities of both claddings remain closely matched . These requirements suggest the need for high-performance optomechanical components and working with index-matching oils or custom-fabricated structures.
Another difficulty is that the mode fields of the LRSPP extend deeply into the claddings; so the loading effect of the higher-index prism on the mode is quickly apparent, perturbing its fields and propagation characteristics and rendering it radiative into the prism, as is readily apparent from the plane wave computations of Wendler and Haupt . If the prism is too far from the metal slab, the coupling efficiency is poor, but if it is too close, loading destroys the LRSPP (cf. Figs. 1 and 7 of ). Under conditions of perfect coupling to the LRSPP, radiative damping equals intrinsic damping (damping without the prism), so its range is halved (i.e., its attenuation is doubled) [40, 48], as is the case for the single-interface SPP in the Otto [18, 37] and Kretschmann–Raether [3, 23] configurations.
Thus, the reflected field in a prism coupling experiment includes a specularly reflected contribution as well as an outcoupled (reradiated) contribution due to the propagating LRSPP. Since the LRSPP has a long propagation length, the outcoupled contribution extends over a large spatial cross section and can interfere with the specular reflection. The contributions can, however, be distinguished if a finite-size input beam is used [41, 42, 45] (as for the single-interface SPP  in the Otto  and Kretschmann–Raether  configurations). For example, Fig. 7(a) shows a finite-width input beam prism coupled to the LRSPP on a Ag film supported by a Pyrex substrate and covered by index-matching oil, and Fig. 7(b) shows the measured intensity profile at the observation plane (, , , , , ) . From Fig. 7(b) two bumps separated by a null are noted, the leftmost one due to specular reflection of the input beam and the rightmost one due to the outcoupled LRSPP. In such a situation, the measured profile of the outcoupled contribution follows the decay of the LRSPP [42, 45]. The importance of considering the finite size of the input beam and of its angular spread when interpreting prism coupling experiments with LRSPPs has also been highlighted in other papers [46, 49, 50, 51, 53].
Barnes and Sambles  excited the LRSPP via prism coupling in a thick Ag slab bounded symmetrically by Langmuir–Blodgett layers (22-tricosenoic acid) about thick at . They achieved a modest increase in propagation length over the corresponding single-interface SPP, due in part to the large thickness of the Ag film, to loss in the Langmuir–Blodgett layers, and to suspected damage caused to one of the Langmuir–Blodgett layers during Ag evaporation. However, it is interesting to note that a multilayer structure incorporating organic claddings could be fabricated with some success.
The LRSPP was also seemingly observed in a thin Ag film bounded on both sides by thin layers of Teflon and excited by using prism coupling, presumably at .
Sarid et al.  computed the field enhancement in a prism coupling arrangement at assuming Ag films 70 and thick bounded by dielectrics of index 1.5. The field enhancement was defined as the ratio of the squared magnitude of the magnetic field in the vicinity of the metal surface to the squared magnitude of the magnetic field in the prism. They reported field enhancements of 40 and 600 at the excitation angle of the mode in the 70 and thick films, respectively, for near perfect coupling into the modes. The mode in the thick case is not long-range and corresponds almost to the single-interface SPP in the Otto geometry  (distinct and modes still occur at this thickness); the mode in the thick case corresponds to the LRSPP. They mentioned computing the field enhancement in the Kretschmann–Raether geometry and observing that they are similar to the Otto case. Thus, they concluded that the field enhancement associated with the excitation of the LRSPP via prism coupling can be 1 order of magnitude larger than that in the Otto and Kretschmann–Raether geometries . Indeed, assuming that Eq. (10) of  holds, then the field enhancement of the LRSPP relative to the single-interface SPP should follow the ratio of their propagation lengths (’s) or the inverse ratio of their attenuations (α’s). LRSPP field enhancements were later shown to exist on corrugated gratings as well (e.g., ).
Lévy et al.  compared theoretically the prism-coupled LRSPP field enhancement to that of the prism-coupled and modes supported by a dielectric slab waveguide (equivalently, a 1D cavity). Results for the structures investigated show that the field enhancement of the and modes along the core center was about greater than the enhancement of the LRSPP fields.
This enhancement [39, 40, 44] might loosely be termed an intensity enhancement, since another definition for the field enhancement found in the literature forms the ratio of the field magnitudes, yielding the square root of the former.
Symmetric periodically corrugated metal slabs supporting LRSPPs have also been explored as an alternative to prism coupling and as structures having rich and interesting properties in their own right. Generally (TM) p-polarized light is used to illuminate a corrugated grating, with the plane of incidence oriented perpendicular to the grooves, and the response in reflection and/or transmission recorded as a function of the angle of incidence (angle scan) or the wavelength of operation (wavelength scan).
Different theoretical approaches have been used to model grating responses (cf. Barnes et al.  for a good discussion), many rooted in Chandezon’s coordinate transformation formalism , whereby corrugated boundaries are mapped to flat surfaces allowing the straightforward application of boundary conditions and transfer matrices for handling many layers. Cotter et al.  improved on the method by using scattering matrices instead of transfer matrices to eliminate numerical instability when modeling thick structures.
Inagaki et al.  reported the excitation of the and modes at in a thin free-standing thick corrugated Ag slab in air, with the measured resonance width of the mode being considerably narrower than that of the mode. Further experiments were conducted by the same group  on free-standing corrugated Ag slabs having thicknesses in the range of , showing splitting of the and modes for a slab thickness of about . The measured phase and attenuation constants of these modes compared qualitatively well with theoretical expectations for a flat slab over the thickness range investigated, but larger phase and attenuation constants were measured throughout. The coupling efficiency was also found to be strongly dependant on the thickness of the slab.
This symmetric free-standing corrugated slab was subsequently studied theoretically by Dutta Gupta et al. , supporting the main conclusions of Inagaki et al. [55, 57]. They also study the effects of damping by outcoupling to free radiation for various corrugation amplitudes.
Cavalcante et al.  studied theoretically the metal slab bounded symmetrically by dielectrics where one of the metal–dielectric interfaces is smooth and the other interface has a sinusoidal profile. They computed the reflectivity and the intensity of the field near the metal slab, both as a function of the angle of incidence for various grating amplitudes. They observed, among other points, an increase in the linewidth of the reflectivity dip for both the and modes, indicating increased attenuation, as the corrugation amplitude increases.
The study of Chen and Simon  also pointed out that an additional loss contribution to the LRSPP in corrugated gratings may be expected that is due to scattering from the grating grooves.
Bryan-Brown et al.  measured the reflection response of symmetrically cladded thin corrugated Ag slabs through the LRSPP excitation angle and deduced the grating groove depth as well as the thickness and optical parameters of the embedded Ag films via comparisons with a theoretical model of the experiment.
Salakhutdinov et al.  investigated sinusoidally corrugated metal slabs where the upper and lower surfaces of a slab are corrugated in phase or phase shifted by π. They computed the perturbed propagation constant of the LRSPP for corrugations of both types in a Cu film, as a function of groove depth and t. They found that the corrugations increase the attenuation of the mode, but much more so for π-shifted corrugations than for in-phase corrugations. They found that increases for the π-shifted corrugations but that it either increases or decreases for in-phase corrugations depending on t. They modeled and demonstrated experimentally anomalous reflection (interference between reflected and outcoupled guided fields) from in-phase corrugations.
Hooper and Sambles  reported a theoretical study pertaining to the excitation and nature of the modes in unsupported (air on both sides) corrugated thin Ag slabs where both surfaces of the slab were corrugated either conformally (identical profiles, in phase) or nonconformally (identical profiles, phase shifted). Sinusoidal profiles with and without the first harmonic in grating wavenumber were explored. The reflection, transmission, and absorption of incident radiation were determined for a variety of structure parameters. Four coupled modes (two -like and two -like) were found in the nonconformal case that incorporates the first harmonic in grating wavenumber, which are due to new symmetries for the charge distribution that are allowed by the grating features. The second harmonic also leads to anticrossings and bandgaps in the dispersion diagram. This grating architecture for small t is capable of an unusually large transmittance over a wide frequency band, almost independently of the angle of incidence, with the transmittance mediated by one of the modes.
Chen et al.  further explored this grating concept experimentally and theoretically by way of a thick Ag slab deposited conformally onto a corrugated substrate, covered with index-matching fluid, and clamped with a pair of prisms (one on each side of the slab). A sinusoidal grating profile with the first harmonic in grating wavenumber was implemented. The origin of anticrossings and bandgaps was discussed in terms of structure parameters and the role of the -like and two -like modes.
Lévêque and Martin  investigated theoretically the excitation of the LRSPP by a Gaussian beam normally incident onto gratings consisting of either periodic rectangular grooves etched into the top surface of an Au slab, or periodic rectangular Au protrusions deposited onto the top surface of an Au slab. The Au slab was assumed to be free standing (bounded by vacuum), thick, and the groove depth or protrusion height in the range of about . The length of the gratings was set to 5 periods, the width of the incident Gaussian beam matched the grating length, and the other grating parameters were varied (period, duty cycle, and groove depth or protrusion height). Lévêque and Martin predicted an optimal coupling efficiency of 33% from the Gaussian beam into LRSPPs propagating in both directions along the slab (i.e., along ) at , using a grooved grating. (They also investigated coupling into the single-interface SPP, using similar grating structures.)
Sellai and Elzain  computed the total p-polarized reflectivity from gratings formed by modulating the thickness of the metal slab, predicting sharp reflection dips at specific operating wavelengths where the incident light is efficiently coupled into the and modes of the structure.
Agarwal  studied a different kind of structure, consisting of a smooth metal slab bounded by dielectrics of the same index but with the top dielectric being of finite thickness and having a sinusoidal corrugation applied to its top surface. He showed that the excitation of the LRSPP under plane wave illumination leads to a field enhancement near the grating surface.
A less studied geometry for exciting corrugated gratings is where the plane of incidence is rotated 90° in the plane of the grating (azimuthally) such that it is parallel to the grating grooves. This case was investigated experimentally and theoretically by Chen et al.  for a thick Ag slab deposited conformally onto a sinusoidally corrugated substrate, covered with index-matching fluid, and clamped with prisms (a structure similar to those investigated in  but without the harmonic in the grating vector). They measured a coupling efficiency of about 24% into the mode and predicted strong coupling into the LRSPP but under different design and excitation conditions.
Korovin  developed a formulation for modeling multilayer corrugated gratings based on a curvilinear coordinate transformation, but, in contrast to Chandezon , he used an established solution to Maxwell’s equations in Cartesian coordinates having eigenvalues that are determined analytically, thus improving on the efficiency and accuracy of the method. He then applied this formulation to model the reflectance from a thick Au film bounded by glass ( thick on one side, semi-infinite on the other) with all interfaces conformally corrugated. Reflectance dips due to coupling to the and modes are apparent from his computations.
2.5. Modal Studies [71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104]
A good number of studies have been reported concerning the modes supported by metal slabs [Fig. 1(b)] and variations thereof. The modes are typically obtained by solving a suitably defined 1D boundary value problem based on the wave equations written in the frequency domain for lossy isotropic media. Necessary boundary conditions are applied between layers, and modes are found as solutions to a characteristic (transcendental) equation much in the same way as for 1D dielectric waveguides. The permittivity of the metal layer is usually obtained from the Drude model of the metal or from measurements. Additional 1D layers (metallic or dielectric) can be included in the analysis, and transfer matrices can be developed for the stack (e.g., ). The solution approach typically does not involve the discretization of spatial variables, and although the transcendental equation is solved numerically, the approach is essentially analytical.
Tomaš and Lenac  derived expressions to estimate the damping properties (lifetime and propagation length) of the LRSPP in thin unsupported metal slabs at long wavelengths.
The study by Stegeman et al. on asymmetric metal slabs  included consideration of leaky modes in addition to the bound modes supported by the structure. The mode nomenclature, or , originated here. A cutoff thickness exists in their asymmetric structures below which the LRSPP no longer propagates as a purely bound mode, corroborating the results of Sarid .
The follow-up study by Burke et al.  reports a comprehensive treatment of the structure and its modes. The bound modes in symmetric structures, and the bound and leaky modes in asymmetric structures, are discussed in detail, and their evolution over good ranges of structure parameters and operating wavelength are given. A physical interpretation is given to the leaky modes.
Burton and Cassidy  presented modal results for the metal slab at for various metals and cladding indices near silica, as a function of slab thickness and asymmetry.
Smith et al.  also investigated the metal slab, emphasizing the power dissipation spectra of a dipole into the modes supported, and modal dispersion.
The excitation of the modes was also considered by Burke et al. , and a suggestion was made to use end-fire coupling, as sketched in Fig. 6, as the means to excite the modes (following , where Stegeman et al. computed high excitation efficiencies for the single-interface SPP). This coupling technique was used by Vaicikauskas  to obtain Fourier-transform infrared spectra of the LRSPP in a multilayer 1D structure.
Stegeman and Burke  explored the bound modes of two coupled thin metal slabs as a function of their separation for a few refractive index values of the intervening dielectric. They found supermodes of the system consisting of symmetrical and asymmetrical couplings of the mode, and of symmetrical and asymmetrical couplings of the mode. The supermodes formed from coupled modes can be long-range. Yoon et al.  computed the dispersion of the symmetrically coupled modes supported by the same structure as a function of the separation between the thin metal slabs for an intervening dielectric that has a higher or lower index compared with the claddings. They found for a high-index intervening dielectric that the dispersion curves insect a common point for all separations between the metal slabs. At the intersection frequency, the effective index of the mode is independent of separation and equal to the index of the intervening dielectric, the propagation length varies linearly with separation, the transverse mode fields have a constant magnitude between the metal slabs (i.e., they are uniform), and the mode fields normalized to the separation are invariant.
Economou  also studied the modes of two and three coupled thin metal slabs, but without loss in the metals. Avrutsky et al.  investigated the supermodes supported by a multilayer stack constructed from (five) alternating thin metal and (six) dielectric slabs and found (among other modes) one long-range supermode consisting of symmetrical couplings of modes. Davis  investigated a similar structure and found two long-range supermodes for their stack (three metal and four dielectric layers), symmetrically and asymmetrically distributed.
Stegeman and Burke  investigated the effects of one air gap on the propagation of the LRSPP. In this study, the gap was positioned along the top surface of the metal film and various parameters of the system were altered. In general, the effects of the air gap were found to be deleterious, causing the LRSPP to become cut off for very thin gaps . Methods for mitigating air gaps include using index-matched fluids to fill them  or carefully controlling the thickness of the gap while compensating for its effects by using a high-index layer deposited on the opposite surface to form an index-matched effective medium .
Lenac and Tomaš [77, 78, 79] studied theoretically the metal–dielectric slab system and found that the modes of the system were hybridized or coupled modes of the individual constitutive structures. They found that the LRSPP was very sensitive to the presence of a thin dielectric layer, leading to increased damping. They also studied the frequency dispersion of the system applying Drude [77, 78, 79] and dielectric  dispersion models to the metal and dielectric layers, respectively. Yang et al.  studied a similar system.
Kou and Tamir  proposed placing a high-index dielectric layer within one of the claddings near the metal film, essentially forming a dielectric waveguide therein that couples with the metal slab. When the structural parameters are properly selected, one of the coupled modes of the system retains the character of the LRSPP but has a longer propagation length. Prism coupling to this structure (among others) was discussed in [49, 51].
Guo and Adato [95, 97] proposed to extend the range of the LRSPP in the metal slab by introducing thin low-index layers on both sides of the metal film with high-index outer claddings bounding the system and operating the LRSPP near cutoff, where it has lower loss. In , they also investigated the effects of placing thin high-index layers alongside the film, finding that they increase both the confinement and attenuation of the LRSPP. Durfee et al.  investigated similar structures and, additionally, the case where only one low-index layer is placed alongside the metal film. Adding low-index layers in this manner [95, 97, 101] allows a thicker metal slab to be used with lower attenuation compared with a corresponding conventional symmetric slab, but care must be taken to operate sufficiently far from cutoff to avoid radiation loss ( of outer claddings). As with the LRSPP in the conventional symmetric structure, consideration to mode size and coupling efficiency must also be given.
Yun et al.  investigated theoretically the modes supported by a buried rectangular dielectric waveguide with a thin metal slab bisecting the structure along its horizontal plane of symmetry. The modes investigated are similar in character to the and modes of the metal slab but perturbed by the high-index dielectric core which provides strong horizontal and vertical confinement. Propagation lengths and mode areas are reported revealing that the -like mode is capable of long-range propagation when the metal slab is thin enough.
Stegeman  studied the metal slab bounded symmetrically by birefringent media , finding that the effective index of the LRSPP was determined primarily by the perpendicular (to the metal slab) refractive index and that birefringence was essentially negligible. Mihalache et al.  investigated the modes of thin metal films on a uniaxial substrate having its optic axis in the plane of the film, covered with an isotropic cladding. They computed mode cutoff thicknesses assuming quartz as one cladding (, ), glass as the other, various metals, and various orientations of the optic axis.
Wendler and Haupt  studied the evolution of the LRSPP as a function of structure asymmetry , showing that the LRSPP cuts off in an asymmetric structure, and that as the LRSPP nears cutoff, in agreement with other studies (e.g., [30, 35, 80]). They proposed operating the LRSPP near cutoff, suggesting that an increase in of up to 3 orders of magnitude over the LRSPP in a symmetric structure is achievable. However, in subsequent work , they considered the excitation of the LRSPP near cutoff by plane waves assuming prism coupling, and they showed that prism loading limited the increase in to a factor of about 3 over the LRSPP in the corresponding prism-loaded symmetric case, and to a factor of about 2 over the LRSPP in the corresponding unloaded symmetric case (Fig. 11 of ). Interestingly, the measurements of Dohi et al. , obtained at the same operating wavelength in a system similar to that modeled , are in good agreement with these prism-loaded computations: an increase in by a factor of about 3 is apparent for the LRSPP near cutoff compared with the LRSPP of their symmetric case (Fig. 1 of ).
Breukelaar and Berini  modeled a section of asymmetric metal slab placed between two butt-coupled symmetric ones , and they computed the insertion loss through the system for end-fire excitation while varying the asymmetry from none through well beyond cutoff of the LRSPP. Radiation spreading through the asymmetric portion was modeled via normal mode decomposition by discretizing the radiation continuum into an appropriate orthonormal basis. The lowest insertion losses were always obtained for no asymmetry . This is due to the fact that the end-fire coupling efficiency into the LRSPP decreases as the mode nears cutoff because its fields expand more deeply into the higher-index cladding. Indeed, at cutoff the mode fields extend infinitely into the higher-index cladding, and the end-fire coupling efficiency is zero. Thus, while it is true that near cutoff , this range extension is not readily accessible owing to difficulties in coupling the mode to sources [48, 94]. At present, the literature suggests no extension , or extension by a factor of about 3 [35, 48], for the LRSPP in asymmetric structures compared with in symmetric ones.
Zervas  presented computations suggesting that a bound LRSPP could be supported by a thin metal slab in a highly asymmetric structure, well beyond the cutoff point (see the curve labeled in region I of Figs. 2 and 3 of ). However, one notes from his computations that the effective index of the mode in the highly asymmetric region is far below the index of one of the claddings ; so the mode is radiative. A subsequent study by Liu et al.  ignored the effective index and this consideration, focusing only on the attenuation constant of the mode.
Tournois and Laude  found negative group velocities for the mode in the lossless metal slab when becomes similar to . For metals this condition occurs near the energy asymptote, as is apparent from the computations of Kliewer and Fuchs (Fig. 3 ).
As mentioned in Subsection 2.1, the confinement and attenuation of the LRSPP vanish together as the thickness of the metal film is reduced, leading to a trade-off between these two fundamental mode properties (this is apparent from many of the modal studies conducted for the slab, e.g., ). Zia et al.  reported computations illustrating this trade-off for the metal slab and the metal clads. Berini  proposed three figures of merit (FoMs) for 1D waveguides, defined as benefit-to-cost ratios, where three different confinement measures were used as the benefit and attenuation was used as the cost. Closed-form expressions of the FoMs were derived for the single-interface SPP. The FoMs and the quality factor (Q) were used to assess and compare the single interface, the metal slab, and the metal clads, implemented by using Ag and , over ’s ranging from the infrared to the SPP energy asymptote. Preferred ’s emerged depending on which FoM was used and thus on how confinement was measured. The largest FoMs were obtained for the LRSPP, but they can also be large for the symmetric mode in the metal clads. Q’s of about 10,000 were found for the LRSPP.
Al-Bader and Imtaar studied the bound , leaky , and bound hybrid  modes of cylindrical structures, comprising a metal film wrapped around a core and surrounded by a cladding. In the case where the cladding is index matched to the core, a long-range mode similar to the mode in the slab is supported . They showed that some of the modes in this structure evolve into those of the corresponding metal slab at large radii.
There are few studies of the effects of roughness on the propagation of the LRSPP, due largely to the arduous task of producing an accurate theoretical treatment. The effects of roughness are that it modifies the LRSPP’s propagation characteristics, outcoupling (scattering) it into free radiation. These effects are of considerable importance to the LRSPP, especially as attempts are made to produce lower-loss waveguides.
Farias and Maradudin  computed the effect of roughness on the propagation length of the mode supported by a rough Ag film of variable thickness and bounded by vacuum at the optical frequencies of and , where is the plasma frequency of Ag. The Ag roughness profiles were taken as identical along both surfaces and characterized by an RMS (root mean squared) deviation of and correlation lengths of 50 and . They found, as expected, that the propagation length is reduced by roughness, but not very significantly for their chosen (typical) roughness parameters. Extrapolating their results down to a Ag thickness, one finds that roughness reduces the propagation length by about 10% and 4% at and , respectively, compared with the corresponding smooth films.
Paulick  constructed a smooth surface model to account for roughness, by introducing an equivalent surface current incorporating one empirical roughness parameter. Results generated by the model were then compared with the experiments of Inagaki et al. [55, 57], yielding encouraging agreement for the dispersion characteristics of the and modes. He also considered enhanced transmission through an otherwise opaque metal film, in the spirit of [19, 22, 43], but using corrugations instead of prisms to couple the and modes with free radiation. He found strong transmission at normal incidence through a corrugated free-standing thick Ag film and demonstrates that the mode mediates the transmission, more so than the mode.
Modeling roughness as sinusoidal corrugation [56, 58, 59, 60] simplifies the theoretical treatment and seems reasonable for generating approximate estimates (or better) of the effects of roughness on the propagation of the LRSPP and other modes. However, much work remains to be done to determine the range of validity, and indeed the appropriateness, of such models for the LRSPP subject to real roughness profiles.
2.7. Islandized Metal Films [107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122]
It is natural to try very thin metal slabs in order to increase the propagation length. However, as the thickness of a metal film decreases below a threshold value, its (as-deposited, unannealed) optical parameters begin to differ [108, 110, 112, 113, 114, 115] from the bulk values [109, 111]. The main cause is that the volume mass density of metal films typically decreases from the bulk as the thickness is reduced beyond a threshold value because of the formation of voids. Further reduction in thickness eventually leads to a discontinuous (islandized) film. The volume mass density can start decreasing at a thickness of , and as it does, decreases while increases. This trend for the permittivity persists as the film turns into a collection of islands if the film is treated as an equivalent continuous plane parallel film and its thickness and permittivity are interpreted as effective optical values. For an islandized metal film in air (or vacuum), as , goes through 0, changes sign, and hits a peak value before approaching , whereas hits a peak value then approaches 0; i.e., as (e.g., ). If unmitigated, the main implications of this for the LRSPP are that its attenuation does not drop as rapidly with thickness as suggested by using the bulk permittivity of the metal, and a thickness is eventually reached where the metal permittivity is such that the LRSPP is no longer supported. Annealing treatments postdeposition have a strong impact on the microstructure and thus can significantly improve the quality of a thin metal film , so the density problem can be mitigated to some extent. Figure 8 shows trends for the resistivity and the relative permittivity of Au at , measured in situ as a function of Au thickness during vacuum deposition (by thermal evaporation at a rate of onto borosilicate glass substrates held at room temperature); here the bulk values (annealed) are achieved for .
Experiments were reported by Yang et al. , where a propagation length of was measured for the long-range mode propagating along an islandized Ag film on quartz covered with index-matching fluid and excited at via prism coupling (the experimental arrangement is sketched in the inset to Fig. 9(a)). Figure 9(a) compares the measured ATR spectrum (crosses) to the fitted theoretical response (solid curve) from which the effective relative permittivity and thickness of the islandized Ag film were extracted [ and —the relative permittivity follows the trends shown in Figs. 8(b), 8(c)]. Figure 9(b) shows the computed distribution of the z-directed component of the Poynting vector at , from which it is noted that essentially no energy propagates within the islandized film, explaining the long range of the mode (the field distribution is otherwise similar to that of the mode in a continuous film).
Takabayashi et al.  studied theoretically long-range modes supported by symmetric ultrathin Ag films at , using the thickness-dependant permittivity measurements for Ag reported in  combined with an effective medium theory. They reported that under certain conditions, s-polarized light can excite a surface wave in this structure. In , Takabayashi et al. investigated the impact of absorbing claddings in such structures.
Long-range surface mode propagation lengths in the range of were measured by Wu et al.  for Au, Cu, Al, and Fe islandized films on a glass substrate, covered with index-matching fluid, and excited at via prism coupling.
Takabayashi et al.  measured long-range surface mode propagation lengths of 300 and for islandized Ag films 4 and thick, respectively, on BK7 glass covered with index-matching fluid and excited at by using prism coupling; these propagation lengths are longer by factors of 75 and 120, respectively, than that of the corresponding single-interface SPP.
Kume et al.  observed both (TM) p- and (TE) s-polarized long-range modes along a composite layer comprising dispersed and isolated Ag spherical nanoparticles about in diameter embedded in an film, as a function of the composite layer thickness, which was varied from 13 to . The volume ratio of Ag to the total volume of the film was about 0.05, and the composite film was bounded slightly asymmetrically by claddings. Modes were excited by using prism coupling at . The longest propagation length measured was for the p-polarized long-range mode in the thick composite. They also investigated the effect of prism loading on their responses by varying the thickness of the intervening layer. Theoretical responses to the measurements using a relative permittivity for the composite layer estimated based on Maxwell–Garnett theory agreed reasonably well.
Wood et al.  studied the long-range mode propagated by islandized thick Ni films on covered with index-matching fluid and excited via prism coupling at . They then attempted to fit a theoretical response to the measurements in order to uniquely determine the optical effective permittivity and thickness of their Ni layer only to find that degenerate fits (i.e., different combinations of , , t produce equally good fits). They concluded generally that the thickness and permittivity of films cannot both be determined from theoretical fits to prism coupling experiments conducted with any long-range mode and that the thickness is required in order to determine the permittivity.
Though long propagation lengths are evidently achievable on islandized metal films, it is noted that fabricating such structures reproducibly is challenging and that an islandized film eliminates applications where a continuous metal is essential. However, the ability to support s-polarized waves in such structures is an interesting and potentially useful attribute.
2.8. Long-Range Surface Exciton Polariton[116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134]
Kovacs [26, 123] investigated an unusual system, consisting of a thin continuous Fe film bounded symmetrically by cryolite claddings. What renders the system unusual is that the relative permittivity of Fe at his chosen operating wavelength is and . Modal computations for semi-infinite claddings revealed that the mode had and only for the case of thin films ( or less). These inequalities no longer held as t was increased and did not hold for the mode at any t. (The corresponding single-interface does not support a purely bound SPP.) Additionally, it was noted for the mode that no energy propagates in the Fe film in directions parallel or antiparallel to that of modal propagation (i.e., ), in contrast to the LRSPP where energy propagates in the metal film along the antiparallel direction. The existence of this mode was confirmed experimentally via prism coupling at in a structure consisting of an Fe film thick bounded by about of cryolite.
Yang et al.  also showed that a thin slab of material satisfying can support a long-range surface mode when cladded symmetrically, but rather surprisingly, that its range increases as increases. Since the condition is readily met in continuous materials near an excitonic resonance, they termed the mode long-range surface exciton polariton (LRSEP). Experimental excitation of the mode was conducted at via prism coupling to a thick V film deposited onto a quartz substrate and covered with an index-matched fluid. (The effective optical permittivity of islandized metal films can also satisfy ; see Subsection 2.7, Fig. 8 and [108, 110, 112, 113, 114]; so the long-range surface wave propagating along them is often termed an LRSEP, although no excitonic resonances are involved.)
In a subsequent paper, Yang et al.  reported useful small-t asymptotic expressions for the long-range surface mode and explored them for various limiting cases of . Slight asymmetry and absorption in the claddings were also considered. An expression for the propagation constant of the LRSPP in the symmetric structure is given, from which it is observed that . They found that long-range surface modes exist in a symmetric structure for almost any value of , including the case and , which does not support a single-interface SPP. Expressions for the cutoff thickness and cutoff asymmetry ( and ) are derived for the LRSPP in slightly asymmetric structures; the expressions exhibit and dependencies, respectively, indicating that a larger t is required at longer wavelengths to maintain guidance given an asymmetry, and that the cladding permittivities must be more closely matched at longer wavelengths for a given t. They also discussed the role of the light line in the case of absorbing structures, indicating, as had been previously noted [72, 80], that it does not clearly separate the radiative and nonradiative regions of the dispersion curve; this was also observed in . Prism coupling experiments were conducted at on thick Pd and thick V films deposited onto a quartz substrate with index-matched fluid forming the other cladding. The Pd ( and V films were observed to support LRSPP and LRSEP waves, respectively.
At about the same time, Prade et al.  produced a study that included consideration of lossless metal slabs bounded by dielectrics where the permittivities satisfy or with , , and similar in magnitude. They showed and modes existing for as long as t remains below a cutoff thickness (recall that the purely bound single-interface SPP is not supported for ). In the asymmetric case the mode has an additional cutoff thickness at smaller t. The mode has a negative group velocity. For the case , only the mode exists, exhibiting a cutoff thickness at small t and evolving into the single-interface SPP at the interface with increasing t.
Yang et al.  applied the virtual mode treatment  to study theoretically prism coupling to the LRSPP and LRSEP. They demonstrated that the perturbation caused by the prism could increase the propagation length of the LRSEP, as opposed to only reducing it as is generally observed for the LRSPP.
Bryan-Brown et al.  excited the LRSEP at along Cr films 33.6, 20.9, 17.8 and thick deposited onto quartz substrates and covered by an index-matched fluid. Flat films excited by using prism coupling and corrugated films excited by grating coupling were investigated.
Bryan-Brown et al.  investigated at the LRSPP in thick corrugated Pd films and the LRSEP in thin corrugated Pd films both excited by grating coupling. The thickness dependence of the Pd permittivity allowed exploration of both types of wave by using the same material.
Giannini et al.  investigated the LRSEP along a thin amorphous Si layer cladded symmetrically by . They investigated a few wavelengths of operation but emphasize the ultraviolet ( and ) where the relative permittivity of Si satisfies . They found that the LRSEP in this structure slightly outperforms the LRSPP in the corresponding Au structure at these wavelengths, by comparing the propagation length, the field extension, and the associated FoM  of each wave. Prism coupling experiments were conducted with broadband light and with a laser source at , on structures comprising a or thick amorphous Si layer cladded by , confirming the excitation of the LRSEP.
2.9. Nonlinear Interactions [135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168]
Sarid et al.  computed the intensity-dependant propagation constant of the and modes in a symmetric structure consisting of a Cu film bounded by InSb at as a function of the Cu thickness. They found that the power required for a phase shift of over a distance of is about 1 order of magnitude lower for the LRSPP at compared with the corresponding single-interface SPP.
Stegeman et al.  computed the cross sections for copropagating and counterpropagating second-harmonic generation, degenerate four-wave mixing, and intensity-dependant phase shifts, driven by LRSPPs. Ag films cladded symmetrically by CdSe or InSb were considered. Ag thicknesses in the range of were found to be optimal for the nonlinear processes considered. Except for copropagating second-harmonic generation, which is limited by the inability to achieve phase matching, they obtain cross sections that are significantly larger than those related to focused Gaussian beams, stating as a reason that a better trade-off between confinement and can be achieved for LRSPPs.
Deck and Sarid  investigated theoretically second-harmonic generation by LRSPPs in a prism coupling arrangement excited at a fundamental of , assuming a nonlinear medium (quartz) as one cladding and an index-matched linear medium as the other. They predicted that the intensity of the outcoupled second-harmonic wave would be 2–4 orders of magnitude greater than in the corresponding Kretschmann–Raether and Otto configurations for the same incident field, and they attributed this improvement to the strongly enhanced fields associated with the LRSPP in this arrangement .
Quail et al.  verified this prediction , measuring an outcoupled second-harmonic signal having an intensity that is 2 orders of magnitude larger than in the Kretschmann–Raether and Otto geometries. The LRSPP was excited at a fundamental of by using prism coupling in a structure consisting of a thick Ag film on an X-cut quartz substrate covered with index-matching fluid, as shown in Fig. 10(a). The outcoupled second-harmonic signal was observed at a 5° offset compared with the angle of the fundamental because of dispersion in the prism. Figure 10(b) compares the measured second-harmonic reflection coefficient R (defined as the ratio of the second-harmonic irradiance to the square of the fundamental irradiance) to that expected for the single-interface SPP. In a subsequent study, Quail and Simon  monitored the transmitted second-harmonic signal generated by the LRSPP in a similar experimental arrangement, noting greater generation efficiency.
Stegeman et al.  computed the efficiency of second-harmonic generation by two counterpropagating LRSPPs excited at a fundamental of by prism coupling in a structure consisting of a thick Ag film on MNA (2-methyl-4-nitroaniline, an organic crystal; see Ref. 15 of ) used as the nonlinear cladding and assuming a linear index-matched material as the other cladding. Second-harmonic generation efficiencies of about were predicted.
Moshrefzadeh et al.  computed the second-harmonic signal generated by an LRSPP propagating along thin ( and ) Ag films bearing a thick nonlinear adlayer and bounded symmetrically by benzene. The LRSPP was excited at a fundamental of by prism coupling. An increase in second-harmonic intensity of , compared with simple reflection off the corresponding metal surface, was predicted for the thick Ag case.
Stegeman and Karaguleff  investigated theoretically degenerate four-wave mixing via LRSPPs excited by prism coupling at , predicting conversion efficiencies into the degenerate fourth wave of about 10%–50% for an Ag film thick on a nonlinear substrate (PTS, bis-(p-toluene sulphonate) of 2,4-hexadiyne-1,6-diol; see Ref. 15 of ) and covered with a linear index-matched material.
Liao et al.  computed the cross sections for second-harmonic generation by counterpropagating waves in dielectric waveguides and Ag slabs supporting LRSPPs. Fundamental ’s of 1.06 and were considered, assuming MNA and CdSe as the nonlinear materials, respectively. The Ag slab was bounded symmetrically by the nonlinear material and the dielectric waveguides incorporated the nonlinear material as either the core or the lower cladding. They find the highest cross sections when the nonlinear medium is used as the core of a dielectric waveguide.
Karaguleff and Stegeman  conducted a similar study for degenerate four-wave mixing, at assuming PTS, and at assuming InSb, as the nonlinear materials. They reached essentially the same conclusion as for second-harmonic generation , which is that the highest cross sections are obtained when the nonlinear medium is used as the core of a dielectric waveguide. They attributed the difference to the better confinement–attenuation trade-off available in the dielectric waveguides. They also noted that the LRSPP cross section in the corresponding prism-coupled geometry  is greater by a factor of about two compared with the freely guided (no prism) LRSPP.
Stegeman and Seaton  investigated theoretically the modes supported by the metal slab bounded on one or both sides by media with intensity-dependant refractive indices . They assumed lossless Cu films of thicknesses 50 and , and InSb (, self-defocusing) as the nonlinear medium and give results at . They found in this case that the modes approach cutoff with increasing mode power. The details of the approach are different depending one whether one or both claddings are nonlinear, but the mode in the thinnest film (i.e., the LRSPP) bounded on both sides by the nonlinear material approaches cutoff most rapidly. Assuming (self-focusing) generates a more complex power dependency, including the observation of a maximum in mode power transmission, mode fields developing maxima away from the film, and new modes existing above minimum power thresholds. Stegeman and Seaton searched for TE-polarized SPP’s in these structures but did not find them. However, subsequent work  revealed their existence for very thin metal films having a small . Ariyasu et al.  expanded on this work [144, 145], considering different combinations for the nonlinear claddings, including attenuation.
Mihalache et al. [153, 157] also found TE-polarized SPPs in the self-focusing case, for slightly asymmetric claddings. Boardman and Twardowski [156, 162] studied the interaction between TE- and TM-polarized modes in similar intensity-dependant nonlinear waveguides and considered to some extent the thin metal slab.
Hickernell and Sarid  developed a theory for prism coupling to the LRSPP in structures having a metal slab on a substrate exhibiting an intensity-dependant refractive index and covered by an index-matched medium. They found that an intensity 2 orders of magnitude lower than that required for the single-interface SPP is needed to observe bistability via the LRSPP.
Agarwal and Dutta Gupta  developed a theory for a multilayer structure on a similar nonlinear substrate and studied bistability with prism-coupled single-interface SPPs (Kretschmann–Raether), and the and LRSPP modes in a thin metal slab. They found that LRSPP bistability occurs at much lower intensity thresholds (by at least 1 order of magnitude).
Nunzi and Ricard  observed optical phase conjugation using single-interface SPPs and attributed the main nonlinear contribution in their experiments to heating from the metal film (i.e., thermal effects) rather than to field-induced nonlinearities in the media. They made the points that all of the energy coupled into the SPP is converted into heat in the metal film and that the adjacent dielectric region is where nonlinear interactions are sought. They commented that thermal aspects would remain important in comparable experiments conducted with the LRSPP even though the latter exhibits less attenuation. Sambles and Innes  raised similar points, emphasizing the need for consideration of thermal issues even when working with the LRSPP.
Yang and Sambles  measured, via prism coupling at , optical power-dependant changes in the excitation curve and angle of LRSEPs in a Pd film thick deposited onto BK7 and bounded on the other side by index-matched fluid. High- and low-power angle scans were conducted, along with power scans at fixed angles. The optical nonlinearities were thermally induced by heating of the Pd film and were deemed to be changes in the index of the matching fluid (thermo-optic effect) and changes in the thickness of the fluid layer (coupling gap).
Quail and Simon  observed from modal computations that the mode at the second harmonic had the same phase constant as the mode at the fundamental for a prescribed Ag thickness in their experimental arrangement. Based on this observation, they demonstrated phase-matched copropagating second-harmonic generation by exciting the mode at the fundamental and measuring the prism outcoupled mode (LRSPP) at the second harmonic.
Fukui et al.  considered the effects of a finite-width beam in prism-coupled second-harmonic generation experiments, computing the spatial profile of the reflected second-harmonic signal for an incident wide square profile (stepped intensity) beam. The profiles show a global maximum for a prescribed metal film thickness.
Building on this theory, Schmidlin and Simon  measured the profile of the reflected second-harmonic signal in the same experimental situation as that reported in , but using a wide output slit (moved via a translation stage) and a wide input beam. They deduced the propagation length of the LRSPP (excited at the fundamental) from the measured reflected second-harmonic spatial profile in the region outside the incident beam width.
Chen and Simon [154, 158] studied experimentally and theoretically second-harmonic generation from LRSPPs excited at a fundamental of in a corrugated Ag film bounded on both sides by quartz, although only one of the quartz claddings was responsible for the generation of the second harmonic owing to the presence of index-matching fluid at the other Ag/quartz interface. They observed enhanced second-harmonic generation due to the LRSPP but conclude that scattering from the grating grooves significantly limits the possible enhancement compared with flat films in a prism-coupled geometry. Phase matching between the fundamental and the second-harmonic LRSPP was also observed in one of the corrugated structures having the appropriate Ag film thickness .
Tzeng and Lue [159, 160] studied theoretically second-harmonic generation via prism-coupled excitation of the LRSPP in Ag , and Ag, Au, Cu, Al  films, bounded by linear dielectrics, where the second harmonic is generated by nonlinearities in the metal (bulk and selvedge regions) and outcoupled by the prism. Electron gas hydrodynamic theory  was used to describe the nonlinear response of electrons in the metal. They found an increase in the second-harmonic generation that is 2 orders of magnitude larger for the LRSPP compared with the corresponding single-interface SPP. They also noted effects due to the thickness of the coupling gap. Lue and Dai  extended the study to include coupled LRSPPs in structures comprising a pair of Ag slabs bounded on all sides by .
Simon et al.  and Wang and Simon  reported the backscattering of the second-harmonic wave from excitation of the LRSPP via prism coupling at a fundamental of in thin Ag films on quartz, bounded on the other side by index-matched fluid. The measured backscattered second-harmonic signal exhibited a peak in the direction of the incident beam (antispecular) with a slight angular offset due to dispersion in the coupling prism.
Li and Zhang  compared theoretically the nonlinear coefficients of the single interface, the metal slab, and the metal clads.
2.10. Biosensors [169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195]
The high surface sensitivity of the single-interface SPP has been heavily exploited in (bio)chemical sensors ever since the initial demonstrations . The conventional and mainstream approach to SPP sensing [169, 170] rests on the Kretschmann–Raether configuration (Fig. 4, case 3), where efficient coupling to the single-interface SPP depends strongly on the thickness and index of an adlayer located at the metal/fluid interface (surface sensing) and on the index of the sensing fluid itself (bulk or refractometric sensing). Various modifications and alternatives to this approach have been proposed over time in the quest to achieve greater sensitivity and lower detection limits. Reviews discussing different SPP sensor architectures and interrogation schemes were published by Homola and co-workers [171, 173] and by Chien and Chen .
Some of the alternative approaches proposed involve using the LRSPP instead of the single-interface SPP to perform the sensing function along metal slabs. The typical approach uses a prism-coupled arrangement such as that depicted in Fig. 4, case 1, with the lower cladding being the sensing fluid (usually an aqueous buffer) contained within a flow cell, and the upper cladding being a dielectric that has a refractive index closely matched to that of the sensing fluid. Popular materials for this cladding layer are Cytop or Teflon, which have a refractive index slightly above and slightly below (respectively) that of de-ionized water, so the sensing structures are typically slightly asymmetric. Au is normally used as the sensing surface given its chemical stability and the availability of good surface chemistries for this metal.
The use of the LRSPP in a prism-coupled arrangement was proposed as a sensor by Matsubara et al. . They computed the angular response at for structures comprising Ag as the metal film, bounded by ethanol on one side, and by on the other side also in contact with the prism. was selected because its index approximately matches that of ethanol. Corresponding Kretschmann–Raether configurations were also modeled. Matsubara et al. measured the angular response of the arrangements in ethanol, demonstrating a narrower response in the case of the LRSPP.
Similar LRSPP structures excited via prism coupling were reported in , where Teflon and were compared as cladding materials and Au was used for the metal film. The approach was demonstrated in bulk sensing experiments.
A similar structure was reported in , where Cytop was selected for the cladding. The approach was demonstrated in bioaffinity sensing experiments on Au spots using an imaging prism-coupled arrangement. An improvement of 20% over the conventional Kretschmann–Raether configuration was claimed.
Slavík and Homola  demonstrated bulk sensing using a Teflon-Au-water structure. The optimal thickness of the Au film and of the Teflon layer were determined from modeling at and selected such that the highest coupling and bulk sensitivity would be achieved. Interrogation was conducted by using an LED centered at and a spectrometer. They reported a bulk sensing detection limit of and stated that this is the best value reported to date for any bulk sensor based on SPPs. The measured spectral position of the LRSPP dip is shown in Fig. 11 as aqueous sensing fluids of different refractive indices are injected into the sensor. The baseline noise signal is also shown on an expanded scale, from which the noise was deduced to be .
Vala et al.  compared a similar Teflon-Au-water structure with a conventional Kretschmann–Raether structure for detecting large analytes such as biotinylated latex beads (through a streptavidin bridge) and bacteria (heat-killed Escherichia coli through antibodies). Cutoff of the LRSPP was noted in the case of latex beads for the test structure having the thinnest Au film. The sensitivity of the LRSPP sensor was found to be about greater than the conventional Kretschmann–Raether configurations for the detection of E.coli.
Slavík and Homola  considered the Teflon-Au-water and -water structures, optimizing the layer thicknesses (, Teflon, Au) such that efficient prism coupling could be obtained into both the and modes of the same structure under either angular or wavelength interrogation. They found for a properly optimized structure that near perfect coupling can be achieved into both modes and verify this experimentally for a Teflon-Au-water structure under wavelength interrogation. They demonstrated bulk sensing and argue that the availability of two sensing modes having different bulk and surface sensitivities is advantageous for removing unwanted perturbations in biosensor applications. Further work on this approach was reported in , including noise and cross-sensitivity analyses and biosensing experimentation with IgE. The approach was compared with a two-channel compensated Kretschmann–Raether configuration.
Hastings et al.  also considered dual-mode ( and ) operation in the Teflon-Au-water structure under wavelength interrogation. They showed that monitoring both resonances allows changes in adlayer thickness and changes in the bulk index of the sensing solution to be monitored separately. Through modeling they found optimal designs that minimize the detection limit for surface and bulk sensing by using a Cramer–Rao lower bound to estimate the smallest detectable shift in the coupling dips. They conducted experiments demonstrating the ability to monitor changes in bulk index and changes in surface coverage, the former by alternating between two buffer solutions and the latter via the biotin-streptavidin system (with biotin immobilized onto the sensor surface). The detection limits for bulk and surface sensing were estimated to be RIU (refractive index units) and , respectively. In  they investigated the same concept but under angular interrogation and found a reduction in cross sensitivity compared with wavelength interrogation.
Enhanced fluorescence from tagged antibodies pumped with LRSPPs was demonstrated by Kasry and Knoll . The LRSPPs were excited by prism coupling in a Teflon-Au-water structure and the tagged antibodies were adsorbed onto a thick polysterene spacer layer deposited on the Au surface. Fluorescence intensities about 22 times larger were observed for the LRSPP pump compared with pumping with the single-interface SPP in corresponding experiments. The enhanced fluorescence intensity was attributed to the enhanced field of the LRSPP . The authors also pointed out that the larger penetration depth into the sensing medium of the LRSPP could lead to further enhancement if a thick matrix were used to bind tagged analyte to the sensor along the third dimension (into the sensing medium). In a subsequent paper, Knoll et al.  reported the binding of a hydrogel matrix to a Cytop-Au structure and then demonstrated sensing using the LRSPP of free prostate specific antigen to its antibody immobilized in the matrix.
Dostálek et al.  compared prism-coupled Cytop-Au-water and Teflon-Au-water structures at . They explored the sensitivity of both structures for bulk and surface sensing and deduced the detection limits from their measured responses and from the width of their reflection dips. Both structures improve on the detection limits of the conventional Kretschmann–Raether configuration, with the Teflon structure being about better for bulk sensing if operated near cutoff, and the Cytop structure being about better for surface sensing. They measured the fluorescence intensity of dye pumped by the LRSPP (following ) as a function of dye distance from the metal film set by a spacer layer consisting of a protein stack or a thin layer of Cytop. Quenching was noted for a distance of , maximum intensity was measured for distances of , and significantly more fluorescence was measured compared with the conventional Kretschmann–Raether configuration. Their setup for conducting fluorescence studies in this manner is sketched in Fig. 12. They also deposited a series of Au films of thickness ranging from to on Cytop and on Teflon and extracted the optical parameters of their films by fitting measured ATR spectra to theoretical responses. They observed some substrate dependence in the extracted optical parameters, probably due to differences in the roughness of the starting surfaces, but otherwise observe the same trends in the relative permittivity with decreasing t as discussed in Subsection 2.7. The bulk optical parameters for Au were achieved in their cases at .
Wang et al.  demonstrated a biosensor capable of detecting aflatoxin in milk using a Cytop-Au structure. The sensor detects LRSPP pumped fluorescence emitted from labeled antibodies (Cy5-GAR) bound to the sensor surface following the inhibition immunoassay format (Au/thiol/BSA- (where BSA is bovine serum albumin) immobilized onto the sensor followed by detection of /Cy5-GaR).
Dostálek et al.  proposed a biosensor using a pair of broadside coupled thin Au slabs separated from each other and from the prism (on one side) by thin layers of Cytop. As discussed in Subsection 2.5 (with regards to ), coupled thin slabs support supermodes consisting of symmetrical and asymmetrical couplings of and modes, and when the metals are thin, the two supermodes involving the mode are long range. In  the symmetric and asymmetric long-range supermodes are excited at different angles of incidence at the same operating wavelength. The modes have different probing depths into the sensing medium; so they exhibit different bulk and surface sensitivities. Thus, by fitting measured ATR spectra involving both supermodes to theoretical responses (computed using the transfer matrix method), the thickness and refractive index of an adlayer can be extracted. Dostálek et al. demonstrated the approach by characterizing hydrogel layers (thickness and index) deposited onto their sensors and by observing the diffusion of BSA into the hydrogels.
Rajan et al.  modeled a sensor consisting of a multimode step-index fiber with the cladding removed over a length, whereupon an LRSPP supporting structure is deposited directly onto the core, the structure consisting of a Teflon layer followed by a thin Au film covered with the sensing fluid. The bulk sensitivity of the sensor was determined under wavelength interrogation as a function of structure parameters. In a subsequent paper, Jha et al.  studied a similar structure with tapered fiber sections added on either side of the sensing region.
Chen et al.  investigated theoretically the temperature stability of prism-coupled LRSPP sensors in a Cytop-Ag-water configuration by assuming a temperature-dependant model for the permittivity of the Ag film. They found that the configuration is thermally stable over a large temperature range , but they ignored the temperature dependence of the prism (SF10), and of Cytop and water (which have large thermo-optic effects).
In the above-described studies, the index symmetry required to ensure propagation of the LRSPP was provided by selecting a cladding material that is closely index-matched to the index of the sensing fluid. An alternative approach involves using a thin high-index layer, which when taken in combination with the low index of the sensing fluid, creates an effective index that closely matches the index of the material on the other side . Such a structure is physically asymmetric but appears symmetric from an effective medium point of view. Sensing using the LRSPP in such a structure has been considered, with coupling provided by a corrugated grating [176, 179], and where one of the claddings comprises two layers, a thin high-index dielectric layer followed by the sensing medium (air or water), such that the effective index of the combination matches the index of the silica cladding on the other side of the metal film.
A similar idea was applied in  to create an LRSPP sensor based on prism coupling, but where a Teflon/ bilayer system was used to match (effectively) the index of the aqueous medium on the other side of an Au film. This approach was demonstrated in bulk sensing experiments.
The LRSPP in Kou and Tamir’s structure  was considered theoretically in  for sensing by using a corrugated grating as the coupling means. Liao et al.  explored theoretically prism-coupled LRSPPs in a multilayer system and commented on the suitability for chemical sensing. (Supporting the metal slab on a 1D finite photonic crystal as discussed in Subsection 2.L  or suspending the metal slab as discussed in Subsection 2.4  and Subsection 3.9  represents a further alternative to satisfying the index symmetry requirement for sensing with the LRSPP.)
2.11. Emission and Molecular Scattering [196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217]
The decay of excited molecular emitters directly into SPPs via the near field [196, 197] has been a subject of intense study . Reports involving the LRSPP in this process (and the like) are reviewed here.
Gruhlke et al.  observed optical emission mediated by the LRSPP and modes in symmetric corrugated Ag films of various thicknesses , whereby molecular emitters adjacent to the Ag film excited directly the modes, which were then outcoupled into free radiation by the corrugation. The dispersion of the modes over a broad spectral range was also measured via this emission process. The LRSPP in these experiments was responsible for most of the observed emission.
In a subsequent publication , Gruhlke and Hall indicated that in their experiments [199, 200] the excited molecules are on one side of the metal film and the emission is observed on the other side. They emphasized that the process provides a dispersive channel through their otherwise (nearly) opaque metal films (as in studies involving prism coupling [19, 22, 43]). Figure 13 shows the measured radiation spectrum emitted from the symmetric structure sketched at the origin of the polar plot. Fluorescent molecules are located in the lower photoresist layer and are pumped at through the backside at normal incidence. The spectrum is normalized to the intensity measured from a sample that contains only the lower (fluorescent) photoresist layer. The emitted intensity is stronger at owing to mediation by the LRSPP. Indeed, the fluorescence reradiated by the LRSPP into this peak is of that radiated isotropically (at the same wavelength) in the absence of the metal layer. They stated that field enhanced fluorescence likely had a role to play.
Leung et al.  computed the decay rate of a dipole placed above asymmetrically and symmetrically cladded (free-standing) corrugated gratings. Among all the cases considered, the strongest enhancement of the rate was obtained for coupling into the LRSPP of a thick free-standing corrugated Ag grating, where the decay rate was about that of a free dipole.
Lenac and Tomaš [202, 203] derived the cross sections for absorption of the and modes by molecules placed above the metal in a symmetric configuration as a function of t. They found decreasing and increasing cross sections with decreasing t for the and modes, respectively. They also found larger cross sections at shorter ’s. The differences were attributed to differences in modal confinement.
Tomaš and Lenac  also derived the cross sections for scattering of the mode into the mode (as well as other modal combinations) by molecules placed above the metal in a symmetric metal slab as a function of t. They found that the cross section decreases as t decreases and estimated that the LRSPP scattering cross section along an Ag film thick is 1–2 orders of magnitude smaller than that for scattering between single-interface SPPs. They also commented that the LRSPP cross section would compete with the LRSPP field enhancement in prism-coupled arrangements , such that surface-enhanced Raman scattering (SERS) mediated by prism-coupled LRSPPs is not likely to be significantly more enhanced than SERS mediated by prism-coupled single-interface SPPs.
In a subsequent paper , Lenac and Tomaš derive the power lost by a molecular dipole to the and modes of a symmetric metal film as a function of t, with and without prism outcoupling. The dipole was located near and below the metal film within the bottom cladding. Their computations without prism outcoupling (freely guided) showed decreasing and increasing power coupled into the and modes, respectively, with decreasing t. For an Ag film thick, about more power is coupled into the mode than into the LRSPP at ; about more power would be coupled into the corresponding single-interface SPP. However, for a dipole located farther from the metal, the LRSPP is preferentially excited, since its fields penetrate more deeply into the cladding than those of the mode. The angular distribution of power emitted into the prism by outcoupling from the and modes, themselves excited by the dipole, was computed, showing that the emission is dominated by the mode rather than the mode. The computations for the mode with decreasing t showed a decrease in the angular width of emission and an increase in the peak of the angular distribution of emitted power (i.e., in ). The peak is about greater for mediation via the LRSPP in an Ag film thick than for a free dipole in the same medium as the claddings. Computations for mediation via the single-interface SPP in the corresponding (optimized) Kretschmann–Raether configuration yield a peak about greater, so mediation via the LRSPP produces an improvement of about . They also estimate a SERS intensity enhancement of about for mediation via the LRSPP versus the single-interface SPP in prism-coupled geometries. Thus, although the coupling of a dipole to the freely guided LRSPP (no prism) and the molecular scattering cross section of freely guided LRSPP’s  are lower than for the corresponding single-interface SPP, the field enhancement associated with prism coupling  compensates sufficiently to produce greater and sharper emission peaks for the LRSPP. Subsequent calculations comparing SERS mediated by single-interface SPPs in the Otto geometry with mediation by the prism-coupled LRSPP also yield an enhancement of about in the peak for the LRSPP .
Barnes and co-workers have addressed emission extraction through metal films in light emitting diodes (e.g., [198, 207, 208, 209, 210, 211]) for different device architectures and from different perspectives, such as the direct coupling of molecular emitters to SPP modes and the conversion of SPPs into free radiation. The most promising approaches to date involve mediation via the coupled modes of an effectively symmetrical thin metal contact [208, 209, 210, 211], where symmetry is achieved by coating the top of the metal contact with a thin dielectric overlayer of a prescribed thickness, which taken in combination with air yields an effective medium that is index matched to the underlying active medium. Outcoupling of the and modes to free radiation is achieved via corrugations in the full structure or in the thin dielectric overlayer. The metal contact in these studies is thick (Ag).
Chiu et al.  compared the measured emission from a two-layer structure consisting of on photoresist on Si with the emission from a four-layer structure consisting of on photoresist on Si. The layers were thick, and the Au layers were thick. Each structure was formed into a 1D rectangular lamellar grating having the purpose of outcoupling the surface plasmons propagating along the structure into free radiation. The structures were pumped by light, and the molecules decay into the modes of the structure including the surface plasmon modes supported therein. Enhanced emission in the range of to was measured from the four-layer structure compared with the two-layer structure and a flat (planar) structure. The enhanced emission was attributed to mediation by LRSPPs. The modes supported by these multilayer structures were not investigated at the operating wavelengths.
Andrew and Barnes  demonstrated the transfer of energy between donor and acceptor molecules (dipoles) through thin Ag films, with the transfer being mediated by the and modes. The Ag film was cladded below by a thick layer of -doped polymethyl methacrylate (PMMA, donor) on , and above by a thick layer of R6G-doped PMMA (acceptor, where R6G is Rhodamine 6G) , as sketched in Fig. 14(a). The structure is essentially symmetric at the wavelengths of interest. Figure 14(b) shows the measured emission spectrum of an -doped layer only, and the emission and absorption spectra of an R6G-doped layer only. The emission spectrum of the -doped layer overlaps strongly with the absorption spectrum of the R6G-doped layer (in the range to ); so, when pumped from the back side (at ), the excited molecules can donate their energy to the R6G molecules through interaction with the and modes of the structure. The R6G molecules then decay spontaneously at , and part of this emission is captured by the optical detection setup aligned with the top of the structure. Figure 14(c) shows the measured spectrum obtained from the structure for the case of a thick Ag film (black curve), along with the spectra of two control samples, each having only one type of dipole—donor only (no R6G, blue curve) and acceptor only (no , red curve). The emission was significantly larger than in the case of the control samples, indicating efficient energy transfer through the Ag film via mediation by the and modes. The mediation accounted for 70% of the total emission.
Okamato et al.  suggested that lasing in the LRSPP in a corrugated grating bounded symmetrically by gain media [4-dicyanomethylene-2-methyl-6- (p-dimethylaminostyryl)-4H-pyran (DCM)-doped , peak ] could be achieved, by pointing out that the loss of the LRSPP in flat Ag films thinner than is lower than the available gain in the medium considered. In their concept, standing LRSPP waves would be amplified by stimulated emission and partially outcoupled by the grating forming an output laser beam.
Winter et al.  investigated further the concept proposed by Okamoto et al. , considering the existence of the mode as well as the LRSPP. Photoluminescence measurements on (effectively) symmetric corrugated structures similar to those reported in , where the Ag film thickness was varied from 20 to , show that a significant amount of power is indeed coupled into both the and modes. They also compute the fraction of total power coupled by a dipole emitting at , and positioned 20 and away from a flat symmetrically cladded Ag film, into the and modes of the structure, as a function of t. For a dipole separation of and for , about 80% of the power is coupled into the mode while about 7% couples into the LRSPP. However, for a dipole separation of and for , about 15% of the total power is coupled into the mode, and about 15% is coupled into the LRSPP. For larger dipole separations, more power would (proportionally) be coupled into the LRSPP than into the mode, a trend consistent with earlier computations , owing to the greater extension of the LRSPP mode fields into the cladding. They then estimate the gain available to each mode by ascribing to each a fraction of the total gain taken to be the same as the fraction of total power coupled in from a dipole at a particular location. Taking into account in this manner the gain lost to the mode, their results suggest that lasing in the LRSPP would be possible, but an Ag film thinner than about would be needed.
Wang and Zhou  considered the prospects for amplification of LRSPPs in a structure consisting of an Au film on Si and covered with a multilayer system consisting of alternating Si and Er:Si nanolayers. Ignoring the mode, they predicted net gain into the LRSPP at for Au films thinner than about .
De Leon and Berini  proposed a model for SPP amplification that accounts for the nonuniform gain distribution of a dipolar gain medium pumped at broadside and placed along a symmetric metal slab. The model takes into account four channels for excited state decay (coupling to , and radiative modes, and coupling to electron–hole pairs in the metal), leading to a position-dependant dipole lifetime. Additionally, the model uses a realistic pump irradiance distribution within the gain medium as computed by using a transfer matrix method. The rate equations for the standard four-level pumping model are then applied locally, with the lifetime and irradiance taking on their position-dependant values, leading to a nonuniform gain distribution. The distribution is then incorporated into a multilayer waveguide mode solver from which mode power gains are computed for the modes of the system. Using this approach, they predict that net amplification of the LRSPP is possible in the visible by using a reasonable concentration of R6G molecules and a reasonable pump irradiance ( at , pulsed) assuming a thick Ag slab bounded by Cytop on one side and the index-matched gain medium on the other.
2.12. Other Studies [218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245]
Other studies of the LRSPP in the metal slab include investigations involving magnetic materials [218, 219, 220, 233, 236], electro-optic materials [221, 226, 229, 230, 235] and photonic crystals [237, 245], studies of polarizing devices [222, 223, 224, 225, 228, 231, 234] and pulse reshaping devices , a study of propagation through gaps in the metal slab , and a study of bending out of plane . Electro-optic and filtering devices involving the LRSPP propagating across an array of metal nanowires have been explored [242, 244]. The LRSPP supported by semiconductor heterostructures has been investigated . Approaches to extend the range of the single-interface SPP by involving a large multimode waveguide [232, 243] or by using an anisotropic photonic crystal  have been reported.
Sarid [218, 219] modeled the modes supported by a metal film bounded by identical magnetic semiconductor claddings subject to a magnetization field in the plane of the metal film and perpendicular to the direction of propagation of the modes (transverse magnetization). He found that the range of the LRSPP decreases with increasing magnetization owing to the destruction of symmetry in the dielectric claddings. Cutoff thicknesses were noted for the LRSPP.
Hickernell and Sarid  investigated LRSPPs in a thin magnetic metal film (Ni) bounded symmetrically by dielectrics under transverse magnetization. They computed the change in the propagation constant of the LRSPP as a function of t due to the magnetization, showing that the change decreases as t decreases. The differential reflectance of a thin Ni film due to the application of a magnetic field was measured by using prism coupling through the excitation angle of the LRSPP. A small magnetically induced modulation was detected.
Sepúlveda et al.  studied the metal slab bounded on both sides by magneto-optic dielectric claddings (yttrium iron garnet, YIG), or bounded on one side by YIG and on the other by an index-matched nonmagnetic dielectric. Both configurations were considered under different magnetization directions (transverse, longitudinal). Nonreciprocity was predicted for both the and the mode in certain cases and was quantified with decreasing t. Their results indicate that nonreciprocal LRSPPs can be obtained, in particular, for transverse magnetizations of opposite direction in the lower and upper claddings for the YIG cladded film, and for both transverse magnetization directions in the case of the YIG/dielectric cladded structure. They also discussed the LRSPP in a ferromagnetic metal film (Co) bounded by dielectrics.
Khurgin  studied a similar system, consisting of a thin metal film bounded on one side by a magneto-optic material (bismuth-doped gadolinium iron garnet, Bi:GdIG) and on the other by an index-matched dielectric. He also concluded that nonreciprocal LRSPPs are possible in this structure under transverse magnetization.
Device applications using electro-optic materials include modulators [221, 226, 230] and a tunable filter , where the LRSPP is excited via prism coupling in a thin metal film bounded on one side by an electro-optic dielectric, and on the prism side by an index-matched dielectric. The electric field is applied to the electro-optic medium via the metal film and an additional electrode deposited onto its other surface. The narrow width of the LRSPP reflection dip improves the performance of the devices compared with the single-interface SPP.
Liu and Xiao  proposed and modeled an electro-optic switch consisting of a lossless metal slab bounded on both sides by an electro-optic material (, BST). Switching occurs by varying electro-optically the phase (beating) between the and the mode of the slab and taking as the outputs of the switch the top cladding in isolation and the bottom in isolation.
Konopsky and Alieva  reported the excitation of LRSPPs on a thin Au film bounded by air on one side and on the other by a finite 1D photonic crystal implemented as a multilayer stack on a BK7 prism. Semi-infinite 1D photonic crystals support Bloch surface waves in the bandgap, having fields that oscillate and decay exponentially into the crystal. Thus the full structure (finite 1D photonic crystal, intermediate dielectric layer, metal slab) was designed such that essentially the same LRSPP field distribution was achieved within the metal and air regions as in the corresponding hypothetical air–metal–air structure. An LRSPP propagation length estimated to be about 2 orders of magnitude longer than that of the single-interface SPP in the corresponding Kretschmann–Raether configuration was measured. In subsequent work , the authors described a similar structure, shown in Fig. 15, but using an thick layer of Pd as the slab. The computed ATR angular reflection spectrum is sketched as the red curve in Fig. 15, showing the intensity dip due to coupling to the LRSPP and fringes due to interference between the outcoupled field and the reflected beam (see also Subsection 2.3). They demonstrated sensing of 3% in a atmosphere by using this structure. A finite 1D photonic crystal provides a practical means to support the metal slab for an arbitrary dielectric bounding its other side (air, water, etc.) [237, 245], but the design is specific to the operating wavelength, metal slab (thickness and index), and index of the bounding medium.
The LRSPP has a role to play in polarizers and in polarizing couplers and splitters [222, 223, 224, 225, 231, 234], principally because of the strong coupling that can be achieved with fibers and dielectric waveguides in a broadside arrangement. In , specific modes of a dual-mode fiber were selectively excited by coupling to the LRSPP in a structure similar to that of Kou and Tamir , and where the LRSPP was excited using prism coupling.
Andaloro et al.  studied theoretically the reshaping of pulses on reflection from prism-coupled metal films where either LRSPPs or single-interface SPPs are excited.
Sidorenko and Martin  investigated the tunneling of SPPs across a perpendicular gap in the metal slab filled with the background dielectric, for symmetrically and asymmetrically cladded structures, as a function of the gap length and metal slab thickness. Results at show that the LRSPP on a thick Au film can tunnel through a long gap while retaining 80% of its field amplitude. Standing wave patterns in front of the gap are caused by reflection of the LRSPP by the gap. Transmission of the mode through the gap is much less efficient and excites the LRSPP on the other side.
Sun  investigated bending along metal slabs that are curved out of plane [radius of curvature in the plane, Fig. 1(b)]. The metal slab was cladded symmetrically by thin dielectrics, and the structure was bounded by air. The dielectric/air interfaces provide additional confinement to the LRSPP along this direction, allowing low-radiation-loss bends. High transmission levels were computed by using the finite difference time domain method for the LRSPP around 90° bends having radii on the order of the free-space wavelength of operation.
Wu et al.  considered a structure where the metal slab was replaced with an array of ( thick) metal nanowires, cladded on one side by an electro-optic polymer and on the other by an index-matched glass substrate. They computed ATR angular spectra for the structure under prism coupling, showing a narrower dip compared with a continuous metal slab of the same thickness, due to the excitation of the LRSPP across the wires. Low-loss surface waves can be supported along such a structure, as in discontinuous and islandized films (Subsection 2.7). The narrower dip leads to improved modulator performance compared with, e.g., . In a subsequent study  they investigated a similar structure, with the electro-optic polymer replaced by glass, used as a notch filter under prism-coupled wavelength interrogation. They reported filter designs having a bandwidth and tunable over about by varying the angle of incidence of the input light.
Plumridge and Phillips  considered theoretically the SPPs and the dielectric modes guided by a multiple quantum well heterostructure. The quantum wells are modeled as a quasi-2D electron gas, where the electrons are free to move in the plane of the well as a gas of almost free electrons, but are confined and restricted to intersubband transitions in the perpendicular plane. The main electric field component of the LRSPP is perpendicular to the wells (i.e., not in the plane of the electron gas); so damping is due principally to intersubband transitions (i.e., not to free-electron scattering), which by design can have energies far from that of the LRSPP, leading to lower propagation loss. Structures were modeled by using an anisotropic permittivity for the quantum wells, based on the Drude model for the permittivity components in the plane, and a Lorentz oscillator model for the permittivity along the perpendicular. The model was applied to a multilayer structure, and the LRSPP (among other modes) was explored in the infrared ( to ). Propagation lengths of many centimeters were predicted.
Lee and Gray  proposed a Kretschmann–Raether type of structure for re-exciting single-interface SPPs, whereby the thin metal film is placed directly onto the core of a large multimode dielectric waveguide. The multimode waveguide effectively traps the reradiated light outcoupled by the prism (and normally lost) and redirects it toward the metal film at the proper angle for re-excitation of the single-interface SPP. Simulations confirm the viability of range extension via this approach. Montgomery and Gray  further explored the concept and its design space through finite difference time domain simulations and modal analyses.
Krokhin et al.  derived the dispersion relation for the SPPs propagated along a thin metal film bounded by vacuum on one side and a strongly anisotropic photonic crystal on the other. They showed for the SPP localized at the metal–crystal interface that its range can be increased by 50% in the infrared, compared with the corresponding isotropic case, by orienting the optical axis of the substrate along the perpendicular to the metal film. The penetration depth of the SPP into the substrate increases with its range.
3. Metal Stripe
3.1. Modes of the Metal Stripe
The thin metal film of thickness t, finite width w, and relative permittivity , bounded by optically semi-infinite dielectrics of relative permittivity and , is sketched in front cross-sectional view in Fig. 16 and is henceforth referred to as the metal stripe. The metal stripe is obtained from the metal slab (Fig. 1(b)) by limiting its width. This increase in dimensionality leads to major changes, principally, the creation of lateral confinement and thus of an enriched mode spectrum, and an LRSPP having a lower attenuation than its counterpart in the corresponding metal slab. Another important change is that modal solutions to Maxwell’s equations must be obtained numerically, increasing considerably the analysis effort, whereas the modes of the slab can be derived analytically rather straightforwardly. Despite this difficulty, the metal stripe can be handled by well-established numerical techniques and by some commercial modeling tools, if appropriate care is taken. Theoretical studies using vectorial formulations of the method of lines (MoL), the finite element method (FEM), and the finite difference method (FDM) have been reported. The effective index method (EIM) has also been shown to approximate reasonably well some of the modes, including the LRSPP.
There are four fundamental modes supported by the metal stripe, labeled , , , and . Given the finite width of the structure, higher-order modes having extrema along x in their field distribution can also be supported. The field component dominates for all modes when ; so the modes are TM in character, but not purely TM because all field components including are always nonzero.
The nomenclature adopted extends that used for the slab and describes the field component of the mode: a and s refer to asymmetric and symmetric, respectively, the first position being associated with the horizontal dimension and the second with the vertical one ; b signifies purely bound (nonradiative), and the superscript counts the number of extrema in the horizontal distribution of , not counting the corner peaks.
The evolution of the modes with dimensions , materials and operating wavelength is complex, especially for asymmetric structures . The difficulty arises because all modes are supermodes (coupled modes) created from the coupling of elemental corner and/or edge (or finite-width interface) modes, with the selection of specific ones depending on the similarity of their phase constants. Since the elemental modes also change with structure parameters and operating wavelength, the supermodes can at times evolve unpredictably. However, trends have been noted over a range of studies pertaining to symmetric structures , as discussed in the following paragraphs.
The evolution of the modes as 0 resembles the evolution of the and modes of the symmetric metal slab in that all modes eventually become partitioned into either lower-attenuation ( and are -like) or higher-attenuation (, are -like) modes, as determined by the distribution (symmetric or asymmetric) of along y. This partitioning is readily observed from Fig. 17 in the case of a wide stripe (compare Fig. 2). Unlike in the metal slab, the modes in the metal stripe are not asymptotic with the single-interface SPP as w and t increase. Instead, the four fundamental modes are asymptotic and degenerate with an elemental corner mode, and the higher-order modes are asymptotic with elemental edge modes, as is also readily observed from Fig. 17.
The mode in the thin symmetric metal stripe is the fundamental long-range mode and, following convention, is termed a LRSPP.
One of the fundamental modes, the mode, evolves smoothly and predictably as w, into the vertically polarized TEM wave of the background. Its and fields evolve from being highly localized near the metal corners, as shown in Figs. 18(a), 18(c), to being spread out over the waveguide cross section, as shown in Figs. 18(b), 18(d). (Figure 18 plots the normalized , where is the complex power density carried by the mode.) The resulting field distribution [Figs. 18(b), 18(d)] can be well matched to Gaussian-like fields, such as those emerging from dielectric waveguides (e.g., single-mode fiber, SMF), leading to efficient end-fire excitation. This modal transformation as is accompanied by a reduction in confinement and attenuation, as noted from Fig. 17 [red curve in Fig. 17(b)], due to reduced field penetration into the metal. The mode in the thin symmetric metal stripe is the fundamental long-range mode and, following convention, is termed a LRSPP.
The LRSPP in wide metal stripes has an attenuation similar to that of the LRSPP in the corresponding metal slab, but narrowing the width can reduce the attenuation further by orders of magnitude. However, as in the metal slab, confinement must also be traded off against attenuation. Furthermore, it is important for the structure to be symmetric in order for the mode to remain purely bound, long range, and well behaved with varying structure parameters and operating wavelengths. The degree of asymmetry that can be tolerated depends principally on the confinement provided by the stripe.
Efficient end-fire excitation of the LRSPP in the metal stripe can be achieved by using a free-space beam in the manner depicted in Fig. 6 or by the polarization-aligned fundamental mode of a SMF, butt coupled directly to the input of the structure as depicted in longitudinal cross-sectional view in Fig. 19(a). Butt coupling can also be achieved with a polarization-maintaining SMF (PM-SMF) or tapered SMF.
One of the other fundamental modes, the mode, evolves in a similar manner as , except that its field develops two extrema along x, and the mode becomes unguided below cutoff dimensions that depend on the materials and operating wavelength. This mode, which is fundamental for large w and t, evolves into the first long-range higher-order mode as . Long-range modes of orders higher than the mode may also exist, originating from the mode families (, odd) and (, even). As observed from Fig. 17, they have cutoff dimensions that increase with mode order m and are long-range only near cutoff.
All modes having the and symmetries increase in attenuation as and do not couple efficiently with Gaussian-like fields in an end-fire arrangement because is asymmetric along y. The and modes are guided for all film dimensions and remain localized near the corners.
In asymmetric structures , true field symmetries exist only with respect to the y axis. Mode fields exhibit symmetriclike or asymmetriclike distributions with respect to the x axis, and localization to either the top or bottom metal–dielectric interfaces. As in the slab, the symmetriclike modes are localized along the interface with the lowest dielectric, while the asymmetriclike modes are localized along the interface with the highest dielectric. The evolution of modes with structure parameters is not obvious, since different elemental modes can merge in and out of a supermode as parameters change. Also, different numbers of extrema may occur along the top and bottom interfaces (they are counted along the interface where the field is localized, and this number is used in the nomenclature). Long-range modes can exist in asymmetric structures but only near cutoff and having perturbed field distributions that compromise excitation in an end-fire arrangement. The long-range modes are localized on the low-index side with fields that extend deeply into the high-index cladding.
3.2. Straight Waveguides [246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282]
The symmetric metal stripe provides confinement in the plane transverse to its longitudinal axis. The stripe can be dimensioned such that the mode propagates with low loss (as the LRSPP) and couples efficiently in an end-fire arrangement to TM-polarized Gaussian-like beams (in fiber or free space) while all long-range high-order modes are cut off and any remaining high-loss modes are excited with very low efficiency.
The majority of the experimental work conducted to date has used butt-coupling to SMF to excite the waveguides; this is in contrast to the metal slab (Section 2), where prism and grating coupling are normally used. As mentioned in Subsection 2.1, end-fire coupling is easier to implement than prism coupling, but all modes (including radiative ones) that overlap to some extent with the input fields will be excited; outputs must therefore be interpreted carefully, especially in cases where the input coupling is inefficient. Fortunately, LRSPP butt-coupling excitation efficiencies greater than 90% are readily achievable, rendering the excitation of any other modes (usually) immaterial.
Work conducted on straight waveguides propagating the LRSPP [246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282] is reviewed in this subsection. Passive integrated structures [283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299, 300, 301, 302] are reviewed in Subsection 3.3, but some of the work also reports results on straight waveguides that are discussed here.
Berini  studied theoretically the four fundamental bound modes of the symmetric metal stripe using the MoL and extended the mode nomenclature of  to identify them. The study was conducted at on a wide Ag stripe for a range of t, with the stripe embedded in an optically infinite dielectric background having . The constitution of modes from elemental corner and edge modes was discussed. The partitioning of modes with decreasing t into -like and -like modes was observed. The evolution of the mode into the LRSPP, including the extension of its fields from the corners into a Gaussian-like distribution, was described. The need to trade off attenuation versus confinement was noted, end-fire coupling was suggested for exciting the mode, and its use for optical signal transmission was proposed. (In , parts (a) and (b) of Fig. 2 were inverted with respect to its caption; this is corrected here as Figs. 18(c), 18(d).)
A subsequent study on symmetric structures  investigated via the MoL the fundamental and higher-order modes as a function of stripe dimensions, background permittivity, and operating wavelength. Figure 17 gives the effective index and normalized attenuation of the first eight modes of the stripe and of the and modes of the corresponding slab, and Fig. 18 shows some mode intensity distributions. Cutoff dimensions were noted for the higher-order long-range modes. It was found that decreasing the stripe dimensions (w and t) and the background permittivity, and operating at longer wavelengths, decreased the attenuation and confinement of the LRSPP. MPAs of were predicted for the LRSPP in the infrared . It was suggested that the stripe operating in the LRSPP could be used to implement integrated structures such as interconnects, splitters, couplers, and interferometers.
Charbonneau et al.  demonstrated propagation of the LRSPP in the metal stripe by measuring output intensity distributions from wide, thick, long Au stripes embedded in . The structures were excited at via end-fire coupling to a standard PM-SMF as sketched in Fig. 19(a). Outputs were captured for various polarization angles of the input light, demonstrating efficient coupling for vertically (TM) polarized light and substantially no coupling for horizontally (TE) polarized light, as shown in Fig. 19(b). The propagation constants of the and modes were also computed via the MoL as a function of t, along with those of the and modes of the corresponding slab, yielding an MPA of for an wide thick stripe.
Theoretical studies using the MoL on the bound modes of asymmetric stripes were reported by Berini [250, 251]. The study in  was conducted at on a wide Ag stripe for a range of t, with the stripe bounded on both sides by optically semi-infinite dielectric claddings having and . None of the fundamental modes, including the and modes, were found to be long-range in this structure. One of the higher-order modes was long range but only near cutoff, and it had a perturbed field distribution. The subsequent study  considered higher-order modes, other widths, and more cases of cladding asymmetry. Modes were found to change character as t was reduced, because elemental modes evolve and merge in or out of supermodes. Long-range modes were again found but only near cutoff and having perturbed field distributions. Other studies involving the (short-range) modes of the asymmetric metal stripe include, for example, [247, 252, 253, 255, 257, 260, 265].
Charbonneau  measured the MPA of the LRSPP at via the cutback method, as well as its output intensity distribution. The waveguides consisted of thick Au stripes, of widths 8, 6 and , embedded in , and the measured MPAs were 2.4, 1.6, and , respectively. The Au stripes in these structures (and in those of ) were deposited on a thick lower cladding thermally oxidized from the underlying Si substrate. The upper cladding consisted of a thick layer of plasma-enhanced chemical vapor deposition covered with index-matching fluid. The claddings in these structures were likely slightly index mismatched and insufficiently thick to allow the modes their full expansion, resulting in higher than expected attenuation. However, the trend of decreasing attenuation with decreasing width was observed, as well as long-range single-mode guidance.
It was surmised in  that narrowing the width of the stripe would lead to reduced polarization sensitivity. This was confirmed by Berini , who modeled via the MoL square cross-section stripes in a homogeneous dielectric background, and showed that the mode had two degenerate orthogonally polarized counterparts in this structure, with their main transverse electric field component polarized along x (TE polarized) and y (TM polarized). The modes were named and , respectively, and modeling revealed that they exhibit the same qualitative evolution with decreasing dimensions as the mode of the rectangular stripe. MPAs in the range of were computed at for stripe dimensions ranging from to with Au used for the stripe and as the background dielectric. Coupling losses to SMF in the range of per facet were computed, with the smallest coupling loss occurring for and largest one occurring for . A trade-off between MPA and coupling loss is necessary, and it was suggested that the case (, coupling loss ) would provide a good compromise.
Nikolajsen et al.  reported insertion loss measurements at for the LRSPP in Au stripes of various lengths, wide and thick, deposited on of benzocyclobutene (BCB) spin coated and cured on Si and covered by of BCB deposited by the same process. They reported MPA values in the range of about and coupling losses to PM-SMF of about . The vertical and horizontal profiles of the output mode were also measured. The MPA of amounts to a increase in propagation length over the corresponding single-interface SPP.
Al-Bader  modeled symmetric Ag stripes in Si at , using the FDM. The width of the stripe was fixed to , and t was varied. The propagation constants of the fundamental modes , , , , and of the first higher-order mode , were computed and compared with the and modes of the corresponding slab. He computes an MPA of about for the LRSPP in stripes thick.
Boltasseva and co-workers [285, 286] report MPA measurements at for wide Au stripes having to and cladded on both sides by of BCB. The lowest MPA measured among their waveguides was for . The coupling loss to SMF was measured for and found to vary from 0.5 to per facet for to , respectively. Mode outputs and profiles were measured for various stripe dimensions. The EIM was proposed for modeling the metal stripe, and propagation constants for the fundamental and first three higher-order long-range modes were computed as a function of w for .
Charbonneau et al.  measured LRSPP mode outputs at for long, thick Au stripes of various widths () on of on Si and covered with index-matched polymer. The mode outputs captured under identical measurement conditions show reduced mode confinement and attenuation with decreasing w.
Zia et al.  applied the EIM to the metal stripe and compared their computed propagation constants for the , , and modes with those of Al-Bader  over his range of t. Good agreement was achieved for the and modes, but not for the mode (the authors reported achieving a better agreement with  for this mode). They also applied the FDM to extend the comparison to other stripe dimensions ( to and ), again achieving good agreement with the EIM for the and modes.
Berini et al.  measured the MPA and coupling efficiency to PM-SMF of the LRSPP at , in waveguides comprising one or many metal stripes, deposited directly (no adhesion layer) onto of on Si, and covered with index-matched polymer. Au stripes 31, 25, and thick, and Ag stripes thick, were used to implement the structures. The lowest MPAs measured among the set of Au and Ag stripes were 0.42 and ( and , respectively). These MPAs are and lower than those of the corresponding Au and Ag single-interface SPPs, respectively. The largest coupling efficiency measured among the set of Au structures was 98%, corresponding to a coupling loss of per facet. Theoretical results were obtained via the MoL for all of the structures characterized. Theory and experiment agreed to within about 5% for all of the 31 and thick Au structures, but a thickness-dependant permittivity had to be assumed in order to achieve agreement to within 12% for the thick Au structures, possibly because films of this thickness are of lower density than the bulk (as discussed in Subsection 2.7).
Charbonneau et al.  characterized the LRSPP in thick, 4 and wide Au stripes on of sputtered on Si and covered by of sputtered (same process). The MPA and coupling efficiency of the LRSPP to SMF were measured at and compared with theoretical expectations (MoL). The MPA was observed to decrease with increasing as expected, and errors between theory and experiment of 4% to 8% were achieved for the wide structure, and of 12% to 17% for the wide structure. The lowest MPA measured was for the wide structure at . All measured coupling efficiencies were greater than 90%, approaching 98% for the wide structure. Horizontal and vertical mode profiles were measured at for both waveguides and compared with far-field diffraction-limited theoretical profiles (computed with the FEM) with near perfect agreement being observed.
Leosson et al.  characterized at narrow ( wide) Au stripes thick, on of BCB and covered with of BCB. The structures were excited via butt coupling to a PM-SMF such that the incident polarization was aligned at 45° to the waveguide axes. Mode outputs and profiles were measured as a function of w through an analyzer (polarizer) aligned either along the vertical (TM) or horizontal (TE) waveguide axes. Figure 20 shows the measured outputs; the stripe aspect ratio is indicated in the white boxes for output pairs a, b, and c. Propagation of both long-range and modes (TM- and TE-polarized, respectively) was observed for the structure having an almost square cross-section , as shown in output pairs c of Fig. 20. The MPA and coupling loss to SMF of the long-range mode were measured as a function of w, yielding about and per facet in the best cases ( and , respectively). The measured MPAs were noted to be larger than the theoretically expected ones , and the discrepancy was attributed to solvable fabrication issues.
Degiron and Smith  modeled symmetric and asymmetric Ag and Au stripes, using a commercial 3D FEM modeling tool. They investigated the effects caused by rounding the corners of the stripe, concluding that the effects on the LRSPP were negligible. Asymmetric stripes were investigated for a few cases of asymmetry and a few t values. They also modeled roughness as randomly distributed metal cylinders, each cylinder having a random height between 0 and and a diameter smaller than its height. They reported that the cylinders cause minor perturbations on the propagation constant of the LRSPP, but act nevertheless as subwavelength scatterers, outcoupling the LRSPP into freely propagating radiation.
Hosseini et al.  modeled symmetric stripes via the FDM with the eigenvalues and eigenvectors of the propagating modes computed by using the Arnoldi method. They compared their results with those reported in , noting good agreement for the and modes.
Boltasseva and Bozhevolnyi  investigated the LRSPP in straight waveguides implemented by using thick Au stripes of width and cladded on both sides by of BCB. MPA and coupling loss measurements were carried out on stripes of different widths over the range with the results following expected trends with w and . The measured MPAs of 10 and wide stripes at were 0.7 and , respectively.
Rao et al.  measured the MPA and mode profile at of the LRSPP in thick Au stripes cladded on both sides by of BCB. The measured MPA of their wide structure was .
Jung et al.  modeled square cross-section Au stripes of dimensions ranging from in BCB at using the FEM. The propagation constant, mode fields and mode size of the first few modes were computed. They found two degenerate LRSPPs polarized along the x and y axes ( and ) for , and they computed for the case . They also discussed field symmetries, mode nomenclature, and the evolution of the modes as the square cross-section changes into a rectangular one and then into a larger square one.
Berini  via the MoL investigated at the effects on the LRSPP of air gaps in various locations adjoining an Au stripe (top, side, and wings) in , finding in general that the gaps were deleterious (as in  for the slab), strongly perturbing its mode fields and causing its MPA and confinement to decrease as the gaps become more invasive such that only nanometric gaps could be tolerated. The FoM [96, 271] decreased with increasing gap size, indicating that confinement decreased more rapidly than attenuation. The mode fields developed strong maxima and localization in the gaps, a feature that could be interesting if high-intensity fields in nanometric air gaps are sought, but only if coupling and radiation losses are essentially irrelevant.
Degiron et al.  enclosed the metal stripe within the core of a dielectric slab and found that the confinement–attenuation trade-off of the LRSPP was favorably altered. They computed that the LRSPP in this structure can propagate over the same distance (or longer) as the conventional one, but with a smaller mode size, using a significantly thinner metal stripe. The LRSPP in this structure is a hybridized SPP–dielectric waveguide mode.
Buckley and Berini  extended the FoMs introduced in  to 2D waveguides and applied them to the LRSPP in symmetric metal stripes, comparing different geometries, metals, and operating wavelengths. Depending on which FoM was used, and hence on how confinement was measured, different preferred designs and operating wavelengths emerged. Each of the metals analyzed showed wavelength regions where their performance was best. They also modeled the LRSPP in the metal stripe embedded within the core of a dielectric slab  as a function of slab thickness, finding increasing attenuation, decreasing mode size, and decreasing effective index with decreasing slab thickness. All of the FoMs were improved for this structure over those of the conventional LRSPP for a good range of slab thickness.
Berini  generated design curves using the MoL for the LRSPP in metal stripes embedded in . Stripe dimensions within the range and were modeled, Au, Ag and Al stripes were compared, the wavelength range covering was considered, and cutoff curves due to index asymmetry were generated. It was found for moderately confining stripes that an index asymmetry of 1 to at most could be tolerated before the onset of radiation.
Kim et al.  measured the attenuation and coupling losses of the LRSPP propagating along Au stripes ( and , ), on of polymer and covered with the same material. The polymers used were UV-curable ZPU12-450 and ZPU12-460 (fluorinated acrylate, Chemoptics), which have refractive indices of about 1.45 and 1.46, respectively. The measurements were obtained at through cutback by using butt-coupled PM-SMFs. They measured decreasing attenuation and increasing coupling losses as the thickness and width of the stripe decreased. They reported an MPA of and total coupling losses of for the smallest stripe (, ) in the lower-index polymer; this value of MPA is lower than that of the corresponding single-interface SPP. The large coupling losses indicate that the mode size is larger than the fiber mode. Lower total coupling losses were measured for wider stripes of the same thickness. They compared their measurements with computations performed by the MoL. In a subsequent paper, Park et al.  reported comparable MPAs (0.14 and ) and coupling losses for similar Au stripes in ZPU 450 at and , respectively.
Salakhutdinov et al.  investigated LRSPP propagation along Au stripes (, ) on a thick BCB lower cladding covered by of BCB, using a scanning near-field optical microscope. The waveguides were excited by butt coupling to PM-SMF over the range of . Salakhutdinov et al. collected aerial scans at various distances above the upper cladding for waveguides having different stripe widths and upper cladding thicknesses. They observed from their scans propagation of the LRSPP and concluded that the scanning near-field optical microscope head was capturing scattered light rather than the evanescent tail of the LRSPP.
Guo and Adato  investigated the LRSPP propagating along a metal stripe bounded on both sides by thin low-index dielectric regions of the same width as the stripe, with thick high-index outer claddings bounding the system. They found via modal analysis (film mode-matching method) that increasing the thickness of the low-index regions decreases the attenuation of the LRSPP, increases its size, and pushes it toward cutoff (as in the 1D case [95, 97]). In a subsequent study on a similar structure , they reported dispersion curves and other computations for the LRSPP as a function of stripe dimensions and as a function of the thickness and index of the thin dielectric regions, including refractive indices that are above and below that of the outer claddings. They investigated the trade-off between confinement and attenuation, concluding that an improvement can be achieved compared with the conventional structure by adding high-index regions, which is consistent with conclusions reached for the metal stripe embedded in the core of a dielectric slab [295, 271]. They also concluded that adding low-index regions does not improve the trade-off compared with the conventional structure.
Jiang et al.  investigated LRSPP propagation along Au stripes (, ) embedded in polymer deposited on of on Si. The polymer used was PFCB (perfluorocyclobutane), which has an index of 1.476 at . The upper and lower claddings had the same thickness, which was varied from 3.5 to ( total thickness) from structure to structure. The thinnest claddings do not allow the LRSPP its full expansion, so in these cases the structure performs as a metal stripe embedded in the core of a dielectric slab . The waveguides were excited at by a butt-coupled PM-SMF and characterized through cutback. They reported MPAs and coupling losses in the range of and per facet. Their lowest MPA of was measured in a wide stripe with thick claddings, and is lower than that of the corresponding single-interface SPP (they mention measuring an MPA of for the same stripe with thicker claddings). They also characterized Au stripes embedded in the center of a thick core of PFCB cladded on each side with of a lower-index polymer . They found for this structure (as for their claddings of limited thickness) that the trade-off between confinement and attenuation is improved, as discussed in [295, 271]. Modal computations were performed by using the FDM, verifying some of their measurements.
Chattopadhyay et al.  proposed an approximate modeling approach for computing the modes of 2D plasmonic waveguides. They demonstrated their approach by application to the symmetric metal stripe and compared their results with those reported in [248, 258, 259].
Park et al.  described a liftoff process for fabricating metal stripes on a thick polymer lower cladding, and fabricated Ag stripes (, ) on of ZPU 450 (Chemoptics) covered by of the same material. They reported MPAs in the range of , and associated mode widths in the range of about , at . The MPA of is lower than that of the corresponding single-interface SPP, but the confinement is extremely weak (mode width ); their more practical structures have MPAs of about .
Berini and Buckley  compared the convergence and accuracy of the MoL (in-house code) and FEM (commercial code) by computing the propagation constant of modes supported by the metal slab, the metal stripe, and the 90° metal corner. A discretization strategy yielding monotonic convergence was demonstrated for both methods, allowing more accurate results to be extrapolated from the convergence histories. Both methods yielded similar anticipated errors for a comparable minimum discretization spacing, with the FEM being slightly more accurate but the MoL requiring less computational effort. Convergence to within an anticipated error of was easily achieved with both methods using a minimum discretization of (except for modes that are highly confined to corners). The percentage difference between the extrapolated results computed from the FEM and MoL convergence histories ranged from 3.75% to 0.012% for the modes of the metal stripe.
Satuby and Orenstein  explored theoretically (finite difference beam propagation method) and experimentally 6 and wide gaps between two thick Au films bounded symmetrically with of BCB at . They claimed to observe TE-polarized long-range modes localized in the gap of each structure, but their modeling of the structure reveals effective indices that are lower than the index of the BCB claddings (so the modes are radiative). Pile et al.  studied the bound modes supported by this structure, finding no evidence of long-range modes and suggesting that the experimentally observed outputs were bulk waves diffracting off metal edges. In a subsequent paper  Satuby and Orenstein observed experimentally propagation at of higher-order LRSPPs in 10 and wide thick Au stripes bounded by similar claddings (they also investigated a structure similar to that of [266, 274]).
3.3. Passive Integrated Structures [283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299, 300, 301, 302]
The attributes of the metal stripe and the mode render them useful as the foundations of an integrated optics technology, and indeed a good deal of work has already been reported toward achieving this goal.
Charbonneau  demonstrated at various passive elements implemented with Au stripes (, ) in ( thick thermal lower cladding, thick PECVD upper cladding covered with index-matched fluid). The elements demonstrated included S bends of varying radii of curvature, sharp angle bends of different angles, Y junctions as mirrored S bends, Mach–Zehnder interferometers (MZIs), and couplers with a variable spacing between the arms. The minimum bending radius was found to be in the range of , and sharp angle bends of 1° produced of radiation loss.
In  Charbonneau et al. measured LRSPP mode outputs at for long passive elements including S bends, Y junctions, and couplers of various designs implemented with Au stripes (, ) on of thermal and covered with index-matched polymer. A minimum bending radius of about was measured.
In a subsequent paper, Charbonneau et al.  characterized at passive elements operating in the LRSPP and implemented with Au stripes (, ) bounded on both sides by of sputtered . S bends, Y junctions, couplers, and MZIs were characterized. A minimum bending radius of about was measured. The experimental results were compared with theoretical results obtained from models constructed for each structure based on decomposition into local modes (including bend modes [284, 290]) and overlap computations. Errors between theory and experiment of about 5% were observed for all of the structures.
Figure 21 collects measured mode outputs at of passive elements operating in the LRSPP and implemented with the metal stripe [283, 287, 288]. Figure 21(a) shows a mosaic of outputs for couplers where the spacing between the parallel stripes is varied from (lower image) to (upper image) in steps of . Figures 21(b), 21(c), 21(d) give outputs for a straight waveguide, an S bend, a Y junction splitter, a MZI, and a sharp angle bend.
Lu and Berini [284, 290] applied the MoL in cylindrical coordinates to model the LRSPP at along curved Ag stripes of various dimensions (, , ) embedded in a dielectric having a refractive index of 2. An absorbing boundary condition was used on the outside (radiating side) of the curve. The 90° insertion and radiation losses, the effective index, mode contours and profiles, and the transition loss to a straight section were computed as a function of the radius of curvature. An optimal radius of curvature minimizing the insertion loss was found, where radiation loss dominates at radii smaller than the optimal radius and propagation loss dominates at radii greater than the optimal radius. The more confining structures led to lower insertion loss bends having a tighter radius of curvature. A stripe having and yielded about of insertion loss for a 90° bend at its optimal radius of curvature of about . Using the same numerical method (MoL in cylindrical coordinates), Berini  reported the computed radiation loss of the LRSPP propagating along 90° curved Au stripes (, , ) in at , and found for moderately confining stripes that radii of curvature in the range of were needed. Berini and Buckley  investigated the convergence and accuracy of this numerical method and found that the absorbing boundary condition used along the radiating side of the curve introduces errors that can further limit the accuracy of the computations. Degiron and Smith  modeled the propagation of a well-confined mode along a curved stripe as a function of radius of curvature by using a commercial 3D FEM modeling tool.
Boltasseva and co-workers [285, 286] investigated passive elements operating in the LRSPP, reporting measurements at for elements implemented with Au stripes in BCB. Multimode long-range guidance was observed in 40 and wide stripes, along with twofold and fourfold images for a centrally aligned excitation. S bends, Y junctions, and couplers of various designs, implemented with or wide stripes were also characterized. They measured a minimum bending radius of about for the wide thick stripe. They also fitted models used in conventional integrated optics to some of their experimental results.
Won et al.  demonstrated broadside couplers of various lengths operating in the LRSPP at , where thick wide Au stripes were deposited one on top of the other, separated by , as sketched in Fig. 22(a). Polymer (ZPU12-470, Chemoptics) was used as the cladding bounding all stripes. The structures exhibited stronger coupling than edge couplers (where the stripes are deposited on the same level—Fig. 21(a)), yielding a shorter length for total coupling. Broadside and edge couplers were also modeled by using the MoL, in support of their experimentation and conclusions. Figure 22(b) plots the normalized power output from both channels (ports) of a series of couplers as a function of interaction length , and Fig. 22(c) shows measured mode outputs for the case .
Boltasseva and Bozhevolnyi  investigated couplers operating in the LRSPP and implemented by using thick Au stripes cladded on both sides by of BCB. Couplers of various interaction lengths having different separations between stripes were designed with and wide stripes. Figure 23(a) shows a sketch of a coupler in top view. Coupling curves were measured, first as a function of interaction length at , then over the range , demonstrating wavelength multiplexing. Figure 23(b) shows the measured outputs for one of their designs (, , ) normalized to the transmittance of a similar straight waveguide, showing E-band light emerging from the direct port and L-band light from the coupled port (the bands are identified as the hashed areas). The structures were excited with a butt-coupled polarization-maintaining photonic crystal fiber carrying either laser light or polarized supercontinuum white light (for the broadband measurements). The measurements were compared with theoretical results generated by using the EIM, with good agreement being achieved.
Kim et al.  modeled the LRSPP propagating along stripes that are curved out of plane (radius of curvature in the plane), using the MoL formulated in cylindrical coordinates. They computed the field distribution and losses (propagation and radiation) at as a function of the radius of curvature for Au stripes of dimensions and embedded in a cladding of index . They found that the radius of curvature should be in the range of to keep radiation losses low.
Liu et al.  used the FEM to model couplers comprised of an Au stripe in positioned alongside a higher-index dielectric core in the same cladding at . They found dimensions allowing efficient coupling between the LRSPP and the TM mode of the dielectric waveguide such that total power coupling can be achieved. They also found that the TE mode of the dielectric waveguide is unaffected by the nearby metal stripe and so propose the device as a polarization splitter. In subsequent papers, Liu et al. [294, 298] investigated theoretically a similar structure consisting of an Au stripe coupled to two dielectric waveguides, one on either side of the stripe, for use as a triple output coupler with a polarization splitting property. In  they explored the design space of the structure, including the effects of offsetting the metal stripe relative to the dielectric waveguides finding that this significantly improves the extinction ratio of the coupler.
Degiron et al.  found via computations with the FEM that the metal stripe embedded at the center of the core of a dielectric slab could support the LRSPP around tighter bends and with lower insertion loss compared with the conventional structure (no dielectric slab). They reported an insertion loss of about at an optimal radius of curvature of for an Au stripe wide and thick embedded in the center of a thick core of BCB bounded on either side by air.
In a subsequent paper , Degiron et al. reported insertion loss measurements at for 90° bends of varying radii operating in the LRSPP, implemented with Au stripes (, , ) embedded in BCB and bounded by on one side and air on the other. A bend was excited via end-fire coupling to a SMF, and its output was redirected by a turning prism, passed through a polarizer, and captured by a detector array mounted to a microscope (the output power was computed by integrating the measured intensity over the output spot). The structures and the experimental setup are sketched in Figs. 24(a), 24(b). Their measured insertion losses, shown in Fig. 24(c) by squares for a thin BCB layer (, , ) and by diamonds for a thick BCB layer (, , ), exhibit a minimum at the optimal radius of curvature [284, 290]. A tighter optimal radius and a lower insertion loss are noted for the thin BCB structure (squares) compared with the thick one (diamonds), which is representative of the conventional metal stripe. Theoretical results computed by using the FEM, plotted as the solid curves in Fig. 24(c), are in very good agreement with the measurements. Figure 24(d) shows measured mode outputs from two bends of radii and in the case of the thin BCB layer (, , for both); radiation loss along the outside of the bend is evident for .
Degiron et al.  investigated couplers consisting of a dielectric core (SU-8) coupled to an Au stripe, both embedded in BCB bounded by on one side and air on the other. They reduced the thickness of the BCB cladding to match the momentum of the LRSPP in the Au stripe to the fundamental TM mode of the SU8 waveguide and showed that this leads to good coupler performance through an analysis of its supermodes (FEM) and through experimental demonstration at . Theoretical projections suggest that 95% power transfer is achievable from the TM mode of the SU8 waveguide to the LRSPP of the Au stripe for light launched into the SU8 waveguide. The structure can be used as an alternative to end-fire excitation of the LRSPP or as a polarization splitter, since TE light launched into the SU8 waveguide is essentially unaffected.
Joo et al.  proposed as an LRSPP waveguide an Au stripe (, ) on a thin low-index polymer core (ZPU12-450, Chemoptics) on a thin Au slab , the structure cladded on both sides by a high-index polymer (ZPU12-470, Chemoptics). Conceptually, the structure builds on coupled broadside stripes  and slabs [74, 100], and the LRSPP investigated is a coupled symmetric supermode of the Au stripe–slab system. Joo et al. investigated the structure experimentally and theoretically, using the FEM. They showed that the structure can be dimensioned such that the LRSPP is readily excited via end-fire coupling to a SMF. They studied the cutoff characteristics of the LRSPP as a function of the thickness and index of the low-index polymer core, showing that the LRSPP remains supported over a broad range of core–cladding index differences, broader than the index difference of the asymmetric metal stripe, but at the expense of greater attenuation. They demonstrated straight waveguides, an S bend, and a Y junction operating in the LRSPP at , commenting that the bending loss may be larger than the symmetric metal stripe.
Xu and Aitchison  investigated experimentally discrete diffraction (beam expansion) across an array of coupled parallel Ag stripes (, , center-to-center stripe separation of ) cladded on both sides by of and operating in the LRSPP. The structure was excited at via free-space end-fire coupling, and the output intensity distribution across the array was measured. The output intensity distribution among the waveguides of the array was found to vary as a function of the angle of incidence (in the plane of the structure, plane of Fig. 16) of the input beam. They found an angle of incidence for which the intensity out of the outermost waveguides of the array is minimized, thus compensating for discrete diffraction. Modeling of the structure was also conducted by using the EIM in combination with a commercial beam propagation tool, as well as a commercial mode solver.
Buckley and Berini  proposed, as an alternative to the metal stripe, a metallodielectric waveguide capable of aggressive bends and long-range propagation. The structure uses a step-index dielectric slab for vertical confinement and a pair of metallic parallel plates disposed alongside the core [Fig. 1(c)] for lateral confinement. They added high-index dielectric plugs near the metal edges to reduce the attenuation. The structure is dimensioned such that all modes that would cause radiation loss in a bend are cut off by the parallel plates. Buckley and Berini modeled the fundamental TE mode of the structure at by using the MoL in cylindrical coordinates and showed that the waveguide can be bent to a vanishing radius of curvature with virtually no loss . Low-loss propagation of the fundamental TE mode along straight sections is also possible by flaring out the waveguide width. In a subsequent paper , Buckley and Berini further explored the design space of this structure and showed that a similar structure used at microwave frequencies, the nonradiative dielectric waveguide, performs comparatively worse when scaled for use at optical wavelengths.
The metal stripe and the LRSPP propagating along it have been explored as an interconnect solution for high-speed optical data transmission over short reaches.
Park and co-workers [276, 304] carried out optical data transmission measurements at via the LRSPP propagating along a straight thick, wide, long Au stripe in polymer (ZPU 450, Chemoptics), reporting a power penalty of relative to back-to-back transmission at a bit error rate (BER) of . Ju et al.  transmitted optical data signals at via the LRSPP in a similar straight waveguide, butt coupled to SMF. The waveguide was long, and its fiber-to-fiber insertion loss was ( MPA, coupling loss per facet). They report a power penalty of relative to back-to-back transmission at a BER of . No significant pulse broadening was measured.
Kim et al.  demonstrated an optical chip-to-chip interconnect operating in the LRSPP, consisting of a straight long array of four parallel Au stripes (, ) embedded in polymer [273, 276], butt coupled to an array of four TM-polarized vertical cavity surface emitting laser (VCSELs) operating at on the transmit side and to an array of four photodetectors on the receive side. The coupling tolerance of the VCSEL array to the array of stripes was investigated, yielding alignment tolerances of about for their design. Optical data was transmitted though the interconnect at a rate of per channel for an aggregate transmission rate of . The BER per channel was estimated at .
Lee et al.  investigated flexible long interconnects operating in the LRSPP, consisting of an thick, wide Ag stripe cladded symmetrically by a pair of thick high-index polymer claddings followed by a pair of thick low-index polymer claddings . The index step between the pair of outer claddings adds vertical confinement [95, 295, 271], allowing sharper bends out of the plane  than in the case of a single thick cladding . The structures were fabricated on a substrate and then detached, forming a flexible ribbonlike structure. Insertion loss measurements at revealed that the out-of-plane bending loss did not change appreciably from the straight value for radii down to . Figure 25(a) shows the interconnect curved out-of-plane at a radius of along with its measured mode output in Fig. 25(b); the mode size is along the horizontal and vertical directions. They also demonstrated operation of the interconnect with a 90° twist applied about the propagation direction, as shown in Figs. 25(c), 25(d).
Kim et al.  reported LRSPP waveguides consisting of Ag ( and ) and Au stripes, ranging in width from 1.5 to , cladded on both sides by of polymer having an index of 1.431 at (ZPU13-430, Chemoptics). The measurements reported at include MPA, mode size, coupling loss to butt-coupled SMF, and alignment tolerance to a TM-polarized VCSEL. The lowest MPA measured among their waveguides was for a thick wide Ag stripe, and the coupling loss of this waveguide to SMF was per facet. Optical transmission measurements at were conducted with the LRSPP propagating along two butt-coupled long Ag stripes for a total length of (, ), aligned to a TM-polarized VCSEL and a photodetector. The BER achieved was estimated at .
Bragg gratings operating in the LRSPP have been implemented by perturbing the metal stripe over a length, either by changing its width symmetrically (step-in-width), its thickness symmetrically (step-in-thickness), or the claddings symmetrically (step-in-index). Corrugating the stripe or adding bumps near the stripe have also been considered.
Jetté-Charbonneau and co-workers [309, 311] explored step-in-width Bragg gratings, constructed by stepping the width of the metal stripe symmetrically about the longitudinal axis over a length , as shown in top view in Fig. 26(a). The unit cell of period Λ is comprised of stripe segments of width and , each of length and , as shown in Fig. 26(b). In  they demonstrated uniform third-order gratings operating in the LRSPP, where the width of an wide Au stripe was stepped to a narrower width of 4, 3, 2 or over a period of . The grating lengths varied from 1 to . The Au stripes were thick, patterned by using contact lithography, deposited onto , and covered with an index-matched polymer. The measured Bragg wavelengths were near , the FWHM reflection bandwidths were in the range of about 0.3 to , and the largest reflectance measured was 66% for a long grating over a bandwidth. The reflectance and bandwidth increased with grating strength. The bandwidth narrowed and the reflectance increased with grating length, then saturated because of attenuation along the structure. The transmittance on and off resonance decreased with grating length for the same reason. A lower transmittance for was evident in the strongest gratings, due to scattering into free radiation, and this asymmetry in the transmittance increased with grating length.
In [309, 315], Jetté-Charbonneau and co-workers characterized similar gratings and proposed and validated a model to represent them. Figure 27 shows the measured transmittance and reflectance spectra of uniform third-order gratings of different length implemented with an Au stripe stepped in width from 8 to in a 50% duty cycle over a period of . The inset to Fig. 27(a) shows a measured LRSPP mode output from one of the gratings, and the inset to Fig. 27(b) shows an atomic force microscope (AFM) scan of a grating section. The largest reflectance measured among the gratings tested was 47.4%, the bandwidths ranged from 0.2 to , and the ’s were near . Jetté-Charbonneau and co-workers modeled the gratings as an equivalent dielectric stack where slices of the stack take on the complex effective index computed via the MoL of the corresponding waveguide cross section. Physical measurements of the grating profile on chip were used to construct an accurate representation of the grating using about 75 slices to model one period. Very good quantitative agreement with the principal features of the measured transmittance and reflectance responses was achieved with this model, except for the transmittance for because of radiation losses not being included in the theory. They find that the effective index perturbation in these gratings was about .
Using this model, Jetté-Charbonneau and co-workers produced theoretical responses for step-in-width gratings of various architectures and duty cycles, including first-order uniform, interleaved, and apodized designs, and third-order uniform designs [309, 318]. An wide thick Au stripe embedded in was assumed as the nominal guiding structure, and the ’s were positioned near . First-order uniform gratings provide the strongest reflectance, 97% over a bandwidth of for gratings a few millimeters long. Third-order uniform gratings provide at most a reflectance of 77% over a bandwidth of . A smaller duty cycle led to a higher reflectance, since the loss of gratings decreases with the duty cycle. Specific ’s were found to be attainable by using interleaving, and apodized designs lowered the sidelobe levels. Off-resonance insertion losses of a few decibels were predicted.
Boltasseva and co-workers [285, 310, 312, 317] investigated step-in-thickness gratings as sketched in Fig. 26(b), which can also be interpreted as the side view of the grating, with and representing the total stripe thickness. In  they demonstrated first-order uniform step-in-thickness Bragg gratings operating in the LRSPP, where the thickness of an wide, thick Au stripe was stepped over a period of to a thickness of by adding of Au on both sides of the stripe. Grating lengths ranged from 20 to . Two metallization–lithography steps were used to make the structures, and the gratings were patterned by using electron-beam lithography. BCB was used as the cladding material. The ’s were near , the bandwidths were about , and the largest reflectance was about 50% for a long grating. A lower transmittance for was noted, which was due to scattering into free radiation. An add–drop filter using gratings and stripes arranged in a zigzag pattern was also demonstrated. (In  they further investigated this zigzag filter, exploring different angles for the zigzag and gratings of different lengths; a bandwidth of was measured for one of the designs.) The gratings were represented as a metal slab having a modulated thickness and were modeled by using the Lippmann–Schwinger equation, which implicitly includes scattering into free radiation. A good correspondence between experimental and theoretical features was noted.
Søndergaard and Bozhevolnyi  investigated these gratings further, via the theoretical model of . They computed the transmission, reflection, scattering (into free radiation), and absorption spectra of three grating structures. They find that scattering occurs over the entire spectrum and that it is greatest for . Longer gratings produce more scattering for , and stronger gratings increase the scattering across the entire spectrum, but especially for . Field plots show scattered fields primarily at small angles from the structure in the forward and backward directions and show that the start, end, and central portions of the gratings may scatter strongly depending on the wavelength. In a subsequent paper , Søndergaard et al. extended their computations to include more grating designs and more comparisons with experiments, noting a good qualitative agreement between both.
Boltasseva and co-workers  investigated a wider design range for these gratings, including thickness steps of 10 to (half added to each side of the stripe), and duty cycles of 0.25 to 0.75. Figure 28 shows the measured reflectance and transmittance spectra of four uniform first-order gratings implemented with an Au stripe (, ) as a function of the thickness step; the inset to Fig. 28(b) shows an AFM scan of a grating section. Measured bandwidths ranged from 5 to , and the largest reflectance was 64% for a long grating. An long grating had 55% reflectance over about . Boltasseva and co-workers estimated the effective index perturbation to be 2 to . They found asymmetry in the transmittance spectrum, and trends with grating strength and length similar to those noted in [311, 312].
Pedersen et al.  fabricated via nanoimprinting wide, thick straight and corrugated Au stripes buried in Mr-I T85 (nanoimprinting resist). The corrugations had a period of , were deep, had different duty cycles, and varied in length from 0.25 to . Broad dips in transmittance at were measured for the corrugated stripes and were found to be deepest for the longest gratings and for a duty cycle of 50%. The measured MPA in straight stripes was about at .
Mu and Huang  investigated theoretically (complex mode-matching method) the propagation of the LRSPP along first-order gratings implemented with a continuous metal slab but where the cladding index is stepped within one period (step-in-index). They considered a few thicknesses of an Au film and index steps that range from 0.02 to 0.06 above that of . They predicted reflectance bandwidths of about and peak reflectances close to 90% for operation at . They predict lower losses than a comparable step-in-thickness grating .
Boltasseva et al.  explored the LRSPP propagating along Au stripes augmented with photonic bandgap structures implemented as a 2D array of Au bumps protruding above and below the stripe. The stripes were thick, and the bumps had a diameter of , a height that varied from 90 to , and were arranged in a 2D triangular lattice having a period in the range from 500 to . Transmission and reflection spectra were obtained for the LRSPP interacting with arrays in various arrangements, including blocks and defects along certain orientations.
Motivated by the fact that the LRSPP is highly sensitive to asymmetry (e.g., ), Breukelaar and co-workers [321, 327] studied the performance of stripes that start out being symmetric and gradually become asymmetric to well beyond the cutoff point of the LRSPP, in order to assess whether cutoff might be suitable for variable attenuation or modulation. Radiation spreading through the structure was modeled via normal mode decomposition, whereby the radiation continuum was represented as a discretized orthonormal basis of radiation modes, and overlaps with end-fire coupled symmetric stripes were used to determine the weighted excitation of the modes. All modes were computed by using the FEM. Radiation loss curves were computed as a function of asymmetry and length for Au stripes in at , showing that extinction ratios beyond were possible over moderate lengths (millimeters) for small index asymmetries in stripes of moderate confinement. Fiber-to-fiber insertion loss measurements were also obtained for a wide, thick, long Au stripe on covered with index-matching fluid, as a function of device temperature controlled by a thermoelectric cooler. Index asymmetries were thermally induced via the large difference in the of the fluid and . Measured extinction curves were compared with the theoretical one, with near perfect quantitative agreement being achieved. The curves were linear over a wide range of asymmetry, and an extinction ratio of was measured for an asymmetry of . The minimum insertion loss occurred on symmetry, as noted earlier in the case of the slab , for the same reasons.
Details of the formulation  were given in , along with results for more cases. was also considered as a cladding material. Waveguide design curves (propagation constant and mode size) and radiation loss curves were computed. Measured mode outputs and radiation loss curves for Au stripes of different lengths on were given and compared with theory with very good quantitative agreement being noted. It was concluded that the initial level of confinement, the length of the stripe, and the index asymmetry were the primary factors influencing the radiation loss and that, other than setting the initial confinement level, the stripe cross section and the cladding materials (e.g., versus ) had little to do with the radiative process.
Using the metal stripe to affect the immediate optical environment, and therefore the propagation characteristics of the LRSPP, through, for example, a material effect in the claddings, has stimulated strong interest; means for connecting electrically to the stripe in an optically noninvasive manner were soon devised. Gagnon and co-workers [322, 323, 330], used arms of the same width and thickness as the waveguide stripe to connect the latter to contact pads, and they isolated the stripe electrically by introducing long gaps as depicted in Fig. 29(a). The effect of the arms and the gaps on the optical insertion loss was found to be negligible.
Gagnon and co-workers described various thermo-optic devices operating in the LRSPP, including variable optical attenuators (VOAs), switches, tunable filters, and tunable external cavity lasers, where the upper and lower claddings in the immediate vicinity of the stripe are heated by passing current along the stripe . They considered structures where both claddings have the same (e.g., positive in the case of glasses or negative in the case of polymers), and where they have ’s of opposite sign (e.g., glass on one side and polymer on the other). The thermo-optic performance of straight stripes, MZIs, Y-junction splitters, couplers, and Bragg gratings was considered. Monitoring the temperature of the stripe through changes in resistivity was also discussed, providing means for controlling the temperature of the chip or for detecting local heating due to the propagation of the LRSPP along it.
Gagnon and co-workers constructed VOAs configured as an electrically connected and straight thick Au stripe [Fig. 29(a)] deposited on of on Si and covered with index-matched polymer [322, 323, 330]. The stripe widths were 4, 6, and , the heated length varied from , and the overall length was . Electrical contact with the pads was achieved by probing. The attenuation mechanism is radiation loss due to thermally induced index asymmetry driven by the passage of current and commensurate heating of the stripe. Figures 29(b), 29(c), 29(d), 29(e) show a sequence of LRSPP outputs for a wide stripe measured at as a function of dissipated electrical power in the stripe . An extinction ratio of about was achieved with a dissipated electrical power of in an wide stripe having a heated length. Figure 29(f) shows the measured extinction ratio for two stripe widths as a function of the dissipated electrical power density . Response times of (fall time) and a few milliseconds (rise time) were measured. The electromigration limit of stripes was estimated experimentally to be by monitoring over a short period of time (minutes) the resistance of stripes subjected to increasing levels of current. The VOAs were driven at a maximum current density equal to about half of this value without any visible or measurable (electrical, optical) degradation in the stripes. The resistivity of stripes was measured as a function of the drive current and temperature, and the simultaneous use of a stripe as a thermal monitor was discussed and demonstrated. Thermal modeling was conducted by using the FDM, revealing that a temperature increase near the stripe of at least was readily achievable with current densities well bellow the electromigration limit. Computed thermal contours plotted over half of an wide stripe are shown in Fig. 29(g); evidently, a very good overlap with the optical mode was achieved.
Nikolajsen et al.  reported similar VOAs, consisting of an wide, thick Au stripe on of BCB on Si covered by of BCB; long gaps were used to isolate the stripe. The VOAs characterized were long and had a heated stripe length in the range of . A extinction ratio was measured at with about of electrical drive power and a response time. Devices were tested over , showing reasonable uniformity . The index of BCB decreases in proportion to its temperature; so as heating is induced and progresses, the effective index of the LRSPP decreases, eventually reaching cutoff.
Leosson et al.  also investigated this VOA concept but using narrower thicker Au stripes in BCB. They measured extinction ratios at of with of electrical drive power in a long device with a long heated section. In a subsequent paper , Leosson et al. reported improved devices implemented with 1 and long electrically contacted square cross-section Au stripes cladded on both sides by of BCB. They reported fiber-to-fiber on-state insertion losses of about and polarization-dependant loss values of about for in the range of . Their best VOA was long and produced about of extinction, with of dissipated electrical power, and a polarization-dependant loss that remained within about over the full attenuation range. Horizontal profiles of the long-range and modes were measured and compared with far-field diffraction-limited theoretical profiles (computed with the FEM) with excellent agreement being achieved.
Park and Song  reported a similar VOA consisting of a wide, thick, long Au stripe on of polymer (ZPU 450 Chemoptics) and covered by of the same polymer. The stripe was not contacted, but an overlying resistive long heater was used, causing an asymmetric thermal distribution along the vertical axis of the device as verified through modeling, and thus a strong index asymmetry that induces mode cutoff and radiation. A measured extinction ratio of was obtained at with of electrical drive power. A response time was measured for a extinction ratio.
Nikolajsen et al.  demonstrated a thermo-optic VOA and a thermo-optic switch based on a MZI and a coupler, respectively, implemented by using wide, thick Au stripes in BCB; long isolation gaps were used. The VOA was long, had a heated stripe length of , and an arm separation of . As shown in Fig. 30(a), deep extinction was measured at for a drive power of . The switch was long and had a coupled stripe separation of . Heating one stripe of the coupler over the interaction length resulted in switching, with a extinction ratio being achieved for a drive power of , as shown in Fig. 30(b). The response time of both devices was .
Bozhevolnyi et al.  exploited the dependence of resistivity on temperature to monitor the power carried by the LRSPP along thick, wide, long Au stripes embedded in BCB. The power absorbed by the stripe due to the propagation of the LRSPP caused its temperature to increase, thus increasing its resistance as monitored in a Wheatstone bridge configuration. The bridge was implemented by using similar stripes to mitigate environmental effects. A linear response was measured for LRSPP powers in the range of about 0.25 to at . A sensitivity of a few microwatts was estimated, and their wavelength response was investigated over the range to .
Berini et al.  explored direct wafer bonding (fusion) as an approach to fabricate metal stripes in cladding materials that are not readily deposited. The approach was applied with Pyrex and z-cut claddings. The structures consisted of a thin Au stripe deposited into a shallow trench etched into one of the claddings, to which another cladding of the same material was directly bonded. In some cases, the claddings were thinned by wafer grinding and polishing to thick layers. The MPA of the LRSPP in an wide thick Au stripe in Pyrex was measured at and found to be about larger than theoretical expectations. Electro-optic LRSPP cutoff measurements were made on -cladded Au stripes, and the results compared favorably with theoretical expectations , indicating that the claddings retained their bulk electro-optic properties.
Mattiussi et al.  described in detail the fabrication of these structures  for the case of z-cut claddings and gave results of physical characterization conducted on completed structures and on structures that underwent intermediate processing steps. The main processing steps consisted of etching a shallow trench on the surface of a wafer via ion beam milling, deposition and liftoff of Au stripes into the trenches, direct bonding to , grinding and polishing thick layers down to , and direct bonding a stack to a Si wafer through intermediate layers of . All processing occurred at low temperature to preserve the embedded Au stripe and mitigate potential problems due to the difference in the coefficient of thermal expansion of and Si. Edge breakaway and point defects in layers were noted after the polishing process, and point defects were noted after the lithography process. Figure 31(a) shows a scanning electron microscope (SEM) cross section of a finished device, and Fig. 31(b) shows an optical microscope image of a device in top view. The inset to Fig. 31(b) shows the measured LRSPP output at of a long device where the Au stripe (, nominally) was buried at about the midpoint of a thick stack.
Jetté-Charbonneau et al.  developed a polishing process to produce high-quality end facets on these heterogeneous devices . The process successfully removed much of the damage caused by dicing and led to a significant reduction in fiber-to-fiber insertion loss, ranging from 4.7 to less loss per waveguide, attributed to less loss per facet. The LRSPP insertion loss was measured at on long structures comprising a 1.2 or wide thick Au stripe buried in the middle of a thick stack of z-cut .
Berini et al.  considered parallel or antiparallel crystal orientations for the z-cut claddings and investigated the implementation of various LRSPP structures theoretically and experimentally. The main theoretical findings were that an MPA of or less was achievable by using Au stripes of reasonable dimensions ( wide, thick); low overlap loss to large and small modes was achievable for the same stripe thickness; radii of curvature in the range of were required for stripes having a moderate confinement, first-order step-in-width gratings having a reflectance of 0.9, a bandwidth of , and electro-optically tunable over a range were achievable; and electro-optically induced index asymmetry of about was sufficient to cut off weakly confined LRSPPs. Measured optical insertion losses for long straight structures varied from 10 to , which is somewhat higher than theoretical expectations. Low-frequency electro-optic mode cutoff measurements produced extinction ratios of and a linear transfer characteristic.
Kim et al.  explored structures where a thin dielectric slab waveguide having a core composed of an electro-optic polymer was modified by depositing a pair of thin Au stripes ( or ) of various separations on top of the bottom cladding. Modes were supported in the gap between the stripes, with vertical confinement provided by the core–cladding index step and horizontal confinement provided by the metal stripes. (This structure is a metal-cladded dielectric waveguide, which differs significantly from those investigated in [266, 274]). A long-range TM mode was observed in a multimode structure. Thermo-optic variable attenuation and electro-optic phase modulation were also observed in some structures.
Nezhad et al.  investigated theoretically the propagation of SPPs along metal structures adjacent to gain media. They computed, using the FEM, that the LRSPP in a thick Ag slab bounded on both sides by semi-infinite InGaAsP could propagate with no loss at if the claddings provide a gain of about . For a wide thick Ag stripe in the same background they found that a gain of about is needed for lossless LRSPP propagation, and they noted that this is about one order of magnitude less than required for the corresponding single-interface SPP. They also computed the gain required from a thin layer, thick, located within one of the InGaAsP claddings, as a function of its proximity to the Ag stripe. For a layer adjoining the stripe they found a required gain in the range of about and pointed out that such gain is available from a stack of quantum wells (for example). They showed that the required gain increases monotonically as the layer moves away from the stripe.
Alam et al.  investigated theoretically at an Ag stripe on a thick barrier followed by five thick quantum wells , each separated by thick barriers on a semi-infinite InP substrate and covered with a finite-thickness superstrate index matched to the barriers. They investigated the effects caused by changing the stripe dimensions and the thickness, index, and loss of the superstrate, finding a good design such that the LRSPP is supported with an MPA of and a good overlap with the quantum wells. For a superstrate with no loss they found a required gain of about for lossless propagation of the LRSPP in their structure. The high-index finite-thickness superstrate taken in combination with air is roughly and effectively index matched to the quantum wells, barriers, and lower-index InP substrate on the other side of the Ag stripe.
Jetté-Charbonneau and Berini  demonstrated an external cavity laser using a TE-polarized InP gain chip coupled to a step-in-width Bragg grating operating in the LRSPP. The gain chip had a high-reflection coating on one end facet, to form one end of the lasing cavity, and a length of PM-SMF butt coupled to its other end facet. The other end of the PM-SMF was polarization aligned and butt coupled to a Bragg grating defining the other end of the lasing cavity. The grating was a long third-order uniform periodic step-in-width design similar to those discussed in [311, 315] (Au on covered with index-matched polymer), providing a 10% reflectance over a bandwidth (FWHM) centered at with of insertion loss on resonance. The onset of lasing occurred at a threshold current of , and the maximum output power was about at . The total cavity length was about (length of the PM-SMF), so lasing occurred in many longitudinal modes, but single-mode lasing is feasible for shorter cavity length. The TM-to-TE polarization extinction ratio of the laser output was .
Ambati et al.  investigated experimentally a metal stripe ( of Au on of Cr, ) cladded on both sides by Er-doped glass. The structure was fabricated via wafer bonding, as in . The optical pump and probe were combined by a fiber multiplexer and injected into the structure via a butt-coupled SMF, thus exciting LRSPPs as copropagating pump and probe waves. The output at the probe wavelength was extracted by using a butt-coupled SMF and a fiber demultiplexer. The fiber-to-fiber insertion loss of an long structure, measured at the probe wavelength without pumping, was . The measured insertion loss decreased as the pump power increased, dropping by for a continuous-wave pump power of . Experiments conducted with pulsed pump and probe signals attribute the loss reduction to stimulated emission into the LRSPP.
Charbonneau et al.  reported surface sensing experiments conducted with the LRSPP propagating along MZIs implemented by using Au stripes on and covered with index-matched oil. BSA was deposited on MZIs via microspotting , either covering the entire MZI or only a portion of one arm. The wavelength response of the MZIs was measured near , and a wavelength shift of was obtained for an MZI having a partially covered arm relative to uncovered and fully covered MZIs, a shift comparable with the response of conventional prism-coupled single-interface SPP sensors under wavelength interrogation . A thick adlayer of BSA accounted for the observed shift, as determined by modeling by using the MoL.
Tencer et al.  developed a technique for depositing and confining small amounts of solution on a flat hydrophilic substrate using the capillarity of a flat high surface free energy guide positioned slightly above. Interfacial forces strongly confined the solution between the guide and the substrate, and the solution followed the movement of the guide over the surface to specific thin Au features thereupon. The thermodynamic background of the method was given, and its application to coat one arm of Au MZIs with BSA was demonstrated. AFM measurements were used to verify BSA coverage on the stripes.
Tencer et al.  conducted an AFM study of Au stripes on covered with BSA adsorbed from phosphate buffer solutions, finding that BSA forms an adlayer of average thickness of . Comparisons with an adsorption model suggested the adsorption of one monolayer followed by denaturation and flattening. The BSA coated stripes had higher roughness measures (RMS, average and peak to valley) than a similar virgin stripe, and the most affected measure was the peak-to-valley roughness, reflecting the nonuniformity of adsorption, possibly due to the formation of aggregates.
Tencer and Berini  described an electrochemical method to selectively desorb thiol self-assembled monolayers from gold surfaces separated by tens of micrometers on a substrate. The potentials of both surfaces were controlled independently with a multichannel potentiostat because the resistance between adjacent surfaces can be much lower in electrolyte than the resistance between one of the surfaces and the counter electrode, and both reductive and oxidative thiol desorption may occur. Desorption potentials were determined for 1-dodecanethiol-based self-assembled monolayers in phosphate buffer, phosphate-buffered saline, and NaOH. Desorption was verified via contact angle measurements. The method should be transferable to electrically contacted MZIs [e.g., Fig. 32(c)] in order to chemically differentiate its arms, with a receptor chemistry applied to one (the sensing arm) and a blocking chemistry applied to the other (the reference arm).
Berini et al.  reported a membrane waveguide capable of propagating LRSPPs in any gaseous or liquid environment, consisting of a large-area free-standing ultrathin dielectric membrane upon which thin metal stripes were deposited, and where the environment (gaseous or liquid) was allowed to surround the structure, in effect becoming the claddings. As long as the membrane was not too invasive optically, the mode propagating along the stripe retained the essential characteristics of the LRSPP with slight localization along the top surface of the stripe opposite the membrane. membranes and ( atoms) thick, clamped around their perimeter to the underlying Si substrate, supporting thick, wide Au stripes, were fabricated and characterized at . The measured MPA of the LRSPP in air and in a liquid having an index close to water was 3.6 and , respectively, which agreed with theory (MoL) to within 7%. The excitation of the LRSPP was accomplished by using prism coupling via the top surface because the waveguide end facets were not accessible. Mode outputs were obtained for a membrane cut in half by using a focused ion beam, as shown in Figs. 32(a), 32(b).
In a subsequent paper, Berini et al.  explored theoretically the design space of the membrane waveguide, considering Au on a thin Cr or Ti adhesion layer as the stripe, or as the membrane, vacuum or as the background, and three operating wavelengths . Good thicknesses were found to be in the range of for the stripe and for a membrane or for a membrane. Computations revealed that the LRSPP in such structures had a moderate MPA of at , good confinement, and a reasonable range extension factor of to 25 when in vacuum or . The adhesion layer was found to modestly impact the MPA of the LRSPP, since the latter develops localization along the opposite metal surface away from the adhesion layer. Structures were fabricated as Au on Cr stripes (25 or thick) on free-standing (20 or thick) large area membranes. Propagation of the LRSPP was observed and characterized in air, in , and in optical fluids of similar index at and . The measurements reported include an optical streak, mode outputs and MPAs. Errors between measured and theoretical (MoL) MPAs were below for all of the structures tested.
Berini et al.  described the fabrication of these membrane waveguides [342, 350] and gave results of physical characterization conducted on completed structures and on structures that underwent intermediate processing steps. The main fabrication steps consisted of timed deep reactive ion etching and TMAH (tetramethylammonium hydroxide) back-side etches to release the membranes, and bilayer liftoff photoresist and metal evaporation processes to define metal stripes and features on the front side. They reported optical, SEM, and AFM images, ellipsometry measurements, and membrane segmentation by focused ion beam. Mechanical probing using a tipless AFM probe revealed that the membranes routinely sustain of point force without shattering. Figure 32(c) shows a microscope image of a fabricated electrically contacted MZI on a membrane.
Charbonneau and Berini  investigated techniques for broadside coupling to LRSPPs and demonstrated them by measuring MPAs along membrane waveguides in air and in fluid. Input coupling was achieved by using a high-index right-angle prism placed in close proximity to the membrane waveguide and launching an input beam such that it struck its base close to the 90° corner (to minimize reradiation of the LRSPP into the prism—see Subsection 2.3). Output coupling using a similar (mirrored) arrangement was investigated. A simple alternative output coupling technique was also devised using a tipless AFM probe placed in contact with the stripe to generate out-of-plane scattering, which was then captured by a nearby multimode fiber. The output coupling was reproducible as the AFM probe and its multimode fiber were moved along the stripe to generate cutback measurements, and so measured MPAs were in very good agreement with theoretical ones. Out-of-plane scattering caused by micrometer-sized particles on the stripe was also investigated. The stability of the experimental setup was assessed and found to be about peak-to-peak over a few minutes at constant temperature by using a reference optical signal.
Charbonneau et al.  proposed a broadside excitation technique using a PM-SMF cleaved at a steep angle and positioned near the top surface of the structure to be excited such that the slow mode of the fiber couples evanescently to the LRSPP propagating on it. The technique is analogous to prism coupling except that it is more easily implemented because there is only one element to manipulate and align (the PM-SMF). They demonstrated the technique at by exciting the LRSPP along membrane waveguides in fluid, observing a streak as the LRSPP propagated, and measuring its MPA by cutback, using an AFM probe with a multimode fiber for output coupling. The measured MPAs agreed very well with the theoretical ones.
Tencer and Berini  described an approach for integrating microfluidic channels and optical input–output coupling means with the membrane waveguide. The solution envisaged consisted of clamping a pair of glass chips with etched channels on either side of the membrane such that the membrane forms a common wall between the channels, and accurately controlling the spacing between the top chip and the metal stripe such that prism coupling of incident light to the LRSPP can occur along its bottom surface. The spacing could be set by a metal ring deposited around the membrane, also serving as a fluidic seal. Tencer and Berini computed the pressure drop required along the channels to drive a laminar flow of as a function of channel cross section and length. They also computed the deflection of the membrane as a function of pressure difference across the membrane, concluding that the difference should be minimized for reliable operation.
Fong et al.  described the fabrication of membrane waveguides using Cytop as the membrane material. Cytop has a low index , slightly above that of , so membranes can be made considerably thicker and possibly more robust than if made with or , while remaining optically not too invasive. The main process steps developed include spin coating and curing thin Cytop layers, depositing Au stripes and features (without an adhesion layer) on Cytop via evaporation and liftoff, and applying a TMAH etch through a back-side mask to release the membranes. Fabricated structures had a released Cytop membrane thickness of about , bearing Au stripes (, ) and components such as MZIs.
In a subsequent paper , Fong et al. investigated the deflection and failure of Cytop membranes subjected to a pressure difference applied across the membranes. The membranes tested had a thickness of about , widths that ranged from to , and lengths of (much greater than their widths). A membrane under test was sealed to a custom chuck having a pinhole, and a pressure difference was applied by changing the volume of air below the membrane by using a manual syringe. The pressure difference was measured by using a manometer, and the membrane deflection was measured by using an interferometric optical profilometer. The Young’s modulus and the residual stress in the membranes were extracted from a simple 1D model fitted to the measured pressure-deflection curves. Cytop membranes having widths less than withstood pressure differences of at least while deflecting by .
As an alternative to membrane structures, Daviau et al.  fabricated Au stripes ( to , ) and integrated components on a thick Cytop lower cladding intended for use with an index-matched aqueous solution used as the upper cladding. The Cytop lower cladding was deposited by using a multilayer spin-coat–cure process, and the Au stripes and features were defined (without an adhesion layer) via electron-beam evaporation and liftoff. Physical characterization, consisting of optical, SEM and AFM inspections and of Cytop index and thickness measurements, was conducted on structures that had undergone intermediate process steps and on finished structures. Dicing was found to produce optical-quality end facets. The LRSPP was excited in finished structures covered with an index-matched glycerol solution ( in de-ionized ) at by butt coupling to a PM-SMF. The MPA of the LRSPP measured through cutback was in very good agreement with theory (FEM), and mode outputs were of very good quality.
Daviau et al.  described a broadside excitation technique using a tapered PM-SMF positioned in physical contact with the top surface of the metal stripe such that the slow mode of the fiber couples through partial modal overlap with the LRSPP propagating thereon. The technique was inspired from  except that it used an off-the-shelf tapered fiber. They demonstrated the technique at by exciting the LRSPP along Au stripes on Cytop  covered with index-matched aqueous buffer and by measuring the MPA via cutback. The measured MPAs agreed well with theory.
Berini  assessed theoretically the potential of various surface plasmon waveguides for bulk and surface (bio)chemical sensing, anticipating their use in a MZI sensor. The sensitivity of a generic MZI implemented with attenuating waveguides (e.g., Fig. 21(e)) was derived, revealing that maximum sensitivity occurs when the sensing length is set equal to the propagation length of the mode used. The MZI sensitivities for bulk and surface sensing were then found to be proportional to the ratio of the corresponding waveguide sensitivity (bulk or surface) to its normalized attenuation, so maximizing a ratio leads to preferred waveguide designs and operating wavelengths. The propagation constant, sensitivities, and ratios were determined for modes in 1D surface plasmon waveguides, such as the SPP in the single-interface [Fig. 1(a)], the mode in the metal clads [Fig. 1(c)] and the mode in three variants of the metal slab [Fig. 1(b)] consisting of a free-standing metal slab, a metal slab on a substrate index-matched to the sensing solution (e.g., [180, 185, 188, 355]), and a metal slab on a thin free-standing dielectric membrane (e.g., ). Au was assumed as the metal, as the sensing solution, and adlayers representative of biochemical matter were placed at all metal/ interfaces. Operating wavelengths spanning the range were considered. It was found that the surface sensitivity in the thin metal claddings was larger than in the single interface, whereas that in the thin metal slab was smaller; but the ratio of surface sensitivity to normalized attenuation was larger in the thin metal slab than in the single interface. Also, the bulk sensitivity in the thin metal clads was larger than in the single interface, whereas that in the thin metal slab was slightly smaller; but the ratio of bulk sensitivity to normalized attenuation was larger in the thin metal slab than in the single interface. It was concluded that the LRSPP in thin metal slabs was very competitive for both bulk and surface sensing using an MZI, given the superior ratios of waveguide sensitivity to normalized attenuation. Preferred wavelengths for surface sensing were found to be near the short-wavelength edge of the Drude region, where detection limits of about were predicted.
4. Prospects for Applications
It is the attributes of a technology that guide it toward applications. But it is competitiveness versus other technologies, on the basis of performance, cost, and size, and other considerations such as market size, time to market, and capital requirements, that determine the level of interest for a particular application. Performing a competitiveness evaluation is nontrivial but necessary for any application that is seriously considered. Doing this exhaustively goes beyond the scope of this paper. But the attributes that make LRSPPs competitive for certain applications are outlined and discussed in the subsections that follow.
4.1. Competitive Context
The single-interface SPP and dielectric waveguide modes are the main (established) competition for the LRSPP. Compared with the single-interface SPP, the LRSPP is less attenuated [28, 29, 30, 248], and because of this exhibits a narrower excitation linewidth and a larger field enhancement in a prism-coupled setup , but is less confined (smaller and larger mode size). The single-interface SPP is TM polarized, as is the LRSPP in the slab and stripe, but the square cross section stripe supports TM and TE LRSPPs [254, 262]. Both are surface waves. The LRSPP derives its competitive edge over the single-interface SPP primarily from its lower attenuation.
Dielectric waveguides support essentially bulk waves, and compared with the LRSPP can provide more confinement for less attenuation and a similar or larger field enhancement in a prism-coupled setup . Confinement in dielectric waveguides comes for free, meaning with essentially no attenuation (in principle and in practice for technologies like on Si). But confinement in any SPP waveguide is always accompanied by the cost of attenuation . Applications must therefore rely on other attributes to achieve competitiveness. The presence of the metal is a fundamental distinction, having implications for fabrication, device architecture and chemical interaction, as highlighted below, and is the main feature from which the LRSPP can derive a competitive edge.
The LRSPP is complementary to the symmetric mode of the metal clads [Fig. 1(c)] and 2D variants thereof, in that the waves capture end points in the confinement–attenuation trade-off . Different applications might be targeted altogether by these waves, but if they do target the same application, then the former (LRSPP) relies on length for a competitive edge and the latter (metal clads) on confinement. Interestingly, both can lead to improved sensor performance compared with the single-interface SPP in the Kretschmann–Raether geometry .
4.2. Nonlinear Optics
The field enhancement and long propagation length of the LRSPP offer improved nonlinear optical interactions over reflection from flat surfaces, over single-interface SPPs in the Otto and Kretschmann–Raether geometries, and over some focused beams in bulk media [136, 137, 138]. However, dielectric waveguides, where the core is used as the nonlinear medium, present stiff competition, since they provide a better confinement–attenuation trade-off than the LRSPP [142, 143]. Also, absorption and heating in the metal may cause thermal nonlinearities to dominate over certain field-based ones , which is troublesome given the large fields needed for nonlinear interactions and thus the typically high intensity of the incident beam. These points explain in part why research on nonlinear optics using LRSPPs (and other SPPs) appears to have waned.
Optical biosensing is dominated by surface plasmons, particularly the single-interface SPP on Au in the Kretschmann–Raether geometry [169, 170], as can be observed from the large number of optical biosensing studies reported using this technology in any given year (e.g., ). The technological reasons for this dominance include good sensor performance (sensitivity, detection limit), the ability to integrate microfluidics, and importantly, the availability of good surface chemistries for Au (e.g., ), which are necessary for interfacing the sensor to the (bio)chemical world. The convergence of all of these requirements is essential in a biosensor and may explain in part why dielectric waveguides have yet to make significant inroads into this application space, even though they offer good sensitivity . Other reasons for the dominance of the single-interface SPP include the fact that this technology was early to market, and user-friendly instruments have been available for about . But the approach is nearing maturity, and further substantive improvement seems difficult to achieve.
The LRSPP is less surface sensitive than the single-interface SPP because the mode overlap with a thin adlayer is lower (e.g., ). This is due to its greater spatial extent and to the fact that it is symmetric; so only half of it interacts with the adlayer. But it is overall sensitivity, detection limit, and the rejection of interference (background fluctuations, nonspecific bindings) that really matter, so LRSPPs may lead to better biosensors via their other attributes.
For instance, the LRSPP in the metal slab exhibits a very narrow resonance (coupling) width in a prism-coupled arrangement; so the location of the resonance can be tracked with higher accuracy, and intensity changes are larger near the resonance, features that may more than compensate for the lower surface sensitivity to yield a larger overall sensitivity . The larger field enhancement  is advantageous for fluorescence-based sensing . The LRSPP in the metal stripe can propagate along integrated optic structures, such as long MZIs , potentially leading to greater overall sensitivity [345, 347]. Also, the ability of the MZI to reference out common interferences and (in principle) nonspecific bindings potentially leads to a lower detection limit. The larger penetration depth of the LRSPP into the sensing medium offers an advantage in that larger analyte (e.g., bacteria) can be sensed . Also, a receptor matrix (e.g., receptor molecules in Dextran) deposited on the metal and extending into the sensing medium increases the binding capacity per unit area, compensating for the lower surface sensitivity . A clear drawback, however, is that index symmetry must be ensured at all times, which leads to increased sensor complexity and/or to tighter constraints on the sensing fluid in terms of composition and temperature. But promising solutions to this problem have been proposed [180, 185, 237, 342, 352], and recent experiments pointing to improved biosensor performance are encouraging.
4.4. Emission, Surface-Enhanced Raman Scattering, and Lasing
The extraction of light from emissive devices through metallic contacts is a problem attracting considerable research effort given its practical importance, especially for organic light emitters, where certain architectures position one or two metallic contacts above and/or below the active region. In such cases, a significant amount of power can be coupled directly from the emitting medium into the SPP modes of these contacts, and unless recovered, represents a loss channel reducing the efficiency of the devices. Both coupled modes of the slab ( and ) have an important role to play in this process and in improving the performance of such devices. Using symmetric (or effectively symmetric) corrugated metal slabs appears to be a promising avenue [200, 204, 211].
An improvement in SERS intensity mediated via the LRSPP is also potentially at hand , although this prediction has apparently not yet been verified experimentally.
Although the gain available from certain media (e.g., semiconductors and dyes) is larger than the attenuation of the LRSPP, it is still unclear whether net amplification and oscillation (lasing) in this mode is possible. The challenge with lasing in the LRSPP is that the mode (slab) or the and modes (stripe) are always supported and more tightly confined, and thus will deplete a portion of the gain region nearest the metal through spontaneous emission and amplified spontaneous emission. Recent computations [214, 217] suggest that these challenges can be overcome in dipolar gain media (dyes), and initial experiments conducted in Er-doped glass are encouraging . If net amplification of the LRSPP were possible, then better biosensors (at least) might result, as well as unforeseen applications. Whether better lasers and amplifiers would be enabled is unknown, but at least a greater degree of polarization is expected from the output .
The basic task of an optical interconnect is to confine and transport electromagnetic energy over a prescribed distance with the lowest possible loss while allowing high-speed optical signaling. High interconnect packing density with low crosstalk is also desirable in certain prospective applications such as on-chip interconnects. Relevant interconnect length scales are meters for board-to-board interconnects, centimeters for chip-to-chip or board level interconnects, millimeters for on-chip global and clocking interconnects, and micrometers for transistor-to-global (on chip) interconnection. As the length scale diminishes, the packing density increases. Except for interconnects requiring very high packing density, the dielectric waveguide (optical fiber or other) seems to satisfy the general requirements very well, especially at longer length scales.
The lowest MPAs reported to date for the LRSPP propagating along a metal stripe of reasonable confinement all cluster around [261, 291, 278, 280] at and . So the LRSPP in the metal stripe might challenge other schemes in certain short-reach (centimeters) medium-density ( waveguides per mm width) applications . Competitive features include single (long-range) mode operation, the ability to transmit electrical signals along the same structure, and the fact that buried metal stripes can be fabricated very inexpensively via, for example, evaporation and liftoff of metal and spin-coated or sputtered claddings. The packing density could be increased relatively easily by stacking waveguides vertically along the perpendicular axis (i.e., along y in Fig. 16). Attaching connectors to the waveguides or aligning them to sources and detectors remains an important part of the cost equation, but the relaxed alignment tolerances of larger-sized LRSPPs could be very helpful in this respect (e.g., [329, 306]). The fact that the LRSPP in the stripe is TM polarized only may or may not be a problem, depending on whether one has control over the transmit portion of the interconnection. Alternatively, polarization-insensitive LRSPPs along square cross-section metal stripes could be used [254, 262].
At the other end of the trade-off scale, the metal clads and 2D variants thereof might be suitable for extremely short very high-density interconnections such as transistor-to-global interconnection [11, 12, 359].
4.6. Integrated Optics
Passive integrated optics is dominated by dielectric waveguides because they provide low fiber-to-fiber insertion loss and polarization independence with adequate confinement. These are crucial requirements for many applications, especially in optical fiber communications. Planar waveguides on Si are particularly well suited.
An attractive feature of integrated optics with the LRSPP is the low cost of fabrication, where only one metal formation cycle (e.g., lithography, metal deposition, liftoff) is required for fabricating all elements. Another attractive point is the ability to define metal stripes and features down to the critical dimension of the litho process. This may be difficult to do with etched core dielectric waveguides because of limitations in producing vertical walls and the need to infill trenches and gaps between cores with upper cladding material.
Accurate fabrication leads to very high-quality Y junctions and highly reflective narrowband step-in-width Bragg gratings . Also, very deep extinction is possible in inteferometric structures such as the MZI because they can be fabricated accurately and the TM-only nature of the structure eliminates unwanted interference occurring with cross-polarized light.
Integrated optics with LRSPPs fits very well with sensing applications, particularly biosensing, where Au waveguides are advantageous given the large user base working on this material and given mature and stable surface chemistries available for Au (as pointed out in Subsection 4.3). Integrated elements might also find applications in interconnects and thermo-optic devices.
4.7. Thermo-Optic and Electro-Optic Devices
The embedded metal slab or stripe is a feature of central importance to thermo-optic or electro-optic devices in that it can also be used as the heating element by passing current through it, or as an electrode for applying electric fields. Combining the heating and LRSPP guiding functions, for example, leads to an excellent overlap between thermal and optical contours [330, 324], and in an optimized design may reduce the electrical power consumption compared with conventional thermo-optics with dielectric waveguides.
5. Concluding Remarks
The LRSPP has appeared in a large number of studies, either as the focal point or peripherally within another context. The studies involving the LRSPP range from the fundamental to the applied and span a very wide landscape including modal studies, prism coupling, field enhancement, corrugated gratings, material characterization, roughness, islandized metal films, nonlinear interactions, biosensors and biosensing, molecular scattering and SERS, transmission through opaque films and emission extraction, fluorescence enhancement, amplification and lasing, waveguides and interconnects, integrated optics (bends, couplers, Y junctions, etc.), Bragg gratings, thermo-, electro- and magneto-optics and associated devices such as VOAs, switches, and modulators.
Despite the breadth and vigor of the research conducted to date with LRSPPs, much remains to be uncovered, and the scope for investigations is broad. Topics involving the LRSPP needing and warranting further study include amplification and lasing, waveguides and integrated optics, corrugated and Bragg gratings, roughness, fluorescence enhancement, emission extraction, SERS, thermo- and electro-optic (bio)chemical sensing, and fabrication techniques. The prospects for applications in biosensing and thermo- or electro-optic devices are bright.
Additional works involving the LRSPP have been reported since this review was written [360, 361, 362, 363, 364, 365, 366, 367, 368]. They include a proposal for a sensor based on the excitation of the LRSPP by a long-period grating defined in the core of a dielectric waveguide placed in proximity and parallel to a thin Au slab ; a proposal for a sensor based on integrating a Bragg grating and a microfluidic channel with a Au stripe coupled to a thin Au slab and operating in the LRSPP as a coupled symmetric supermode of the Au stripe–slab system (building on ); a scheme for tracking the LRSPP coupling angle in a prism-coupled sensor where light outscattered in transmission is detected and used to control an angle-scanning piezomirror ; a study of the LRSPP propagating along a stripe embedded in a Si slab bounded by air at mid- to far-infrared wavelengths ; an investigation of the decay of dipoles into the and (LRSPP) modes of a conformally corrugated Ag grating, including gratings where the mode is suppressed ; the quantum mechanical description of LRSPP excitation in a prism-coupled metal slab ; a demonstration of LRSPP-assisted transmission of nonclassical light (squeezed vacuum states) along an Au stripe buried in BCB ; a proposal for a Schottky contact detector consisting of a symmetric stripe buried in Si capable of detecting LRSPPs at infrared wavelengths (below the bandgap of Si) ; and an investigation of SPP and LRSPP amplification in planar metallic structures incorporating dipolar gain media (R6G in solution) .
Appendix A: Alphabetical List of Acronyms
- 1D One dimensional
- 2D Two dimensional
- AFM Atomic force microscope
- ATR Attenuated total reflection
- BER Bit error rate
- BCB Benzocyclobutene
- BSA Bovine serum albumin
- EIM Effective index method
- DCM 4-Dicyanomethylene-2-methyl-6- (p-dimethylaminostyryl)-4H-pyran
- FDM Finite difference method
- FEM Finite element method
- FoM Figure of merit
- FWHM Full width at half-maximum
- LRSEP Long-range surface exciton polariton
- LRSPP Long-range surface plasmon polariton
- MoL Method of lines
- MPA Mode power attenuation
- MNA 2-Methyl-4-nitroaniline
- MZI Mach–Zehnder interferometer
- PFCB Perfluorocyclobutane
- PMMA Polymethyl methacrylate
- PM-SMF Polarization-maintaining single-mode fiber
- PTS bis-(p-Toluene sulphonate
- Q Quality factor
- RMS Root mean squared
- R6G Rhodamine 6G
- SEM Scanning electron microscope
- SERS Surface-enhanced Raman scattering
- SMF Single mode fiber
- SPP Surface plasmon polariton
- TE Transverse electric
- TEM Transverse electromagnetic
- TM Transverse magnetic
- TMAH Tetramethylammonium hydroxide
- VCSEL Vertical-cavity surface-emitting laser
- VOA Variable optical attenuator
- YIG Yttrium iron garnet
Tables and Figures
1. A. D. Boardman, ed., Electromagnetic Surface Modes (Wiley, 1982).
2. V. M. Agranovich and D. L. Mills, eds., Surface Polaritons: Electromagnetic Waves at Surfaces and Interfaces (North Holland (1982).
3. H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings (Springer, 1988).
4. J. R. Sambles, G. W. Bradbery, and F. Yang, “Optical excitation of surface plasmons: an introduction,” Contemp. Phys. 32, 173–183 (1991). [CrossRef]
5. K. Welford, “Surface plasmon-polaritons and their uses,” Opt. Quantum Electron. 23, 1–27 (1991). [CrossRef]
6. W. L. Barnes, “Surface plasmon-polariton length scales: a route to sub-wavelength optics,” J. Opt. A, Pure Appl. Opt. 8, S87–S93 (2006). [CrossRef]
7. S. A. Maier, Plasmonics: Fundamentals and Applications (Springer, 2007).
9. A. V. Zayats, I. I. Smolyaninov, and A. A. Maradudin, “Nano-optics of surface plasmon polaritons,” Phys. Rep. 408, 131–314 (2005). [CrossRef]
10. S. A. Maier and H. A. Atwater, “Plasmonics: localization and guiding of electromagnetic energy in metal/dielectric structures,” J. Appl. Phys. 98, 011101 (2005). [CrossRef]
13. T. W. Ebbesen, C. Genet, and S. I. Bozhevolnyi, “Surface plasmon circuitry,” Phys. Today 61, 44–50 (2008). [CrossRef]
14. A. Degiron, P. Berini, and D. R. Smith, “Guiding light with long-range plasmons,” Opt. Photonics News 19, 28–34 (2008). [CrossRef]
15. R. H. Ritchie, “Plasma losses by fast electrons in thin films,” Phys. Rev. 106, 874–881 (1957). [CrossRef]
16. T. Turbadar, “Complete absorption of light by thin metal films,” Proc. Phys. Soc. London 73, 40–44 (1959). [CrossRef]
17. K. L. Kliewer and R. Fuchs, “Collective electronic motion in a metallic slab,” Phys. Rev. 153, 498–512 (1967). [CrossRef]
18. A. Otto, “Excitation of nonradiative surface plasma waves in silver by the method of frustrated total reflection,” Z. Phys. 216, 398–410 (1968). [CrossRef]
19. A. Otto, “Excitation by light of and surface plasma waves in thin metal layers,” Z. Phys. 219, 227–233 (1969). [CrossRef]
20. E. N. Economou, “Surface plasmons in thin films,” Phys. Rev. 182, 539–554 (1969). [CrossRef]
21. P. K. Tien, R. Ulrich, and R. J. Martin, “Modes of propagating light waves in thin deposited semiconductor films,” Appl. Phys. Lett. 14, 291–294 (1969). [CrossRef]
22. A. Salwén and L. Stensland, “Spectral filtering possibilities of surface plasma oscillations in thin metal films,” Opt. Commun. 2, 9–13 (1970). [CrossRef]
23. E. Kretschmann, “Die Bestimmung optischer Konstanten von Metallen durch Anregung von Oberfliichenplasmaschwingungen,” Z. Phys. 241, 313–324 (1971). [CrossRef]
24. E. T. Arakawa, M. W. Williams, R. N. Hamm, and R. H. Ritchie, “Effect of damping on surface plasmon dispersion,” Phys. Rev. Lett. 31, 1127–1129 (1973). [CrossRef]
25. F. Abelès and T. Lopez-Rios, “Decoupled optical excitation of surface plasmons at the two surfaces of a thin film,” Opt. Commun. 11, 89–92 (1974). [CrossRef]
26. G. J. Kovacs, “Optical excitation of resonant electromagnetic oscillations in thin films,” Ph.D. Thesis, (University of Toronto, 1977).
27. G. J. Kovacs and G. D. Scott, “Attenuated total reflection angular spectra of a Ag film bounded by dielectric slabs,” Can. J. Phys. 56, 1235–1247 (1978). [CrossRef]
28. G. J. Kovacs, “Optical excitation of surface plasma waves in an indium film bounded by dielectric layers,” Thin Solid Films 60, 33–44 (1979). [CrossRef]
29. M. Fukui, V. C. Y. So, and R. Normandin, “Lifetimes of surface plasmons in thin silver films,” Phys. Status Solidi B 91, K61–K64 (1979). [CrossRef]
30. D. Sarid, “Long-range surface-plasma waves on very thin metal films,” Phys. Rev. Lett. 47, 1927–1930 (1981). [CrossRef]
31. D. Sarid, “Long-range surface-plasma waves on very thin metal films,” (erratum) Phys. Rev. Lett. 48, p. 446 (1982). [CrossRef]
32. Y. Kuwamura, M. Fukui, and O. Tada, “Experimental observation of long-range surface plasmon polaritons,” J. Phys. Soc. Jpn. 52, 2350–2355 (1983). [CrossRef]
35. H. Dohi, Y. Kuwamura, M. Fukui, and O. Tada, “Long-range surface plasmon polaritons in metal films bounded by similar-refractive-index materials,” J. Phys. Soc. Jpn. 53, 2828–2832 (1984). [CrossRef]
37. A. Otto, “The surface polariton resonance in attenuated total reflection,” in Polaritons, E. Burstein and F. De Martini, eds. (Pergamon, 1974).
38. W. F. Chen, G. Ritchie, and E. Burstein, “Excitation of surface electromagnetic waves in attenuated total-reflection prism configurations,” Phys. Rev. Lett. 37, 993–997 (1976). [CrossRef]
40. D. Sarid, R. T. Deck, A. E. Craig, R. K. Hickernell, R. S. Jameson, and J. J. Fasano, “Optical field enhancement by long-range surface-plasma waves,” Appl. Opt. 21, 3993–3995 (1982). [CrossRef] [PubMed]
42. H. Dohi, S. Tago, M. Fukui, and O. Tada, “Spatial dependence of reflected light intensity in ATR geometry: long-range surface plasmon polariton case,” Solid State Commun. 55, 1023–1026 (1985). [CrossRef]
44. Y. Lévy, Y. Zhang, and J. C. Loulergue, “Optical field enhancement comparison between long-range surface plasma waves, and waves induced by resonant cavity,” Opt. Commun. 56, 155–160 (1985). [CrossRef]
46. M. Fukui, S. Tago, and K. Oda, “Characteristics of long-range surface plasmon polaritons excited by fundamental Gaussian beam,” J. Phys. Soc. Jpn. 55, 973–980 (1986). [CrossRef]
47. W. L. Barnes and J. R. Sambles, “Antisymmetric coupled surface-plamon-polaritons in an LB/metal/LB structure,” Opt. Commun. 60, 117–122 (1986). [CrossRef]
48. L. Wendler and R. Haupt, “An improved virtual mode theory of ATR experiments on surface polaritons: application to long-range surface plasmon-polaritons in asymmetric layer structures,” Phys. Status Solidi B 143, 131–147 (1987). [CrossRef]
50. G. Li and S. R. Seshadri, “Excitation of long-range surface polaritons in silver films by a finite-width light beam,” J. Appl. Phys. 67, 555–557 (1989). [CrossRef]
52. T. Sterkenburgh and H. Franke, “Observation of a long range surface mode with a polymer/silver/polymer multilayer,” J. Appl. Phys. 81, 1011–1013 (1997). [CrossRef]
53. M. Mansuripur, A. R. Zakharian, and J. V. Moloney, “Surface plasmon polaritons on metallic surfaces,” Opt. Photonics News 18(4), 44–49 (2007). [CrossRef]
54. J. Chandezon, M. T. Dupuis, G. Cornet, and D. Maystre, “Multicoated gratings: a differential formalism applicable in the entire optical region,” J. Opt. Soc. Am. 72, 839–846 (1982). [CrossRef]
55. T. Inagaki, M. Motosuga, E. T. Arakawa, and J. P. Goudonnet, “Coupled surface plasmons excited by photons in a free-standing thin silver film,” Phys. Rev. B 31, 2548–2550 (1985). [CrossRef]
56. G. S. Agarwal, “Electromagnetic scattering, local field enhancements, and long-range surface plasmons in layered structures with roughness,” Phys. Rev. B 31, 3534–3539 (1985). [CrossRef]
57. T. Inagaki, M. Motosuga, E. T. Arakawa, and J. P. Goudonnet, “Coupled surface plasmons in periodically corrugated thin silver films,” Phys. Rev. B 32, 6238–6245 (1985). [CrossRef]
58. S. Dutta Gupta, G. V. Varada, and G. S. Agarwal, “Surface plasmons in two-sided corrugated thin films,” Phys. Rev. B 36, 6331–6335 (1987). [CrossRef]
59. M. G. Cavalcante, G. A. Farias, and A. A. Maradudin, “Light scattering by thin films with a grating surface,” J. Opt. Soc. Am. B 4, 1372–1378 (1987). [CrossRef]
60. Z. Chen and H. J. Simon, “Attenuated total reflectance from a layered silver grating with coupled surface waves,” J. Opt. Soc. Am. B 5, 1396–1400 (1988). [CrossRef]
61. G. P. Bryan-Brown, S. J. Elston, and J. R. Sambles, “Coupled surface plasmons on silver coated gratings,” Opt. Commun. 82, 1–5 (1991). [CrossRef]
62. N. P. K. Cotter, T. W. Preist, and J. R. Sambles, “Scattering-matrix approach to multilayer diffraction,” J. Opt. Soc. Am. A 12, 1097–1103 (1995). [CrossRef]
63. W. L. Barnes, T. W. Preist, S. C. Kitson, and J. R. Sambles, “Physical origin of photonic energy gaps in the propagation of surface plasmons on gratings,” Phys. Rev. B 54, 6227–6244 (1996). [CrossRef]
64. I. F. Salakhutdinov, V. A. Sychugov, A. V. Tishchenko, B. A. Usievich, O. Parriaux, and F. A. Pudonin, “Anomalous light reflection at the surface of a corrugated thin metal film,” IEEE J. Quantum Electron. 34, 1054–1060 (1998). [CrossRef]
65. I. R. Hooper and J. R. Sambles, “Coupled surface plasmon polaritons on thin metal slabs corrugated on both surfaces,” Phys. Rev. B 70, 045421 (2004). [CrossRef]
66. G. Lévêque and O. J. F. Martin, “Optimization of finite diffraction gratings for the excitation of surface plasmons,” J. Appl. Phys. 100, 124301 (2006). [CrossRef]
67. Z. Chen, I. R. Hooper, and J. R. Sambles, “Coupled surface plasmons on thin silver gratings,” J. Opt. A, Pure Appl. Opt. 10, 015007 (2008). [CrossRef]
68. A. V. Korovin, “Improved method for computing of light–matter interaction in multilayer corrugated structures,” J. Opt. Soc. Am. A 25, 394–399 (2008). [CrossRef]
69. Z. Chen, I. R. Hooper, and J. R. Sambles, “Strongly coupled surface plasmons on thin shallow metallic gratings,” Phys. Rev. B 77, 161405(R) (2008). [CrossRef]
70. A. Sellai and M. Elzain, “Characteristics of a dielectric-metal–dielectric plasmonic waveguide,” Physica E (Amsterdam) 41, 106–109 (2008). [CrossRef]
71. M. S. Tomaš and Z. Lenac, “Thickness dependence of the surface-polariton relaxation rates in a crystal slab,” Solid State Commun. 44, 937–939 (1982). [CrossRef]
74. G. I. Stegeman and J. J. Burke, “Long-range surface plasmons in electrode structures,” Appl. Phys. Lett. 43, 221–223 (1983). [CrossRef]
75. G. I. Stegeman and J. J. Burke, “Effects of gaps on long range surface plasmon polaritons,” J. Appl. Phys. 54, 4841–4843 (1983). [CrossRef]
77. M. S. Tomaš and Z. Lenac, “Long-range surface polaritons in a supported thin metallic slab,” Solid State Commun. 50, 915–918 (1984). [CrossRef]
78. M. S. Tomaš and Z. Lenac, “Coupled surface polariton with guided wave polariton modes in asymmetric metal clad dielectric waveguides,” Opt. Commun. 55, 267–270 (1985). [CrossRef]
79. Z. Lenac and M. S. Tomaš, “Attenuation of long-range surface polaritons in a thin metallic slab with a dielectric coating,” Surf. Sci. 154, 639–657 (1985). [CrossRef]
80. J. J. Burke, G. I. Stegeman, and T. Tamir, “Surface-polariton-like waves guided by thin, lossy metal films,” Phys. Rev. B 33, 5186–5201 (1986). [CrossRef]
81. L. Wendler and R. Haupt, “Long-range surface plasmon-polaritons in asymmetric layer structures,” J. Appl. Phys. 59, 3289–3291 (1986). [CrossRef]
84. F. A. Burton and S. A. Cassidy, “A complete description of the dispersion relation for thin metal film plasmon-polaritons,” J. Lightwave Technol. 8, 1843–1849 (1990). [CrossRef]
86. S. J. Al-Bader and M. Imtaar, “Azimuthally uniform surface-plasma modes in thin metallic cylindrical shells,” IEEE J. Quantum Electron. 28, 525–533 (1992). [CrossRef]
87. S. J. Al-Bader and M. Imtaar, “TM-polarized surface plasma modes on metal-coated dielectric cylinders,” J. Lightwave Technol. 10, 865–872 (1992). [CrossRef]
88. S. J. Al-Bader and M. Imtaar, “Optical fiber hybrid-surface plasmon polaritons,” J. Opt. Soc. Am. B 10, 83–88 (1993). [CrossRef]
89. D. Mihalache, D.-M. Baboiu, M. Ciumac, L. Torner, and J. P. Torres, “Hybrid surface plasmon polaritons guided by ultrathin metal films,” Opt. Quantum Electron. 26, 847–863 (1994). [CrossRef]
90. P. Tournois and V. Laude, “Negative group velocities in metal-film optical waveguides,” Opt. Commun. 137, 41–45 (1997). [CrossRef]
91. C. Chen, P. Berini, D. Feng, V. I. Tolstikhin, S. Tanev, and V. P. Tzolov, “Efficient and accurate numerical analysis of multilayer planar lossy and active optical waveguides in anisotropic media,” Opt. Express 7, 260–272 (2000). [CrossRef] [PubMed]
92. R. Zia, M. D. Selker, P. B. Catrysse, and M. L. Brongersma, “Geometries and materials for subwavelength surface plasmon modes,” J. Opt. Soc. Am. A 21, 2442–2446 (2004). [CrossRef]
93. V. Vaicikauskas, “Fourier transform infrared analysis of long-range surface polaritons excited by the end-fire method,” Thin Solid Films 493, 288–292 (2005). [CrossRef]
94. I. Breukelaar and P. Berini, “Long range surface plasmon-polariton mode cutoff and radiation in 1D waveguides,” J. Opt. Soc. Am. A 23, 1971–1977 (2006). [CrossRef]
98. I. Avrutsky, I. Salakhutdinov, J. Elser, and V. Podolskiy, “Highly confined optical modes in nanoscale metal–dielectric multilayers,” Phys. Rev. B 75, 241402(R) (2007). [CrossRef]
99. F. Liu, Y. Rao, Y.-D. Huang, W. Zhang, and J.-D. Peng, “Abnormal cutoff thickness of long-range surface plasmon polariton modes guided by thin metal films,” Chin. Phys. Lett. 24, 3462–3465 (2007). [CrossRef]
101. C. G. Durfee, T. E. Furtak, R. T. Collins, and R. E. Hollingsworth, “Metal-oxide-semiconductor-compatible ultra-long-range surface plasmon modes,” J. Appl. Phys. 103, 113106 (2008). [CrossRef]
102. L. H. Smith, M. C. Taylor, I. R. Hooper, and W. L. Barnes, “Field profiles of coupled surface plasmon-polaritons,” J. Mod. Opt. 55, 2929–2943 (2008). [CrossRef]
103. T. J. Davis, “Surface plasmon modes in multi-layer thin-films,” Opt. Commun. 282, 135–140 (2009). [CrossRef]
104. B. Yun, G. Hu, and Y. Cui, “Bound modes analysis of symmetric dielectric loaded surface plasmon-polariton waveguides,” Opt. Express 17, 3610–3618 (2009). [CrossRef]
105. G. A. Farias and A. A. Maradudin, “Effect of surface roughness on the attenuation of surface polaritons on metal films,” Phys. Rev. B 27, 7093–7106 (1983). [CrossRef]
106. T. C. Paulick, “Nonradiating electromagnetic oscillations of a thin metal film with rough surfaces,” J. Appl. Phys. 64, 1384–1397 (1988). [CrossRef]
107. M.-L. Thèye, “Investigation of the optical properties of Au by means of thin semitranparent films,” Phys. Rev. B 2, 3060–3078 (1970). [CrossRef]
108. C. Reale, “Optical constants of vacuum deposited thin metal films in the near infrared,” Infrared Phys. 10, 175–181 (1970). [CrossRef]
109. P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 6, 4370–4379 (1972). [CrossRef]
110. P. Gadenne and G. Vuye, “In situ determination of the optical and electrical properties of thin films during their deposition,” J. Phys. E 10, 733–736 (1977). [CrossRef]
111. E. D. Palik, ed., Handbook of Optical Constants of Solids (Academic, 1985).
112. T. Inagaki, J. P. Goudonnet, P. Royer, and E. T. Arakawa, “Optical properties of silver island films in the attenuated-total-reflection geometry,” Appl. Opt. 25, 3635–3639 (1986). [CrossRef] [PubMed]
113. M. Fukui and K. Oda, “Studies on metal film growth through instantaneously observed attenuated total reflection spectra,” Appl. Surf. Sci. 33/34, 882–889 (1988). [CrossRef]
114. M. Yano, M. Fukui, M. Haragichi, and Y. Shintani, “In situ and real-time observation of optical constants of metal films during growth,” Surf. Sci. 227, 129–137 (1990). [CrossRef]
115. W.-J. Lee, J.-E. Kim, H. Y. Park, S. Park, M.-S. Kim, J. T. Kim, and J. J. Ju, “Optical constants of evaporated gold films measured by surface plasmon resonance at telecommunication wavelengths,” J. Appl. Phys. 103, 073713 (2008). [CrossRef]
117. M. Takabayashi, H. Shiba, M. Haraguchi, and M. Fukui, “Studies on surface polaritons in ultrathin films sandwiched by identical dielectrics,” J. Phys. Soc. Jpn. 61, 2550–2556 (1992). [CrossRef]
118. M. Takabayashi, H. Shiba, M. Haraguchi, and M. Fukui, “Surface optic waves in ultrathin films bounded by identical slightly lossy materials,” J. Phys. Soc. Jpn. 62, 2719–2724 (1993). [CrossRef]
119. Z.-C. Wu, E. T. Arakawa, T. Inagaki, T. Thundat, and L. J. Schowalter, “Experimental observations of a long-range surface mode in metal island films,” Phys. Rev. B 49, 7782–7785 (1994). [CrossRef]
120. E. L. Wood, J. R. Sambles, F. A. Pudonin, and V. Yakolev, “Degenerate long range surface modes, supported on thin nickel films,” Opt. Commun. 132, 212–216 (1996). [CrossRef]
121. M. Takabayashi, M. Haraguchi, and M. Fukui, “Propagation length of long-range surface optic waves in islandized silver films,” J. Mod. Opt. 44, 119–125 (1997). [CrossRef]
122. T. Kume, T. Kitigawa, S. Hayashi, and K. Yamamoto, “Long-range surface modes supported by composite thin films,” Surf. Sci. 395, 23–29 (1998). [CrossRef]
123. G. J. Kovacs, “Surface polariton in the ATR angular spectra of a thin iron film bounded by dielectric layers,” J. Opt. Soc. Am. 68, 1325–1332 (1978). [CrossRef]
125. F. Yang, J. R. Sambles, and G. W. Bradberry, “Long-range surface modes supported by thin films,” Phys. Rev. B 44, 5855–5872 (1991). [CrossRef]
126. F. Yang, J. R. Sambles, and G. W. Bradberry, “Prism coupling to long-range coupled-surface modes,” J. Mod. Opt. 38, 707–717 (1991). [CrossRef]
127. G. P. Bryan-Brown, F. Yang, G. W. Bradberry, and J. R. Sambles, “Prism and grating coupling to long-range coupled-surface exciton-polaritons,” J. Opt. Soc. Am. B 8, 765–769 (1991). [CrossRef]
128. B. Prade, J. Y. Vinet, and A. Mysyrowicz, “Guided optical waves in planar heterostructures with negative dielectric constant,” Phys. Rev. B 44, 13556–13572 (1991). [CrossRef]
129. G. P. Bryan-Brown, S. J. Elston, and J. R. Sambles, “Coupled surface excitations on palladium-coated gratings,” Thin Solid Films 207, 252–257 (1992). [CrossRef]
130. R. J. Crook, F. Yang, and J. R. Sambles, “Long-range optical modes supported by a strongly absorbing thin organic film,” J. Opt. Soc. Am. B 10, 237–243 (1993). [CrossRef]
131. R. J. Crook, F. Yang, and J. R. Sambles, “An investigation of p- and s-polarized long-range optical modes supported by a strongly absorbing normal uniaxial thin organic film,” J. Mod. Opt. 40, 243–253 (1993). [CrossRef]
132. M. Takabayashi, M. Haraguchi, and M. Fukui, “Propagation length of guided waves in lossy Si film sandwiched by identical dielectrics,” J. Opt. Soc. Am. B 12, 2406–2411 (1995). [CrossRef]
133. F. Yang and J. R. Sambles, “Determination of the optical permittivity and thickness of absorbing films using long range modes,” J. Mod. Opt. 44, 1155–1163 (1997). [CrossRef]
135. D. Sarid, R. T. Deck, and J. J. Fasano, “Enhanced nonlinearity of the propagation constant of a long-range surface-plasma wave,” J. Opt. Soc. Am. 72, 1345–1347 (1982). [CrossRef]
136. G. I. Stegeman, J. J. Burke, and D. G. Hall, “Nonlinear optics of long range surface plasmons,” Appl. Phys. Lett. 41, 906–908 (1982). [CrossRef]
137. R. T. Deck and D. Sarid, “Enhancement of second-harmonic generation by coupling to long-range surface plasmons,” J. Opt. Soc. Am. 72, 1613–1617 (1982). [CrossRef]
138. J. C. Quail, J. G. Rako, H. J. Simon, and R. T. Deck, “Optical second-harmonic generation with long-range surface plasmons,” Phys. Rev. Lett. 50, 1987–1989 (1983). [CrossRef]
139. G. I. Stegeman, C. Liao, and C. Karaguleff, “Second harmonic generation by oppositely traveling long range surface polaritons,” Opt. Commun. 46, 253–256 (1983). [CrossRef]
140. R. Moshrefzadeh, R. Fortenberry, C. Karaguleff, G. I. Stegeman, N. E. Van Wijck, and W. M. Hetherington III, “Second harmonic generation by monolayers using long range surface plasmon excitation,” Opt. Commun. 46, 257–259 (1983). [CrossRef]
141. G. I. Stegeman and C. Karaguleff, “Degenerate four wave mixing with long range surface plasmons in attenuated total reflection geometries,” J. Appl. Phys. 54, 4853–4855 (1983). [CrossRef]
142. C. Liao, P. Bundman, and G. I. Stegeman, “Second harmonic generation with surface guided waves in signal processing geometries,” J. Appl. Phys. 54, 6213–6217 (1983). [CrossRef]
143. C. Karaguleff and G. I. Stegeman, “Degenerate four-wave mixing with surface guided waves,” IEEE J. Quantum Electron. QE-20, 716–722 (1984). [CrossRef]
145. G. I. Stegeman, J. D. Valera, and C. T. Seaton, “Nonlinear s-polarized surface plasmon polaritons,” Solid State Commun. 53, 293–297 (1984). [CrossRef]
146. J. C. Quail and H. J. Simon, “Optical second-harmonic generation in transmission with long-range surface plasmons,” J. Opt. Soc. Am. B 1, 317–319 (1984). [CrossRef]
147. J. C. Quail and H. J. Simon, “Second-harmonic generation with phase-matched long-range and short-range surface plasmons,” J. Appl. Phys. 56, 2589–2591 (1984). [CrossRef]
148. J. M. Nunzi and D. Ricard, “Optical phase conjugation and related experiments with surface plasma waves,” Appl. Phys. B 35, 209––216 (1984). [CrossRef]
149. M. Fukui, S. Tago, H. Dohi, and O. Tada, “Optical second harmonic generation by long-range surface plasmon polaritons excited by a finite 1-D beam,” Appl. Opt. 24, 1220–1223 (1985). [CrossRef] [PubMed]
150. J. Ariyasu, C. T. Seaton, G. I. Stegeman, A. A. Maradudin, and R. F. Wallis, “Nonlinear surface polaritons guided by metal films,” J. Appl. Phys. 58, 2460–2466 (1985). [CrossRef]
151. R. K. Hickernell and D. Sarid, “Optical bistability using prism-coupled, long-range surface plasmons,” J. Opt. Soc. Am. B 3, 1059–1069 (1986). [CrossRef]
152. G. S. Agarwal and S. Dutta Gupta, “Exact results on optical bistability with surface plasmons in layered media,” Phys. Rev. B 34, 5239–5243 (1986). [CrossRef]
153. D. Mihalache, D. Mazilu, and F. Lederer, “Nonlinear TE-polarized surface plasmon polaritons guided by metal films,” Opt. Commun. 59, 391–394 (1986). [CrossRef]
155. J. R. Sambles and R. A. Innes, “A comment on nonlinear optics using surface plasmon-polaritons,” J. Mod. Opt. 35, 791–797 (1988). [CrossRef]
156. A. D. Boardman and T. Twardowski, “Theory of nonlinear interaction between TE and TM waves,” J. Opt. Soc. Am. B 5, 523–528 (1988). [CrossRef]
157. D. Mihalache, D. Mazilu, M. Bertolotti, and C. Sibila, “Exact solution for transverse electric polarized nonlinear guided waves in saturable media,” J. Mod. Opt. 35, 1017–1027 (1988). [CrossRef]
158. H. J. Simon and Z. Chen, “Optical second-harmonic generation with grating-coupled surface plasmons from a quartz-silver-quartz grating structure,” Phys. Rev. B 39, 3077–3085 (1989). [CrossRef]
159. C.-C. Tzeng and J. T. Lue, “Optical second-harmonic generation from silver films with long-range surface-plasmon excitation,” Surf. Sci. 216, 579–586 (1989). [CrossRef]
162. A. D. Boardman and T. Twardowski, “Transverse-electric and transverse-magnetic waves in nonlinear isotropic waveguides,” Phys. Rev. B 39, 2481–2492 (1989). [CrossRef]
163. H. J. Simon, Y. Wang, L.-B. Zhou, and Z. Chen, “Coherent backscattering of optical second-harmonic generation with long-range surface plasmons,” Opt. Lett. 17, 1268–1270 (1992). [CrossRef] [PubMed]
164. F. Yang and J. R. Sambles, “Experimental observation of thermo-optic nonlinearity with a long-range surface wave supported by a very thin palladium film,” J. Mod. Opt. 39, 485–498 (1992). [CrossRef]
165. Y. Wang and H. J. Simon, “Coherent backscattering of optical second-harmonic generation in silver films,” Phys. Rev. B 47, 13695–13699 (1993). [CrossRef]
166. J.-T. Lue and C. Dai, “Nonlinear optical generation from multiple layers of metal–dielectric films,” Phys. Rev. B 47, 13653–13658 (1993). [CrossRef]
167. J. E. Sipe, V. C. Y. So, M. Fukui, and G. I. Stegeman, “Analysis of second harmonic generation at metal surfaces,” Phys Rev B 21, 4389–4402 (1980). [CrossRef]
168. Y. Li and X. Zhang, “SPM of nonlinear surface plasmon waveguides,” Opt. Commun. 281, 5009–5013 (2008). [CrossRef]
169. B. Liedberg, C. Nylander, and I. Lundstrom, “Surface plasmon resonance for gas detection and biosensing,” Sens. Actuators 4, 299–304 (1983). [CrossRef]
171. J. Homola, S. S. Yee, and G. Gauglitz, “Surface plasmon resonance sensors: review,” Sens. Actuators B 54, 3–15 (1999). [CrossRef]
172. F.-C. Chien and S.-J. Chen, “A sensitivity comparison of optical biosensors based on four different surface plasmon resonance modes,” Biosens. Bioelectron. 20, 633–642 (2004). [PubMed]
173. J. Homola, “Surface plasmon resonance sensors for detection of chemical and biological species,” Chem. Rev. (Washington, D.C.) 108, 462–493 (2008). [CrossRef]
175. S. Glasberg, A. Sharon, D. Rosenblatt, and A. A. Friesem, “Long-range surface plasmon resonances in grating-waveguide structures,” Appl. Phys. Lett. 70, 1210–1212 (1997). [CrossRef]
176. N. M. Lyndin, I. F. Salakhutdinov, V. A. Sychugov, B. A. Usievich, F. A. Pudonin, and O. Parriaux, “Long-range surface plasmons in asymmetric layered metal–dielectric structures,” Sens. Actuators B 54, 37–42 (1999). [CrossRef]
177. C.-H. Liao, C.-M. Lee, Y.-T. Cheng, J.-S. Shyu, and W.-K. Su, “Resonant properties of long-range surface plasmons in an arbitrary multilayer structure,” Jpn. J. Appl. Phys., Part 1 38, 5938–5944 (1999). [CrossRef]
178. G. G. Nenninger, P. Tobiška, J. Homola, and S. S. Yee, “Long-range surface plasmons for high-resolution surface plasmon resonance sensors,” Sens. Actuators B 74, 145–151 (2001). [CrossRef]
179. F. Pigeon, I. F. Salakhutdinov, and A. V. Tishchenko, “Identity of long-range surface plasmons along asymmetric structures and their potential for refractometric sensors,” J. Appl. Phys. 90, 852–859 (2001). [CrossRef]
181. A. Kasry and W. Knoll, “Long range surface plasmon fluorescence spectroscopy,” Appl. Phys. Lett. 89, 101106 (2006). [CrossRef]
183. R. Slavík and J. Homola, “Simultaneous excitation of long and short range surface plasmons in an asymmetric structure,” Opt. Commun. 259, 507–512 (2006). [CrossRef]
184. R. Slavík, J. Homola, and H. Vaisocherová, “Advanced biosensing using simultaneous excitation of short and long range surface plasmons,” Meas. Sci. Technol. 17, 932–938 (2006). [CrossRef]
185. R. Slavík and J. Homola, “Ultrahigh resolution long range surface plasmon-based sensor,” Sens. Actuators B 123, 10–12 (2007). [CrossRef]
186. A. Rajan, K. Sharma, and B. D. Gupta, “Fibre optic sensor based on long-range surface plasmon resonance: a theoretical analysis,” J. Opt. A, Pure Appl. Opt. 9, 682–687 (2007). [CrossRef]
187. C.-W. Chen, C.-H. Lin, H.-P. Chiang, Y.-C. Liu, P. T. Leung, and W. S. Tse, “Temperature dependence of the sensitivity of a long-range surface plasmon optical sensor,” Appl. Phys. A 89, 377–380 (2007). [CrossRef]
188. J. Dostálek, A. Kasry, and W. Knoll, “Long range surface plasmons for observation of biomolecular binding events at metallic surfaces,” Plasmonics 2, 97–106 (2007). [CrossRef]
189. J. T. Hastings, J. Guo, P. D. Keathley, P. B. Kumaresh, Y. Wei, S. Law, and L. G. Bachas, “Optimal self-referenced sensing using long- and short-range surface plasmons,” Opt. Express 15, 17661–17672 (2007). [CrossRef] [PubMed]
191. R. Jha, R. K. Verma, and B. D. Gupta, “Surface plasmon resonance-based tapered fiber optic sensor: sensitivity enhancement by introducing a Teflon layer between core and metal layer,” Plasmonics 3, 151–156 (2008). [CrossRef]
192. W. Knoll, A. Kasry, F. Yu, Y. Wang, A. Brunsen, and J. Dostálek, “New concepts with surface plasmons and nano-biointerfaces,” J. Nonlinear Opt. Phys. Mater. 17, 121–129 (2008). [CrossRef]
193. Y. Wang, J. Dostálek, and W. Knoll, “Long range surface plasmon-enhanced fluorescence spectroscopy for the detection of aflatoxin in milk,” Biosens. Bioelectron. 24, 2264–2267 (2009). [CrossRef]
194. J. Dostálek, R. F. Roskamp, and W. Knoll, “Coupled long range surface plasmons for the investigation of thin films and interfaces,” Sens. Actuators B 139, 9–12 (2009). [CrossRef]
195. M. Vala, S. Etheridge, J. A. Roach, and J. Homola, “Long-range surface plasmons for sensitive detection of bacterial analytes,” Sens. Actuators B 139, 59–63 (2009). [CrossRef]
197. I. Pockrand, A. Brillante, and D. Möbius, “Nonradiative decay of excited molecules near a metal surface,” Chem. Phys. Lett. 69, 499–504 (1980). [CrossRef]
198. W. L. Barnes, “Fluorescence near interfaces: the role of photonic mode density,” J. Mod. Opt. 45, 661–669 (1998). [CrossRef]
200. R. W. Gruhlke and D. G. Hall, “Transmission of molecular fluorescence through a thin metal film by surface plasmons,” Appl. Phys. Lett. 53, 1041–1042 (1988). [CrossRef]
201. P. T. Leung, Y. S. Kim, and T. F. George, “Decay of molecules at corrugated thin metal films,” Phys. Rev. B 39, 9888–9893 (1989). [CrossRef]
202. Z. Lenac and M. S. Tomaš, “Absorption of surface polaritons by molecules near the surface of a metallic slab,” Solid State Commun. 61, 261–265 (1987). [CrossRef]
203. Z. Lenac and M. S. Tomaš, “Absorption of surface polaritons by molecules near the surface of a metallic slab,” (erratum), Solid State Commun. 62, i (1987).
204. M. S. Tomaš and Z. Lenac, “Scattering of surface polaritons by molecules near a metallic slab surface,” Surf. Sci. 189/190, 543–547 (1987). [CrossRef]
205. Z. Lenac and M. S. Tomaš, “Enhanced molecular fluorescence mediated by long-range surface plasmons,” Surf. Sci. 215, 299–318 (1989). [CrossRef]
206. M. S. Tomaš and Z. Lenac, “Enhanced Raman scattering in four-layered ATR configurations,” Surf. Sci. 251/252, 310–313 (1991). [CrossRef]
207. W. L. Barnes, “Electromagnetic crystals for surface plasmon polaritons and the extraction of light from emissive devices,” J. Lightwave Technol. 17, 2170–2182 (1999). [CrossRef]
209. S. Wedge, J. A. E. Wasey, W. L. Barnes, and I. Sage, “Coupled surface plasmon-polariton mediated photoluminescence from a top-emitting organic light-emitting structure,” Appl. Phys. Lett. 85, 182–184 (2004). [CrossRef]
210. S. Wedge, S. H. Garrett, I. Sage, and W. L. Barnes, “Photoluminescence emission through thin metal films via coupled surface plasmon-polaritons,” J. Mod. Opt. 52, 833–843 (2005). [CrossRef]
211. S. Wedge, A. Giannattasio, and W. L. Barnes, “Surface plasmon-polariton mediated emission of light from top-emitting organic light-emitting diode type structures,” Org. Electron. 8, 136–147 (2007). [CrossRef]
213. T. Okamoto, F. H’Dhili, and S. Kawata, “Towards plasmonic band gap laser,” Appl. Phys. Lett. 85, 3968–3970 (2004). [CrossRef]
214. G. Winter, S. Wedge, and W. L. Barnes, “Can lasing at visible wavelengths be achieved using the low-loss long-range surface plasmon-polariton mode?,” New J. Phys. 8, 125 (2006). [CrossRef]
215. Y. Wang and Z. Zhou, “Silicon amplifier based on surface-plasmon-polariton enhancement,” Appl. Phys. Lett. 91, 053504 (2007). [CrossRef]
216. N.-F. Chiu, C.-W. Lin, J.-H. Lee, C.-H. Kuan, K.-C. Wu, and C.-K. Lee, “Enhanced luminescence of organic/metal nanostructure for grating coupler active long-range surface plasmonic device,” Appl. Phys. Lett. 91, 083114 (2007). [CrossRef]
217. I. De Leon and P. Berini, “Theory of surface plasmon-polariton amplification in planar structures incorporating dipolar gain media,” Phys. Rev. B 78, 161401(R) (2008). [CrossRef]
218. D. Sarid, “Enhanced surface-magnetoplasma interactions in a semiconductor,” Phys. Rev. B 29, 2344–2346 (1984). [CrossRef]
219. D. Sarid, “Enhanced magnetic interaction of surface magnetoplasmon polaritons,” IEEE J. Quantum Electron. QE-20, 943–948 (1984). [CrossRef]
222. W. Johnstone, G. Stewart, T. Hart, and B. Culshaw, “Surface plasmon polaritons in thin metal films and their role in fibre optic polarizing devices,” J. Lightwave Technol. 8, 538–544 (1990). [CrossRef]
224. M. N. Zervas, “Surface plasmon-polariton fibre-optic polarizers using thin-nickel films,” IEEE Photon. Technol. Lett. 2, 253–256 (1990). [CrossRef]
225. M. N. Zervas, “Surface plasmon-polariton fier-optic polarizers using thin chromium films,” IEEE Photon. Technol. Lett. 2, 597–599 (1990). [CrossRef]
226. Y. Wang and H. J. Simon, “Electrooptic reflection with surface plasmons,” Opt. Quantum Electron. 25, S925–S933 (1993). [CrossRef]
228. S. Barcelos, M. N. Zervas, and P. St. J. Russell, “Selective excitation of fiber-modes using surface plasmons,” IEEE Photon. Technol. Lett. 7, 1051–1053 (1995). [CrossRef]
229. P. J. Kajenski, “Tunable optical filter using long-range surface plasmons,” Opt. Eng. 36, 1537–1541 (1997). [CrossRef]
230. C.-H. Liao and J.-S. Shyu, “Resonance modulations of surface plasmons in arbitrary modified Otto and Sarid structures,” Jpn. J. Appl. Phys., Part 1 39, 6334–6338 (2000). [CrossRef]
231. K. T. Kim, J. W. Moon, and J. P. Mah, “Side-polished polarization-maintaining fibre coupler with a thin metal interlay and its application as a polarization splitter,” Jpn. J. Appl. Phys., Part 1 43, 598–601 (2004). [CrossRef]
232. T.-W. Lee and S. K. Gray, “Regenerated surface plasmon polaritons,” Appl. Phys. Lett. 86, 141105 (2005). [CrossRef]
233. B. Sepúlveda, L. M. Lechuga, and G. Armelles, “Magnetooptic effects in surface-plasmon-polaritons slab waveguides,” J. Lightwave Technol. 24, 945–955 (2006). [CrossRef]
234. D. Kumar, V. K. Sharma, and K. N. Tripathi, “Design and fabrication of multilayer metal-clad dielectric surface plasmon waveguide polarizers,” Opt. Eng. 45, 054601 (2006). [CrossRef]
235. S. W. Liu and M. Xiao, “Electro-optic switch in ferroelectric thin films mediated by surface plasmons,” Appl. Phys. Lett. 88, 143512 (2006). [CrossRef]
236. J. B. Khurgin, “Optical isolating action in surface plasmon polaritons,” Appl. Phys. Lett. 89, 251115 (2006). [CrossRef]
237. V. N. Konopsky and E. V. Alieva, “Long-range propagation of plasmon polaritons in a thin metal film on a one-dimensional photonic crystal surface,” Phys. Rev. Lett. 97, 253904 (2006). [CrossRef]
239. A. A. Krokhin, A. Neogi, and D. McNeil, “Long-range propagation of surface plasmons in a thin metallic film deposited on an anisotropic photonic crystal,” Phys. Rev. B 75, 235420 (2007). [CrossRef]
240. J. Plumridge and C. Phillips, “Ultralong-range plasmonic waveguides using quasi-two-dimensional metallic layers,” Phys. Rev. B 76, 075326 (2007). [CrossRef]
241. Z. Sun, “Vertical dielectric-sandwiched thin metal layer for compact, low-loss long range surface plasmon waveguiding,” Appl. Phys. Lett. 91, 111112 (2007). [CrossRef]
243. J. M. Montgomery and S. K. Gray, “Enhancing surface plasmon polariton propagation lengths via coupling to asymmetric waveguide structures,” Phys. Rev. B 77, 125407 (2008). [CrossRef]
244. Z. Wu, J. W. Haus, Q. Zhan, and R. L. Nelson, “Plasmonic notch filter design based on long-range surface plasmon excitation along metal grating,” Plasmonics 3, 103–108 (2008). [CrossRef]
246. P. Berini, “Plasmon-polariton modes guided by a metal film of finite width,” Opt. Lett. 24, 1011–1013 (1999). [CrossRef]
247. J.-C. Weeber, A. Dereux, C. Girard, J. R. Krenn, and J.-P. Goudonnet, “Plasmon polaritons of metallic nanowires for controlling submicron propagation of light,” Phys. Rev. B 60, 9061–9068 (1999). [CrossRef]
248. P. Berini, “Plasmon-polariton waves guided by thin lossy metal films of finite width: bound modes of symmetric structures,” Phys. Rev. B 61, 10484–10503 (2000). [CrossRef]
249. R. Charbonneau, P. Berini, E. Berolo, and E. Lisicka-Skrzek, “Experimental observation of plasmon-polariton waves supported by a thin metal film of finite width,” Opt. Lett. 25, 844–846 (2000). [CrossRef]
251. P. Berini, “Plasmon-polariton waves guided by thin lossy metal films of finite width: bound modes of asymmetric structures,” Phys. Rev. B 63, 125417 (2001). [CrossRef]
252. B. Lamprecht, J. R. Krenn, G. Schider, H. Ditlbacher, M. Salerno, N. Felidj, A. Leitner, and F. R. Aussenegg, “Surface plasmon propagation in microscale metal stripes,” Appl. Phys. Lett. 79, 51–53 (2001). [CrossRef]
253. J.-C. Weeber, J. R. Krenn, A. Dereux, B. Lamprecht, Y. Lacroute, and J.-P. Goudonnet, “Near-field observation of surface plasmon polariton propagation on thin metal stripes,” Phys. Rev. B 64, 045411 (2001). [CrossRef]
254. P. Berini, “Optical waveguide structures,” U.S. patent 6,741,782 (May 25, 2004).
255. J.-C. Weeber, Y. Lacroute, and A. Dereux, “Optical near-field distributions of surface plasmon waveguide modes,” Phys. Rev. B 68, 115401 (2003). [CrossRef]
256. R. Nikolajsen, K. Leosson, I. Salakhutdinov, and S. I. Bozhevolnyi, “Polymer-based surface-plasmon-polariton stripe waveguides at telecommunication wavelengths,” Appl. Phys. Lett. 82, 668–670 (2003). [CrossRef]
257. J. R. Krenn and J.-C. Weeber, “Surface plasmon polaritons in metal stripes and wires,” Philos. Trans. R. Soc. London, Ser. A 362, 739–756 (2004). [CrossRef]
258. S. J. Al-Bader, “Optical transmission on metallic wires—fundamental modes,” IEEE J. Quantum Electron. 40, 325–329 (2004). [CrossRef]
260. R. Zia, M. D. Selker, and M. L. Brongersma, “Leaky and bound modes of surface plasmon waveguides,” Phys. Rev. B 71, 165431 (2005). [CrossRef]
261. P. Berini, R. Charbonneau, N. Lahoud, and G. Mattiussi, “Characterisation of long-range surface plasmon-polariton waveguides,” J. Appl. Phys. 98, 043109 (2005). [CrossRef]
262. K. Leosson, T. Nikolajsen, A. Boltasseva, and S. I. Bozhevolnyi, “Long-range surface plasmon polariton nanowire waveguides for device applications,” Opt. Express 14, 314–319 (2006). [CrossRef] [PubMed]
265. R. Zia, J. A. Schuller, and M. L. Brongersma, “Near-field characterization of guided polariton propagation and cutoff in surface plasmon waveguides,” Phys. Rev. B 74, 165415 (2006). [CrossRef]
267. Y. Rao, F. Liu, Y.-D. Huang, O. Dai, W. Zhang, and J.-D. Peng, “Long-range surface plasmon polaritons guided by a thin metal stripe,” Chin. Phys. Lett. 24, 1626–1628 (2007). [CrossRef]
268. Y. Satuby and M. Orenstein, “Surface plasmon polariton waveguiding: from multimode stripe to a slot geometry,” Appl. Phys. Lett. 90, 251104 (2007). [CrossRef]
269. J. Jung, T. Søndergaard, and S. I. Bozhevolnyi, “Theoretical analysis of surface plasmon-polariton waveguides for long-range polarization-independent waveguiding,” Phys. Rev. B 76, 035434 (2007). [CrossRef]
270. P. Berini, “Air gaps in metal stripe waveguides supporting long-range surface plasmon-polaritons,” J. Appl. Phys. 102, 033112 (2007). [CrossRef]
272. P. Berini, “Long-range surface plasmon-polariton waveguides in silica,” J. Appl. Phys. 102, 053105 (2007). [CrossRef]
273. J. T. Kim, S. Park, J. J. Ju, S. K. Park, M.-S. Kim, and M.-H. Lee, “Low-loss polymer-based long-range surface plasmon-polariton waveguide,” IEEE Photon. Technol. Lett. 19, 1374–1376 (2007). [CrossRef]
275. I. Salakhutdinov, J. S. Thakur, and K. Leosson, “Characterization of long-range surface plasmon-polariton in stripe waveguides using scanning near-field optical microscopy,” J. Appl. Phys. 102, 123110 (2007). [CrossRef]
276. S. Park, M.-S. Kim, J. T. Kim, S. K. Park, J. J. Ju, and M.-H. Lee, “Long range surface plasmon polariton waveguides at 1.31 and wavelengths,” Opt. Commun. 281, 2057–2061 (2008). [CrossRef]
278. J. Jiang, C. L. Callender, S. Jacob, J. P. Noad, S. Chen, J. Ballato, and D. W. Smith Jr., “Long-range surface plasmon polariton waveguides embedded in fluorinated polymer,” Appl. Opt. 47, 3892–3900 (2008). [CrossRef] [PubMed]
279. S. Chattopadhyay and P. K. Saha, “Study of bound surface plasmon modes on subwavelength structures of rectangular cross section using a simple analytical technique,” J. Opt. A, Pure Appl. Opt. 10, 095201 (2008). [CrossRef]
280. S. Park, J. J. Ju, J. T. Kim, M.-S. Kim, S. K. Park, J.-M. Lee, W.-J. Lee, and M.-H. Lee, “Sub-dB/cm propagation loss in silver stripe waveguides,” Opt. Express 17, 697–702 (2009). [CrossRef] [PubMed]
281. R. Adato and J. Guo, “Modification of dispersion, localization, and attenuation of thin metal stripe symmetric surface plasmon-polariton modes by thin dielectric layers,” J. Appl. Phys. 105, 034306 (2009). [CrossRef]
282. P. Berini and R. Buckley, “On the convergence and accuracy of numerical mode computations of surface plasmon waveguides,” J. Comput. Theor. Nanosci. 6, 2040–2053 (2009). [CrossRef]
283. R. Charbonneau, “Demonstration of a passive integrated optics technology based on plasmons,” M.A.Sc. thesis (University of Ottawa, 2001).
284. J. Lu, “Modelling optical waveguide bends and applications to plasmon-polariton waveguides,” M.A.Sc. thesis (University of Ottawa, 2003).
285. A. Boltasseva, “Integrated-optics components utilizing long-range surface plasmon polaritons” Ph.D. thesis (Technical University of Denmark, 2004).
286. A. Boltasseva, T. Nikolajsen, K. Leosson, K. Kjaer, M. S. Larsen, and S. I. Bozhevolnyi, “Integrated optical components utilizing long-range surface plasmon polaritons,” J. Lightwave Technol. 23, 413–422 (2005). [CrossRef]
287. R. Charbonneau, N. Lahoud, G. Mattiussi, and P. Berini, “Demonstration of integrated optics elements based on long-ranging surface plasmon polaritons,” Opt. Express 13, 977–984 (2005). [CrossRef] [PubMed]
288. R. Charbonneau, C. Scales, I. Breukelaar, S. Fafard, N. Lahoud, G. Mattiussi, and P. Berini, “Passive integrated optics elements based on long-ranging surface plasmon polaritons,” J. Lightwave Technol. 24, 477–494 (2006). [CrossRef]
289. H. S. Won, K. C. Kim, S. H. Song, C.-H. Oh, P. S. Kim, S. Park, and S. I. Kim, “Vertical coupling of long-range surface plasmon polaritons,” Appl. Phys. Lett. 88, 011110 (2006). [CrossRef]
291. A. Boltasseva and S. I. Bozhevolnyi, “Directional couplers using long-range surface plasmon polariton waveguides,” IEEE J. Sel. Top. Quantum Electron. 12, 1233–1241 (2006). [CrossRef]
292. W.-K. Kim, W.-S. Yang, H.-M. Lee, H.-Y. Lee, M.-H. Lee, and W.-J. Jung, “Leaky modes of curved long-range surface plasmon-polariton waveguide,” Opt. Express 14, 13043–13049 (2006). [CrossRef] [PubMed]
293. F. Liu, Y. Rao, Y. Huang, W. Zhang, and J. Peng, “Coupling between long range surface plasmon polariton mode and dielectric waveguide mode,” Appl. Phys. Lett. 90, 141101 (2007). [CrossRef]
294. F. Liu, Y. Rao, X. Tang, R. Wan, Y. Huang, W. Zhang, and J. Peng, “Hybrid three-arm coupler with long range surface plasmon polariton and dielectric waveguides,” Appl. Phys. Lett. 90, 241120 (2007). [CrossRef]
295. A. Degiron, C. Dellagiacoma, J. C. McIlhargey, G. Shvets, O. J. F. Martin, and D. R. Smith, “Simulations of hybrid long-range plasmon modes with application to 90° bends,” Opt. Lett. 32, 2354–2356 (2007). [CrossRef] [PubMed]
296. A. Degiron, S.-Y. Cho, C. Harrison, N. M. Jokerst, C. Dellagiacoma, O. J. F. Martin, and D. R. Smith, “Experimental comparison between conventional and hybrid long-range surface plasmon waveguide bends,” Phys. Rev. A 77, 021804(R) (2008). [CrossRef]
297. Y. H. Joo, M. J. Jung, J. Yoon, S. H. Song, H. S. Won, S. Park, and J. J. Ju, “Long-range surface plasmon polaritons on asymmetric double-electrode structures,” Appl. Phys. Lett. 92, 161103 (2008). [CrossRef]
298. F. Liu, R. Wan, Y. Rao, Y. Zheng, Y. Huang, W. Zhang, and J. Peng, “Hybrid three-arm coupler consisted of long range surface plasmon polariton and dielectric waveguides,” J. Lightwave Technol. 26, 3872–3882 (2008). [CrossRef]
299. A. Degiron, S.-Y. Cho, T. Tyler, N. M. Jokerst, and D. R. Smith, “Directional coupling between dielectric and long-range plasmon waveguides,” New J. Phys. 11, 015002 (2009). [CrossRef]
302. R. Buckley and P. Berini, “Radiation suppressing metallo-dielectric optical waveguides,” J. Lightwave Technol. 27, 2800–2808 (2009). [CrossRef]
303. J. J. Ju, S. Park, M.-S. Kim, J. T. Kim, S. K. Park, Y. J. Park, and M.-H. Lee, “ light signal transmission in long-range surface plasmon waveguides,” Appl. Phys. Lett. 91, 171117 (2007). [CrossRef]
304. J. J. Ju, M.-S. Kim, S. Park, J. T. Kim, S. K. Park, and M.-H. Lee, “ optical signal transmission via long-range surface plasmon polariton waveguide,” ETRI J. 29, 808–810 (2007). [CrossRef]
305. J. J. Ju, S. Park, M.-S. Kim, J. T. Kim, S. K. Park, Y. J. Park, and M.-H. Lee, “Polymer-based long-range surface plasmon polariton waveguides for optical signal transmission applications,” J. Lightwave Technol. 26, 1510–1518 (2008). [CrossRef]
306. J. T. Kim, J. J. Ju, S. Park, M.-S. Kim, S. K. Park, and M.-H. Lee, “Chip-to-chip optical interconnect using gold long-range surface plasmon polariton waveguides,” Opt. Express 16, 13133–13138 (2008). [CrossRef] [PubMed]
307. J.-M. Lee, S. Park, M.-S. Kim, S. K. Park, J. T. Kim, J.-S. Choe, W.-J. Lee, M.-H. Lee, and J. J. Ju, “Low bending loss metal waveguide embedded in a free-standing multilayered polymer film,” Opt. Express 17, 228–234 (2009). [CrossRef] [PubMed]
308. J. T. Kim, J. J. Ju, S. Park, S. K. Park, M.-S. Kim, J.-M. Lee, J.-S. Choe, M.-H. Lee, and S.-Y. Shin, “Silver stripe optical waveguide for chip-to-chip optical interconnections,” IEEE Photon. Technol. Lett. 21, 902–904 (2009). [CrossRef]
309. S. Jetté, “A study of Bragg gratings based on plasmon-polariton wave guides,” M.A.Sc. thesis (University of Ottawa, 2003).
310. A. Boltasseva, S. I. Bozhevolnyi, T. Søndergaard, T. Nikolajsen, and K. Leosson, “Compact Z-add-drop wavelength filters for long-range surface plasmon polaritons,” Opt. Express 13, 4237–4243 (2005). [CrossRef] [PubMed]
311. S. Jetté-Charbonneau, R. Charbonneau, N. Lahoud, G. Mattiussi, and P. Berini, “Demonstration of Bragg gratings based on long-ranging surface plasmon polariton waveguides,” Opt. Express 13, 4674–4682 (2005). [CrossRef] [PubMed]
312. S. I. Bozhevolnyi, A. Boltasseva, T. Søndergaard, T. Nikolajsen, and K. Leosson, “Photonic bandgap structures for long-range surface plasmon polaritons” Opt. Commun. 250, 328–333 (2005). [CrossRef]
313. A. Boltasseva, T. Søndergaard, T. Nikolajsen, K. Leosson, S. I. Bozhevolnyi, and J. M. Hvam, “Propagation of long-range surface plasmon polaritons in photonic crystals” J. Opt. Soc. Am. B 22, 2027–2037 (2005). [CrossRef]
314. T. Søndergaard and S. I. Bozhevolnyi, “Out-of-plane scattering properties of long-range surface plasmon polariton gratings,” Phys. Status Solidi B 242, 3064–3069 (2005). [CrossRef]
315. S. Jetté-Charbonneau, R. Charbonneau, N. Lahoud, G. Mattiussi, and P. Berini, “Bragg gratings based on long-range surface plasmon-polariton waveguides: comparison of theory and experiment,” IEEE J. Quantum Electron. 41, 1480–1491 (2005). [CrossRef]
316. T. Søndergaard, S. I. Bozhevolnyi, and A. Boltasseva, “Theoretical analysis of ridge gratings for long-range-surface plasmon polaritons,” Phys. Rev. B 73, 045320 (2006). [CrossRef]
317. A. Boltasseva, S. I. Bozhevolnyi, T. Nikolajsen, and K. Leosson, “Compact Bragg gratings for long-range surface plasmon polaritons,” J. Lightwave Technol. 24, 912–918 (2006). [CrossRef]
318. S. Jetté-Charbonneau and P. Berini, “Theoretical performance of Bragg gratings based on long-range surface plasmon-polariton waveguides,” J. Opt. Soc. Am. A 23, 1757–1767 (2006). [CrossRef]
319. R. H. Pedersen, A. Boltasseva, D. M. Johansen, T. Nielsen, K. B. Jørgensen, K. Leosson, J. E. Østergaard, and A. Kristensen, “Nanoimprinted reflecting gratings for long-range surface plasmon polaritons,” Microelectron. Eng. 84, 895–898 (2007). [CrossRef]
320. J.-W. Mu and W.-P. Huang, “A low-loss surface plasmonic Bragg grating,” J. Lightwave Technol. 27, 436–439 (2009). [CrossRef]
321. I. G. Breukelaar, “Surface plasmon-polaritons in thin metal strips and slabs: Waveguiding and mode cutoff,” M.A.Sc. thesis (University of Ottawa, 2004).
322. G. Gagnon, “Thermo-optic variable optical attenuators using plasmon-polariton waveguides,” M.A.Sc. thesis (University of Ottawa, 2004).
323. P. Berini, G. Gagnon, and S. Jetté-Charbonneau, “Thermo-optic plasmon-polariton devices,” U.S. Patent 7,043,134 (May 9, 2006).
324. T. Nikolajsen, K. Leosson, and S. I. Bozhevolnyi, “Surface plasmon polariton based modulators and switches operating at telecom wavelengths,” Appl. Phys. Lett. 85, 5833–5835 (2004). [CrossRef]
325. T. Nikolajsen, K. Leosson, and S. I. Bozhevolnyi, “In-line extinction modulator based on long-range surface plasmon polaritons,” Opt. Commun. 244, 455–459 (2005). [CrossRef]
326. S. I. Bozhevolnyi, T. Nikolajsen, and K. Leosson, “Integrated power monitor for long-range surface plasmon polaritons,” Opt. Commun. 255, 51–56 (2005). [CrossRef]
327. I. Breukelaar, R. Charbonneau, and P. Berini, “Long range surface plasmon-polariton mode cutoff and radiation,” Appl. Phys. Lett. 88, 051119 (2006). [CrossRef]
328. S. Park and S. H. Song, “Polymeric variable optical attenuator based on long range surface plasmon polaritons,” Electron. Lett. 42, 402–404 (2006). [CrossRef]
329. I. Breukelaar, R. Charbonneau, and P. Berini, “Long-range surface plasmon-polariton mode cutoff and radiation in embedded strip waveguides,” J. Appl. Phys. 100, 043104 (2006). [CrossRef]
330. G. Gagnon, N. Lahoud, G. Mattiussi, and P. Berini, “Thermally activated variable attenuation of long-range surface plasmon-polariton waves,” J. Lightwave Technol. 24, 4391–4402 (2006). [CrossRef]
332. P. Berini, G. Mattiussi, N. Lahoud, and R. Charbonneau, “Wafer-bonded surface plasmon waveguides,” Appl. Phys. Lett. 90, 061108 (2007). [CrossRef]
333. S.-K. Kim, Y.-C. Hung, W. Yuan, D.-H. Cha, K. Geary, H. R. Fetterman, S. M. Ha, Q. Pei, J. Luo, X.-H. Zhou, A. K.-Y. Jen, M.-S. Kwon, and W. H. Steier, “Metal-slotted polymer optical waveguide device,” Appl. Phys. Lett. 90, 243507 (2007). [CrossRef]
334. P. Berini, R. Charbonneau, S. Jetté-Charbonneau, N. Lahoud, and G. Mattiussi, “Long-range surface plasmon-polariton waveguides and devices in lithium niobate,” J. Appl. Phys. 101, 113114 (2007). [CrossRef]
335. G. Mattiussi, N. Lahoud, R. Charbonneau, and P. Berini, “Fabrication of long-range surface plasmon-polariton waveguides in lithium niobate on silicon,” J. Vac. Sci. Technol. A 25, 692–700 (2007). [CrossRef]
336. S. Jetté-Charbonneau, N. Lahoud, R. Charbonneau, G. Mattiussi, and P. Berini, “End-facet polishing of surface plasmon waveguides in lithium niobate,” IEEE Trans. Adv. Packag. 31, 479–483 (2008). [CrossRef]
339. S. Jetté-Charbonneau and P. Berini, “External cavity laser using a long-range surface plasmon grating as a distributed Bragg reflector,” Appl. Phys. Lett. 91, 181114 (2007). [CrossRef]
341. M. Tencer, R. Charbonneau, and P. Berini, “Confinement and manipulation of droplets on hydrophilic surfaces using a flat high surface energy guide,” Lab Chip 7, 483–489 (2007). [CrossRef] [PubMed]
343. M. Tencer, R. Charbonneau, N. Lahoud, and P. Berini, “AFM study of BSA adlayers on Au stripes,” Appl. Surf. Sci. 253, 9209–9214 (2007). [CrossRef]
344. R. Charbonneau, E. Lisicka-Shrzek, and P. Berini, “Broadside coupling to long-range surface plasmons using an angle-cleaved optical fibre,” Appl. Phys. Lett. 92, 101102 (2008). [CrossRef]
345. R. Charbonneau, M. Tencer, N. Lahoud, and P. Berini, “Demonstration of surface sensing using long-range surface plasmon waveguides on silica,” Sens. Actuators B 134, 455–461 (2008). [CrossRef]
346. R. Charbonneau and P. Berini, “Broadside coupling to long-range surface plasmons in metal stripes using prisms, particles and an atomic force microscope probe,” Rev. Sci. Instrum. 79, 073106 (2008). [CrossRef]
347. P. Berini, “Bulk and surface sensitivities of surface plasmon waveguides,” New J. Phys. 10, 105010 (2008). [CrossRef]
348. P. Berini, N. Lahoud, and R. Charbonneau, “Fabrication of surface plasmon waveguides and integrated components on ultra-thin free-standing membranes,” J. Vac. Sci. Technol. A 26, 1383–1391 (2008). [CrossRef]
350. P. Berini, R. Charbonneau, and N. Lahoud, “Long-range surface plasmons along membrane-supported metal stripes,” IEEE J. Sel. Top. Quantum Electron. 14, 1479–1479 (2008). [CrossRef]
351. M. Tencer and P. Berini, “Design of microfluidic channels separated by an ultra-thin free standing dielectric membrane,” Microfluid. Nanofluid. 6, 17–26 (2009). [CrossRef]
352. R. Daviau, E. Lisicka-Skrzek, R. N. Tait, and P. Berini, “Broadside excitation of surface plasmon waveguides on Cytop,” Appl. Phys. Lett. 94, 091114 (2009). [CrossRef]
353. N. Fong, P. Berini, and R. N. Tait, “Fabrication of surface plasmon waveguides on thin Cytop membranes,” J. Vac. Sci. Technol. A 27, 614–619 (2009). [CrossRef]
354. N. Fong, P. Berini, and R. N. Tait, “Mechanical properties of thin free-standing Cytop membranes,” submitted to J. Microelectromech. Syst..
355. R. Daviau, A. Khan, E. Lisicka-Skrzek, R. N. Tait, and P. Berini, “Fabrication of surface plasmon waveguides and integrated components on Cytop,” submitted to Microelectron. Eng..
357. J. C. Love, L. A. Estroff, J. K. Kriebel, R. G. Nuzzo, and G. M. Whitesides, “Self-assembled monolayers of thiolates on metals as a form of nanotechnology,” Chem. Rev. (Washington, D.C.) 105, 1103–1169 (2005). [CrossRef]
358. W. Lukosz, “Principles and sensitivities of integrated optical and surface plasmon sensors for direct affinity sensing and immunosensing,” Biosens. Bioelectron. 6, 215–225 (1991). [CrossRef]
359. J. A. Conway, S. Sahni, and T. Szkopek, “Plasmonic interconnects versus conventional interconnects: a comparison of latency, cross-talk and energy costs,” Opt. Express 15, 4474–4484 (2007). [CrossRef] [PubMed]
361. T. Okamoto, J. Simonen, and S. Kawata, “Plasmonic crystal for efficient energy transfer from fluorescent molecules to long-range surface plasmons,” Opt. Express 17, 8294–8301 (2009). [CrossRef] [PubMed]
362. D. Ballester, M. S. Tame, C. Lee, J. Lee, and M. S. Kim, “Long-range surface-plasmon-polariton excitation at the quantum level,” Phys. Rev. A 79, 053845 (2009). [CrossRef]
364. Y. H. Joo, S. H. Song, and R. Magnusson, “Long-range surface plasmon-polariton waveguide sensors with a Bragg grating in the asymmetric double-electrode structure,” Opt. Express 17, 10606–10611 (2009). [CrossRef] [PubMed]
365. A. Huck, S. Smolka, P. Lodahl, A. S. Sørensen, A. Boltasseva, J. Janousek, and U. L. Andersen, “Demonstration of quadrature-squeezed surface plasmons in a gold waveguide,” Phys. Rev. Lett. 102, 246802 (2009). [CrossRef] [PubMed]
366. V. N. Konopsky, D. V. Basmanov, E. V. Alieva, D. I. Dolgy, E. D. Olshansky, S. K. Sekatskii, and G. Dietler, “Registration of long-range surface plasmon resonance by angle-scanning feedback and its implementation for optical hydrogen sensing,” New J. Phys. 11, 063049 (2009). [CrossRef]
367. C. Scales, I. Breukelaar, and P. Berini, “Surface-plasmon Schottky contact detector based on a symmetric metal stripe in silicon,” submitted to Opt. Lett..
368. I. De Leon and P. Berini, “Modeling surface plasmon-polariton gain in planar metallic structures,” submitted to Opt. Express.