## Abstract

We show that in plasmonic systems, exact loss compensation can be achieved with the help of spasers pumped over a wide range of pumping values both below and above the spasing threshold. We demonstrate that the difference between spaser operation below and above the spasing threshold vanishes, when the spaser is synchronized by an external field. As the spasing threshold loses its significance, a new pumping threshold, the threshold of loss compensation, arises. Below this threshold, which is smaller than the spasing threshold, compensation is impossible at any frequency of the external field.

© 2013 OSA

## 1. Introduction

The problem of loss compensation in plasmonic systems, such as plasmonic composites, metamaterials, and plasmonic transmittance lines, has become increasingly pressing as these systems gain attention owing to their prospects for attractive applications [1–10]. These applications cannot be realized due to the high level of Joule losses in metal components of plasmonic composites. To overcome this problem in metamaterials, it has been suggested to incorporate gain media (atoms, molecules or quantum dots) into these systems [1–7]. The general goal of such compensation is to construct a gain plasmonic material in which the electromagnetic response mimics the response of an ordinary composite without loss. In other words, such a material should be characterized by a dielectric permittivity. In particular, the wave generated by an external harmonic field in this metamaterial should have the same frequency as the external field.

The goal of loss compensation in metamaterials can be achieved with the help of spasers [7, 11–15]. Without external field a spaser operating above the pumping threshold (the spasing regime) is an autonomic (self-oscillating) system exhibiting undamped harmonic oscillations. These oscillations are characterized by their own frequency ${\omega}_{a}$ and the amplitude [16–18]. An external harmonic field with frequency $\nu $ can synchronize the spaser forcing it to operate at frequency $\nu $ [12]. However, in this case, an active metamaterial cannot be described by a dielectric function because the synchronization is possible only if the amplitude of the external field exceeds a threshold value ${E}_{synch}\left({\Delta}_{E}\right)$, which depends on the frequency detuning ${\Delta}_{E}=\nu -{\omega}_{a}$. The domain in which $E({\Delta}_{E})>{E}_{synch}({\Delta}_{E})$ is called the Arnold tongue [19]. Synchronization inside the Arnold tongue known from the theory of self-sustained oscillations was also observed for lasers [20]. In addition, for the field not much stronger than ${E}_{synch}({\Delta}_{E})$, the response of a synchronized spaser depends only weakly on the external field [13].

Generally, inside the Arnold tongue the dipole moment of the spaser has a nonzero imaginary part. This means that the spaser operates either as a gain or loss inclusion. Exact compensation occurs when the amplitude of the wave travelling in an active metamaterial is equal to a special value ${E}_{com}({\Delta}_{E})$ which depends on the frequency detuning and pump [13]. Interestingly, inside the Arnold tongue, the amplitude of the travelling wave reaches ${E}_{com}({\Delta}_{E})$ automatically. If the amplitude of the travelling wave is greater than ${E}_{com}({\Delta}_{E})$, the energy is transferred from the wave to spasers and the amplitude of the travelling wave drops. In the opposite case, the energy is transferred from the spasers to the wave and the wave amplitude grows. Thus, eventually the wave travels with the amplitude independent on the external field [14].

Below the pumping threshold, *D _{th}*, spasers are always synchronized by an external harmonic field, and therefore the system can be characterized by an effective dielectric permittivity. In this regime, spasers seem to be good candidates for loss compensation if not for the apparent energy shortage. Indeed, self-oscillations of a spaser above the pumping threshold are due to the energy delivered by pumping which cannot be smaller than Joule losses. An external field breaks this balance causing additional losses in the plasmonic nanoparticle (NP). Since below the threshold, the pumping does not compensate for losses even in the absence of an external field, it may seem that in order to achieve compensation, one should pump the driven spaser above the threshold. These arguments qualitatively agree with the result obtained in [8] where it was shown that loss compensation occurs simultaneously with the start of spasing. On the other hand, in the absence of the interaction between plasmonic and active media [1], as well as in the linear regime [7, 11], below-threshold compensation was predicted.

