Design optimization of spasers considering the degeneracy of excited plasmon modes

We model spaser as an n-level quantum system and study a spasing geometry comprising of a metal nanosphere resonantly coupled to a semiconductor quantum dot (QD). The localized surface plasmons are assumed to be generated at the nanosphere due to the energy relaxation of the optically excited electron-hole pairs inside the QD. We analyze the total system, which is formed by hybridizing spaser’s electronic and plasmonic subsystems, using the density matrix formalism, and then derive an analytic expression for the plasmon excitation rate. Here, the QD with three nondegenerate states interacts with a single plasmon mode of arbitrary degeneracy with respect to angular momentum projection. The derived expression is analyzed, in order to optimize the performance of a spaser operating at the triple-degenerate dipole mode by appropriately choosing the geometric parameters of the spaser. Our method is applicable to different resonator geometries and may prove useful in the design of QD-powered spasers. © 2013 Optical Society of America OCIS codes: (250.5403) Plasmonics; (140.3430) Laser theory; (230.0230) Optical devices; (250.5590) Quantum-well, -wire and -dot devices. References and links 1. M. Premaratne and G. P. Agrawal, Light Propagation in Gain Media: Optical Amplifiers (Cambridge University, 2011). 2. D. K. Gramotnev and S. I. Bozhevolnyi, “Plasmonics beyond the diffraction limit,” Nature Photon. 4, 83–91 (2010). 3. S. A. Maier and H. A. Atwater, “Plasmonics: Localization and guiding of electromagnetic energy in metal/dielectric structures,” J. Appl. Phys. 98, 011101 (2005). 4. S. Maier, Plasmonics: Fundamentals and Applications (Springer, 2007). 5. 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, 027402 (2003). 6. R. F. Oulton, “Plasmonics: Loss and gain,” Nature Photon. 6, 219–221 (2012). 7. M. I. Stockman, “Spasers explained,” Nature Photon. 2, 327–329 (2008). 8. J. Seidel, S. Grafström, and L. Eng, “Stimulated emission of surface plasmons at the interface between a silver film and an optically pumped dye solution,” Phys. Rev. Lett. 94, 177401 (2005). 9. M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature 460, 1110–1112 (2009). 10. R. A. Flynn, C. S. Kim, I. Vurgaftman, M. Kim, J. R. Meyer, A. J. Mäkinen, K. Bussmann, L. Cheng, F. S. Choa, and J. P. Long, “A room-temperature semiconductor spaser operating near 1.5 μm,” Opt. Express 19, 8954–8961 (2011). 11. N. Zheludev, S. Prosvirnin, N. Papasimakis, and V. Fedotov, “Lasing spaser,” Nature Photon. 2, 351–354 (2008). #190208 $15.00 USD Received 8 May 2013; revised 12 Jun 2013; accepted 12 Jun 2013; published 19 Jun 2013 (C) 2013 OSA 1 July 2013 | Vol. 21, No. 13 | DOI:10.1364/OE.21.015335 | OPTICS EXPRESS 15335 12. S. W. Chang, C. Y. A. Ni, and S. L. Chuang, “Theory for bowtie plasmonic nanolasers,” Opt. Express 16, 10580– 10595 (2008). 13. A. Lisyansky, I. Nechepurenko, A. Dorofeenko, A. Vinogradov, and A. Pukhov, “Channel spaser: Coherent excitation of one-dimensional plasmons from quantum dots located along a linear channel,” Phys. Rev. B 84, 153409 (2011). 14. M. Grundmann, J. Christen, N. N. Ledentsov, J. Böhrer, D. Bimberg, S. S. Ruvimov, P. Werner, U. Richter, U. Gösele, J. Heydenreich, V. M. Ustinov, A. Y. Egorov, A. E. Zhukov, P. S. Kop’ev, and Z. I. Alferov, “Ultranarrow luminescence lines from single quantum dots,” Phys. Rev. Lett. 74, 4043–4046 (1995). 15. S. Mukamel, Principles of Nonlinear Optical Spectroscopy, Oxford series on optical sciences (Oxford University, 1999). 16. A. V. Baranov, A. V. Fedorov, I. D. Rukhlenko, and Y. Masumoto, “Intraband carrier relaxation in quantum dots embedded in doped heterostructures,” Phys. Rev. B 68, 205318 (2003). 17. A. V. Fedorov, A. V. Baranov, I. D. Rukhlenko, T. S. Perova, and K. Berwick, “Quantum dot energy relaxation mediated by plasmon emission in doped covalent semiconductor heterostructures,” Phys. Rev. B 76, 045332 (2007). 18. M. I. Stockman, “The spaser as a nanoscale quantum generator and ultrafast amplifier,” J. Opt. 12, 024004 (2010). 19. M. I. Stockman, “Nanoplasmonics: past, present, and glimpse into future,” Opt. Express 19, 22029–22106 (2011). 20. I. E. Protsenko, A. V. Uskov, O. A. Zaimidoroga, V. N. Samoilov, and E. P. O‘reilly, “Dipole nanolaser,” Phys. Rev. A 71, 063812 (2005). 21. J. B. Khurgin and G. Sun, “Injection pumped single mode surface plasmon generators: threshold, linewidth, and coherence,” Opt. Express 20, 15309–15325 (2012). 22. J. B. Khurgin and G. Sun, “How small can nano be in a nanolaser?” Nanophotonics 1, 3–8 (2012). 23. E. Andrianov, A. Pukhov, A. Dorofeenko, A. Vinogradov, and A. Lisyansky, “Dipole response of spaser on an external optical wave,” Opt. Lett. 36, 4302–4304 (2011). 24. E. Andrianov, A. Pukhov, A. Dorofeenko, A. Vinogradov, and A. Lisyansky, “Forced synchronization of spaser by an external optical wave,” Opt. Express 19, 24849–24857 (2011). 25. E. S. Andrianov, A. A. Pukhov, A. V. Dorofeenko, A. P. Vinogradov, and A. A. Lisyansky, “Rabi oscillations in spasers during nonradiative plasmon excitation,” Phys. Rev. B 85, 035405 (2012). 26. A. Ridolfo, O. Di Stefano, N. Fina, R. Saija, and S. Savasta, “Quantum plasmonics with quantum dot-metal nanoparticle molecules: influence of the fano effect on photon statistics,” Phys. Rev. Lett. 105, 263601 (2010). 27. A. Rosenthal and T. Ghannam, “Dipole nanolasers: A study of their quantum properties,” Phys. Rev. A 79, 043824 (2009). 28. A. K. Sarychev and G. Tartakovsky, “Magnetic plasmonic metamaterials in actively pumped host medium and plasmonic nanolaser,” Phys. Rev. B 75, 085436 (2007). 29. 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, 19785–19798 (2008). 30. 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, 127401 (2010). 31. A. Fang, T. Koschny, and C. M. Soukoulis, “Lasing in metamaterial nanostructures,” J. Opt. 12, 024013 (2010). 32. D. Sarid and W. Challener, Modern Introduction to Surface Plasmons: Theory, Mathematica Modeling, and Applications (Cambridge University, 2010). 33. P. K. Jain, K. S. Lee, I. H. El-Sayed, and M. A. El-Sayed, “Calculated absorption and scattering properties of gold nanoparticles of different size, shape, and composition: Applications in biological imaging and biomedicine,” J. Phys. Chem. B 110, 7238–7248 (2006). 34. K. L. Kelly, E. Coronado, L. L. Zhao, and G. C. Schatz, “The optical properties of metal nanoparticles: The influence of size, shape, and dielectric environment,” J. Phys. Chem. B 107, 668–677 (2003). 35. A. L. Aden and M. Kerker, “Scattering of electromagnetic waves from two concentric spheres,” J. Appl. Phys. 22, 1242 –1246 (1951). 36. C. Bohren and D. Huffman, Absorption and scattering of light by small particles (Wiley, 1983). 37. E. Sondheimer, “The mean free path of electrons in metals,” Adv. Phys. 1, 1–42 (1952). 38. M. I. Stockman, “Nanoplasmonics: The physics behind the applications,” Phys. Today 64, 39–44 (2011). 39. J. B. Khurgin, G. Sun, and R. Soref, “Practical limits of absorption enhancement near metal nanoparticles,” Appl. Phys. Lett. 94, 071103–071103 (2009). 40. J. B. Khurgin, G. Sun, and R. Soref, “Electroluminescence efficiency enhancement using metal nanoparticles,” Appl. Phys. Lett. 93, 021120–021120 (2008). 41. I. D. Rukhlenko, D. Handapangoda, M. Premaratne, A. V. Fedorov, A. V. Baranov, and C. Jagadish, “Spontaneous emission of guided polaritons by quantum dot coupled to metallic nanowire: Beyond the dipole approximation,” Opt. Express 17, 17570–17581 (2009). 42. L. Landau, E. M. Lifshitz, and L. P. Pitaevskii, Course of Theoretical Physics Vol 8: Electrodynamics of Continuous Media (Elsevier, 2004). 43. M. Ventra, S. Evoy, and J. Heflin, Introduction to Nanoscale Science and Technology, Nanostructure Science and Technology (Springer, 2004). #190208 $15.00 USD Received 8 May 2013; revised 12 Jun 2013; accepted 12 Jun 2013; published 19 Jun 2013 (C) 2013 OSA 1 July 2013 | Vol. 21, No. 13 | DOI:10.1364/OE.21.015335 | OPTICS EXPRESS 15336 44. A. Fedorov, A. Baranov, and Y. Masumoto, “Coherent control of optical-phonon-assisted resonance secondary emission in semiconductor quantum dots,” Opt. Spectrosc. 93, 52–60 (2002). 45. I. Rukhlenko, A. Fedorov, A. Baymuratov, and M. Premaratne, “Theory of quasi-elastic secondary emission from a quantum dot in the regime of vibrational resonance,” Opt. Express 19, 15459–15482 (2011). 46. A. Fedorov and I. Rukhlenko, “Study of electronic dynamics of quantum dots using resonant photoluminescence technique,” Opt. Spectrosc. 100, 716–723 (2006). 47. A. Ansel’m, Introduction to Semiconductor Theory (Mir, 1981). 48. D. Bimberg, R. Blachnik, P. Dean, T. Grave, G. Harbeke, K. Hübner, U. Kaufmann, W. Kress, O. Madelung et al., Physics of Group IV Elements and III-V Compounds / Physik der Elemente der IV. Gruppe und der III-V Verbindungen, v. 17 (Springer, 1981). 49. U. Fano, “Description of states in quantum mechanics by density matrix and operator techniques,” Rev. Mod. Phys. 29, 74–93 (1957). 50. K. Blum, Density Matrix Theory and Applications (Springer, 2010). 51. F. Wang and Y. R. Shen, “General properties of local plasmons in metal nanostructures,” Phys. Rev. Lett. 97, 206806 (2006). 52. G. Sun, J. B. Khurgin, and C. Yang, “Impact of high-order surface plasmon modes of metal nanoparticles on enhancement of optical emission,” Appl. Phys. Lett. 95, 171103–171103 (2009). 53. J. Lim, A. Eggeman, F. Lanni, R. D. Tilton, and S. A. Majetich, “Synthesis and single-particle optical detection of low-poly


