A Novel Metal Nanoparticles-Graphene Nanodisks-Quantum Dots Hybrid-System-Based Spaser

Active nanoplasmonics have recently led to the emergence of many promising applications. One of them is the spaser (surface plasmons amplification by stimulated emission of radiation) that has been shown to generate coherent and intense fields of selected surface plasmon modes that are strongly localized in the nanoscale. We propose a novel nanospaser composed of a metal nanoparticles-graphene nanodisks hybrid plasmonic system as its resonator and a quantum dots cascade stack as its gain medium. We derive the plasmonic fields induced by pulsed excitation through the use of the effective medium theory. Based on the density matrix approach and by solving the Lindblad quantum master equation, we analyze the ultrafast dynamics of the spaser associated with coherent amplified plasmonic fields. The intensity of the plasmonic field is significantly affected by the width of the metallic contact and the time duration of the laser pulse used to launch the surface plasmons. The proposed nanospaser shows an extremely low spasing threshold and operates in the mid-infrared region that has received much attention due to its wide biomedical, chemical and telecommunication applications.


Introduction
Over the last decade, plasmonic nanosources have attracted significant attention as coherent and intense near-field generators of nanolocalized optical fields, eliminating the need for a coupling mechanism between photons and surface plasmons (SPs), and paving the way to many promising applications [1,2]. The surface plasmon waves, resulting from collective oscillations of free electrons near the surface of metals, propagate along the interface between metal and dielectric and decay in the direction perpendicular to the interface. The propagation length of surface plasmon polaritons is limited by the plasmonic losses that depends on the dielectric properties of metals [3]. For a metal bounded by an ideal dielectric, the losses are caused by free electron scattering and absorption through interband transitions at significantly shorter wavelengths [4].
Because the losses are often a serious limitation for practical applications of nanoplasmonics, many efforts have been devoted to proposing loss compensation and amplification approaches by introducing a gain medium in the dielectric surrounding the metal. Specifically, Seidel et al. reported the first experiment demonstrating amplification of surface plasmons on a flat silver film surrounded by dye molecules [5]. This work was followed by many theoretical and experimental investigations [6][7][8][9][10][11][12][13] that have paved the way to surface plasmon amplification by stimulated emission of radiation, i.e., spaser, introduced by Bergman and Stockman [14] and demonstrated experimentally with various plasmonic resonators, gain media and geometries [15][16][17][18][19]. Moreover, several theoretical approaches have been established to explain the spasing quantum mechanically, using two-level and three-level models for the gain medium [20,21]. Richter et al. have constructed a numerical approach to account for a large number of identical chromophores with an arbitrary number of energy levels non-perturbatively, based on density matrix theory through the use of the Tavis-Cummings model that has been introduced to describe the collective behavior of multiple atomic dipoles interacting with electromagnetic radiation [22].
Due to their ability to confine the optical energy near the surface, the spaser can utilize plasmonic components as its resonators to support the plasmonic modes and externally excited population-inverted gain media in order to provide the energy for spasing modes. The demonstration of spasing involves a resonant energy transfer from optical transitions in the gain medium to plasmon excitations in the metal, as well as a stimulated emission of surface plasmons due to the high local fields created by plasmons that excite the gain medium and stimulate more emission of selected plasmonic modes, leading to the required amplification [23].
The performance of a spaser is limited by the relaxation rates of the surface plasmons and the gain medium. To overcome the former limitation, plasmons in graphene provide a suitable alternative to those of metals due to the high mobility of its charge carriers, leading to tight confinement and a relatively long propagation distance, as well as the tunability of graphene's plasmons via electrostatic gating [24][25][26][27]. A spaser formed by doped graphene nanoribbons surrounded by semiconductor quantum dots (QDs) has been theoretically proposed [11]. This spaser has been shown to support a wide frequency generation region from terahertz to infrared, small plasmon damping and a low pumping threshold. Moreover, a tunable spherical graphene spaser that supports localized surface plasmon modes has been proposed. It was found that the spasing could occur when the quality factors of some localized modes becomes larger than some critical values given in terms of the Fermi energy of graphene [13]. Additionally, Apalkov et al. have proposed a novel nanospaser made of a graphene nanopatch and a quantum well cascade emitter [12]. With this spaser, optical fields have been generated, exhibiting a high nanolocalization and coherent generation of SPs in the graphene nanopatch.
However, the quantum well cascade suffers from nonradiative relaxation of electrons in the upper radiative state, due to thermally activated electron-longitudinal optical phonon scattering [28]. Interestingly, due to the discrete nature of their energy levels, it is possible to greatly increase the lifetime of the upper levels of QDs via a phonon bottleneck that suppresses electron-longitudinal optical phonon scattering [29]. Thus, a spaser that utilizes a quantum dots cascade emitter as its gain medium could be characterized by a low spasing threshold [30]. Moreover, based on selection rules, the optical transitions in QDs cascade emitters are only allowed along the direction of the QDs growth [31]. This can significantly simplify the calculations of a spaser.
Recently, it has been shown that a metal nanoparticles-graphene nanodisks-quantum dots hybrid system can support the ultrafast energy transfer between excitons and plasmons [32]. Specifically, within the near field approximation, having a relatively large size of metal nanoparticles (MNPs) of polarizability comparable to that of highly doped graphene nanodisks (GNDs) that support plasmons that are resonant with excitons in the QD, leads to a controllable and ultrafast energy transfer within the system [32]. Therefore, it is expected that by using an MNP-GND hybrid system as a resonator, one can enhance the performance of a spaser and exercise more control on its performance characteristics, since the plasmons in graphene can be launched and controlled effectively with resonant metal nanoantennas [32][33][34].
Clearly, there is continuous interest in developing novel spaser devices with lower spasing thresholds, a simplified but robust structure, high amplification and stable performance. In the present work, we propose a novel nanospaser composed of an MNP-GND hybrid plasmonic system as its resonator and a quantum dots cascade as its gain medium to take advantage of the tunability of long-lived plasmons in graphene [27] and the possibility to control the graphene plasmons by plasmons in noble metals [33], as well as the low spasing threshold with a quantum dots cascade emitter. The proposed spaser operates in the mid-infrared region that has wide chemical, biomedical and telecommunication applications [35]. The properties of the spaser as a plasmonic amplifier will be investigated using the density matrix theory with a quantized plasmonic field in the transient regime.

Theoretical Formalism
We aim to study the properties of the proposed spaser depicted in Figure 1 that utilizes the MNP-GND hybrid system as its resonator and a gallium nitride (GaN) QDs cascade stack as its gain medium. The GaN QDs stack is characterized by mechanical and thermal stability, as well as low sensitivity to ionizing radiation [36]. Moreover, GaN QDs cascade terahertz emitters, modeled for one period shown in Figure 1, where the bold numbers denote the size of the QDs that are sandwiched between Al 0.18 Ga 0.82 N layers, have demonstrated a relatively long lifetime of the upper level, leading to enhanced population inversion that supplies the spasing mode [37]. Let us first find the electric and magnetic field components of the plasmonic excitations near the GNDs arrays, taking into account the presence of the MNPs lattice. To this end, consider a two-dimensional periodic lattice of silver nanospheres of width W, on top of a periodic lattice of GNDs, located at z = 0, encompassing the |x| < W/2 region and centered at x = 0. Consider a transform limited pulse propagating along the z-direction, with its linear polarization along the interface between the GNDs layer and the dielectric (i.e., along the x-axis). The launching pulse induces the fields E ind and B ind that obey the Maxwell's equations, and due to the symmetry in the y-direction, the fields have transverse magnetic polarization (TM) with E x , E z and B y nonzero field components. Thus, by using Fourier transformation, one can write the field induced along the x-direction as: where q is the wavenumber of the SPs waves propagating along the interface. In the following, the dependence on time will be dropped in order to simplify the notation and focus on the spatial dependence. Using Maxwell's equations in the spectral representation, one can find: where k 2 q = q 2 − ω 2 c 2 , and is the dielectric constant of the embedding medium. Note thatẼ(ω) represents the envelope of the pulse in its angular frequency domain. J x (q) is the Fourier component of the current density in the plane z = 0, which can be written as: where Θ is the Heaviside step function and E(x, 0) is the total electric field at z = 0 given as the sum of the external and induced fields. σ GNDs and σ MNPs are the effective conductivity of GNDs and MNPs arrays, respectively. To model the optical response of nanoparticles arrays, it is fairly convenient to employ the effective medium approximation (EMA) that is only valid when the inner structure length scale is much smaller than the incident wavelength. Within EMA, the system comprising of interacting nanoparticles (NPs) is replaced by a homogeneous layer with an effective dielectric constant representing the overall response of the modeled system [38]. The effective dielectric constant for nanodisks and nanospheres arrays has been derived analytically by Genov et al. through considering the RLC model of NPs arrays, in which the negative permittivity of plasmonic material and the positive permittivity of the dielectric are represented by inductance R-L and capacitance C, respectively [39]. Based on this model for a plasmonic material, graphene (G) or metal (M), of a dielectric constant written as M, represents the loss parameter that is given as the ratio between the relaxation rate of plasmons γ M,G and the incident frequency ω, the effective dielectric constants of the periodic arrays of GNDs and MNPs embedded in a dielectric of d are given as: is the packing density of arrays given as a ratio between the diameter of particles and the interparticle distances, and with the approximation of plasma frequency (ω p ) M,G >> ω, M,G as (ω p /ω) 2 M,G . Substituting from Equation (2) into Equation (1), through the use of the Fourier transformation of J x (x), leads to [40]: Combining Equations (5) and (3) leads to the following equation for E x (x, 0) in all space [40]: where ξ q = 1 + (iσ GNDs k q /2ω 0 ) is the dielectric function of the GNDs layer. Clearly, E x (x , 0) involved in Equation (6), represents the field inside the metallic contact of width W, which can be written as a Fourier expansion: Thus, Equation (6) can be reduced to a system of algebraic equations to be solved numerically in a matrix form running the sum up to a judicially chosen upper value n = N. After determining the coefficients A n , the electric field induced along the x-direction near the graphene layer due to the inhomogeneity introduced by the metallic contact is obtained in terms of the dimensionless quantities u = qW and k u = k q W as [41]: where ζ σ = (σ MNPs /σ GNDs ) − 1, and It is clear that the parameter a contains a dependence on the ratio of the width of the metallic contact to the incident wavelength. Thus, a can be considered as a measure of raterdated effects. Specifically, the case of small a, i.e., u > a, corresponds to an evanescent wave confined near the graphene layer. Equation (8) gives the electric field component of the induced field along the interface. The transverse component can be found by using the following Maxwell's equations valid for TM polarization: Using Equation (2), and substituting from Equation (10b) into Equation (10a) while inserting J x given through Equation (3), we obtain the two-dimensional electric field induced along the z-direction due to the metallic contact: The corresponding transverse component of the magnetic field can be obtained using Equation (10a): The propagation of surface plasmons is governed by the dispersion relation derived by using the condition of the discontinuity of B y at z = 0. For a relatively large σ GNDs , i.e., 1−ξ u ξ u = −1, the discontinuity of B y is reduced to the following dispersion relation: where σ is the total conductivity of the GNDs and MNPs arrays at z = 0 and 1 ( 2 ) is the electric permittivitiy of the dielectric medium in z > 0 (z < 0). Note that the frequency of the surface plasmon polaritons for the system can be determined by solving for the pole of the dispersion relation that is given in terms of geometrical features and optical properties of the system.
To study the properties of the proposed MNPs-GNDs-QDs hybrid-system-based spaser, we focus on the z-component of the electric field, which describes the confined plasmonic field in the direction perpendicular to the interface, since the dipole transitions of the gain medium are only allowed along this direction [31]. Let the wavenumber-dependent amplitude in the field equation be: u can also be written in terms of the surface plasmon polaritons energyhω u where the field can be quantized by using the Brillouin expression for the field mean energy in the dispersive medium [42]: where µ and (ω, r) are the magnetic permeability and electric permittivity of the dispersive medium, respectively. The latter can be defined as: Solving Equation (15) yields the z-component of the quantized plasmonic field of mode ω u in the z < 0 region of the dielectric medium (medium 2) where the gain medium is placed: whereâ † u andâ u are the raising and lowering operators of the SPs, respectively. |z| represents the penetration length of plasmonic field that can be adjusted to be the distance between the GNDs layer and the active medium. The area of the GNDs layer is denoted by S 0 , and du/dω u included in Equation (17) can be obtained by using the dispersion relation, Equation (13): where for GNDs layer of thickness t and Fermi energy E F : The model of the proposed mid-infrared MNPs-GNDs-QDs hybrid-system-based spaser is shown in Figure 2. Following injection current pumping I in , the electrons injected into the upper sub-band (level 3) of the GaN QDs cascade emitter undergo transitions to the lower level (level 2) stimulated by resonant surface plasmons. Subsequently, these electrons relax non-radiatively to the lowest energy level (level 1) to eventually tunnel via the chirped superlattice to the adjacent QD with rate γ out [31]. In the presence of a SPs mode resonant with the dipole transition energyhω 32 of the GaN QDs cascade gain medium, energy transfer will occur through the coupling between the dipole transitions of the gain medium and plasmon excitations, leading to the amplification of this mode of SPs. Note, that the number of periods of the QDs cascade emitter that can be included in the energy transfer depends on the penetration length of the plasmonic fields.The dynamics of the system can be investigated by solving the Lindblad master equation through the use of the following Hamiltonian for the u mode: where Ω 32 is the coupling parameter illustrated in Figure 2 and defined as d 32 E uz /h. d 32 is the dipole moment of transition |3 ↔ |2 . Since the dipole transition of the GaN QDs cascade is along the z-direction, we consider only the z-component of the plasmonic fields. The two first terms of the Hamiltonian refer to the plasmonic resonator with mode u and the gain medium, respectively. The third term describes the energy transfer between the optical transitions in the gain medium and plasmon excitations. Specifically, |3 2|â u corresponds to the absorption of a SP in order to make the resonant transition from |2 to 3|. On the other hand,â † u |2 3| represents the emission of a SP to make the transition from |3 to 2|. The Liouvillian of the system that describes the decay channels is given by [43]: where γ out is the rate of tunneling given in terms of the dimension of QD (L QD ), and its effective mass (m QD ) as [37]: L b and m b in Equation (21) are the thickness of the barrier and its effective mass, respectively. The height of the potential barrier is denoted by H b . With the above consideration, the density matrix equations of the system are given as [44]: where a u = a 0u e −iω sp t , ρ 32 = ρ 032 e −iω sp t , ρ 31 = ρ 031 e −i(ω sp −iγ 21 )t and ρ 21 = ρ 021 e −γ 21 t , with the slowly varying amplitudes a 0u , ρ 032 , ρ 031 and ρ 021 . γ sp is the damping rate of the SPs. Note that the rate of the injection current, i.e., I in /e, where e is the charge of an electron, affects the coherence and population terms that are related to the upper level [31]. γ ij are the damping rates corresponding to the relaxation times (τ ij ) shown in Figure 2.

Analysis of the Confined Plasmonic Fields near the GNDs Layer
This paper aims to study the optical properties of a spaser consisting of a GNDs lattice excited by means of an MNPs one as its resonator and a GaN QDs cascade emitter as its gain medium. Thus, since the QDs emitter has optical transitions only along the direction of the QDs growth, i.e., the z-direction, we focus on the z-component of the electric field given by Equation (11), that describes the confined plasmonic fields near the GNDs layer. Figure 3 shows the confined plasmonic field over a few tens of nanometers, corresponding to the penetration length of the plasmonic field, where the sum in Equation (11) is running up to N = 20. Clearly, with a 10 fs-sech launching pulse of time dependence given in terms of the time duration of the pulse ∆τ as [45], We get a relatively strong confined plasmonic field of MV/m strength that is significantly enhanced because of the ultrashort pulsed excitation. Moreover, the confined plasmonic fields are enhanced for a relatively small ratio between the width of the metallic contact and the wavelength of the incident radiation, as well as a small dielectric constant of the medium where the two lattices are embedded. This can be attributed to the corresponding limited retardation effects that are responsible for the decay of the confined surface plasmons as an electromagnetic radiation [46]. Additionally, we examine the effect of the packing density of the GNDs and Ag NPs lattices on the confinement of the plasmonic fields near the GNDs layer, as illustrated in Figure 4. We observe that a relatively strong confined plasmonic field is obtained when the packing density of the Ag NPs lattice is larger than that of the GNDs one. The confined plasmonic field decreases as the difference of packing density between the two lattices decreases. We observe a significant decrease of the confined plasmonic field for the case of δ GNDs > δ AgNPs . This is reasonable, since the plasmonic fields are enhanced for a large δ AgNPs /δ GNDs corresponding to a highly conductive stripe. Note, that the conductivities of the two lattices are increased for a large packing density [39].  Figure 5 shows the effect of the thickness of the GNDs (t) and their doping levels, which are defined by the Fermi energy (E F ), on the confined plasmonic fields. Clearly, the latter are enhanced for a relatively small thickness of GNDs with large doping levels that lead to a relatively large bulk plasma frequency for graphene, ω 2 p,G = (e 2 E F )/(πh 2 0 t). It can be seen that the penetration length is not affected by the thickness of the GNDs and their doping levels. Additionally, the confined plasmonic field near the GNDs layer is sensitive to the thickness of the GNDs more than to the doping level.

Analysis of the Dynamics of the Proposed MNPs-GNDs-QDs Hybrid System-Based Spaser
Unlike the conventional lasers and amplifiers in quantum electronics, spasers have an inherent feedback provided by plasmonic structures that typically cannot be removed [2]. Thus, a spaser will develop the accumulation of a large number of coherent SPs, leading to a continuous wave (CW) regime where the gain compensates exactly for the losses with zero net amplification. One way to operate a spaser as a plasmonic amplifier is to consider the transient regime, recalling the fact that the establishment of the CW regime requires a relatively long time [20]. Thus, we consider the transient amplification of SPs during hundreds of femtoseconds, a time interval smaller than the relaxation times of the gain medium and SPs. The transient dynamics of the spaser can be directly obtained by the numerical solution of the density matrix equations (Equation (22)).
At first, we solve numerically for the pole of the dispersion relation of SPs for the planar structure given by Equation (13). The parameters of the system are adjusted to get the SPs mode energy to be resonant with that of the GaN QDs cascade emitter. Specifically, for a GNDs monolayer of a thickness of 0.34 nm, Fermi energy of 0.7 eV and mobility of 10 4 cm 2 /V s with the GNDs lattice (Ag NPs lattice) embedded in a dielectric medium of dielectric constant 3 (12) with packing density 100 (150), we get an SPs mode energy of 41.5 meV, which is resonant with the laser transition of GaN QD terahertz cascade laser [37]. The GaN QD is modeled as a three-level system of cascade configuration, where the energy spacing between the two lower levels is resonant with longitudinal optical phonon scattering, leading to the fast depopulation of the lower state of the laser transition, which makes this quantum emitter attractive as a gain medium. The relaxation rates of the gain medium are set to be γ 31 = 7.14 × 10 10 s −1 , γ 32 = 7.7 × 10 10 s −1 and γ 21 = 2.2 × 10 12 s −1 [37], whereas that of the SPs corresponding to E F = 0.7 eV is γ sp = 7.12 × 10 11 s −1 . The rate of injection is set to be 10 10 s −1 and the rate of tunneling is calculated with a GaN QD of dimension L QD = 4 nm and an effective mass of 0.2m 0 with an AlGaN barrier of L b = 3 nm and m b = 0.3m 0 , correspondingly [47]. The dipole moment of the gain medium is taken as 2.5e nm [47] and the area of the GNDs layer is set to be 75 µm 2 .
The plasmonic components are launched by a 1ps-sech pulse that is shorter than the relaxation times of the gain medium, so that the slowly varying approximation is valid [48] and our calculations for the energy transfer between the plasmonic excitations and the optical transitions in the gain medium are simplified. The dynamics of energy transfer between the SPs and the optical transitions of the gain medium are shown in Figure 6 by numerically solving for the time-dependent density matrix equations (Equation (22)) with initial conditions ρ 33 (0) = a * 0u (0) = 1; that is, electrons are initially in the upper level (level 3) and a resonant SP mode is present to induce the transition to the lower level (level 2). The penetration length of the plasmonic field is set to be equivalent to one period of the GaN QDs cascade emitter, i.e., z = 17 nm [37].
It can be seen that transient amplification of the SPs mode of energy 41.5 meV is obtained over hundreds of femtoseconds before they are damped depending on the relaxation rates of the SPs and the gain medium. It is remarkable that at the instant the number of SPs is maximized, the population inversion reaches its minimum, implying that the energy of optical transition is completely transferred to the SPs to be amplified. Interestingly, the amplification of SPs is significantly affected by the width of the metallic contact, emphasizing the important role of the metal in the energy transfer as shown in [32]. Specifically, for a large width of the metallic contact, the amplification of SPs is decreased and the time corresponding to the maximum number of SPs is displaced to later values. This can be attributed to the retardation effects that dominate as the width of the metallic contact approaches the wavelength of the incident radiation, leading to the decay of the SPs [49].  Figure 7 shows the effect of the time duration of the normalized sech-pulse used for launching SPs on the dynamics of the SPs mode of energy 41.5 meV. It can be seen that the amplification of SPs is enhanced for a small time duration of the launching pulse that is associated with a relatively large intensity. This enhancement is due to the large local field enhancement induced near the surface of plasmonic materials with pulsed excitation [50]. Interestingly, the transient amplification occurs over the hundreds of femtoseconds, a time duration smaller than the relaxation time of the upper level, so that the spontaneous emission noise can be avoided. Moreover, the spasing threshold of our proposed spaser consisting of a GNDs-MNPs plasmonic system and a GaN QDs cascade emitter that are separated by a distance l can be easily estimated by a comparison between losses and gains: h 2 |Ω 32 | 2 e −2ql ≥ Γ SP Γ 32 . Interestingly, the spasing threshold of the proposed spaser is smaller than that demonstrated for the graphene spaser proposed by Apalkov et al. [12], due to the relatively small inter-sub-band-polarization relaxation width of the gain medium (Γ 32 = 50 µeV ) and the relatively low damping width of the SPs (Γ SP = 470 µeV). Specifically, the spasing condition is given in terms of the time relaxation of plasmons by [12]: where q is the wave number of plasmons and Γ 32 is the relaxation width of the QD. The distance between the graphene and the QD is denoted by l. k 0 is the effective Fermi wave vector that can be easily calculated by knowing the energy of plasmons (hω q ), which is chosen to be at resonance with the QD transition frequency in order to facilitate the energy transfer between plasmons in graphene and excitons in QD, which is essential for spasing. Therefore, the operating frequency of our proposed spaser is limited by the transition frequency of the QD. In fact, our focus is on enhancing the amplification and not on introducing broad tunability. For the following parameters: Γ 32 = 50 µeV, E F = 1 eV, l = 17 nm, k 0 = 1.1 nm −1 ,hω q = 41.5 meV, d 32 = 2.5 enm, we get an extremely short minimum relaxation time of graphene plasmons, i.e., τ min = 435, which corresponds to a quality of factor Q min = (ω q τ min )/2 = 0.013. The obtained quality factor is lower than the one obtained by Apalkov et al. [12], implying that the spasing in our proposed hybrid nano-system is relatively easy to achieve. Clearly, the transient regime represents the region where the spaser can act as a plasmonic amplifier with a gain exceeding the losses, achieving an amplification of a selected mode of SPs that can be controlled by the width of metallic contact and the time duration of the launching pulse, leading to an enhancement of plasmonics.

Conclusions
In this paper, we have studied the properties of the proposed MNPs-GNDs-QDs hybrid-systembased spaser. Firstly, we investigated the confined plasmonic fields near the interface of the GNDs lattice with a metallic contact, represented by an Ag NPs lattice, induced by the launching pulse that has its linear polarization along the interface between the GNDs layer and the dielectric environment. We have found that the plasmonic field can be enhanced for a small width of the metallic contact and a small dielectric constant of the surrounding medium, because then the retardation effects are limited. Moreover, the confined plasmonic field is enhanced for a small thickness of the GNDs and large doping levels that lead to a relatively large density of charge carriers. A large packing density in the MNPs lattice compared to that of the GNDs can also significantly enhance the confined plasmonic fields.
Subsequently, we have studied the dynamics of the proposed spaser with a GaN QDs cascade emitter chosen to be the gain medium. A novel MNPs-GNDs-QDs hybrid-system-based spaser is obtained with a transient amplification of the SPs mode, providing a plasmonic amplifier. The number of SPs that determines the intensity of the obtained plasmonic field is significantly affected by the width of the metallic contact, emphasizing the important role of the MNPs lattice in the spasing process. In particular, the presence of an MNPs lattice provides more options for controlling and enhancing the amplification of the plasmonic fields.
Interestingly, the various factors in our scheme, which was identified to affect the amplification of the plasmonic fields, act cooperatively and cumulatively. They are within reach of present-day materials technology, and various combinations of them, within the theoretically identified limits, could result in plasmonic amplifiers with a low spasing threshold and significant amplification. This feature adds a degree of flexibility in the realization of the herein proposed plasmonic amplifier. Therefore, the proposed controllable SPs amplification scheme could prove useful for various interesting applications, such as stand-alone components or gain sections integrated with other plasmonic elements to compensate for losses and improve the performance of plasmonic devices.