Nano plasmon polariton modes of a wedge cross section metal waveguide

Optical plasmon-polariton modes confined in both transverse dimensions to significantly less than a wavelength are exhibited in open waveguides structured as sharp metal wedges. The analysis reveals two distinctive modes corresponding to a localized mode on the wedge point and surface mode propagation on the abruptly bent interface. These predictions are accompanied by unique field distributions and dispersion characteristics.


Introduction
Sub wavelength confinement of guided waves is highly advantageous for downscaling optical components to the nano-scale, such as for inter-chip optical wiring and nano-systems interconnects. Here we present modal analysis of a metallic waveguide with a cross-section shaped as a narrow wedge, exhibiting modal size hundred times smaller than the wavelength in each transverse dimension.
Surface Plasmon-Polariton (SPP) is a guided wave confined in 1D to a sub-wavelength cross-section, though infinite in the other transverse dimension [1]. Confinement in both the transverse dimensions is achievable using the plasmonic slow waves and specific topologies, analogous to "standard" dielectric waveguides -namely metallic nano-wires, which were discussed elsewhere [2,3]. However the surface guiding of SPP allows also for open contour topologies, relaxing the waveguide transversal dimension constraint.
While a metallic surface supports 1D confined SPP, a modification in the curvature of a metallic surface, e.g., a wedge shaped metallic structure, is expected to support Plasmon-Polariton confined in 2D. The enhanced field of the guided wave in the vicinity of the wedge tip may be employed to enhance matter-field interactions. Here we analyze the wedge waveguide in closed form to understand the role of parameters (e.g. wedge angle and dielectric constants), the dispersion relations and the modal confinement.
Full fledged analytic solutions for optical modes of metallic wedge waveguides were never reported, and the most detailed analysis employed a scalar wave-equation assuming electrostatic regime for sharp [4] and smooth [5,6] topologies. Although simplifying the analysis, this eliminates the mere issue of the full vectorial nature of the hybrid modes and their structural dispersion while discarding the retardation effects. The electrostatic solution obtained in [4] is an asymptotical approximation to the one obtained here.
Recently, localized plasmons propagating on metal wedge were demonstrated using FDTD simulation and in experiments [7]. The analysis presented here provides understanding of the structure of the modes and their parametric dependence, among which proposing the existence of a critical wedge angle similar to the one reported for modes of different symmetry, and should assist in enhanced design of such novel nano-waveguides.
Worth mentioning is a related complementary problem of propagation along smooth channels carved in metal, analyzed using a modified moments method [8]. However, the basis set of the functions used there is not suitable for sharp wedge configuration and thus not appropriate for the problem at hand. Additionally, resolving the field profile in [8] involves a tedious solution of a set of integral equations. In the current work we derive closed form solutions for sharp-wedge cross-section metal waveguides, exhibiting a substantial field concentration near the wedge point and study the interesting dispersion characteristics of the guided modes. This analysis can be trivially enhanced to the complementary structures of vgrooves in metal, recently studied via simulations and experimentally [9,10].

Analysis of wedge plasmonic modes
The analysis of the metal wedge waveguide surrounded by a dielectric is performed in cylindrical coordinates system {r,θ,z} assuming indefinitely long wedge, both in the radial and axial directions, as shown in Fig. 1(a). Though actual wedges are finite, the plasmonic fields are located around the wedge tip, which validate this assumption. Since plasmonic modes are slow waves, they exponentially decay into both air and metal, in the radial direction as well as azimuthally. In cylindrical coordinate system the radial solutions of the wave equation are the modified Bessel functions of the second kind, denoted by K. Satisfying the boundary conditions at each interface point is impossible for a single Modified Bessel function at each medium, thus the solution is a series of K-functions.
The order of the K-functions may be chosen as real or imaginary, i.e. set of basis function of the form ~K v (qr)exp(ivθ) or ~K iv (qr)exp (-vθ). q is the radial momentum, and v denoted the K-function order. Taking the order to be imaginary is more suitable since it implies azimuthal hyperbolic dependence, rather than harmonic, which better describes the azimuthal decay of the plasmonic solutions away from the interface. Moreover, the {K iv (qr)} v is a complete set (in contrary to {K v (qr)} v ) with basis functions being square integrable and having vanishing power at infinity (decaying faster than r -1 ). Using this set for the solution at each media: F(r,θ)=∫dv f v K iv (qr)exp(-vθ), the expansion coefficients, as a function of the order, f v , define the v-spectrum of the solution. The general v-spectrum is continuous on the indefinite interval, [0,∞). However the compactness of this expansion allows for a sufficient representation on a limited interval, [0,v max ]. The criteria for sufficiently high v max were the fulfillment of the boundary conditions to 1% accuracy. As the order of a basis function is increased, the azimuthal decay is enhanced while the radial one is decreased. The confinement of the plasmonic fields around the wedge tip is translated to effectively bounding the v-spectrum. This could be inferred also from the Kontorovich-Lebedev transform [11], having K-functions with imaginary order as basis set.
Truncating the v-spectrum components above v max may cause errors mostly at the interfaces. Thus, constructing and testing the solution, according to the boundary conditions, on the interface assures of the solution accuracy on the entire domain.