Introduction
The emerging era of nanoplasmonics is expected to improve the speed and efficiency of optical devices, by allowing miniaturization beyond the diffraction limit using surface plasmons in circuits [1][2][3].Despite their advantages in miniaturization, surface plasmons which excite at metal-dielectric interfaces are highly dissipative and vanish in a few wavelengths when traveling [1,4].Therefore, the energy must be transferred from an external source to the surface plasmon wave to sustain its existence in nanoplasmonic circuits [1].The spaser, which is the nanoplasmonic counterpart of a conventional laser is the prime device for generating these surface plasmon waves and amplifying them during propagation [5,6].
Bergman and Stockman proposed the theory of the spaser [5] by claiming that a nanosystem with an excited active medium can transfer energy nonradiatively to a closely located plasmonic resonator and excite localized fields inside it.This energy transfer takes place due to the interaction between the resonator and active medium through the near field.It was stated that, electronic transitions in active medium are stimulated by the surface plasmons which already exist in the system resulting in a multiplication of the plasmon population, very much resembling the coherent feedback strategy used in conventional laser cavities [5,7].
Since its theoretical formulation, many experimental efforts have been carried out to fabricate a spaser.Seidel et al. [8] proved the possibility of the stimulated emission of surface plasmon by amplifying surface plasmons at the interface between a flat continuous silver film and a liquid containing organic dye molecules.The first demonstration of a spaser was done by Noginov et al. [9] using a gold nanosphere with a 7 nm radius, surrounded by the active medium, which is a 15 nm thick silica shell containing dye molecules.When the dye is optically pumped, it releases energy to the gold nanosphere causing excitation of localized surface plasmons.Another realization of a spaser was reported by Flynn et al. [10] who sandwiched a gold-film plasmonic waveguide between optically pumped InGaAs quantum wells.
Several other structures have also been proposed for spasers, including a V-shaped metallic nanoparticle attached to quantum dots (QDs) [5], a two-dimensional array of split-ring shape plasmonic resonators supported by a substrate acting as the active medium [11], a bowtieshaped metallic structure, in which QDs are placed in the bowtie gap and multiple quantum wells are located in the substrate [12], and a metal groove with QDs placed at its bottom [13].Even though most of these spaser designs employ QDs in the active medium, it is possible to have optically pumped rare-earth ions, dyes or bulk electrical injection as the excitable gain element [1].QDs are widely used in many lasing setups.As the charge carriers in a QD are confined in all three dimensions to a very small size, its density of states almost resembles a delta function, mimicking an atomistic behavior with a well defined spectral response.They also promise a better stability over temperature variations [14][15][16][17].
Although ample research has been done on spasers, it is seen that results depend on the model.Furthermore, most of these studies very much focus on the role of the active medium.In this paper, instead of only considering the states of the active medium, we analyze the electronic and plasmonic states of the whole spaser as a single quantum system.Degeneracy of plasmon modes is also considered.This is a more general approach to describe the spaser quantum mechanically compared to available literature.We rigorously describe our model using a simple spaser geometry comprising of a spherical metal nanoparticle and a QD.A major importance of the spaser is its ability to be effectively used in nanoplasmonic circuits, to generate surface plasmon waves.This paper is based on studying the capability of a spaser to excite these surface plasmons and finding all the parameters affecting the excitation rate.We also consider dissipations to the environment and pay attention on all the material properties influencing the spaser operation.The presenting work also provides clear design guidelines for a spherical spaser.These guidelines allow one to select its operating wavelength, possible values for the size parameters of the spaser components, and optimum QD placement.
It is important to compare our work with the related literature.Many theoretical studies have been performed on interaction between the gain elements and plasmonic systems with either localized or propagating surface plasmons.In order to determine spaser characteristics such as threshold spasing condition, spasing frequency and spaser linewidth, Stockman [5,18,19] analyzed the interaction between a metal nanoparticle and closely located chromophores which are modeled as two level systems (TLS).Introducing a 'dipole nanolaser', Protsenko et al. [20] derived the equations of motion for spasers.Khurgin et al. have also discussed the threshold spasing condition, linewidth and coherence of a spaser in [21,22].Considering a simple spaser formed by a two level QD attached to a metal nanoparticle, Andrianov et al. have examined some attributes of spasers such as dipole response to an optical wave [23], external forced synchronization [24] and occurrence of Rabi oscillations [25].Among the other investigations on radiation emitted from nanoparticle-TLS ensembles, [26] which inspect the resultant photon statistics, and [27] which presents the spectral properties of these photons are conspicuous.The theory of spaser also suggests that propagating surface plasmons can be amplified by placing a gain medium in the vicinity [1].This idea has been exploited to compensate the losses in metamaterials.A few of such utilizations are described in [28,29] where plasmon-TLS interactions and [30,31] in which plasmon-four level system interactions are analyzed.
Among the quantum mechanical models of the spaser, the work of Stockman [5,18,19] and Protsenko et al. [20] are noteworthy.However, they only analyze the quantum states of the active medium assuming it as a TLS.We improve this model by considering the whole spaser as a single n-level quantum system with a 3-level active medium.Stockman considers a third level [18] but assumes that its population is negligible and electron-hole pairs rapidly relax to the second level.In our model, there is a finite third level population and electron-holes pairs relax to the second level by interacting with the environment at a defined rate.Having a three level system, a designer has more control on choosing a suitable pump frequency such that it doesn't overlap with any surface plasmon resonance frequency of the resonator.In addition, we consider dissipations through the generalized master equation in the form of interaction with the system's bath.Here, all the major dissipation rates at each level are well represented.Since our intention is to find all the parameters affecting the plasmon excitation rate, it is required to consider all the forms of losses.
In our analysis, we first find the electric field and the eigenfrequencies of the localized surface plasmon modes in the resonator, then characteristics of the active medium, followed by the derivation of the spaser Hamiltonian, which is needed to analyze the spaser kinetics using the density matrix theory.Then we obtain an expression for the plasmon excitation rate of the spaser, and study how the geometric parameters affect the operation of a dipole spaser.

Spaser model
A spaser consists of a plasmonic resonator (or cavity), which supports surface plasmon modes, and an active medium, which amplifies the surface plasmons [5]. Figure 1 shows a schematic diagram of our spaser model where the resonator is a metal nanosphere surrounded by dielectric shell, and a QD is embedded in the shell as the active medium.Spasing occurs as a consequence of the nonradiative energy transfer from the QD to the nanosphere, exciting localized surface plasmon modes inside it.
Fig. 1.Spaser model under consideration.R 1 is the nanosphere radius, R 2 is the radius of the dielectric shell's outer boundary.A QD is located at the position r 0 with respect to the nanosphere's center.The parameters ε 1 , ε 2 and ε 3 denote the permittivities of the nanosphere, dielectric shell, and ambient, respectively.
The spherical core-shell structure is one of the most studied geometries for surface plasmon resonance [32][33][34], because it is one of the few geometries where Maxwell's equations can be analytically solved [35,36].Suppose the inner radius of the metal nanosphere is R 1 and it is smaller than the skin depth of the metal [37], which is usually around 25 nm for noble metals.The latter assumption is necessary to allow the localized surface plasmon modes to penetrate everywhere within the nanosphere and excite coherent electron cloud oscillations [18,38].Hence, the nanosphere radius is very smaller compared to the wavelength of incident light.Further, we assume that R 1 is greater than v F /ω (where v F is Fermi velocity and ω is the surface plasmon frequency), which is about 1 nm for noble metals, to avoid the effects of Landau damping [38].The outer radius of the shell, R 2 must be chosen such that the shell thickness is large enough to the QD be entirely embedded within the dielectric (i.e.R 2 − R 1 > QD diameter).
When the nanosphere radius and the shell thickness are fixed, the resonator supports a series of plasmon modes (dipole, quadrupole, etc.) with unique energies [32,33].These plasmon modes may overlap with the QD in different degrees, and hence experience gain in different amounts.The plasmon mode receiving the highest gain survives and becomes the dominant spaser mode [9,38].The gain to the spasing mode is received as a result of the electronic transitions in the QD causing recombination of an electron-hole pair.To ensure continuous spasing, we facilitate population inversion of the two energy levels pertaining to these transitions using a suitable pump source.Here we opt for optical pumping for the sake of design simplicity.Even though QD is the target of the pumping field, it may also excite the resonator in the vicinity [39].To minimize this extraneous effect, we intentionally select the optical pumping frequency with considerable detuning from the surface plasmon resonances of the resonator.Therefore, the excited field is much weaker compared with the spaser mode which overlaps with QD emission lines.However, we have not made any restrictive assumption here because one may completely avoid this situation by adopting a different pumping mechanism such as direct electric injection to QD as suggested in [40].