In this paper, within the framework of the single-particle model, we demonstrate that the exact compensation of Joule losses can be achieved with spasers in both below and above pumping threshold regimes. We show that the synchronizing external field destroys self-oscillations of spasers transforming them into a nonlinear oscillator with a new pumping threshold ${D}_{comp}\left(E\right)<{D}_{th}$. When the pumping rate, *D*_{0}, exceeds *D _{comp}*, there is a line of exact loss compensation, $E={E}_{comp}({\omega}_{a}-{\omega}_{E},{D}_{0})$.

## 2. Equations of motions for a free spaser and a spaser driven by external optical wave

We consider a simplest model of spaser as a two-level system (TLS) of size ${r}_{TLS}$ placed near a plasmonic spherical NP of size ${r}_{NP}$ [16, 17]. The energy from the pumped TLS is non-radiatively transferred to the NP exciting surface plasmons (SPs). At the frequency of the SP resonance, the dynamics of the NP polarization is governed by the oscillator equation:

The quantization of this oscillator can be performed in an ordinary way by introducing the Bose operators ${\widehat{\tilde{a}}}^{\u2020}(t)$ and $\widehat{\tilde{a}}(t)$for the creation and annihilation of the dipole SP [21, 22]:The corresponding Hamiltonian is:The quantum description of a TLS is done via the transition operator $\widehat{\tilde{\sigma}}=|g\u3009\u3008e|$ between ground $|g\u3009$ and excited $|e\u3009$ states of the TLS, so that the operator for the dipole moment of the TLS is represented as

The commutation relations for operators $\widehat{\tilde{a}}(t)$ and $\widehat{\tilde{\sigma}}(t)$ are standard: $[\widehat{\tilde{a}}(t),{\widehat{\tilde{a}}}^{+}(t)]=1$ and $[{\widehat{\sigma}}^{\u2020},\widehat{\sigma}]=\text{}\widehat{D}$, where the operator $\widehat{D}={\widehat{n}}_{e}-{\widehat{n}}_{g}$ describes the population inversion of the ground ${\widehat{n}}_{g}=|g\u3009\u3008g|$ and excited states ${\widehat{n}}_{e}=|e\u3009\u3008e|$, ${\widehat{n}}_{g}+{\widehat{n}}_{e}=\widehat{1}$, of the TLS.

We describe the dynamics of the free spaser by the model Hamiltonian [16, 17, 24]

where the operator $\widehat{V}=-{\widehat{d}}_{NP}{\widehat{E}}_{TLS}$ is responsible for the dipole-dipole interaction between the TLS and the NP. Taking into account that ${\widehat{E}}_{TLS}=-{\widehat{\mu}}_{TLS}{r}^{-3}+3({\widehat{\mu}}_{TLS}\cdot r)r{r}^{-5}$, we obtain*r*

_{0}is the distance between the TLS and the NP, and ${e}_{r}=r/r$ is the unitary vector. The last term in the Hamiltonian (9) is responsible for relaxation and pumping processes.

Assuming that ${\omega}_{TLS}$ is close to the frequency of the plasmonic resonance, ${\omega}_{SP}\approx {\omega}_{TLS}$, we can use the approximation of the rotating wave [25] by looking for the solutions in the form $\widehat{\tilde{a}}(t)\equiv \widehat{a}(t)\mathrm{exp}(-i{\omega}_{a}t)$ and $\widehat{\tilde{\sigma}}(t)\equiv \widehat{\sigma}(t)\mathrm{exp}(-i{\omega}_{a}t)$, where $\widehat{a}(t)$, $\widehat{\sigma}(t)$ are slow varying in time operators and ${\omega}_{a}$ is the autonomous frequency of the spaser which we seek. Disregarding fast-oscillating terms proportional to $\mathrm{exp}(\pm 2i{\omega}_{a}t)$, the interaction operator $\widehat{V}$may be written the form of the Jaynes–Cummings Hamiltonian [21]:

Using Hamiltonian (9) and the commutation relations for operators $\widehat{a}(t)$ and $\widehat{\sigma}(t)$ we obtain the Heisenberg equations of motion for the operators $\widehat{a}(t)$, $\widehat{\sigma}(t)$, and $\widehat{D}(t)$ [17]

