Abstract
We report simulations of electrically pumped waveguide emitters in which the emissive layer contains silicon nanoclusters and erbium ions. Plasmonic coupling to metallic or semi-metallic overlayers provides enhancement of the radiative rate of erbium ions, enabling high quantum efficiency emission. Using 2D and 3D finite difference time domain (FDTD) simulations we show that up to 75% of the light emitted from the active layer can be coupled into a nanowire silicon rib waveguide. Our results suggest that such devices, which can readily be fabricated using CMOS processing techniques, pave the way for viable waveguide optical sources to be realized in silicon photonics.
©2011 Optical Society of America
1. Introduction
The lack of an efficient silicon-based optical source is the most pressing challenge facing silicon photonics. Silicon is an indirect band gap semiconductor and suffers from strong nonradiative decay processes—as a result it is a very inefficient light source. Nevertheless, a wide range of passive optical devices such as filters, modulators and detectors have been demonstrated on silicon optical platforms using silicon-on-insulator (SOI) technologies [1]. There are many close-to-market applications driving this research, the leaders being optical interconnect technology, telecommunications (particularly the enablement of fiber-to-the-home) and lab-on-a-chip devices. Further into the future, photonics is well placed to deliver quantum computing, optical data storage, and a range of novel biophotonics technologies. The immediate and pressing requirement to push photonics technology into the mainstream is for a monolithic silicon light source [2].
It is well known that luminescent active materials placed within microstructures that support dense optical modes exhibit enhanced spontaneous emission rates due to Purcell effects [3,4]. This phenomenon has already been exploited to create efficient Light Emitting Diodes (LEDs), using InGaN/GaN quantum dots [5,6] and organic compounds [7,8]. The luminescent active material in these devices was placed between electrode layers in the form of sub-100nm thick films, increasing decay rates by emitting into Surface Plasmon Polariton (SPP) modes [9–12].
We have previously investigated electroluminescence from erbium-doped silicon-rich silica (SiOx) [13,14], a material that has been of considerable interest since the demonstration of sensitization of erbium luminescence by excitation transfer from silicon nanoclusters [15]. This is a very appealing material for the fabrication of light emitting devices as it has the potential to overcome the limitations of silicon as an optical material while retaining CMOS processing compatibility. Recently, due to Purcell enhancement, a five-fold increase in emission from erbium doped silicon nanoclusters were observed in silicon slot waveguides [16]. Theoretical investigations on CMOS compatible metal-dielectric-metal (MDM), metal-dielectric-semiconductor (MDS) waveguiding structures have demonstrated significant rate enhancements even at wavelengths around 1535nm, well away from the plasmonic resonances of the metallic layers [17–20].
Moreover, a further Purcell enhancement of the radiative emission rate of erbium ions is produced by their proximity to silicon nanoclusters in the oxide matrix [21]. Therefore, the possibility of fabricating a silicon-based electrically pumped waveguide optical source with a high quantum efficiency is an appealing one [22]. Here we present an analysis of such structures based on the metal-dielectric-semiconductor-dielectric (MDSD) geometry [19].
The device structure is illustrated in Fig. 1 . A silicon rib waveguide sits on top of an oxide layer that provides optical and electrical isolation from the underlying silicon substrate. On top of this is a thin layer of erbium-doped silicon-rich silica (Er:SiOx), on top of which is an electrode, which may be metallic (e.g., silver) or semi-metallic (e.g., indium tin oxide, ITO). Electrical contact to the electrode and silicon layers is provided by metallization, which is shown only schematically on the figure. The structure may readily be fabricated using silicon-on-insulator (SOI) wafers, allowing the waveguide layer to be single crystal silicon. This allows the waveguide losses to be minimized.
2. Enhancement of emission rates
We begin our analysis by investigating how the decay rates of emitters, specifically erbium in the SiOx system, are modified when embedded in MOS structures. A cross-section of an MDSD structure such as that in Fig. 1 is shown in Fig. 2 .
The changes in erbium decay rates can be evaluated using the dipole approximation derived in [4,12,23]. We consider a dipole located at the origin of the coordinates shown in Fig. 2—that is, at the center of the dielectric layer. The normalized decay rate of the dipole under investigation, with respect to a dipole located far from any interface, is given by
Whereis the internal quantum efficiency of the dipole, the classically dissipated power, andthe in-plane wavevector. Considering the magnitude of the wavevector in the dielectric (), the out of plane wavevector,, the dipole orientation vector () and the modification to the electric field at the dipole (), the integrand in Eq. (1), can be evaluated to bewhereIn our notation, superscripts and refer to the parallel polarized magnetic field (TM) and the parallel polarized electric field (TE). Subscriptsand denote the perpendicular and parallel orientations of the dipole, respectively, and subscriptsare used to indicate the order of the dielectric stack (air-metal-dielectric-semiconductor-dielectric) as shown in Fig. 1. andrefer to the reflection coefficients calculated by applying Fresnel’s equations towards the slab components in and directions from the dipole. We consider the perpendicular and the parallel orientations for the dipoles separately to obtain the decay rates and. Assuming that the dipoles are isotropic, these rates are combined to give the average normalized decay rate [12]:
Figure 3(a) and (b) show the power spectrum (the result of Eq. (2) and the normalized decay rates of a dipole located in the middle of a 20nm thick dielectric layer in MDM, MDS and MDSD waveguide stacks. In all cases we have used silver [24] as the metal and silicon as the semiconductor material. Where necessary, the dielectric constant of the medium with the emitters is taken to be 2.25, approximating that of erbium doped silicon-rich silica [25], and the dielectric constant of the buried oxide layer (box) is taken to be 2.074 (silica).
According to Fig. 3, the rate enhancement obtained using the MDSD or MDS structure is considerably lower than that for an MDM waveguide. However, when used as an on-chip integrated device, a MDSD waveguide structure will have certain advantages over a MDM waveguide. As an example, reduced skin depth of the MDSD modes will reduce absorption losses [26,27]. Furthermore, as we demonstrate in the next section, the MDSD modes can be designed to have a very high mode overlap with low loss, high-index-contrast silicon nanowire waveguides [28], facilitating extremely efficient light coupling between the plasmonic and the photonic modes.
We next calculate the rate enhancement for a distribution of dipoles concentrated in the dielectric layer of an MDSD structure. We consider the dipole distribution along the axis to be, then the total luminescence intensity is written as a function of and time () as
may then be fitted with a stretched exponential function to obtain the averaged decay rates for the ensemble of emitters [29,30]. In Fig. 4 we plot the calculated decay rates for a MDSD structure with different dielectric layer (active layer) thicknesses and a continuous distribution of emitters. The active layer is separated from the top Ag layer (150nm thick) by a 10nm spacer layer to avoid light coupling to lossy surface waves [4]. The calculated rate enhancements at a wavelength of 1535nm vary from ~13 for =1 to ~7 for =0.53, the latter being close to reported quantum efficiencies from silicon nanocluster based active material [31].By integrating over the relevant poles of the power spectrum obtained from the emitter distribution (Fig. 5(a) ), the power coupled into different decay channels (TE and TM modes, lossy waves, or radiation) as a function of the active layer width of the MDSD structure can be calculated. The results of such calculations are shown in Fig. 5(b). Significantly, 85% of the light emitted from the dipoles is coupled in to the TM mode for an active layer thickness of 30nm.
3. Light coupling from active layer to silicon waveguide
3.1. 2D simulations
Initially, we used 2-dimensional FDTD simulations using a freely available software package (MEEP) [32] to investigate light coupling from the above analyzed MDSD waveguide to a low loss silicon nano-wire waveguide [28] for long range propagation. A diagram of the simulation space and the field intensity plot is shown in Fig. 6 .
A line source placed across the height of the active layer was used to simulate the distributed emitters (erbium ions). Transmission and reflection were measured by line flux meters placed across the waveguides at either side of the device edge. The flux meter length was determined by separately evaluating the modal area incorporating more than 90% of the light propagating in infinitely long MDSD/silicon waveguide modes. A 2nm×2nm grid resolution was used in our simulations, which reduced the relative error by ~1.5% in transmission from a grid size of 4nm×4nm at a wavelength of 1535nm. The reduction in relative error obtained by further enhancing the resolution up to 1nm×1nm was less than 0.5%.
Figure 7 shows the calculated transmission, reflection and losses at the device edge as functions of active layer thickness and Si waveguide height. The length of the MDSD section was kept constant at 500nm.
As observed from Fig. 7(a), the efficiency of coupling into the silicon waveguide reduces as the active layer thickness increases. Figure 7(b) was obtained for a constant active layer thickness of 30nm. The optimum height of the silicon waveguide was found to be ~230nm, at which point 88% of the 1535nm light in the MDSD waveguide is coupled into the silicon waveguide. The transmission, reflection and loss spectra near optical wavelengths for these device dimensions are shown in Fig. 7(c). It was also observed that increasing the MDSD device length improves the transmission efficiency, as this allows the light emitted in the active layer to be fully coupled to the MDSD modes before reaching the edge of the device.
3.2. 3D simulations
Considering the values obtained from the 2-D simulations, comparable coupling efficiencies (~89% transmission) were obtained for an optimized 3D device structure. The MDSD length of the simulated device was 600nm, coupling into a silicon nanowire waveguide with a width of 450nm and a height of 260nm (Fig. 8(a) ).
The electric field profile (Fig. 8(b)) of the TM mode for the optimized 3D structure was obtained using the method of Finite Elements using a commercially available software package. The effective index () of the MDSD guided mode was found to be 2.49, and the attenuation coefficient was 5.57x10−2µm−1. Therefore, a 600nm long device will only experience an attenuation of 0.145dB due to the absorption of the metallic contact.
4. Conclusion
An MDSD structure with a silver top contact, an active layer thickness of 30nm, silicon waveguide height of 260nm and a width of 450nm increased the decay rates of an ensemble of emitting erbium ions (η=0.53), by a factor of seven. 85% of the light emitted from the active layer was coupled in to the TM mode of the MDSD structure. For the above device dimensions, the MDSD TM mode has an of 2.49 and an attenuation coefficient of 5.57x10−2 µm−1, indicating very low losses. FDTD simulations indicate that 88% coupling efficiencies at a wavelength of 1535nm are achievable between a 600nm long MDSD waveguide of the said dimensions and a silicon waveguide. Thus, 75% of the light generated in the active layer by electrical pumping can couple into the silicon waveguide beyond the electrical contacts. This can be a promising solution for silicon based emission technologies.
References and links
1. M. Lipson, “Guiding, modulating, and emitting light on silicon-challenges and opportunities,” J. Lightwave Technol. 23(12), 4222–4238 (2005). [CrossRef]
2. L. Pavesi, “Will silicon be the photonic material of the third millenium?” J. Phys. Condens. Matter 15(26), R1169–R1196 (2003). [CrossRef]
3. E. Purcell, “Spontaneous emission probabilities at radio frequencies,” Phys. Rev. 69, 681 (1946).
4. G. Ford and W. Weber, “Electromagnetic interactions of molecules with metal surfaces,” Phys. Rep. 113(4), 195–287 (1984). [CrossRef]
5. K. Okamoto, I. Niki, A. Scherer, Y. Narukawa, T. Mukai, and Y. Kawakami, “Surface plasmon enhanced spontaneous emission rate of InGaN/ GaN quantum wells probed by time-resolved photoluminescence spectroscopy,” Appl. Phys. Lett. 87(7), 071102 (2005). [CrossRef]
6. J. Vuckovic, M. Loncar, and A. Scherer, “Surface plasmon enhanced light-emitting diode,” IEEE J. Quantum Electron. 36(10), 1131–1144 (2000). [CrossRef]
7. K. Saxena, V. Jain, and D. S. Mehta, “A review on the light extraction techniques in organic electroluminescent devices,” Opt. Mater. 32(1), 221–233 (2009). [CrossRef]
8. M. Ramuz, L. Burgi, R. Stanley, and C. Winnewisser, “Coupling light from an organic light emitting diode (OLED) into a single-mode waveguide: Toward monolithically integrated optical sensors,” J. Appl. Phys. 105(8), 084508 (2009). [CrossRef]
9. A. V. Zayats, I. I. Smolyaninov, and A. A. Maradudin, “Nano-optics of surface plasmon polaritons,” Phys. Rep. 408(3-4), 131–314 (2005). [CrossRef]
10. E. Ozbay, “Plasmonics: merging photonics and electronics at nanoscale dimensions,” Science 311(5758), 189–193 (2006). [CrossRef] [PubMed]
11. S. Wedge, J. 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(2), 182 (2004). [CrossRef]
12. W. L. Barnes, “Fluorescence near interfaces: the role of photonic mode density,” J. Mod. Opt. 45(4), 661–699 (1998). [CrossRef]
13. O. Jambois, Y. Berencen, K. Hijazi, M. Wojdak, A. J. Kenyon, F. Gourbilleau, R. Rizk, and B. Garrido, “Current transport and electroluminescence mechanisms in thin SiO2 films containing Si nanocluster-sensitized erbium ions,” J. Appl. Phys. 106(6), 063526 (2009). [CrossRef]
14. O. Jambois, F. Gourbilleau, A. J. Kenyon, J. Montserrat, R. Rizk, and B. Garrido, “Towards population inversion of electrically pumped Er ions sensitized by Si nanoclusters,” Opt. Express 18(3), 2230–2235 (2010). [CrossRef] [PubMed]
15. A. J. Kenyon, P. F. Trwoga, M. Federighi, and C. W. Pitt, “Optical properties of PECVD erbium doped silicon-rich silica—Evidence for energy transfer between silicon microclusters and erbium ions,” J. Phys. Condens. Matter 6(21), L319–L324 (1994). [CrossRef]
16. M. Galli, D. Gerace, A. Politi, M. Liscidini, M. Patrini, L. Andreani, A. Canino, M. Miritello, R. L. Savio, A. Irrera, and F. Priolo, “Direct evidence of light confinement and emission enhancement in active silicon-on-insulator slot waveguides,” Appl. Phys. Lett. 89(24), 241114 (2006). [CrossRef]
17. Y. Jun, R. Kekatpure, J. White, and M. Brongersma, “Nonresonant enhancement of spontaneous emission in metal-dielectric-metal plasmon waveguide structures,” Phys. Rev. B 78(15), 153111 (2008). [CrossRef]
18. Y. C. Jun, R. M. Briggs, H. A. Atwater, and M. L. Brongersma, “Broadband enhancement of light emission in silicon slot waveguides,” Opt. Express 17(9), 7479–7490 (2009). [CrossRef] [PubMed]
19. A. Hryciw, Y. C. Jun, and M. L. Brongersma, “Plasmon-enhanced emission from optically-doped MOS light sources,” Opt. Express 17(1), 185–192 (2009). [CrossRef] [PubMed]
20. A. Hryciw, Y. C. Jun, and M. L. Brongersma, “Plasmonics: electrifying plasmonics on silicon,” Nat. Mater. 9(1), 3–4 (2010). [CrossRef] [PubMed]
21. P. Horak, W. H. Loh, and A. J. Kenyon, “Modification of the Er3+ radiative lifetime from proximity to silicon nanoclusters in silicon-rich silicon oxide,” Opt. Express 17(2), 906–911 (2009). [CrossRef] [PubMed]
22. R. J. Walters, R. V. A. van Loon, I. Brunets, J. Schmitz, and A. Polman, “A silicon-based electrical source of surface plasmon polaritons,” Nat. Mater. 9(1), 21–25 (2010). [CrossRef] [PubMed]
23. R. Chance, A. Prock, and R. Silbey, “Lifetime of an emitting molecule near a partially reflecting surface,” J. Chem. Phys. 60(7), 2744 (1974). [CrossRef]
24. A. D. Rakic, A. B. Djurisic, J. M. Elazar, and M. L. Majewski, “Optical properties of metallic films for vertical-cavity optoelectronic devices,” Appl. Opt. 37(22), 5271–5283 (1998). [CrossRef] [PubMed]
25. N. Daldosso, D. Navarro-Urrios, M. Melchiorri, L. Pavesi, C. Sada, F. Gourbilleau, and R. Rizk, “Refractive index dependence of the absorption and emission cross sections at 1.54μm of Er coupled to Si nanoclusters,” Appl. Phys. Lett. 88(16), 161901 (2006). [CrossRef]
26. J. Dionne, L. Sweatlock, H. Atwater, and A. Polman, “Plasmon slot waveguides: towards chip-scale propagation with subwavelength-scale localization,” Phys. Rev. B 73(3), 035407 (2006). [CrossRef]
27. A. V. Krasavin and A. V. Zayats, “Silicon-based plasmonic waveguides,” Opt. Express 18(11), 11791–11799 (2010). [CrossRef] [PubMed]
28. Y. A. Vlasov and S. J. McNab, “Losses in single-mode silicon-on-insulator strip waveguides and bends,” Opt. Express 12(8), 1622–1631 (2004). [CrossRef] [PubMed]
29. J. Bao, N. Yu, F. Capasso, T. Mates, M. Troccoli, and A. Belyanin, “Controlled modification of erbium lifetime in silicon dioxide with metallic overlayers,” Appl. Phys. Lett. 91(13), 131103 (2007). [CrossRef]
30. P. Worthing, R. Amos, and W. Barnes, “Modification of the spontaneous emission rate of Eu3+ ions embedded within a dielectric layer above a silver mirror,” Phys. Rev. A 59(1), 865–872 (1999). [CrossRef]
31. R. Walters, J. Kalkman, A. Polman, H. Atwater, and M. de Dood, “Photoluminescence quantum efficiency of dense silicon nanocrystal ensembles in SiO2,” Phys. Rev. B 73(13), 132302 (2006). [CrossRef]
32. A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Comput. Phys. Commun. 181(3), 687–702 (2010). [CrossRef]