Localized surface plasmon modes of the resonator
To find the electric field and eigenfrequencies of the localized surface plasmon modes, we assume that QD's presence does not perturb the electric field in the system because it is very small compared to the nanosphere and therefore the resultant permittivity change of the shell is insignificant.We adopt a standard spherical coordinate system (r, θ , ϕ ) where the origin coincides with the center of the spherical shell (see Fig. 1).The supported electromagnetic modes can be found by solving the vector Helmholtz equation for the spherical shell following the method of Debye potentials [35,36].The solution gives a series of localized surface plasmon modes denoted by the angular momentum number l p and its projection m p , where ω l p is the angular frequency of the plasmon mode l p .The stationary electric field is given by the expression where l (kr) is one of the spherical Bessel functions j l (kr), y l (kr), h l (kr) for ν = 1, 2, 3, P m l are associated Legendre polynomials, k is the wavenumber, and a l , d l , g l and f l are constants assuring the field's continuity at boundaries.
Plasmon modes expressed by Eq. ( 1) possess a unique energy determined by hω l p , which can be found by solving the dispersion relation of the resonator.In order to obtain the dispersion relation, we first write a set of equations expressing the constants a l , d l , g l and f l in Eq. ( 2) by equating the tangential components of the field at the metal-dielectric and dielectric-ambient interfaces.This results in a homogeneous system of linear equations in which we apply the condition for a nontrivial solution and obtain the following dispersion relation for the lth mode surface plasmon resonance: where Here the Reccati-Bessel functions are ψ l (x) = x j l (x), ξ l (x) = xh (1) l (x) and χ l (x) = −xy l (x); k m , k d , and k 0 are the wavenumbers of the electromagnetic field in the metal, dielectric, and ambient, respectively.
According to Eq. ( 3), the energy of the plasmon mode, hω l is a function of R 1 , R 2 , ε 1 , ε 2 , and ε 3 , which implies that only the size parameters of the resonator determine the energy of the spaser mode when the spaser materials are chosen.Furthermore, this nontrivial solution of the boundary conditions allows us to express the coefficients a l , d l , g l in terms of f l as To find an expression for f l as a normalization constant, we follow the procedure of secondary quantization for dispersive media [41,42], and equate the total energy of the electric field of the plasmon mode to hω l , 1 4π Here V implies the volume integral over three dimensional Euclidean space.We can note that owing to this equality, f l implicitly depends on the spaser shell parameters R 1 and R 2 .