*c*-numbers. In this case, the Hermitian conjugation turns into the complex conjugation [17, 28-30]. Note that the quantity $D(t)$ that describes the difference in populations of excited and ground states of the TLS is a real quantity because the corresponding operator is Hermitian.

The system of Eqs. (12)-(14) has stationary solutions, which depend on the pumping level ${D}_{0}$. For ${D}_{0}$ smaller than the threshold value,

The Hamiltonian of a spaser driven by an external field of optical wave, which is assumed to be classical, ${E}_{OW}(t)=E\mathrm{cos}\nu t$, may be written in the form (see for details [12, 13, 31, 32]):

*H*is given by Eq. (9), ${\Omega}_{1}=-{\mu}_{NP}E/\hslash $ and ${\Omega}_{2}=-{\mu}_{TLS}E/\hslash $ are the coupling constants of the external field interaction with the NP and the TLS, respectively.

As above, the equations of motion for slow amplitudes $\widehat{a}$, $\widehat{\sigma}$, and $\widehat{D}$ can be obtained as:

## 3. Loss compensation

Since the energy flows through a spaser, it can be considered as an open system. The flow starts at pumping, which causes the population inversion of the TLS, then the nonradiative transition of the TLS excites SPs at the NP, and the Joule losses of SPs at the NP finalize the flow. In the absence of an external field, the energy of pumping is consumed by SP excitations. At fixed pumping, an increase of the SP amplitude is limited by Joule losses at the NP. The self-oscillating state of a spaser (spasing) occurs at exact compensation of losses with pumping [16, 27]. If losses exceed the energy supplied by pumping, the stationary amplitudes of oscillations are equal to zero. The maximum value of pumping, below which there is no spasing, is referred to as the threshold pumping.

The external field performs work on the dipole moments of the TLS and the NP. Thus, subjecting a spaser to an external field leads to additional channels of energy flow; namely, the energy flows from the field to the TLS and to the NP. These energy flows may close up via the interaction of the TLS with the NP and interfere with each other and with the primary energy flow from pumping to the TLS and then to Joule losses at the NP. Rather complicated dynamics of these flows results in non-zero oscillations of the below-threshold spaser, as shown in Fig. 1(a).

For exact loss compensation, the work performed by the field on the spaser should be equal to zero. The time averaged work performed by the external field **E** on the spaser is $~\left({{\mu}^{\u2033}}_{TLS}+{{\mu}^{\u2033}}_{NP}\right)\cdot {E}^{*}$ [33]. Since$\mu \parallel E$, for non-zero dipole moments, this expression turns to zero when $\left({{\mu}^{\u2033}}_{TLS}+{{\mu}^{\u2033}}_{NP}\right)=0$ or the phase difference between the field and spaser dipole oscillations is equal to *π*. If the phase difference is greater than *π*, the wave is amplified (negative work). If the phase is smaller than *π*, the wave attenuates (positive work). In other words, for exact loss compensation the sum of the imaginary parts of the dipole moments of the TLS and the NP should be equal to zero. It has been shown [13], that above the spasing threshold, the imaginary part of the dipole moment of a synchronized spaser is equal to zero on the compensation curve $E={E}_{comp}\left({\Delta}_{E},{D}_{0}\right)$ (see Fig. 2). For the fields below this curve, the energy pumped into the system exceeds losses; for the fields above the compensation curve, the system becomes lossy. Note, that exact compensation can occur only for the fields smaller than${E}_{\mathrm{max}}({D}_{0})=\mathrm{max}\{E={E}_{comp}\left({\Delta}_{E},{D}_{0}\right)\}$.

The possibility of loss compensation below the threshold can be illustrated in a simple limiting case of a vanishing interaction between the TLS and the NP. In this case, we deal with an heterogeneous system in which plasmonic and gain ingredients do not interact. The pumping threshold depends on the coupling constant ${\Omega}_{R}$ (see Eq. (15)). For ${\Omega}_{R}\to 0$ the threshold tends to infinity. Such a spaser cannot spase. In such a TLS-NP couple, the TLS and the NP are independent and play opposite roles when they interact with the external wave. While the TLS amplifies the wave’s field, losses in the NP weaken it. In the limit ${\Omega}_{R}=0$, Eqs. (18)-(20) are reduced to:

In this case, the contradiction with the energy shortage is resolved if we notice that the pumping energy should not compensate for Joule losses in the NP caused by the TLS field. The energy of pumping is transferred to the field by the TLS. At the same time the NP absorbs the field energy. If Eqs. (27) and (28) are satisfied, then the total energy transfer to the system is zero. Note that since the pumping threshold in this toy model (${\Omega}_{R}=0$) is infinity, loss compensation occurs below the pumping threshold.

In a sense, the spaser with ${\Omega}_{R}=0$ is similar to the system suggested in [1] in which gain and plasmonic media are confined to different layers of one-dimensional photonic crystal. Thus, there is no direct near field interaction between the media and, as a consequence, no spasing. The necessary condition for lasing arises when the energy delivered by pumping exceeds loss at plasmonic layers.

Having ${\Omega}_{R}=0$ is sufficient but not necessary for loss compensation below the spasing threshold. As it is shown in Fig. 3(b), in the case ${\Omega}_{R}\ne 0$ below the spasing threshold, there are frequencies for which the spaser’s dipole moment is zero and the energy is not transferred to or from the system. There are also frequencies for which the imaginary part of the NP dipole moment is negative, so that the NP releases energy to the wave. The reason of such unexpected behavior is interference of energy fluxes similar to the case of the Fano resonance [34].

The analogy with the Fano resonance clearly manifests itself for low pumping (${D}_{0}~-1$) when there is still no loss compensation (Fig. 3a). In this case, the spaser response to the external field is practically linear and we can consider the system as two coupled resonators. The first resonance with a low *Q*-factor is the SP resonance at the NP, the other is the high-*Q* resonant transition of the TLS. The response of such a system on the external force has the shape of the Fano curve [35]

*q*describes asymmetry of the line.

In the case of the spaser with a low level of pumping, ${D}_{0}~-1$, at ${\tau}_{a}\ll {\tau}_{\sigma}$, the frequency dependence of the imaginary part of the total dipole moment of the spaser can be obtained as

*q*is equal to ${\Omega}_{R}\left|{D}_{0}\right|{\mu}_{TLS}/{\mu}_{NP}$, which strongly depends on the interaction of the NP and the TLS with the external field.

As one can see in Fig. 3(a), for low pumping, $\mathrm{Im}{\mu}_{tot}>0$, so that the system is lossy for any frequency. At the same time, the minimum of losses is significantly lower than that for an isolated NP. When pumping increases, the resonant line still resembles the Fano resonance line but the minimum value of $\mathrm{Im}{\mu}_{tot}$ becomes negative (Fig. 3(b)), so that for a range of frequencies the spaser releases energy. The compensation curve, $E={E}_{comp}\left(\Delta ,{D}_{0}\right)$, also exists in this case (see Fig. 4). As in the case of the above-threshold spaser, this curve lies below some value of the electric field ${E}_{\mathrm{max}}\left({D}_{0}\right)$. The main difference between spasers in the above and below threshold pumping regimes is that in the latter case, the Arnold tongue occupies the whole half-plane, i.e., the spaser is always synchronized. Below the threshold, the synchronized spaser is not a self-oscillating system but it is rather a non-linear oscillator. Thus, again we arrive at the contradiction of the shortage of the pumping energy [8, 36, 37].

This contradiction is resolved, if one notices that for the fields in which loss compensation exists ($E<{E}_{\mathrm{max}}\left({D}_{0}\right)$) the absolute values of the dipole moments of the above-threshold spaser is greater than that for the below-threshold spaser (see Fig. 1). As a result, below the pumping threshold, losses of the TLS field in the NP are smaller than that above the threshold. Thus, the energy of pumping is sufficient to compensate for these losses as well as for additional Joule losses due to the external field. If $E>{E}_{\mathrm{max}}\left({D}_{0}\right)$, the dipole moments of a spaser above and below the threshold are nearly the same. As a consequence, for these fields pumping greater than the threshold is needed for loss compensation. Therefore, ${E}_{\mathrm{max}}\left({D}_{0}\right)$ is a critical amplitude of the external field, for which the dipole moments of the below-threshold spaser become comparable with those of the *free* above-threshold spaser. If the external field exceeds ${E}_{\mathrm{max}}\left({D}_{0}\right)$, losses in the NP of the field generated by the TLS exceed the energy supplied by pumping and the energy of pumping below the threshold becomes insufficient for loss compensation.