Modal characteristics
Applying the analysis described in the previous section, the symmetry permits for two sets of solutions: odd and even about the axis of symmetry (θ=α). Solving the problem in terms of tangential field components (i.e. E z and H z ), four independent boundary conditions apply for each r>0 on the interface θ=0: The ~ denotes the fields at the interfaces, M and D denote the fields in metal and dielectric respectively. f i (i=1...4) coefficients are frequency and β dependent (β -the modal propagation constant), given explicitly in the Appendix. From Eq.
(1) it is apparent that only hybrid modes are possible. Out of the four symmetries we study in this text the solutions having the E(even)-H(odd) symmetry, exited with symmetrical electrical field. The general solution is: Where the radial momentum: q D,M =(β 2 -k 0 2 ε D,M ) 0.5 , k 0 is the free-space wave number, and ε is the relative dielectric constant. For H z o a,b,cosh are replaced by c,d,sinh in Eq. 2, v which denotes the K-function imaginary order, is real and continuous values. Since each of the basic field functions in both sides of the interface has a different r-dependence, the boundary conditions for all r values dictate a superposition of wave solutions. This is a major complication absent in the electrostatic analysis, for which one assumes substantial β-values resulting with similar basis functions on both sides of interface (i.e. q M =q D ).
Approximating the v-spectrum integral as N-sized discrete series, projecting the four boundary conditions onto the base function of the dielectric {K iv (q D r)}, using orthogonality of K iv (qr)/r, and characteristics of the K-function integral, yields 4N algebraic equations: Expressions for G and H functions are given at the Appendix. The solution of this algebraic equation set is self-consistent for a vanishing determinant -resulting in the dispersion relation of the guided modes and the field v-spectrum, from which the field distribution is extracted. Typically the v-spectrum peak value is observed around v=0 and v=3/α for the fields in the dielectric and metal respectively and the spectrum is rapidly decaying with v. Broadest span of v-spectrum is required for the solution on interfaces.
To elucidate the basic plasmonic mode characteristics, the metal losses are not taken into account in the followings. It can be done due to the relative smallness of loss in respect to metal dispersion and for short propagation distance, as is commonly done. We use the lossless Drude model for the metal permittivity, with plasma wavelength of λ P =137nm. The same calculation scheme may support the case of complex permittivity, namely complex ε M and β.
The analysis is exemplified for a 36 degrees gold wedge waveguide (2α=36 0 ) surrounded by air. The dispersion relations are depicted in fig. 1(b), denoted by two bold curves. Both dispersion curves are reaching asymptotically the single surface SPP frequency (ω SPP =ω P /(1+ε D ) 0.5 ) as β tends to infinity -typical for plasmonic waveguides. The single surface SPP dispersion curve is superimposed in Fig. 1(b) (dashed). The upper branch (frequency-wise) exhibits negative group velocity, which vanishes for infinite β values. The lower dispersion branch is divided by a minimum frequency point (ω min =0.69ω p ) into negative and positive group velocity intervals. For a frequency above the cutoff frequency (ω min ), in the interval (ω min <ω<ω SPP ), two propagation constants are simultaneously supported, with opposite group velocity signs. Similar phenomenon for the same symmetry is shown, e.g., for specific set up of a metal slab [12]. The positive group index at high β values is due to the single surface SPP characteristics (power always predominantly in the dielectric), while negative group index at lower β is related to higher power in the metal due to enhanced interface coupling. Similar phenomena can be seen for a metal wire configuration [12] The frequency of the upper branch decreases monotonically with β, accompanied by tighter localization around the wedge point. The E z field component at pt. (I) is depicted in Fig. 2, and drops off symmetrically for all transverse directions -into the metal and dielectric. At the same effective index value, the lower dispersion branch exhibits different field characteristics, as demonstrated at Fig. 3. Comparison of the E z component reveals that a main lobe is still localized on the wedge point but secondary lobes appear as well on the metal interfaces. Decreasing the modal frequency, the number of side lobes as well as the respective modal cross section increase. The distribution of the pointing vector directed along z (S z ) at dispersion points I and II are depicted in Fig 4(a) and 4(b), accordingly. Especially at point II, although the power propagates predominantly in the dielectric, more than 90% of the mode's intensity and power (S z ) are guided within a cross section of 0.01% of free-space wavelengthwell under the "diffraction limit". In Fig. 5 it is apparent that also for the lower branch -the mode becomes more localized towards the wedge point as β is increased.
The distinct characteristics of the two dispersion branches are associated with the 'dual role' that the wedge plays. The one dimensional discontinuity formed at the wedge tip serves as indefinitely thin plasmonic waveguide "core". On the other hand, a metal wedge is comprised of two metal surfaces, each serves as a single surface SPP mode waveguide, which are coupled at their merging point. These two distinct modes of propagation supported by the wedge configuration are manifested by the two dispersion curves. The propagation associated with wedge tip 1D discontinuity is relegated to the upper branch, while the lower curve is more of the two coupled surface modes.   Similarly to slab modes, the field associated with the lower curve (depicted in Fig. 3 for pt. II) has both tangential component (E z ) and component normal to the interface (E θ ). For the normal component the two surface modes coupled at the tip are evident, being located on the interfaces. However, the dominant component (in terms of maximal field amplitude) is actually E r , which is the transversal parallel component. This component, which is mainly located around the tip, is absent for slab waveguide, and emerges from the hybrid configuration. Since the field related to the upper branch is mainly guided by the wedge tip (rather than the wedge faces) its field components are exhibiting reduced azimuthal dependence, as illustrated in Fig. 2  surface SPP. However SPPs localized on the two wedge facets are exceedingly coupled in the vicinity of the wedge point, which splits the dispersion curve according to the mode symmetry. For E z e the dispersion curve is pushed to lower frequency in respect to the single surface SPP curve. This symmetry, of even tangential E-field, is related to the odd potential mode of metal slab (equivalent to odd transversal H-field) [13]. For low β-values the localization notion is weakening, and the lower branch resembles the upper one for extremely low β-values. The two different asymptotic behaviors of the lower branch results in the minimal frequency point on the dispersion curve at β=5k 0 .
Comparing of the upper branch of Fig 1(b) with the electrostatic analysis [4], is exhibiting a qualitative similarity for large β values: modal frequency is above the asymptotic single surface SPP frequency and decreasing monotonically with β. This is related to the mode localized on the 1D discontinuity. Equipped with the dispersion characteristics, sensitivity to structural parameters can be conjectured. As the surrounding dielectric constant increases, the dispersion curves are pulled down, as ω spp decreases. The cutoff frequency should increase with the wedge angle, as the lower dispersion branch may diverge less from the single surface SPP curve, due to reduced coupling between the modes near the wedge point. A cutoff angle is thus expected at each given wavelength, above which a solution does not exist. This resembles the cutoff angle that was observed in some specific FDTD simulations for different symmetry modes (E e H e ) [7] and also for the complementary structure of V-grooved channel polaritons [10]. It is noteworthy that the modes of the symmetry reported here are have higher effective index (i.e., enhanced plasmonic effect) and tighter confinement, than those reported in [7].

Conclusion
We analyzed in detail nano-scale modes propagating on a metal waveguide with a wedge cross-section. Two branches are resolved, one is related to the wedge point localized mode (wire like) and the second to the waves propagating on the interfaces crossing at the wedge point. This notion is used to explain the curve characteristics, as well to validate with previous asymptotic studies and observations. The cross-section of the modes is of a substantial subwavelength dimension, with most of the guided power propagating within a cross section with each dimension ~1% of the wavelength. Including the metal losses, especially for frequencies near ω P , yields rapidly attenuated modes. Gain-assisted propagation [14] may mitigate this difficulty. We acknowledge the Israeli Ministry of Science for support.

APPENDIX
The boundary conditions coefficients: