Longwave plasmonics on doped silicon and silicides

The realization of plasmo-electronic integrated circuits in a silicon chip will be enabled by two new plasmonic materials that are proposed and modeled in this article. The first is ion-implanted Si (n-type or p-type) at the surface of an intrinsic Si chip. The second is a thin-layer silicide such as Pd2Si, NiSi, PtSi2 WSi2 or CoSi2 formed at the Si chip surface. For doping concentrations of 10 cm and 10 cm, our dispersion calculations show that bound surface plasmon polaritons will propagate with low loss on stripe-shaped plasmonic waveguides over the 10 to 55 μm and 2.8 to 15 μm wavelength ranges, respectively. For Pd2Si/Si plasmonic waveguides, the wavelength range of 0.5 to 7.5 μm is useful and here the propagation lengths are 1 to 2300 μm. For both doped and silicided guides, the SPP mode field extends much more into the air above the stripe than it does into the conductive stripe material. © 2008 Optical Society of America OCIS codes: (240.6680) Surface plasmons; (250.5403) Plasmonics; (310.2790) Guided waves; (160.3918) Metamaterials; (160.3130) Integrated optics materials. References and links 1. M. V. Bashevoy, F. Jonsson, A. V. Krasavin, and N. I. Zheludev, Y. Chen, M. I. Stockman, “Generation of traveling surface plasmon waves by free-electron impact,” Nano Lett. 6, 1113 (2006). 2. J. T. van Wijngaarden, E. Verhagen, A. Polman, C. E. Ross, H. J. Lezec, and H. A. Atwater, “Direct imaging of propagation and damping of near-resonance surface plasmon polaritons using cathodoluminescence spectroscopy,” Appl. Phys. Lett. 88, 221111 (2006). 3. S. Fan, R. A. Höpfel, E. Vass, and E. Gornik, “Thermal excitation of two dimensional plasma oscillations,” Phys. Rev. Lett. 49, 1667 (1982). 4. R. A. Soref, “The opto electronic integrated circuit,” invited lead chapter in Silicon Photonics—State of the Art, G. T. Reed, eds., (John Wiley & Sons, Chichester, UK, to be published 2008). 5 R. A. Soref, “Towards silicon-based longwave integrated optoelectronics (LIO),” paper 6898-5 (invited) SPIE Photonics West, San Jose, CA, 21 January 2008 . 6. S. J. Bozhevolnyi and J. Jung, “Scaling for gap plasmon based waveguides”, Optics Express, 16, 2676-2684 (2008). 7. M. Amiotti, G. Guizzetti, F. Marabelli, A. Piaggi, V. N. Antonov, Vl. N. Antonov, O. Jepson, O. K. Anderson, A. Borghesi, F. Nava, V. V. Nemoshkalenko, R. Madar and A. Roualt, “Optical properties of Pd2Si,” Phys. Rev. B 45, 13285 (1992). 8. H. Raether, “Surface plasmon oscillations and their applications,” in Physics of Thin Films, G. Hass, M. Francombe, and R. W. Hoffman, eds., (Academic Press, NY, 1977). 9. E. D. Palik, ed., Handbook of Optical Constants of Solids, (Academic Press, 1985) Vol. 1.


Introduction
Plasmonics refers to the technology in which surface plasmons are launched, transmitted, modulated, detected, etc, on sub-wavelength conducting structures formed on or within dielectrics or semiconductors.Plasmon optics is a subset of this technology in which the surface plasma polaritons (SPPs) are generated and detected optically.However, the SPPs can also be launched "electronically" by an electron beam [1,2] or by a MOSFET [3].The latter case is most interesting because, as we described earlier [4], it will lead to silicon-based plasmo-electronic integrated circuits (PEICs) in which MOSFETs on the perimeter of a planar "plasmonic network" serve to launch and sense SPPs that are utilized by the plasmonic signalprocessing circuits.This will be a CMOS-compatible, manufacturable IC.
A plasmon-optic circuit is also quite valuable because its sub-wavelength features facilitate the new technology of nano-photonics.In fact, plasmon optics (alongside photoniccrystal technology) is a key contributor to a "Moore's Law for photonics" which indicates that the on-chip packing density of integrated photonic components will increase year-by-year thanks to photonic research and development.In the long term, we can expect that plasmon optics will complement the plasmo-electronics discussed above; in other words, we shall eventually witness plasmo-optoelectronic integrated circuits (POEICs), a seamless integration of three techniques.
The prior art of plasmon optics utilized metal/dielectric and metal/semiconductor composites with their wavelengths generally in the visible or near infrared.In this paper, we focus upon non-metalized silicon substrates.We propose and analyze two new plasmonic materials that are probably more compatible with optoelectronics than are silver or gold films.These materials are highly doped-silicon embedded in undoped silicon, and a silicide film (e.g.PtSi 2 , NiSi, Pd 2 Si, WSi 2 , CoSi 2 ) created in the undoped silicon surface.Both materials are "naturals" for fabrication in a silicon CMOS foundry.The doped-Si and silicided-Si SPP circuits interface seamlessly with their FET partners.Patterned doping can be accomplished in a well-known standard way by "shallow" implantation of ions into Si, or by in-diffusion through a mask.Formation of silicides on Si by metal deposition, thermal reaction, and conversion-to-silicide is also well known in the industry.
The wavelengths for prior-art metal-delectric-metal SPP based waveguides are in the near infrared, but as we shall show in this article, the doped-Si and silicided-Si offer low-loss performance not in the NIR, but rather in the mid-wave, longwave and very longwave infrared.Our generic term "longwave" is shorthand for these three spectral regions [5].The longwave SPP circuit or "network" will perform signal processing functions on the Si chip.If the footprint of the network is sufficiently small, if the processing is fast, and if the power consumption is quite low, then the wavelength and frequency of the plasmons are immaterial.We anticipate that the longwave SPP circuit will meet the size-speed-power requirements, making the longwaves "viable".
The organization of this paper is as follows: we first propose striped and slotted structures for channel waveguides, then we compile the relevant surface plasmon properties calculated from the complex permittivity.These properties include the dispersion curves with their bound and free modes, the propagation losses within feasible longwave regions, and the mode profiles with estimated penetration depths into the air and the conductive medium.

Waveguide structures
Two structures that are beneficial for plasmonic circuits are the simple stripe and the parallelstripe "slotted" plasmonic waveguides [6].These are illustrated in Fig. 1 for the cases of doped Si and silicide-on-Si.The vertically stacked double-plasmonic waveguide (not shown
here) will also work quite well.It consists of two plane-parallel conductive stripes (doped or silicided) that are separated by a dielectric layer.We begin our analysis of waveguiding by examining the permittivity data for p-type Si, n-type Si, Pd 2 Si and a prior-art reference material (gold).

Complex permittivity data
The SPP properties that can be calculated from the permittivity include the complex wavevector, the propagation length L for the SPP intensity along the surface, the SPP spatial confinement as determined by the distance L air that the field extends into air above the waveguide surface and by the field penetration L cond into the doped or silicided conductor.As a concrete example of a silicide where at least some empirical information exists [7] for the complex permittivity ε(ω) we choose Pd 2 Si.Doped Si supposes P or B dopants at levels of 10 20 and 10 21 cm -3 .Empirical ε(ω) data for either materials system is presently incomplete in the longwave regime, so we use instead theoretical extrapolation based on the Drude model.
To determine the longwave permittivity from the available near-IR measurements, the Drude model gives where ε ∞ is the permittivity for ω >> ω p , the plasma frequency is ω p , and τ is the electron relaxation time.For doped Si, τ is determined from the empirical drift mobility μ = eτ/m*, where m* is the orientationally averaged electron effective mass m* = 0.272 m 0 , ε ∞ = 11.7, and ω p 2 = Ne 2 /ε 0 m*, with N the electron concentration per cm 3 .The mobility at the two concentrations considered here is ~50 cm 2 /V-s.For Pd 2 Si, and an orientationally averaged plasma frequency of 3.1 eV, a relaxation time of 1/τ = 0.03 eV [7], and ε ∞ = 1 were used.SPPs are supported at frequencies wherever ε′ is negative.For Pd 2 Si, unlike gold, the probability of SPP generation by electron beams, given by Im[-1/(1-ε)], extends deep into the infrared [8].Longwave SPPs are also generated efficiently in highly doped Si.Little is known about the infrared permittivity of the novel Si-compatible narrow-gap single-crystal GeSn films, but their intrinsic conductivity approaches 10 20 cm -3 with Drude model parameters expected to be comparable to Ge.Thus this alloy also has SPP promise.
Figures 2(a) and 2(b) present the real and imaginary parts of the complex permittivity of our four materials.The Fig. 2(a) plot of Re[ε(ω)], or ε′(ω) in Eq. ( 1), shows that SPPs are supported on Si for frequencies ω corresponding to optical wavelengths longer than 7 μm for 10 20 cm -3 doping and longer than 2 μm for 10 21 cm -3 doping.For Pd 2 Si and gold, the high energy curves are empirical values [7,9] and the low energy sections are our calculated Drude-model values.For the silicide, SPPs can be sustained throughout the visible and infrared.Figure 2(b) presents Im[ε(ω)], or ε″(ω) in Eq.For the silicide, ε″ is significantly larger than for gold in the visible and NIR, suggesting higher loss, but the small Pd 2 Si damping parameter γ = 1/τ of only 0.03 eV [7] results in a calculated Drude curve that is significantly lower than for gold in the mid-and far-infrared.The large discontinuity between empirical-and Drude-based permittivity ε″ values emphasizes the lack of experimental infrared data for the silicide.