As pumping decreases, the compensation curve monotonically shrinks toward the line $E=0$ disappearing at some level of pumping (Fig. 5). The dependence ${E}_{\mathrm{max}}\left({D}_{0}\right)$, shown in Fig. 6, is characterized by the new pumping threshold, ${D}_{comp}$, below which no compensation is possible. As shown in Fig. 7, ${D}_{comp}$ is never greater than the spasing threshold.

Above and below the threshold, the compensation line originates from the points $E=0,$ $\Delta =0$ and $E=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\Delta >0$, respectively (see Fig. 8). The pumping ${D}_{0}={D}_{th}$ is the smallest pumping at which compensation at zero frequency detuning is possible. This is the case considered in [8, 38]. For ${D}_{comp}\le {D}_{0}<{D}_{th}$, compensation can only be achieved for $\Delta >0$.

## 4. Conclusion and discussion

Synchronizing a spaser to an external field leads to the destruction of the spaser as a self-oscillating system and to its transformation into an active nonlinear driven oscillator. Indeed, after synchronization, the qualitative differences between below- and above-threshold spasers disappear; the remaining differences are merely quantitative (i.e., the ratio of dipole moments is of the order of the ratio of pumping values).

The threshold pumping ${D}_{th}$ of a free spaser loses its significance. In particular, the exact compensation of losses by a spaser is realized over a wide range of pumping values above and below ${D}_{th}$. A new pumping threshold, the threshold of compensation ${D}_{comp}$, arises. Below this threshold compensation is impossible at any frequency of the external field.

Only in the absence of the external field and at zero frequency detuning, the pumping threshold may coincide with ${D}_{comp}$. This coincides with conclusion of [8].

Below the threshold, loss compensation by a spaser is only possible if the frequency of the external field is greater than the transition frequency for the TLS. In this connection, it is interesting to consider the results of [7], in which it is shown numerically that below-threshold compensation is achieved at frequencies below the transition frequency. Nevertheless, if we take into account the Lorentz shift of the resonant frequency $\Delta \omega \sim \sqrt{4\pi N{e}^{2}/3m}$ (see, e.g., [39]) appearing due to the difference between local and average fields (in [7], the concentration of active molecules is $N\sim 6\cdot {10}^{18}c{m}^{-3}$), we obtain $\Delta \omega \sim 3\cdot {10}^{13}{s}^{-1}$. This is in good agreement with our conclusion that below-threshold compensation is possible for positive detuning only. Thus, there is no contradiction between the results of [7] and [8].

Recently, in [9, 10], it was demonstrated that in injection pumped spasers, the threshold current density may reach unsustainable values. On this basis, concerns of possibility of loss compensation spasers were raised. As we discussed in [40], for dye-based spasers pumped optically, which we consider in the current paper, the situation is different. For our system, both the pump power and the electric field strength for the pumping wave have reasonable values of ${10}^{-5}\text{\hspace{0.17em}}W$ and 100 *V/m*, respectively.

In conclusion, when pumping is below the spasing threshold, spasers may be used for compensation for Joule losses over a range of frequencies once the necessary pumping exceeds a new compensation threshold ${D}_{comp}$, which is smaller than the spasing threshold.

## Acknowledgments

This work was supported by RFBR Grants Nos. 12-02-01093, 13-02-00407, by a PSC-CUNY grant, and by Dynasty Foundation.

## References and links

**1. **S. Anantha Ramakrishna and J. Pendry, “Removal of absorption and increase in resolution in a near-field lens via optical gain,” Phys. Rev. B **67**(20), 201101 (2003). [CrossRef]

**2. **M. A. Noginov, G. Zhu, M. Bahoura, J. Adegoke, C. Small, B. A. Ritzo, V. P. Drachev, and V. M. Shalaev, “The effect of gain and absorption on surface plasmons in metal nanoparticles,” Appl. Phys. B **86**(3), 455–460 (2007). [CrossRef]