Active medium
Characteristics of the active medium determine the strength of its electron-hole pairs' interaction with the plasmon modes in the resonator and the amount of amplification that each plasmon mode receives.The characteristics of the QD includes its wavefunction and the energy levels of the electron-hole pairs.Electronic transitions between these levels fuel spasing and we assume that, as in the case of lasers, those transitions representing recombination of electron-hole pairs would excite the plasmon modes in the resonator.Normally, an excitable QD like this can be effectively described using three energy levels to account for pumping and stimulated transitions.We denote the ground level of the QD by state vector |0 e ⟩ and two excited levels by |1 e ⟩ and |2 e ⟩.To assign energies to these three levels, we have to take the quantum confinement effects into account.To do this, we approximate the QD by a sphere with a radius R q , confined by an infinite potential barrier (i.e.potential V (r) = 0 for r < R q and it is infinite otherwise).The resulting time-independent Schrodinger equation is separable in spherical coordinates and possesses a solution similar to the hydrogen atom and given by [43] ψ n q ,l q ,m q (r, θ , ϕ ) = where n q , l q , m q are principal, azimuthal and magnetic quantum numbers describing the states, ξ n q l q is the n q th root of the spherical Bessel function of the first kind (i.e j(ξ n q l q ) = 0), and Y m q l q are spherical harmonics.The corresponding eigenvalue of the wave function ψ n q ,l q ,m q gives the excitable energy levels of an electron-hole pair: is the reduced mass of an electron-hole pair, m * e and m * h are effective masses of an electron and a hole, and E g is the bandgap of the QD material.These energy levels denoted by E n q ,l q are 2l q + 1 times degenerate.For a transition from an initial state |s i ⟩ with quantum numbers s i = (n i , l i , m i ) to a final state s f ⟩ with quantum numbers s f = (n f , l f , m f ), the absorbed energy from the system will be ∆E s f ,s i = E n f ,l f − E n i ,l i .Energy will be released in case this quantity is negative.
Having this knowledge, we map three QD energy levels.For |1 e ⟩ |2 e ⟩, we select two lowest energy levels [i.e (n q , l q ) = (1, 0), (1, 1)], assuming that probability of populating the higher energy levels are very small, and |0 e ⟩ is mapped to the ground level.This mapping also enables us to calculate the QD radius, R q required for an efficient energy transfer to the resonator where the energy received by the spaser modem, hω l p matches the energy released by the QD, −∆E s 0e ,s i , giving the resonance QD radius: Although achieving a perfectly matched resonance may be difficult in practice, having a closer value to the quantity in Eq. ( 6) will be adequate for spasing.This analysis of QD reveals how its geometrical parameters affect the spaser operation.γ pl is the plasmon decay rate and ξ i j denotes the rate of transition from the state | j s ⟩ to |i s ⟩ due to the interaction with bath.

Spaser kinetics
In the previous section we analyzed the characteristics of the isolated electronic and plasmonic subsystems.Here, we let them to interact with each other and create the functioning spaser system, as shown in Fig. 2. In the electronic subsystem, we have three states denoted by |0 e ⟩, |1 e ⟩, and |2 e ⟩, as described in Section 2.2.We assume that active medium strongly interacts only with the plasmon mode l p , which is the spaser mode.Here we note that this assumption is not valid for higher modes when frequency spacing between them become smaller than the QD emission linewidth.Let us define the state 0 pl ⟩ with zero plasmons, and 2l p + 1 states 1 (l p ,m p ) pl ⟩ with one plasmon, as the plasmonic subsystem.Then we amalgamate these two subsystems to makeup the total system shown in Fig. 2(c) with n = 2l p + 4 product states defined as The product states |1 s ⟩ , |2 s ⟩ and |3 s ⟩ are associated with zero plasmons and |4 s ⟩ , |5 s ⟩ , ..., |n s ⟩ possess one plasmon of the spaser mode.The |1 s ⟩ → |3 s ⟩ transition, which is the excitation of ground electron-hole pairs to the highest energy level in our model, occurs due to the electronhole pairs' interactions with the pump light which we analyze classically.Transitions from the state |2 s ⟩ to one of the states |4 s ⟩ , |5 s ⟩ , ..., |n s ⟩ is the driving force for the phenomena of spasing because they excite plasmon modes in the resonator.Some transitions may occur from the state | j s ⟩ to |i s ⟩ due to the interaction with the bath.They can be considered as dissipations.Having this model, we analyze the kinetics of the n-level system by first constructing its Hamiltonian, and then deriving the density matrix equations to find the corresponding active state populations.

Hamiltonian of the spaser
Hamiltonian H of the spaser contains the non interacting electronic and plasmonic Hamiltonians, H e and H pl , and Hamiltonian, H i , of the interacting subsystems: where The interacting Hamiltonian H i can be decomposed to represent the interactions between the electron-hole pairs and pump light as H e,L and the interactions between electron-hole pairs and surface plasmons as H e,pl : where Here V s f ,s i is the matrix element for the transition |i⟩ → | f ⟩, i, f = {0 e , 1 e , 2 e }, g l p = √ ω l p /2ε 2 V n , where V n is the normalization volume, φ(t) and ω L are the envelope function and frequency of the pump light, and c.c. represents the complex conjugate [44][45][46].Relaxation process from the state |2 e ⟩ to |1 e ⟩ is taken into account through the relaxation constant ξ 23 in spaser kinetics.The matrix element V l p ,m p s f ,s i corresponding to plasmon's interaction with electron-hole pairs can be written as where |u c ⟩ and |u v ⟩ are the Bloch functions, |s i ⟩ and s f ⟩ are the envelope wavefunctions of the initial and final electronic states characterized by the sets of quantum numbers s i and s f and r is the displacement vector of the electron-hole pairs.Assuming the equality E l p m p = E l p m p e l p m p , we may write The first matrix element in Eq. ( 12) can be expressed through the Kane's parameter P [47,48] as, ⟨u c | e l p m p .r |u v ⟩ = √ 2P/E g and the second matrix element, where the integration is evaluated over the quantum dot's volume V QD , r = r 0 +r ′ is the electron position inside the quantum dot, and r 0 is radius vector of the quantum dot's center.E l p m p (r) can be derived from Eq. (2) and ψ s f is taken from Eq. ( 5).In case QD is very small compared to the nanosphere, it is reasonable to assume that E l p m p (r) is approximately constant over the QD's volume, therefore the integral in Eq. ( 13) can be simplified to ϒ l p ,m p s f ,s i = E l p m p (r 0 )δ s f ,s i .Hence, from Eq. ( 12), the matrix element for the spaser mode's interaction with electron-hole pairs can be given by Since this quantity determines the contribution to the total Hamiltonian by the interaction of the spaser mode with the active medium, it also defines the strength of spasing.If the QD is moved away from the nanosphere this matrix element becomes very much smaller.Therefore, the closer the QD stronger the spasing.When calculating the V s 2e ,s 0e , electric field term E l p m p in Eq. ( 11) should be replaced by the electric field caused by pump light.

