Plasmonics offers a unique opportunity to break the diffraction limit of light and bring photonic devices to the nanoscale. As the most prominent example, an integrated nanolaser is a key to truly nanoscale photonic circuits required for optical communication, sensing applications and high-density data storage. Here, we develop a concept of an electrically driven subwavelength surface-plasmon-polariton nanolaser, which is based on a novel amplification scheme, with all linear dimensions smaller than the operational free-space wavelength λ and a mode volume of under λ3/30. The proposed pumping approach is based on a double-heterostructure tunneling Schottky barrier diode and gives the possibility to reduce the physical size of the device and ensure in-plane emission so that the nanolaser output can be naturally coupled to a plasmonic or nanophotonic waveguide circuitry. With the high energy efficiency (8% at 300 K and 37% at 150 K), the output power of up to 100 μW and the ability to operate at room temperature, the proposed surface plasmon polariton nanolaser opens up new avenues in diverse application areas, ranging from ultrawideband optical communication on a chip to low-power nonlinear photonics, coherent nanospectroscopy, and single-molecule biosensing.
Nanophotonics is ubiquitous in many applications ranging from optical interconnects and sensing to high-density data storage and information processing , , , , , . Particularly, the plasmonic approach offers an opportunity to deal with photonic signals at the nanoscale via coupling to the free-electron oscillations in a metal. It can provide ultracompact components for optical interconnects such as modulators, photodetectors, waveguides and incoherent and coherent nanoscale optical sources , , , , , , . The implementation of the latter is a great challenge since a significant amount of the surface plasmon field is concentrated in the metal that results in high Joule losses. Nevertheless, the idea of SPASER (from surface plasmon amplification by stimulated emission radiation) ,  was proposed inspiring further research for novel truly nanoscale optical sources with efficient pumping schemes. While incoherent sources of surface plasmon polaritons (SPPs) can be relatively easily integrated in a nanophotonic circuitry , , , the situation is more complicated in the case of coherent SPP sources and amplifiers, since to reach the lasing threshold one has to compensate not only Joule but also radiation losses, which can be quite large due to a small mode volume of the plasmonic cavity , . The threshold gain may become very high, posing a serious challenge for the creation of the required population inversion. Although full compensation of the SPP losses at the nanoscale has been demonstrated , , , , , , it was achieved under optical pumping, which is hardly compatible with on-chip integration , , , , the target area of prospective applications. In this regard, an electrical pumping scheme is undeniably more desired and must be developed for practical application of plasmonic components , , , .
The possibility to implement electrical pumping in plasmonic devices is strongly restrained by a very limited choice of low-loss plasmonic metals. Gold, silver and copper form rectifying Schottky contacts to direct-bandgap semiconductors , while low-resistance ohmic contacts are needed for efficient electron and hole injection and creation of the population inversion in the gain medium , . Therefore, one has to use ohmic contacts made of highly absorptive materials (such as titanium, chromium, palladium, and different alloys , , , , ), place it at a large distance from the active region of the structure and rely on the refractive index contrast between semiconductor materials to localize the optical mode and prevent its overlap with the ohmic contact , . On the other hand, many proposed nanolaser designs based on low-absorptive structures (see, for example, Refs. , , ) can hardly be pumped electrically due to the high-resistance of metal-semiconductor contacts, which does not allow to inject a sufficient number of electrons and holes into the gain medium for the creation of the population inversion and full compensation of the high optical losses of the lasing mode. The high-resistance contact limitation greatly complicates the design of electrically pumped nanolasers and increases their physical size compared to the optically pumped counterparts. To overcome this problem, it was proposed to use the minority carrier injection properties of Schottky contacts and inject electrons into the gain semiconductor medium directly from the plasmonic metal, avoiding the commonly used intermediate semiconductor electron injection layer , . However, there is only one binary III–V semiconductor material (InAs) that forms Schottky contacts to plasmonic metals with a barrier height greater than the bandgap energy of the semiconductor , , which is required for the creation of the population inversion , . This limits the operation wavelength to around 3 µm and the operation temperature to below 150 K . In addition, the inability of the Schottky barrier to block the majority carriers under forward bias leads to a quite high leakage current even in the presence of a single heterostructure , which greatly limits the energy efficiency of the amplification scheme.
Here, we propose an efficient electrically driven single-mode deep-subwavelength coherent optical source featuring a unique combination of strong nanophotonic guiding and a novel amplification scheme based on a double-heterostructure tunneling Schottky barrier diode, which gives the possibility to naturally implement electrical pumping by using the SPP supporting interface as a low-resistance n-type tunneling contact. This approach allows dramatic reduction of the physical size of the device, simplifies its design, ensures in-plane emission and increases the operating temperature to above 300 K. Contrary to the electrically pumped metal coated nanopillar lasers , , , , which radiate photons into free space in a dipole-like fashion, the proposed SPP source is intrinsically integrable with nanophotonic circuitry. Unlike nanolasers based on localized optical modes, which are difficult to pump electrically , the subwavelength size of the proposed electrical SPP source is achieved with a propagating plasmonic mode, and the device can be referred as a SPPASER (surface plasmon polariton amplification by stimulated emission of radiation) or an SPP spaser . The proposed approach enables the SPP to be efficiently out-coupled to a bus plasmonic or photonic waveguide without losing energy through the emission into the surrounding space, which is particularly important in large-scale integrated nanophotonic circuits, where performance and scalability are limited by the crosstalk noise.
2 Results and discussion
2.1 Implementation of electrical pumping
The proposed amplification scheme based on the Au/n+-InAs0.4P0.6/In0.72Ga0.28As/p+-Al0.29In0.71As double-heterostructure tunneling Schottky barrier diode is shown in Figure 1 (see Supplementary Information Section 1 for details). It gives the possibility to use the low-loss SPP supporting interface between gold and the 100-nm-thick electron injection n+-InAsP layer as a low-resistance tunneling contact to the n+-InAsP electron injection layer with nearly ideal ohmic characteristics. Therefore, the material gain in InGaAs is provided in the near-field proximity to the gold contact, exactly where the plasmonic modes are localized, which gives a crucial advantage over traditional electrical pumpings schemes of plasmonic and photonic nanolasers based on highly absorptive ohmic contacts , , , , . The n-doped InGaAs layer sandwiched between the p+-AlInAs and n+-InAsP layers with bandgap energies exceeding that of InGaAs acts as an active region of the semiconductor structure. Under high forward bias, electrons and holes are injected into InGaAs from the InAsP and AlInAs sides, respectively, while high potential barriers for electrons at the InGaAs/AlInAs heterojunction and for holes at the InAsP/InGaAs heterojunction confine excess carriers to the active InGaAs region, which gives the double-heterostructure tunneling Schottky barrier diode amplification scheme a significant advantage over the amplification schemes based on Schottky diodes, which suffer from high leakage currents , , . Although the Au/n-InP Schottky barrier is as high as 0.5 eV , the Fermi level in Au/InAs contacts is anomalously pinned to the conduction band of InAs , . Therefore, a wide range of InAsP ternary alloys form contacts to gold with a quite small barrier height. In particular, the height of the Au/n-InAs0.4P0.6 barrier does not exceed 0.2 eV , . The barrier width is determined by the density of donors and is equal to only 13 nm at a moderate doping level of . Electrons can easily tunnel through this barrier thanks to their extremely small effective mass (0.057me0) (see Supplementary Information Section 2). The evaluated specific contact resistance of the tunneling Au/n+-InAs0.4P0.6 contact is as low as 1.2 × 10−7 Ohm cm2 (Methods). The voltage drop across the contact does not exceed 7 mV even at a pump current density of 50 kA/cm2 (Figure 1c), which ensures nearly ideal ohmic characteristics of the Au/n-InAs0.4P0.6 contact and efficient injection of electrons into the InAs0.4P0.6 layer, which can be further injected into the active InGaAs region. In comparison, InP, which would be used as an InGaAs lattice-matched material for the laser operation at 1.55 µm wavelength, with the same donor concentration of 3 × 1018 cm−3 forms a rectifying contact to gold (see Supplementary Information Section 2). Its contact resistance is two orders of magnitude higher than that of the Au/n+-InAs0.4P0.6 contact, which results in high heat dissipation, which does not allow to use the Au/InP contact for electrical pumping of metal-semiconductor nanolasers at room temperature. In other words, the performance of the proposed amplification scheme is much better than that based on compound semiconductors lattice-matched to InP.
2.2 Design of the nanolaser cavity
An essential component of the coherent SPP source is a cavity. Here, we use a subwavelength ring-resonator design based on a T-shaped waveguiding geometry (Figure 1). The InAsP/InGaAs/AlInAs layer stack on AlInAs is ring-shaped and surrounded by silicon oxide, while a gold superstrate completes the cavity structure and forms an SPP supporting interface. The ring width w = 300 nm is designed to be about 7 times smaller than the source operating wavelength to provide high lateral confinement of plasmonic modes propagating along the metal-semiconductor interface of the ring resonator. The ring width can be easily decreased below λ/10, but this procedure does not noticeably improve the mode localization or even deteriorates it, whereas the bending losses are substantially magnified due to the decrease in the mode effective index. The ring height of 1 µm is smaller than the light wavelengths in the targeted spectral region, compelling parasitic photonic modes to leak into the semiconductor substrate with a high refractive index. At the same time, the height is large enough to prevent strongly localized SPP modes from undesirable radiation losses. The resulting aspect ratio between the ring height H and the ring width w is 3.3:1 that should not present any problem for the device fabrication . The thickness of the electron injection InAsP layer is chosen to be 50 nm, which is significantly greater than the depletion layer thickness of the Au/n+-InAsP contact (13 nm) and allows efficient electron injection into the active InGaAs region. At the same time, the InGaAs layer is set to be 350 nm thick to provide more than 74% confinement of the plasmonic modes to the active InGaAs region to reduce the threshold gain (Methods). The ring radius R = 0.85 μm is selected so that one of the fundamental plasmonic modes of the ring resonator is close to the maximum of the gain spectrum of InGaAs at both room and low temperatures. Finally, we note that the proposed pumping scheme is universal and can also be used with other types of resonators.
Figure 2 shows the device mode structure at room temperature, obtained using the three-dimensional eigenmode simulations (Methods). Five types of modes (, , , , and ) with the quality factors (Q-factors) greater than 10 are supported by the resonator within the spontaneous emission bandwidth of InGaAs. Here, the Q-factor is defined as a ratio of the energy of the SPP stored in the resonator to the energy dissipated per radian of the oscillation. and modes having two maxima in the vertical direction are pushed out of the cross section of the ring waveguide and radiate into the substrate. Photonic modes are much better localized in the ring waveguide, but their mode effective indices are significantly lower than the refractive index of AlInAs, which leads to strong coupling to substrate radiation modes. As a result, the Q-factors of these modes are less than 50, corresponding to modal losses of more than 3000 cm−1, which are too high to be compensated by InGaAs. It takes modes out of the competition for the gain since the highly localized plasmonic modes experience 5 times smaller losses and show a significant advantage in mode confinement to the active InGaAs region. Owing to the submicrometer radius of the ring resonator, the spectral spacing between the modes is of the order of the material gain bandwidth of InGaAs, which establishes favorable conditions for the single mode stimulated emission into the eighth order mode ( mode) with a frequency of 0.637 eV (λ = 1.95 µm in free space). For the chosen ring radius of 850 nm, the resonator mode volume is only 0.033λ3, where λ is the free-space wavelength, and can be further reduced by decreasing the ring radius (Methods). However, it will result in the increase of emission into free space, while ideally, this process needs to be suppressed. Thus, the chosen ring radius ensures the ideal trade-off between the device footprint and performance. Finally, to achieve an integrated coherent SPP source, the electrically driven plasmonic resonator is coupled to a bus plasmonic waveguide of the same type, which serves as a device output (Figure 3a).
2.3 SPP nanolaser operation
To analyze the device operation, we have developed a self-consistent steady-state model describing the carrier transport within the semiconductor, which comprises the Poisson equation connecting the electrostatic potential and the charge carrier densities, the drift-diffusion current equations, the electron and hole continuity equations and the rate equations for the resonator modes (Methods). These optoelectronic processes are simulated with the finite difference method combined in a self-consistent way with the finite element analysis of the excited optical modes, and hence the simulations take into account the Purcell effect and coupling to a bus waveguide in the rate equations (for details, see Supplementary Information Sections 5 and 6).
Figure 3 shows the results of the self-consistent optoelectronic simulations. Under a negative voltage applied to the top gold contact (forward bias), electrons and holes are injected into the active InGaAs region of the ring resonator and recombine there both radiatively (spontaneous and stimulated emission) and non-radiatively (Auger recombination) (Figure 3c). At small forward bias (not shown in Figure 3c), the recombination for the spontaneous emission exceeds the Auger recombination, but the high threshold optical gain requires high concentrations of nonequilibrium electrons and holes, and non-radiative processes begin to dominate over radiative ones as the bias increases. Eventually, more than 90% of the total current goes into the Auger recombination just below the threshold. The situation, however, changes dramatically above the threshold, when recombination due to stimulated emission into the spasing mode surpasses the other recombination processes (Figure 3c).
Despite the fact that the Q-factors of , , , and are of the same order of magnitude (Figure 2а), internal absorption of InGaAs due to band-to-band and intervalence band transitions leads to a significant difference in the modal loss at zero bias (see Supplementary Information Section 4). The mode being closer to the bandgap energy of InGaAs suffers less absorption and could be the first to reach the threshold, but the material gain spectrum of InGaAs changes as the electron and hole concentrations increase and its maximum is blue-shifted, which gives the mode an advantage over the other modes (Figure 3f). At a current density of 8.4 kA/cm2, the output power of the mode becomes smaller than that of the mode and does not change appreciably henceforward (Figure 3b). The and modes show behavior very similar to the mode but do not reach the threshold.
In the light-current characteristic of the mode at a current density of about 30 kA/cm2 there is a well-pronounced amplified spontaneous (ASE) kink (the current at which the derivative of the slope of the curve of the output power versus injection current in a log-log scale changes its sign ). At this current, the recombination rate for the stimulated emission exceeds that for the spontaneous emission into all the resonator modes and into free space (Figure 3c,g). The pronounced ASE kink is explained by the fact that the spontaneous emission spectrum of bulk InGaAs is fairly broad (Figure 2b), which, combined with a moderate Purcell factor in the considered waveguiding geometry (see Supplementary Information Subsection 5.3), results primarily in emission into free space. The spontaneous emission coupling factor β ranges from 0.1 at small bias to 0.06 at high forward bias above the threshold. At J = 30 kA/cm2, the quality factor of the mode is equal to 4700 corresponding to a sub-nanometer linewidth, which is extremely difficult to achieve at the nanoscale .
Due to the high β-factors, the emission properties of nanolasers are remarkably different from that of macroscopic lasers , , the ASE kink by itself does not guarantee that the output radiation is coherent right above it , while coherence is the most important and distinct feature of laser radiation. Also, it has been shown that the ASE kink in the light-current characteristic cannot be a measure of the threshold , . To find the spasing threshold of the mode, one should use the definition of the threshold based on the degree of the second-order coherence of the output radiation . The calculations show that the threshold current density is equal to 47 kA/cm2, which is slightly higher than the current density of the ASE kink (Methods). The output power of the mode at the threshold is 29 μW, while the output powers of the and modes are 100 and 300 times lower, respectively. The powers of , , , and modes are even smaller than that of the spurious mode. Thus, single-mode spasing is achieved.
At low temperatures, the device performance is greatly improved due to the decrease in the non-radiative recombination rate, since the threshold gain is achieved at a significantly smaller carrier concentration, and the Auger recombination coefficient in InGaAs is reduced from 3.8 × 10−28 cm6 s−1 at 300 K to 1.2 × 10−28 cm6 s−1 at 150 K and to about 2 × 10−29 cm6 s−1 at 77 K . In addition, the electron-phonon scattering rate in gold γph decreases as the temperature decreases , . Accordingly, the imaginary part of the permittivity of the 100-nm-thick gold film at a frequency of 0.65 eV is reduced by a factor of 1.35 and 1.6 and equals and at 150 K and 77 K, respectively (see Supplementary Information Subsection 3.2). This leads to the increase of the Q-factor of the mode to 260 at 150 K and 290 at 77 K, compared to Q = 210 at room temperature. The temperature decrease induces broadening of the semiconductor bandgap and, consequently, the corresponding shift of the material gain spectrum (Figure 4a). As a result, at 150 K, the mode with a frequency of 0.582 eV is now outside the amplification bandwidth of InGaAs. The Auger recombination rate dramatically decreases, and, below the threshold, most of the current contributes to spontaneous emission. The spontaneous emission spectrum is narrower than at 300 K, and the β-factor is twice larger (from 0.16 to 0.20 depending on the bias voltage). The current at which the ASE kink is observed, is equal to JASE = 1.1 kA/cm2. However, the light-current characteristic demonstrates a very smooth transition from spontaneous to stimulated emission, which is typical for high-β lasers , . At the same time, the stimulated emission rate exceeds the rate of the spontaneous emission into the spasing mode at Jst-sp = 660 A/cm2, which is 45 times lower than at room temperature. Interestingly, the threshold current does not decrease dramatically as temperature decreases, since the coherence of the nanolaser radiation is determined by the fluctuations in the free carrier densities in the gain medium, which cannot be reduced by orders of magnitude . At 150 K, the threshold current density is equal to Jth = 26 kA/cm2. At 77 K, the light-current characteristic is even smoother than at 150 K and the amplified spontaneous emission kink in the curve for the mode is shifted to lower values of the injection current (JASE = 115 A/cm2). The Auger recombination rate becomes negligibly small and does not have any impact on the output characteristics. The β-factor increases to 0.36. The rate of stimulated emission exceeds the rate of spontaneous emission at Jst-sp = 290 A/cm2, which is only twice lower than at 150 K. At the same time, the threshold current density is almost the same as at 150 K and is equal to 20 kA/cm2, which is typical for nearly thresholdless nanolasers .
2.4 Efficiency of the SPP ring nanolaser
Power consumption plays a decisive role in the viability of nanolasers . Generally, high Joule losses in a plasmonic resonator present a serious challenge for obtaining a large efficiency of electro-optical conversion since a significant portion of the generated SPP quanta is absorbed by the metal rather than emitted into the waveguide. In addition, at room temperature, due to the high threshold carrier concentration, non-radiative Auger recombination still gives a significant contribution to the total current above the threshold (Figure 3c). Nevertheless, in the proposed scheme, the energy efficiency (the ratio of the output power to the input electrical power) above the threshold is between 7% and 9% depending on the pump current, which is higher than for the state-of-the-art electrically pumped photonic-crystal nanocavity lasers , . Decreasing the operating temperature dramatically improves characteristics of the proposed SPP nanolaser. At 150 K, the electro-optical conversion efficiency exceeds 16%, which is comparable with the best semiconductor lasers integrated on a chip , . At cryogenic temperatures, the efficiency is improved to about 20%. This is mainly the result of the fundamental relation between spontaneous and stimulated emission into the spasing mode.
The energy efficiency can be further enhanced by optimizing the geometry of the SPP nanolaser. Particularly, by changing the separation distance d between the ring resonator and the bus waveguide, one can significantly increase the ratio of the output power at the bus waveguide to the power flow in the ring. As the separation distance d is decreased to 100 nm, the efficiency of electro-optical conversion at 150 K rises to 37% and is even higher at lower temperatures. With the further decrease of the separation distance, the coupling efficiency can be even higher . However, at room temperature, this does not provide any real advantages, since a higher coupling coefficient corresponds to higher modal loss, which substantially increases the threshold current.
The scheme proposed here for SPP amplification in deep-subwavelength plasmonic waveguides with integrated electrical pumping, based on a double heterostructure tunneling Schottky barrier diode, provides a realistic design of an electrically driven subwavelength coherent SPP source with dimensions orders of magnitude smaller than the typical size of in-plane coherent optical sources integrated on a chip , , , which are hardly compatible with nanophotonic circuitry. At the same time, unlike electrically pumped metal-coated nanolasers , , , the presented device does not merely emit light into free space in a dipole-like fashion, but efficiently couples the radiation into an in-plane plasmonic waveguide. It can also be directly coupled to a bus photonic semiconductor waveguide, which provides the compatibility with a standard low-loss photonic circuitry. The latter, however, comes at a price of smaller output power due to the increased scattering losses underlined by the spatial modal mismatch between the plasmonic (ring) and the photonic waveguides (see Supplementary Information Section 7). The gold/semiconductor contact plays a dual role of the SPP supporting interface and the electron injection electrode, which eliminates the need for any additional electrical contacts, minimizing the device footprint and greatly facilitating the design of large-scale circuits. The vertical injection scheme, large band-edge offsets at the heterojunctions, and very low resistivity of the tunneling Au/n+-InAsP Schottky contact provide favorable conditions for efficient carrier injection into the active layer and minimize the leakage current, which is limited to the trap-assisted recombination at the interfaces. In contrast to most of light-emitting nanophotonic devices, the surface recombination does not play a significant role in the lattice matched InAs0.4P0.6/In0.72Ga0.28As/Al0.29In0.71As system thanks to a very small surface recombination velocity (2 × 103 cm/s)  combined with the low surface-to-volume ratio of the active region. At room temperature, the threshold current is determined by the non-radiative Auger recombination (for details on the effect of surface recombination, see Supplementary Information Section 7). Above the threshold, the contribution of surface recombination to the total current does not exceed 1.5%. In this aspect, the proposed configuration of the active region and the proper material choice provide a remarkable advantage over the tensile strained multiple quantum wells, which, in principle, can also serve as a gain medium , but are characterized by much greater area of the interfaces, substantially higher surface recombination velocities due to the increased density of surface defects and poorer confinement of the SPP field to the active region. Here, we should note that poor-quality interfaces may give rise to the surface density of trapping centers and increase the surface recombination velocity up to ∼1 × 104 cm/s. Nevertheless, our simulations demonstrate that even in this case, at room temperature, the contribution of surface recombination will not exceed 3 kA/cm2 at the threshold, which is only about 10% of the total current.
The presented SPP source demonstrates remarkable characteristics, which can be further improved by structural and material optimization. Under room temperature conditions, the threshold is achieved at a moderate current density of 47 kA/cm2, ensuring continuous-wave operation. The output power of up to 100 μW at the energy efficiency of 8% is high enough to be used in future on-chip optical interconnects . This makes the proposed SPP nanolaser superior to electrically driven photonic-crystal nanocavity lasers, which are much more complicated in fabrication and are characterized by a significantly smaller output power despite the lower threshold currents , . As temperature decreases, the device performance is dramatically enhanced. At temperatures below 150 K, the energy efficiency increases to above 17% and can be further improved to more than 35% by decreasing the separation distance between the output waveguide and the resonator to 100 nm, which makes it comparable with the energy efficiency of the best external semiconductor lasers. The proposed pumping scheme provides a strong basis for the design and development of electrically pumped active plasmonic waveguides that ensure lossless SPP guiding at the nanoscale, subdiffraction-limited electrically-driven SPP nanolasers for on-chip optical interconnects or large-scale nanophotonic and plasmonic circuits. Using heterogeneous integration (either hybrid integration based die-to-wafer bonding  or direct monolithic integration based on heteroepitaxy , ), these nanolasers can also be integrated with silicon waveguides, providing highly efficient electrically-driven on-chip coherent sources for various nanophotonic applications.
4.1 Self-consistent optoelectronic simulations
The proposed SPP source is simulated using the developed comprehensive self-consistent steady-state hybrid electronic/photonic model, which comprises six first order nonlinear differential equations describing the carrier behavior within the semiconductor (see Equation (S35) in Supplementary Information) and M integral rate equations for M optical modes of the plasmonic resonator (see Equation (S51) in Supplementary Information), coupled with the eigenmode FEM optical simulations (see Supplementary Information Subsection 6.2). The carrier transport equations include the Poisson equation (treated as two first order differential equations) connecting the electrostatic potential and the carrier densities, the drift-diffusion current equations and the electron and hole continuity equations. The carrier continuity equations cover processes of carrier generation and recombination, namely non-radiative Auger recombination, spontaneous emission into free space and into the resonator modes enhanced by the Purcell effect and stimulated emission. The latter connects the differential equation for carrier transport with the integral rate equations for the resonator modes. At the same time, the rate equations include the material gain and spontaneous emission into the resonator modes, which are obtained by integrating the transition probabilities between electron states in the conduction and valence bands over the states’ energies. This introduces implicitly the electron and hole densities into M integral equations for the resonator modes.
Providing 18 boundary conditions at the material interfaces and using the electromagnetic field distributions of the resonator modes from the eigenmode optical simulations and the Purcell factor distribution from the 3D frequency domain FEM simulations (see Supplementary Information Subsection 5.3), the SPP source is self-consistently modeled using the finite difference method. The details on the optoelectronic model and simulations are described in Supplementary Information Section 6.
4.2 Electron transport through the tunneling Schottky contact
Tunneling and an over-barrier current across the metal-semiconductor contact is calculated using the exact solutions of the Schrödinger equation. The current density is given by
where q is the electron charge, V is the voltage across the contact, and are the effective electron masses in the semiconductor and metal, respectively, is the Fermi distribution function, υfm is the Fermi velocity in the metal, is density of states in the metal at the Fermi energy, and is the barrier transmission coefficient calculated from an electron wave function (for details, see Supplementary Information Section 2).
4.3 Mode structure simulations
The set of modes supported by the plasmonic resonator is found via the 3D finite element method (FEM) numerical simulations using COMSOL Multiphysics software. The outer domain boundaries are placed at a large distance from the structure to ensure their position in the radiation zone, so that the radiative components of the fields escaped from the resonator are efficiently absorbed using perfectly matched layers. For each resonator mode, a complex eigenfrequency is determined, the real part of which defines the actual resonant frequency, while the imaginary part is responsible for losses. From these values, the resonator quality factor and the modal loss are calculated. Here, υg is the group velocity of the corresponding optical mode of the waveguide of the ring resonator, which is found in the designated 2D FEM eigenmode simulations.
The dielectric function of gold used in the simulations is calculated taking into account electron-phonon, electron-electron, grain boundary and surface scattering, and verified by the experimental data. The temperature dependence of the dielectric function was determined considering the proper temperature variations of each scattering component. The dielectric functions of the compound semiconductors used in the simulations are calculated using the Clausius–Mossotti relation and their temperature dependences are obtained using an empirical model by Herve and Vandamme.
For details on the eigenmode simulations and the dielectric constants evaluation, see Supplementary Information Section 3.
4.4 Quality factor
The quality factor of the resonator is defined as a ratio of the energy stored in the resonator to the energy loss per radian of the oscillation . To evaluate the Q-factor under electrical pumping, we use the results from 3D eigenmode simulations of the passive structure and self-consistent optoelectronic simulations, which give the modal loss of each mode of the resonator under applied bias (for details, see Supplementary Information Section 6).
4.5 Mode volume
The mode volume of the resonator mode is calculated using the data retrieved from the 3D eigenmode simulations 
where is the complex amplitude of the electric field of the optical mode and the integration in the numerator is performed over the entire domain occupied by the mode.
4.6 Mode confinement to the active region
The confinement factor of the mode to the active region of the waveguide of the plasmonic ring resonator is determined as :
where is the energy density of the optical mode, is the refractive index of In0.72Ga0.28As, and is the mode group velocity. The integral in the numerator is taken over the active In0.72Ga0.28As region, while the integration in the denominator is performed over the whole 2D domain occupied by the mode.
4.7 Purcell factor
The Purcell factor, defining the change in the spontaneous decay rate of an emitter located in the active region of the waveguide of the plasmonic ring resonator from its value in the corresponding uniform semiconductor, is calculated as a ratio of the power flows emitted by a dipolar emitter placed in the respective electromagnetic environments:
where is the dipole moment, and and are the total electric fields at the dipole position. To find the electric field values, a dedicated set of 3D frequency domain numerical simulations of the dipole emission was performed. Varying the dipole position and averaging the results over all its orientations, the full Purcell factor map of the active region is obtained (for details, see Supplementary Information Subsection 5.3).
4.8 Threshold current
The spasing threshold was found using its definition through the second-order coherence properties of the emission as the transition from the emission in the thermal state of light to the emission in the coherent state of light. The threshold current density is obtained by solving the following transcendental equation :
where is the volume of the InGaAs active region, is the modal gain of the spasing mode in the ring resonator, is the averaged electron density in the InGaAs active region, is the angular frequency of the spasing mode, is the quality factor of the spasing mode that includes losses due to coupling to a bus waveguide and q is the elementary charge.
Funding source: Ministry of Science and Higher Education of the Russian Federation
Award Identifier / Grant number: 0714-2020-0002
Funding source: Russian Science Foundation
Award Identifier / Grant number: 17-79-20421
Funding source: Russian Foundation for Basic Research
Award Identifier / Grant number: 19-57-12008
Funding source: Engineering and Physical Sciences Research Council
Award Identifier / Grant number: Unassigned
Author contribution: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.
Research funding: D.Y.F. acknowledges support from the Russian Foundation for Basic Research (19-57-12008). A.V.K. and A.V.Z. acknowledge support from the Engineering and Physical Sciences Research Council (U.K.). A.V.A. acknowledges support from the Ministry of Science and Higher Education of the Russian Federation (0714-2020-0002). The development of the software for optoelectronic simulations was supported by the Russian Science Foundation (17-79-20421).
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
 K. C. Y. Huang, M.-K. Seo, T. Sarmiento, Y. Huo, J. S. Harris, and M. L. Brongersma, “Electrically driven subwavelength optical nanocircuits,” Nat. Photonics, vol. 8, pp. 244–249, 2014, https://doi.org/10.1038/nphoton.2014.2.Search in Google Scholar
 P. Wang, A. V. Krasavin, M. E. Nasir, W. Dickson, and A. V. Zayats, “Reactive tunnel junctions in electrically driven plasmonic nanorod metamaterials,” Nat. Nanotechnol., vol. 13, pp. 159–165, 2018, https://doi.org/10.1038/s41565-017-0017-7.Search in Google Scholar
 B. C. Stipe, T. C. Strand, C. C. Poon, et al., “Magnetic recording at 1.5 Pb m−2 using an integrated plasmonic antenna,” Nat. Photonics, vol. 4, pp. 484–488, 2010, https://doi.org/10.1038/nphoton.2010.90.Search in Google Scholar
 P. Wang, M. E. Nasir, A. V. Krasavin, W. Dickson, and A. V. Zayats, “Optoelectronic synapses based on hot-electron-induced chemical processes,” Nano Lett., vol. 20, pp. 1536–1541, 2020, https://doi.org/10.1021/acs.nanolett.9b03871.Search in Google Scholar
 W. Heni, Y. Fedoryshyn, B. Baeuerle, et al., “Plasmonic IQ modulators with attojoule per bit electrical energy consumption,” Nat. Commun., vol. 10, no. 8, p. 1694, 2019, https://doi.org/10.1038/s41467-019-09724-7.Search in Google Scholar
 A. V. Krasavin and A. V. Zayats, “Active nanophotonic circuitry based on dielectric-loaded plasmonic waveguides,” Adv. Opt. Mater., vol. 3, pp. 1662–1690, 2015, https://doi.org/10.1002/adom.201500329.Search in Google Scholar
 C. Haffner, D. Chelladurai, Y. Fedoryshyn, et al., “Low-loss plasmon-assisted electro-optic modulator,” Nature, vol. 556, pp. 483–486, 2018, https://doi.org/10.1038/s41586-018-0031-4.Search in Google Scholar
 J. E. Muench, A. Ruocco, M. A. Giambra, et al., “Waveguide-integrated, plasmonic enhanced graphene photodetectors,” Nano Lett., vol. 19, pp. 7632–7644, 2019, https://doi.org/10.1021/acs.nanolett.9b02238.Search in Google Scholar
 A. V. Krasavin and A. V. Zayats, “Benchmarking system-level performance of passive and active plasmonic components: integrated circuit approach,” Proc. IEEE, vol. 104, pp. 2338–2348, 2016, https://doi.org/10.1109/jproc.2016.2603118.Search in Google Scholar
 A. V. Krasavin, P. Wang, M. E. Nasir, Y. L. Jiang, and A. V. Zayats, “Tunneling-induced broadband and tunable optical emission from plasmonic nanorod metamaterials,” Nanophotonics, vol. 9, pp. 427–434, 2020, https://doi.org/10.1515/nanoph-2019-0411.Search in Google Scholar
 C. P. T. McPolin, J. S. Bouillard, S. Vilain, et al., “Integrated plasmonic circuitry on a vertical-cavity surface-emitting semiconductor laser platform,” Nat. Commun., vol. 7, no. 8, p. 12409, 2016, https://doi.org/10.1038/ncomms12409.Search in Google Scholar
 D. J. Bergman and M. I. Stockman, “Surface plasmon amplification by stimulated emission of radiation: quantum generation of coherent surface plasmons in nanosystems,” Phys. Rev. Lett., vol. 90, p. 027402, 2003, https://doi.org/10.1103/physrevlett.90.027402.Search in Google Scholar
 S. I. Azzam, A. V. Kildishev, R. M. Ma, et al., “Ten years of spasers and plasmonic nanolasers,” Light Sci. Appl., vol. 9, no. 21, p. 90, 2020, https://doi.org/10.1038/s41377-020-0319-7.Search in Google Scholar
 R. J. Walters, R. V. van Loon, I. Brunets, J. Schmitz, and A. Polman, “A silicon-based electrical source of surface plasmon polaritons,” Nat. Mater., vol. 9, pp. 21–25, 2010, https://doi.org/10.1038/nmat2595.Search in Google Scholar
 S. Kena-Cohen, P. N. Stavrinou, D. D. Bradley, and S. A. Maier, “Confined surface plasmon-polariton amplifiers,” Nano Lett., vol. 13, pp. 1323–1329, 2013, https://doi.org/10.1021/nl400134v.Search in Google Scholar
 C.-J. Lee, H. Yeh, F. Cheng, et al., “Low-threshold plasmonic lasers on a single-crystalline epitaxial silver platform at telecom wavelength,” ACS Photonics, vol. 4, pp. 1431–1439, 2017, https://doi.org/10.1021/acsphotonics.7b00184.Search in Google Scholar
 R. M. Ma, R. F. Oulton, V. J. Sorger, G. Bartal, and X. Zhang, “Room-temperature sub-diffraction-limited plasmon laser by total internal reflection,” Nat. Mater., vol. 10, pp. 110–113, 2011, https://doi.org/10.1038/nmat2919.Search in Google Scholar
 Y. L. Ho, J. K. Clark, A. Syazwan, A. Kamal, and J. J. Delaunay, “On-chip monolithically fabricated plasmonic-waveguide nanolaser,” Nano Lett., vol. 18, pp. 7769–7776, 2018, https://doi.org/10.1021/acs.nanolett.8b03531.Search in Google Scholar
 C. Xu, W. E. Hayenga, M. Khajavikhan, and P. Likamwa, “Measuring the frequency response of optically pumped metal-clad nanolasers,” Opt. Express, vol. 27, pp. 21834–21843, 2019, https://doi.org/10.1364/oe.27.021834.Search in Google Scholar
 A. W. Fang, H. Park, O. Cohen, R. Jones, M. J. Paniccia, and J. E. Bowers, “Electrically pumped hybrid algainas-silicon evanescent laser,” Opt. Express, vol. 14, pp. 9203–9210, 2006, https://doi.org/10.1364/oe.14.009203.Search in Google Scholar
 M. T. Hill, Y.-S. Oei, B. Smalbrugge, et al., “Lasing in metallic-coated nanocavities,” Nat. Photonics, vol. 1, pp. 589–594, 2007, https://doi.org/10.1038/nphoton.2007.171.Search in Google Scholar
 B. Romeira and A. Fiore, “Physical limits of nanoleds and nanolasers for optical communications,” Proc. IEEE, vol. 108, pp. 735–748, 2020, https://doi.org/10.1109/jproc.2019.2912293.Search in Google Scholar
 F. Vallini, Q. Gu, M. Kats, Y. Fainman, and N. C. Frateschi, “Carrier saturation in multiple quantum well metallo-dielectric semiconductor nanolaser: is bulk material a better choice for gain media?,” Opt. Express, vol. 21, pp. 25985–25998, 2013, https://doi.org/10.1364/oe.21.025985.Search in Google Scholar
 D. Y. Fedyanin, A. V. Krasavin, A. V. Arsenin, and A. V. Zayats, “Surface plasmon polariton amplification upon electrical injection in highly integrated plasmonic circuits,” Nano Lett., vol. 12, pp. 2459–2463, 2012, https://doi.org/10.1021/nl300540x.Search in Google Scholar
 S. Tiwari and D. J. Frank, “Empirical fit to band discontinuities and barrier heights in III–V alloy systems,” Appl. Phys. Lett., vol. 60, pp. 630–632, 1992, https://doi.org/10.1063/1.106575.Search in Google Scholar
 M. T. Hill, M. Marell, E. S. P. Leong, et al., “Lasing in metal-insulator-metal sub-wavelength plasmonic waveguides,” Opt. Express, vol. 17, pp. 11107–11112, 2009, https://doi.org/10.1364/oe.17.011107.Search in Google Scholar
 G. Stareev, H. Künzel, and G. Dortmann, “A controllable mechanism of forming extremely low‐resistance nonalloyed ohmic contacts to group III‐V compound semiconductors,” J. Appl. Phys., vol. 74, pp. 7344–7356, 1993, https://doi.org/10.1063/1.355002.Search in Google Scholar
 K. Ding, M. T. Hill, Z. C. Liu, L. J. Yin, P. J. van Veldhoven, and C. Z. Ning, “Record performance of electrical injection sub-wavelength metallic-cavity semiconductor lasers at room temperature,” Opt. Express, vol. 21, pp. 4728–4733, 2013, https://doi.org/10.1364/oe.21.004728.Search in Google Scholar
 K. Liu, N. Li, D. K. Sadana, and V. J. Sorger, “Integrated nanocavity plasmon light sources for on-chip optical interconnects,” ACS Photonics, vol. 3, pp. 233–242, 2016, https://doi.org/10.1021/acsphotonics.5b00476.Search in Google Scholar
 D. B. Li and C. Z. Ning, “Interplay of various loss mechanisms and ultimate size limit of a surface plasmon polariton semiconductor nanolaser,” Opt. Express, vol. 20, pp. 16348–16357, 2012, https://doi.org/10.1364/oe.20.016348.Search in Google Scholar
 S. Wang, H. Z. Chen, and R. M. Ma, “High performance plasmonic nanolasers with external quantum efficiency exceeding 10,” Nano Lett., vol. 18, pp. 7942–7948, 2018, https://doi.org/10.1021/acs.nanolett.8b03890.Search in Google Scholar
 A. A. Vyshnevyy and D. Y. Fedyanin, “Self-heating and cooling of active plasmonic waveguides,” ACS Photonics, vol. 3, pp. 51–57, 2016, https://doi.org/10.1021/acsphotonics.5b00449.Search in Google Scholar
 J. H. Lee, M. Khajavikhan, A. Simic, et al., “Electrically pumped sub-wavelength metallo-dielectric pedestal pillar lasers,” Opt. Express, vol. 19, pp. 21524–21531, 2011, https://doi.org/10.1364/oe.19.021524.Search in Google Scholar
 J. B. Khurgin and G. Sun, “Comparative analysis of spasers, vertical-cavity surface-emitting lasers and surface-plasmon-emitting diodes,” Nat. Photonics, vol. 8, pp. 468–473, 2014, https://doi.org/10.1038/nphoton.2014.94.Search in Google Scholar
 D. Costantini, A. Bousseksou, M. Fevrier, B. Dagens, and R. Colombelli, “Loss and gain measurements of tensile-strained quantum well diode lasers for plasmonic devices at telecom wavelengths,” IEEE J. Quantum Electron., vol. 48, pp. 73–78, 2012, https://doi.org/10.1109/jqe.2011.2178818.Search in Google Scholar
 J. N. Walpole and K. W. Nill, “Capacitance‐voltage characteristics of metal barriers on p PbTe and p InAs: effects of the inversion layer,” J. Appl. Phys., vol. 42, pp. 5609–5617, 1971, https://doi.org/10.1063/1.1659990.Search in Google Scholar
 K. Kajiyama, Y. Mizushima, and S. Sakata, “Schottky barrier height of n‐InxGa1−xAs diodes,” Appl. Phys. Lett., vol. 23, pp. 458–459, 1973, https://doi.org/10.1063/1.1654957.Search in Google Scholar
 R. Dylewicz, R. M. De La Rue, R. Wasielewski, P. Mazur, G. Mezősi, and A. C. Bryce, “Fabrication of submicron-sized features in InP/InGaAsP/AlGaInAs quantum well heterostructures by optimized inductively coupled plasma etching with Cl2/Ar/N2 chemistry,” J. Vac. Sci. Technol. B, vol. 28, pp. 882–890, 2010, https://doi.org/10.1116/1.3466811.Search in Google Scholar
 A. A. Vyshnevyy and D. Y. Fedyanin, “Lasing threshold of thresholdless and non-thresholdless metal-semiconductor nanolasers,” Opt. Express, vol. 26, pp. 33473–33483, 2018, https://doi.org/10.1364/oe.26.033473.Search in Google Scholar
 W. W. Chow, F. Jahnke, and C. Gies, “Emission properties of nanolasers during the transition to lasing,” Light Sci. Appl., vol. 3, no. 8, p. e201, 2014, https://doi.org/10.1038/lsa.2014.82.Search in Google Scholar
 T. H. Gfroerer, L. P. Priestley, M. F. Fairley, and M. W. Wanlass, “Temperature dependence of nonradiative recombination in low-band gap InxGa1−xAs/InAsyP1−y double heterostructures grown on inp substrates,” J. Appl. Phys., vol. 94, pp. 1738–1743, 2003, https://doi.org/10.1063/1.1586468.Search in Google Scholar
 J. S. Bouillard, W. Dickson, D. P. O’Connor, G. A. Wurtz, and A. V. Zayats, “Low-temperature plasmonics of metallic nanostructures,” Nano Lett., vol. 12, pp. 1561–1565, 2012, https://doi.org/10.1021/nl204420s.Search in Google Scholar
 K. Takeda, T. Sato, A. Shinya, et al., “Few-fj/bit data transmissions using directly modulated lambda-scale embedded active region photonic-crystal lasers,” Nat. Photonics, vol. 7, pp. 569–575, 2013, https://doi.org/10.1038/nphoton.2013.110.Search in Google Scholar
 K. Takeda, T. Sato, T. Fujii, et al., “Heterogeneously integrated photonic-crystal lasers on silicon for on/off chip optical interconnects,” Opt. Express, vol. 23, pp. 702–708, 2015, https://doi.org/10.1364/oe.23.000702.Search in Google Scholar
 B. R. Koch, E. J. Norberg, B. Kim, et al., “Integrated silicon photonic laser sources for telecom and datacom,” Optical Fiber Communication Conference/National Fiber Optic Engineers Conference 2013, OSA Technical Digest (online) (Optical Society of America, 2013), paper PDP5C.8, 2013, https://doi.org/10.1364/ofc.2013.pdp5c.8.Search in Google Scholar
 M. J. R. Heck and J. E. Bowers, “Energy efficient and energy proportional optical interconnects for multi-core processors: driving the need for on-chip sources,” IEEE J. Sel. Top. Quantum Electron., vol. 20, pp. 332–343, 2014, https://doi.org/10.1109/jstqe.2013.2293271.Search in Google Scholar
 S. Chen, W. Li, J. Wu, et al., “Electrically pumped continuous-wave III–V quantum dot lasers on silicon,” Nat. Photonics, vol. 10, pp. 307-311, 2016.10.1038/nphoton.2016.21Search in Google Scholar
 S. Manipatruni, M. Lipson, and I. A. Young, “Device scaling considerations for nanophotonic cmos global interconnects,” IEEE J. Sel. Top. Quantum Electron., vol. 19, p. 8200109, 2013, https://doi.org/10.1109/jstqe.2013.2239262.Search in Google Scholar
 E. Tournié, L. Cerutti, J.-B. Rodriguez, H. Liu, J. Wu, and S. Chen, “Metamorphic III–V semiconductor lasers grown on silicon,” MRS Bulletin, vol. 41, pp. 218–223, 2016, https://doi.org/10.1557/mrs.2016.24.Search in Google Scholar
 Z. Wang, B. Tian, M. Pantouvaki, et al., “Room-temperature inp distributed feedback laser array directly grown on silicon,” Nat. Photonics, vol. 9, pp. 837–842, 2015, https://doi.org/10.1038/nphoton.2015.199.Search in Google Scholar
 K. Srinivasan, M. Borselli, and O. Painter, “Cavity Q, mode volume, and lasing threshold in small diameter algaas microdisks with embedded quantum dots,” Opt. Express, vol. 14, pp. 1094–1105, 2006, https://doi.org/10.1364/oe.14.001094.Search in Google Scholar
 T. D. Visser, H. Blok, B. Demeulenaere, and D. Lenstra, “Confinement factors and gain in optical amplifiers,” IEEE J. Quantum Electron., vol. 33, pp. 1763–1766, 1997, https://doi.org/10.1109/3.631280.Search in Google Scholar
The online version of this article offers supplementary material (https://doi.org/10.1515/nanoph-2020-0157).
© 2020 Dmitry Yu. Fedyanin et al., published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.