**3. **A. K. Popov and V. M. Shalaev, “Compensating losses in negative-index metamaterials by optical parametric amplification,” Opt. Lett. **31**(14), 2169–2171 (2006). [CrossRef] [PubMed]

**4. **M. A. Noginov, V. A. Podolskiy, G. Zhu, M. Mayy, M. Bahoura, J. A. Adegoke, B. A. Ritzo, and K. Reynolds, “Compensation of loss in propagating surface plasmon polariton by gain in adjacent dielectric medium,” Opt. Express **16**(2), 1385–1392 (2008). [CrossRef] [PubMed]

**5. **A. N. Lagarkov, A. K. Sarychev, V. N. Kissel, and G. Tartakovsky, “Superresolution and enhancement in metamaterials,” Phys. Usp. **52**(9), 959–967 (2009). [CrossRef]

**6. **A. Fang, T. Koschny, and C. M. Soukoulis, “Lasing in metamaterial nanostructures,” J. Opt. **12**(2), 024013 (2010). [CrossRef]

**7. **S. Wuestner, A. Pusch, K. L. Tsakmakidis, J. M. Hamm, and O. Hess, “Overcoming losses with gain in a negative refractive index metamaterial,” Phys. Rev. Lett. **105**(12), 127401 (2010). [CrossRef] [PubMed]

**8. **M. I. Stockman, “Spaser action, loss compensation, and stability in plasmonic systems with gain,” Phys. Rev. Lett. **106**(15), 156802 (2011). [CrossRef] [PubMed]

**9. **J. B. Khurgin and G. Sun, “Injection pumped single mode surface plasmon generators: threshold, linewidth, and coherence,” Opt. Express **20**(14), 15309–15325 (2012). [CrossRef] [PubMed]

**10. **J. B. Khurgin and G. Sun, “Practicality of compensating the loss in the plasmonic waveguides using semiconductor gain medium,” Appl. Phys. Lett. **100**(1), 011105 (2012). [CrossRef]

**11. **M. Wegener, J. L. García-Pomar, C. M. Soukoulis, N. Meinzer, M. Ruther, and S. Linden, “Toy model for plasmonic metamaterial resonances coupled to two-level system gain,” Opt. Express **16**(24), 19785–19798 (2008). [CrossRef] [PubMed]

**12. **E. S. Andrianov, A. A. Pukhov, A. V. Dorofeenko, A. P. Vinogradov, and A. A. Lisyansky, “Forced synchronization of spaser by an external optical wave,” Opt. Express **19**(25), 24849–24857 (2011). [CrossRef] [PubMed]

**13. **E. S. Andrianov, A. A. Pukhov, A. V. Dorofeenko, A. P. Vinogradov, and A. A. Lisyansky, “Dipole response of spaser on an external optical wave,” Opt. Lett. **36**(21), 4302–4304 (2011). [CrossRef] [PubMed]

**14. **A. P. Vinogradov, E. S. Andrianov, A. A. Pukhov, A. V. Dorofeenko, and A. A. Lisyansky, “Quantum plasmonics of metamaterials: loss compensation using spasers,” Physics-Uspekhi **55**(10), 1046–1053 (2012). [CrossRef]

**15. **A. P. Vinogradov, E. S. Andrianov, A. A. Pukhov, A. V. Dorofeenko, and A. A. Lisyansky, “Using spaser-effect for loss compensation in metamaterials,” AIP Conf. Proc. **1475**, 92–94 (2012). [CrossRef]

**16. **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. **90**(2), 027402 (2003). [CrossRef] [PubMed]

**17. **I. E. Protsenko, A. V. Uskov, O. A. Zaimidoroga, V. N. Samoilov, and E. P. O’Reilly, “Dipole nanolaser,” Phys. Rev. A **71**(6), 063812 (2005). [CrossRef]

**18. **M. I. Stockman, “The spaser as a nanoscale quantum generator and ultrafast amplifier,” J. Opt. **12**(2), 024004 (2010). [CrossRef]