Plasmon excitation rate of the spaser
Having the Hamiltonian calculated, we analyze the n states system comprises of the product states |1 s ⟩ , |2 s ⟩ , ..., |n s ⟩, using the density matrix formalism.We define the populations of those states by ρ 11 , ρ 22 , ..., ρ nn and assume that the system has a short-term memory [49] and coupled to a bath which is the reservoir for system's dissipations.Then the relaxation superoperator, which is added to the commuted Hamiltonian and density operator to incorporate dissipations, consists of a set of constants that define the relaxation kinetics of diagonal and off-diagonal elements of the reduced density matrix [44,49].Using the Markov and secular approximations [50], master equation for the system can be given by where γ µ µ is the population relaxation rate of the state |µ s ⟩, γ µν = (γ µ µ + γ νν )/2 + γµν is the coherence relaxation rate between the states |µ s ⟩ and |ν s ⟩, γµν is the pure dephasing rate, and ξ νκ is the transition rate from state |κ s ⟩ to state |ν s ⟩ due to interaction with the bath.We assume that the lifetime of the ground state |1 s ⟩ is very large, by setting γ 11 = 0.The parameters γ 44 , γ 55 , ..., γ nn define the dissipation of degenerated states of the plasmon mode denoted by l p , m p quantum numbers.However, since the energy and dielectric properties are common for all the degenerated plasmon states, we assume that all these plasmon dissipation constants are equal and can be denoted by γ pl .Populations ρ 44 , ρ 55 , ..., ρ nn represent the excitation rates for the plasmon modes l p , m p = (l p , −l p ), (l p , −l p + 1), ...,(l p , l p ) respectively.The sum of the plasmon populations of l p th mode, R l p gives the number of plasmons excited in the spaser at a given time.Hence, this quantity can be referred as the effective plasmon excitation rate of the spaser.We solve the system of partial differential equations given in the Eq. ( 15) for the continuous wave (CW) operation assuming that φ(t) = 1 and obtain an expression for the total plasmon excitation rate: where ∆ L3 = ω L − ω 3 , ∆ 2 j = ω 2 − ω j when ω 2 , ω 3 , ω j and ω L are the energies of the states |2 s ⟩,|3 s ⟩,| j s ⟩ and pump light respectively.As all the degenerate plasmon modes have the same energy, ω j = ω l p ∀ j and therefore ∆ 2 j = ∆ 2p .We can assume that detuning of the energy of the pump light with the energy of the state |3 s ⟩, ∆ L3 << γ 13 and hence, γ 13 γ 2 13 +∆ 2

L3
≃ 1 γ 13 .Further, we also assume that all the degenerate states of the plasmon modes decay equally and therefore γ j j = γ pl and γ 2 j = γ 2p ∀ j.With these simplifications, Eq. ( 16) can be further simplified to The system interacting with the bath introduces various relaxations to the electronic subsystem.We assume that these relaxation rates and the decay rate of the dominant spaser mode depend only on the materials of the spaser [38].Since the nanosphere is much smaller than the wavelength, we can neglect the radiation losses [51] and only consider the nonradiative decay terms.However, this is not an overly restrictive assumption because in case radiation losses are not negligible, the decay rate γ pl of the spaser mode can also incorporate the resultant of both radia- tive and nonradiative decay rates.In addition, decay rates could potentially be different for other plasmon modes [52], which are not taken into account assuming that they weakly overlap with the QD emission spectrum.As discussed in Section 2.1, energy of the spaser mode ω l p depends on resonator's size parameters (i.e.R 1 and R 2 ) when the materials of the spaser components are chosen.Therefore, it can be observed from Eq. ( 17) that the total plasmon excitation rate of the spaser, R l p mainly depends on the matrix elements for the electron-hole pair-plasmon interaction, V l p ,m p s 1e ,s 0e for each degenerate state and the detuning ∆ 2p because we assume that matrix element for the pump light-QD interaction V s 2e ,s 0e is constant under CW operation.In the case of exact resonance (i.e.R q is given by Eq. 6), ∆ 2p = 0 and we achieve the highest plasmon excitation rate.
For fixed materials and size parameters R 1 , R 2 , there is a unique spaser mode energy hω l p and R q may be chosen according to ω l p to achieve a higher plasmon excitation rate.Then the position of the QD plays a major role deciding the amount of amplification of the plasmons as the electric field of the spaser mode changes with the location according to Eq. (2).To analyze these factors, we investigate the spaser's behavior with respect to spaser size parameters and the QD's location in the following section taking a dipole spaser as a case study.