Bound SPPs for plasmonic waveguides
We next calculate the ω-k dispersion curves for our materials; the dependence of the SPP energy as a function of the real SPP wave vector k′.The complex SPP wavevector for the vacuum-"metal" interface is determined from the complex permittivity ε of the metal, according to where ω is the angular frequency, and c is the speed of light.The results are presented in Fig. 3(a) showing the SPP dispersion curve ω vs the real wavevector k′ for doped silicon.The curves for both impurity concentrations fall below the light line ω = c k at low wavevector, but rather than saturating at ω p /√2 as for ideal free-electron metals, there is an inflection point where the curve jumps to the upper polariton branch.The above-lightline branch refers to radiative SPPs which are not really useful on our chip, whereas the below-lightline branches describe the bound SPPs desired in our waveguided circuit.Judging from Fig. 3(a), the lower limit of the useful wavelength range ought to be where the SPP branch crosses over the lightline.However, the actual lower wavelength limits will be The usual distinct polariton branches appear if τ is artificially increased by about two orders of magnitude such that ωτ becomes >> 1.Unless that condition is fulfilled, the imaginary part of the permittivity is significant [Eq.( 1)], suggesting high loss.For the doping concentrations considered, ωτ is only of order unity in the infrared, therefore the Fig. 2(a) S-shaped SPP dispersion is a feature of a fairly small relaxation time (fairly large dissipation).Decreasing ε ∞ similarly makes the curve more polariton-like.Figure 3(b) compares the dispersion curves for Pd 2 Si and gold using empirical permittivity data [7,9].A slight S-shaped kink in the SPP dispersion curve exists for gold at ~500 nm optical wavelength, similar to the infrared situation for doped semiconductors, because of dstate absorption loss in the visible.The dispersion curve for Pd 2 Si generally hugs the light line throughout the spectrum, although this curve is actually "microscopically below" the light line, allowing bound SPP propagation over a wide energy range.Experimentally, we have obtained evidence of SPPs excited on Pd 2 Si at long waves using a grating coupler and a CO 2 laser, results that will be published elsewhere

Mode propagation loss
The characteristic propagation length L for SPP intensity is L -1 = 2 Im[k] and we turn to Eq. ( 2) along with Fig. 2 to determine L as a function of SPP energy hν. Figure 4 presents and compares our predicted SPP propagation lengths for the four plasmonic materials.The SPP wavelengths 2π/k′ are nearly the same as the optical wavelengths λ at the two SPP frequencies since the dispersion curve is close to the light line.At the SPP short-wavelength limits of 8 μm and 2.5 μm, Fig. 4 reveals that L is comparable to λ, a critically damped situation where one can hardly speak of a group velocity or the propagation of information.Therefore we shall demand that L is twice λ to make the lower limit practical.Doing that, we then find from Fig. 4 that the useful minimum wavelengths are 3 μm for 10 21 cm -3 and 9 μm for 10 20 cm -3 .
For the silicide, propagation lengths are always less than those for gold.Nevertheless they are still very long in the long-wave IR and THz.Applying the L = 2λ low criterion shows that silicide is suitable for propagating information well into the visible.Also, L exceeds λ by factors of order 1000 in the longwave infrared, showing that SPP's in silicide are underdamped in that region.Our L calculations with L = 2λ low predict minimum useful freespace wavelengths of 0.5, 0.5, 2.8 and 10 μm, respectively, for gold, Pd 2 Si, doped Si at 10 21 , and doped Si at 10 20 .The calculated ~1 cm propagation length for gold at 0.1 eV has been experimentally verified by us, as will be described elsewhere.

Mode profile
Assuming an air cladding above the stripe-shaped plasmonic silicon waveguide, we next ask how the SPP field fringes out into the air above and into the conductive material below the Si surface.The penetration length of the Z-propagating wave is defined as the lateral dimension X at which the field strength decreases to 1/e of its initial value.The SPP fields extend above the surface by an amount L air and penetrate into the conductor by an amount L cond with L air = (c/ω)/Re[√(−1/(1 + ε))] and L cond = (c/ω)/Re[√(−ε 2 /(1 + ε))] where ε is given in Fig. 2. For SPPs propagating in the Z direction, Fig. 5 presents results of our X-penetration calculations for fields into the media on either side of the surface.The Fig. 5 inset shows these penetrations schematically.Generally, the mode fringing into air is much larger than the conductor penetration.The penetration depth into the doped silicon is < 0.6 μm.Uniform doping of the silicon surface to depths of 1 μm is readily achieved by diffusion or ionimplantation.Thus, surface-doped silicon can behave as bulk material as far as SPPs are concerned.This is a very important observation for making SPP waveguides by patterned diffusion or implantation on future PEIC devices.For silicide, the mode penetration into air grows very large at the longest wavelengths.By contrast, the mode penetration into Pd 2 Si is < 70 nm at all wavelengths, which means that a very thin silicide layer, such as 100 nm, will suffice as an SPP host, and that is advantageous for devices.
different criterion based upon the propagation lengths discussed below.