**19. **A. Pikovsky, M. Rosenblum, and J. Kurths, *Synchronization. A Universal Concept in Nonlinear Sciences* (Cambridge University Press, 2001).

**20. **A. E. Siegman, *Lasers* (University Science Books, Mill Valley, California, 1986).

**21. **M. O. Scully and M. S. Zubairy, *Quantum Optics* (Cambridge University Press, 1997).

**22. **V. P. Bykov, *Laser Electrodynamics* (Cambridge Scholars, 2008).

**23. **F. Wang and Y. R. Shen, “General properties of local plasmons in metal nanostructures,” Phys. Rev. Lett. **97**(20), 206806 (2006). [CrossRef] [PubMed]

**24. **A. K. Sarychev and G. Tartakovsky, “Magnetic plasmonic metamaterials in actively pumped host medium and plasmonic nanolaser,” Phys. Rev. B **75**(8), 085436 (2007). [CrossRef]

**25. **R. H. Pantell and H. E. Puthoff, *Fundamentals of quantum electronics* (Wiley, 1969).

**26. **L. Novotny and B. Hecht, *Principles of Nano-Optics* (Cambridge University Press, 2006), p. 558.

**27. **E. S. Andrianov, A. A. Pukhov, A. V. Dorofeenko, A. P. Vinogradov, and A. A. Lisyansky, “Rabi oscillations in spasers during non-radiative plasmon excitation,” Phys. Rev. B **85**(3), 035405 (2012). [CrossRef]

**28. **A. N. Lagarkov, A. K. Sarychev, V. N. Kissel, and G. Tartakovsky, “Superresolution and enhancement in metamaterials,” Phys. Usp. **52**(9), 959–967 (2009). [CrossRef]

**29. **A. S. Rosenthal and T. Ghannam, “Dipole nanolasers: A study of their quantum properties,” Phys. Rev. A **79**(4), 043824 (2009). [CrossRef]

**30. **M. I. Stockman, “The spaser as a nanoscale quantum generator and ultrafast amplifier,” J. Opt. **12**(2), 024004 (2010). [CrossRef]

**31. **I. E. Protsenko and E. P. O’Reilly, “Dipole lasing phase transitions in media with singularities in polarizabilities,” Phys. Rev. A **74**(3), 033815 (2006). [CrossRef]

**32. **E. S. Andrianov, A. A. Pukhov, A. V. Dorofeenko, A. P. Vinogradov, and A. A. Lisyansky, “Stationary behavior of a chain of interacting spasers,” Phys. Rev. B **85**(16), 165419 (2012). [CrossRef]

**33. **L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii, *Electrodynamics of Continuous Media* (Butterworth-Heinemann, 1995), Vol. 8, p. 460.

**34. **B. Luk’yanchuk, N. I. Zheludev, S. A. Maier, N. J. Halas, P. Nordlander, H. Giessen, and C. T. Chong, “The Fano resonance in plasmonic nanostructures and metamaterials,” Nat. Mater. **9**(9), 707–715 (2010). [CrossRef] [PubMed]

**35. **U. Fano, “Some theoretical considerations on anomalous diffraction gratings,” Phys. Rev. **50**(6), 573 (1936). [CrossRef]

**36. **J. B. Pendry and S. A. Maier, “Comment on ‘Spaser action, loss compensation, and stability in plasmonic systems with gain’,” Phys. Rev. Lett. **107**(25), 259703, discussion 259704 (2011). [CrossRef] [PubMed]

**37. **M. I. Stockman, “Stockman Replies,” Phys. Rev. Lett. **107**(25), 259702 (2011). [CrossRef]

**38. **M. I. Stockman, “Loss compensation by gain and spasing,” Philos Trans A Math Phys Eng Sci **369**(1950), 3510–3524 (2011). [CrossRef] [PubMed]

**39. **I. E. Tamm, *Fundamentals of the theory of electricity* (Mir, Moscow, 1979).

**40. **D. G. Baranov, E. S. Andrianov, A. P. Vinogradov, and A. A. Lisyansky, “Exactly solvable toy model for surface plasmon amplification by stimulated emission of radiation,” Opt. Express **21**(9), 10779–10791 (2013). [CrossRef] [PubMed]