Case study: A dipole spaser
Let us consider a spaser whose dipole mode (l p = 1) is amplified by the active medium.We construct this spaser using gold for the nanosphere, and coating silica (SiO 2 ) over it, making up a dielectric shell.We embed a CdSe QD in this shell.Once the materials for spaser components are chosen, we are ready to investigate the dipole spaser's operation according to its geometric parameters.Here, we use the analytical results obtained in Sections 2 and 3 on the resonator's plasmonic properties and the characteristics of the resonator-QD interactions.
The knowledge of the frequency dependence of the permittivities are required to calculate the electric field of the spaser mode and the plasmon decay rate.For silica, we assume a nondispersive permittivity, ε 2 = 2.15 in contrast to gold in which we assume a frequency dependent permittivity.Furthermore, since our gold nanosphere is very small, we have to consider the size dependency modification of the permittivity as well [53,54].In order to incorporate this size effect, we use the model [54], , where ω p is the bulk plasma fre- quency of gold, Γ is the electron collision frequency with damping by Γ = γ b + v F /2R 1 , γ b is the bulk electron collision frequency of gold, v F is the Fermi velocity, and ε ∞ is the con- tribution from the interband transitions obtained by fitting ε 1 (ω, R 1 )| R 1 →∞ to the experimental data published by Johnson and Christy [55] for bulk material.Also we assume ω p = 1.36×10 16  s −1 , γ b = 3.33 × 10 13 s −1 , v F = 1.4 × 10 6 m/s, and ε ∞ = 9.84 [54,56].Spaser's outer boundary is assumed to be in free space and hence ε 3 = 1.
Using these permittivities, we solve the dispersion relation given in Eq. ( 3) for dipole mode and then plot energy of the spaser mode as a function of the nanosphere radius and shell thickness as shown in Fig. 3(a).The contour plot shows that there can be many (R 1 , R 2 ) pairs which can result in the same spaser mode energy.For example, if the spaser mode energy is 2.385 eV (i.e.equivalent wavelength is approximately 520 nm), it traces a curve on the contour plot as marked.It is important to note that, although we refer energies in electron volts for convenience, they should be converted to SI units when substituted to equations.
To find the appropriate QD radius to match the operating point, we substitute the corresponding spaser mode energy in Eq. ( 6).For the previously set operating point, 2.385 eV, the substitution gives resonant QD radius to be 2.4 nm, when the QD material parameters are E g = 1.74 eV, m * e = 0.13, and m * h = 0.45 [57,58].Calculating the energy levels of the QD with this radius, and considering the fact that dipole mode is triply degenerate with m p = −1, 0 and 1, we redraw the total system diagram given in Fig. 2(c) for this dipole spaser to attain the system illustrated in Fig. 4. The corresponding electronic levels |0 e ⟩ , |1 e ⟩ and |2 e ⟩ of the QD in this system posses energies 0, 2.385 and 3.059 eV respectively.The resulting total system has n = 6 states denoted by Evaluating this 6 state system of the dipole spaser for the highest plasmon excitation rate gives the optimum size parameters.Just opting an arbitrary R 1 , R 2 pair on the preferred spaser mode energy curve may not result in the highest plasmon generation.To investigate the behavior of our dipole spaser, we follow the results derived in Section 3 by substituting l p = 1.As we discussed there, the total plasmon excitation rate of the dipole spaser is represented by the sum of the populations of the states |4 s ⟩ , |5 s ⟩, and |6 s ⟩, R 1 = ρ 44 + ρ 55 + ρ 66 .We use the expression for R l p given in Eq. (17) to study the plasmon excitation rate of the dipole spaser.This equation contains some relaxation constants, γ 13 , γ 22 , γ 33 , γ 2p , and ξ 23 , which depend on the environment and used materials, but do not depend on spaser's geometry.Therefore, we keep them constant and compare the plasmon excitation rate by investigating the normalized plasmon population in the relevant plots, as our objective is to study how the spaser's geometrical parameters result in relatively high or low plasmon populations.While calculating the matrix elements, for the sake of simplicity, we assume that QD is located on the nanosphere's dipole axis in θ = 0 direction.
The only relaxation rate influenced by spaser's geometry is the decay rate of surface plasmons, γ pl is a function of the spaser mode energy given by γ [19].Since the nanosphere is very much smaller than the wavelength, here we consider only nonradiative decay assuming that its radiation loss is negligible [51].By substituting the model for nanosphere's permittivity in this expression , the plasmon decay rate can be simplified to . Since both the quantities Γ and ω l p depend on size parameters R 1 and R 2 , γ pl also depends on them for a given material.Continuing with these simplifications, we plot the plasmon excitation rate, R 1 of the dipole spaser with respect to nanosphere radius and shell thickness, as shown in Fig. 3(b), when the QD's location is fixed to the middle of the dielectric shell and R q is fixed to the resonant radius given by Eq. ( 6).It can be noted that plasmon excita-  tion rate is higher for smaller nanosphere radii and shell thicknesses.It monotonically decreases when the total volume increases.This plot helps us to figure out the (R 1 , R 2 ) pair that gives the highest plasmon excitation rate for a preferred spaser mode energy marked on Fig. 3(a).Hence, a designer has the freedom to optimize the spaser geometry by tuning either nanosphere radius or dielectric shell thickness.
Let us continue with the previous example and draw the curve corresponds to the spaser mode energy 2.385 eV on Fig. 3(b).Based on this curve, we can observe that the plasmon excitation rate varies with geometric parameters, even for the same spaser mode energy.Hence, we can select parameters for the optimum design by picking the R 1 , R 2 pair which results in the highest plasmon excitation rate.In our example, parameters (R 1 , R 2 ) = (6 nm, 11 nm) offer the highest spaser mode amplification for the 2.385 eV energy curve.Here we emphasize that, although we obtained this result for a single QD dipole spaser, a better spaser configuration may consists of many QDs to achieve a higher gain.If this dipole spaser consists of N QDs, then N ∝ (R 3 2 − R 3 1 )/R 3 q , if the QDs are uniformly distributed inside the dielectric shell.In that case, having a bigger shell volume to accommodate more QDs may result in a higher plasmon excitation rates because all QDs contribute the total amplification.However, each QD does not contribute evenly in many QDs case as the location of the QD and its size play vital roles in deciding the plasmon excitation rate.
In order to examine how the location of the QD affects the plasmon excitation rate of the spaser, we fix R 1 and R 2 , and plot R 1 for the case of resonant QD radius (i.e.R q is also fixed).
Here we vary QD's location r 0 within the dielectric shell from the innermost to the outermost position with respect to the nanosphere's center.Such plots for four different R 1 , R 2 pairs are shown in Fig. 3(d).It can be observed from these plots that the plasmon excitation rate rapidly decreases when the QD is moved away from the nanosphere.This happens because the interactions between plasmon modes and electron-hole pairs in QD gets weaker towards the shell's outer boundary as the matrix element for interactions, V l p ,m p 1e,0e decreases in the radial direction according to Eqs. ( 2), ( 13) and ( 14).The decreasing rate is higher for smaller nanosphere radii.Especially, when there are many QDs, their contribution to the spasing mode won't be uniform.Hence, when designing a multiple QD spaser, we must set the optimum size parameters such that a large number of QDs are closely located to the nanosphere center.
Thus far, we have only considered the case where the QD radius is tuned according to Eq. ( 6) such that its emission energy exactly equals to the energy of the spaser mode.To investigate the influence of the QD radius on spasing, we vary R q within an interval of 1-3 nm and plot the plasmon excitation rate, as shown in Fig. 3(c), for different shell thicknesses keeping the nanosphere radius fixed to 10 nm.According to this plot, highest plasmon excitation rate is observable in the case of exact resonance and it rapidly decreases when the R q deviates from the resonant QD radius.For example, if the the QD radius deviates much (i.e. about 0.5 nm) from the resonant value, the resultant plasmon excitation rate will tend to zero.
The plots in Figs.3(a)-3(d), clearly explicate how the geometrical parameters of spaser design can be optimized to attain an elevated plasmon excitation rate.In addition, we need to discuss the threshold power required from optical pumping which is a condition of spasing.Since we use a normalized electric field in this analytical treatment, it is important to calculate the required threshold gain, denoted by g th .Threshold gain can be found by applying the condition of population inversion: ρ 22 ≥ ρ 11 .It can be shown that it does not depend on spaser geometry, and is a function of the dielectric constants of the spaser materials and the spaser mode frequency given by [19] Threshold gain is a function of spaser mode frequency and nanosphere's permittivity.Since spaser mode frequency is a function of R 1 and R 2 and nanosphere's permittivity is a function of R 1 , threshold gain also becomes a function of these two geometrical parameters as we plot in Fig. 3(e).According to the graph, the gain does not vary much with the size parameters but the required threshold gain for spasing is little higher when the nanosphere is smaller.We mark the 2.385 eV curve, which we used in the previous example, on threshold pump power plot and it provides an idea of choosing the correct geometrical parameters.It should be noted that, the expression for threshold gain is derived assuming that stimulated emission dominates in the spaser.However, there can be cases where spontaneous emission dominates spaser kinetics when the number of plasmons is one.In such situations, the semiclassical concept of threshold is not applicable and plasmon population in the resonator increases linearly with the pump rate [59].In addition, extending this analysis to find the required pumping intensity will be useful in practical implementations.The developed model can also be used to study the polarization of resultant plasmons, which has been discussed in [20], as it allows us to consider the case where QD dipole moment is neither parallel nor orthogonal to the line connecting the centers of nanosphere and QD.However, in this study, we only focus on the optimization of spaser geometry to enhance plasmon generation based on the introduced quantum mechanical model.The case study which we performed throughout this section provides us design guidelines on choosing the optimum size parameters of spaser components and where to place QDs to achieve a maximum plasmon excitation rate.In the process of designing a spaser, one might choose its operating wavelength as the first step, then pick possible values for the size parameters R 1 , R 2 and R q by examining the resulting in plasmon excitation rate.Then, it is necessary to figure out where to embed the QD.Placing the QD closer to the nanosphere will give a higher plasmon excitation rate.However, if multiple QDs are to be placed, then one need to consider about the precise placement because not all QDs can be placed close to the nanosphere's boundary.
Although we studied a spaser realization with spherical structure, the derivations done in Section 3 is valid for a spaser with any geometry.It might not be possible to obtain the electric field and the energies of plasmon modes analytically for most of the other resonator geometries as we did in Section 2.1.But numerical solutions for the field can be admitted into the expressions derived in Section 3 to analyze the spaser kinetics and investigate the plasmon excitation rate which characterize its performance.

Conclusions
We have theoretically modeled spaser as an n-level quantum system formed by amalgamating spaser's electronic and plasmonic subsystems.With this model, we studied a simple spaser geometry consists of a metal nanosphere resonantly coupled to a QD.The energy transfer between the QD and the nanosphere accompanied the relaxation of electron-hole pairs resonantly generated inside the QD by a continuous-wave laser.By employing the density matrix theory, we analytically found the excitation rate of surface plasmons where three nondegenerate electron-hole pair states are coupled to a single plasmon mode of arbitrary angular momentum.The obtained expression was then examined numerically for the special case of a spaser operating at the triple-degenerate dipole mode.It was shown that the plasmon excitation rate can be significantly enhanced by appropriately choosing the geometric parameters of the spaser.

Fig. 2 .
Fig.2.(a) Electronic subsystem of the spaser representing the states of the QD's electronhole pairs.(b) Plasmonic subsystem of the spaser representing the degenerate localized surface plasmon modes in the resonator.(c) Total system of the spaser with n = 2l p + 4 levels.γ pl is the plasmon decay rate and ξ i j denotes the rate of transition from the state | j s ⟩ to |i s ⟩ due to the interaction with bath.
hω 1 e |1 e ⟩ ⟨1 e | + hω 2 e |2 e ⟩ ⟨2 e | , and H pl = ∑ l p m p =−l p hω l p b † l p ,m p b l p ,m p .Here b † l p ,m p ,b l p ,m p are the creation and annihilation operators of the surface plasmons corresponding to the quantum numbers l p and m p .
of the QD, r (nm)

0 - 1 Fig. 3 .
Fig.3.Plots of (a) energy of the spaser mode and (b) normalized plasmon excitation rate with respect to nanosphere radius and shell thickness when QD's location is fixed to the middle of the dielectric shell, (c) plasmon excitation rate with respect to QD radius and shell thickness when the nanosphere radius is 10 nm and QD's location is fixed to the middle of the dielectric shell, (d) plasmon excitation rate with respect to QD's location for different R 1 , R 2 pairs, (e) threshold gain with respect to nanosphere radius and shell thickness.