Multipolar, time-dynamical model for the loss compensation and lasing of a spherical plasmonic nanoparticle spaser immersed in an active gain medium

Abstract

The plasmonic response of a metal nanoparticle in the presence of surrounding gain elements is studied, using a space and time-dependent model, which integrates a quantum formalism to describe the gain and a classical treatment for the metal. Our model fully takes into account the influence of the system geometry (nanosphere) and offers for the first time, the possibility to describe the temporal evolution of the fields and the coupling among the multipolar modes of the particle. We calculate the lasing threshold value for all multipoles of the spaser, and demonstrate that the dipolar one is lowest. The onset of the lasing instability, in the linear regime, is then studied both with and without external field forcing. We also study the behaviour of the system below the lasing threshold, with the external field, demonstrating the existence of an amplification regime where the nanoparticle’s plasmon is strongly enhanced as the threshold is approached. Finally, a qualitative discussion is provided on later, non-linear stages of the dynamics and the approach to the steady-state of the spaser; in particular, it is shown that, for the considered geometry, the spasing is necessarily multi-modal and multipolar modes are always activated.

Introduction

Metallic nanostructures have since long attracted interest in different fields of nanotechnology, since they can sustain localized plasmon resonances which act as nanoscale concentrators of optical fields. In the presence of optical gain (active medium) near or within the nanostructures, new possibilities open up with the ability to further control, amplify and tune these localized surface plasmon resonances^{1,2,3,4,5,6,7,8}. The interest of coupling metallic nanostructures with active media culminates in the concept of the “spaser” (Surface Plasmon Amplification by Stimulated Emission of Radiation, the plasmonic equivalent of a laser), which widened their potential applicability to nanoscale lithography, probing, microscopy, optoelectronics and more^9,10. In order to effectively shape optical energy to the nanoscale according to our needs, it is of high importance to dispose of accurate theoretical descriptions taking fully into account not only the effects of the metal to active medium coupling and the associated dynamics, but also, on equal footing, the effects of the nanostructure geometry which governs the spatial structure of the spaser field. In this report, we propose a semi-classical model describing the coupling of the temporal quantum dynamics of the gain elements with a homogeneous, spherical nanoparticle (NP), taking into account not only the lowest, but also all other higher electromagnetic modes of the NP as well.

The past few years have seen a number of works on the optical properties of gain-assisted nanoparticles (NPs) with various geometries: spheres, core-shell, multiple core-shell, ellipsoids…, using more or less refined analytical descriptions, or numerical simulation tools like the finite-element method^{3,11,12,13,14,15,16}. All these works were carried out in stationary regimes for the NP and the gain (assuming flat, or Lorentzian, emission lineshapes for the latter). While such a simplification is acceptable far from the lasing threshold and it has been successfully used to describe experimental results^17,18,19, it has to fail near and above it, making the obtained predictions inconsistent: gain saturation or complex time-dependent effects may rise, begging for a more elaborate dynamical and non-linear description of the system.

Time-dependent approaches were also developed but are more scarce, starting with the seminal work of Bergman and Stockman introducing the spaser²⁰, and then in subsequent works^{4,21,22,23,24}. Nevertheless, often these approaches are incomplete: either space-dependent modes were not considered for the sake of generality, or they were averaged out through homogeneization, or only the dipolar mode was considered, leaving out all the effects related to the specific geometry of the actual nanoresonator used to generate the spasing. A few purely numerical works have integrated all time and space-dependent effects using two or four-level population dynamics for gain carriers (see ref. 2 and references therein) locally coupled to Maxwell equations, but they were only concerned with so-called “fishnet” geometries, unrelated to NP-based spasers. Moreover, powerful as they are, full-wave simulations do not necessarily allow for a deep understanding of the mechanisms at work, and a complementary model-based approach is undoubtedly useful. Finally, ref. 25 provides an analytical, time-dependent study of a NP in a gain core-metal shell configuration, where only the dipolar mode can exist.

In this article, we consider a single, homogeneous, metallic nanosphere (NP) of radius a (as illustrated in Fig. 1), immersed in a gain medium consisting of a dielectric host in which active elements are dispersed (e.g., a laser dye or quantum dot solution). This is the simplest imaginable geometry for spasing, but we will show that it gives rise to rich effects nonetheless. We emphasise that the aim of this work is not to study and optimize a given design for experimental spaser realizations, but rather to lay some understanding of basic physical effects that may arise in spaser systems across various geometries. Still, the situation of homogeneous plasmonic spheres immersed in a gain medium is of experimental relevance to some existing experimental studies^26,27.

Figure 1

The system under study: a single, homogeneous spherical nanoparticle made of metal (in light gray), immersed in an active host medium made of many active dipolar elements (red dots) dispersed into a background dielectric host.

The model we present is combining together aspects treated separately in the literature, namely: (i) taking into account the appropriate electromagnetic response of the NP, as induced by its specific geometry and inclusive of all multipolar fields; (ii) including the population dynamics of quantum emitters, in order to explore the time-dependent behavior of the NP-cum-emitters system; (iii) integrating interaction with a sea of emitters around the NP, which is coupled to the spatially inhomogeneous local field created by the NP.

The system will be studied in conditions where it is either excited by a probe field linearly polarized along the z-axis (with unit vector ): , or it is left isolated (E₀ = 0). Note that E₀ refers to the value of the exciting field taken inside the active medium. We refer to the region inside the spherical inclusion as region 1 and to the region outside as region 2, with respective electric fields E₁ and E₂ and respective polarizations P₁ and P₂. The system is taken in the quasi-static limit: where λ denotes the exciting probe wavelength, so that E₀ is assumed to be uniform over the NP.

Metallic Nanoparticle and Active Medium Description

We start by describing how the field E₁(r, t), where r is the spatial coordinate with origin at the particle center and t is time, acts on the electrons in the metallic nanoparticle. We describe this interaction using a free-electron model²⁸:

where the polarization P₁ = n_eed is produced by the displacement d of the electron cloud with respect to the equilibrium position (n_e is the electron density, e is the electron charge); γ is the ionic collisions friction coefficient. We also used the plasma frequency with m_e the electron mass. Assuming e^−iωt harmonic notation for fields, Eq. 1 classically gives back the Drude formula for the metal’s relative permittivity:

In the following, we will focus on the metal, silver, using experimental data from ref. 29 for numerical values of ε_m. Values for ω_pl and γ are found by fitting Eq. 2 against these data (also adding a constant offset ε_∞ to account for non-Drude contributions seen in the experimental measurements). We find: ħω_pl = 9.6 eV and ħγ = 0.0228 eV.

We now consider how the field E₂ interacts with the active medium outside the nanoparticle, which is described as a continuum of two-level emitters in a thermal bath using the the optical Bloch equations and the density matrix formalism^4,30,31:

Here ρ_ij is the i, j element of the density matrix; is an effective relaxation rate, with W the phenomenological pumping rate, and and τ₂ the time constants associated with energy (spontaneous emission) and phase relaxation processes due to the interaction with the thermostat; ω₂₁ is the transition frequency between levels 1 and 2; N = ρ₂₂ − ρ₁₁ is the population inversion, where ρ₂₂ and ρ₁₁ are the diagonal matrix density elements, and finally . The term μ · E₂ accounts for the non-radiative coupling to the metal, with μ the transition dipole moment, while the radiative contribution is taken into account in Eq. 4 through the τ₁ relaxation time. Importantly, note that all quantities in Eqs 3 and 4 are taken implicitly as both space and time-dependent (functions of r and t).

With χ_b the susceptibility of the background (passive) dielectric host, and n the number density for the active elements (e.g., dye molecules), the space and time-dependent polarization of region 2 is obtained by averaging over the dipole moments of the gain elements, assuming they are randomly oriented with respect to the field E₂ (due to Brownian motion in the host solution):

where Ω denotes solid angle.

We have defined the new quantity Π as

Physically, the role of Π is to keep track of the part of the polarization in medium 2 contributed to by the emitters. Integrating the system of equations 3 and 4 over solid angles, we obtain:

We next use the rotating waves approximation, assuming that the frequency of the probe field E₀ is near-resonant and assuming an e^−iωt harmonic form for all time-dependent quantities. The system to solve becomes:

together with

(Now Π, E_1,2 and P_1,2 from here on denote only the slowly-varying enveloppes of previously introduced physical quantities.) It is worth noting that the population inversion N = N(r, θ, t) as it appears in this system of equations is, in the general case, non-uniform due to the non-uniformity of the source term appearing in the r.h.s. of Eq. (10).

We can now in particular calculate the permittivity ε_h of the active medium under pumping (without particle), when excited by a uniform field E₂ = E₀: using Eq. (12), and the steady state-value of Π extracted from Eq. (9), one finds a Lorentzian lineshape centered on ω = ω₂₁,

where ε_b = 1 + 4πχ_b is the permittivity of the background medium, is the maximum level of gain and Δ = 2/τ₂ is the emission linewidth. Such Lorentzian emission lineshapes have been widely used in classical steady-state models for plasmon-gain interaction^3,11.

A metallic nanoparticle such as the one studied here naturally has a dipolar localized plasmon resonance frequency ω₀, which can be obtained from the permittivities ε_m and ε_h just calculated through the standard Fröhlich condition^11,32:

We will from here assume that the gain lineshape is centered exactly at the plasmon frequency, i.e., that ω₂₁ = ω₀. Note that, while this ensures a better efficiency of coupling between gain and metal, it is by no means a requirement for the effects described below to occur.

Multipolar Mode Equations and Boundary Conditions

We now look into describing the complete system composed of the NP immersed in the active medium. We start by assuming that the probe field amplitude E₀, as well as E₂, are sufficiently small that the right-hand-side term of equation 10, which measures the rate at which the population inversion of the gain elements is depleted, will remain negligible. This means that the population inversion, according to Eq. (10), remains spatially uniform: N(r, θ, t) = N₀(t). (Situations where this approximation does not hold anymore will be discussed in the last Section of this Report.)

Due to the spherical geometry of the problem, we introduce spherical harmonics, i.e., multipolar polarization modes of the particle and electromagnetic fields. Introducing the potentials ϕ_1,2 and ψ_1,2 such that , and , one can look for solutions of the Laplace equations for potentials in terms of a superposition of multipolar modes which are obtained by expanding the angular dependency on the Legendre polynomials . Using spherical coordinates centered on the NP, with r = |r| as the radial distance, and θ as the polar angle, taking into account that the potentials should be regular at r = 0 and that for , the field has to reconnect to the external, uniform field E₀, the following expressions for ϕ_1,2 and ψ_1,2 are obtained:

Here, δ stands for the Kronecker symbol; , are the mode amplitudes for the electrical fields E_1,2; and are the mode amplitudes for the polarizations P₁ and Π. The quantity Π₀ is defined as the uniform component of Π, reflecting the polarization response (in the z-direction) arising from only the active elements in the absence of nanoparticle, when submitted to the probe field E₀.

We can now insert these expansions in system (9–11), to obtain the full set of time-dynamical evolution equations for Π₀, N₀, and as functions of the probe field E₀:

where , and the dielectric permittivity ε_m of the metal from Eq. 2 were used.

Finally, it is necessary to take boundary conditions at the NP surface (r = a) into account. The continuity condition for the tangential electrical field writes as , and has been already implemented to obtain Eq. (22). Continuity of the normal electrical displacement links to the variables Π₀, and , closing the system in the following way:

The set of Eqs (19, 20, 21, 22, 23) represents the first complete treatment of a spherical NP immersed in an active gain medium, taking into account the full array of multipolar modes. We are now in a position to consider spasing and amplification effects in the optical response of this system.

Spasing Threshold and Onset of Instability Under Zero External Field

In a similar way to classical lasing, spasing occurs as a self-oscillation of the electromagnetic field inside and outside the NP, in the absence of an external driving (probe field E₀ = 0), when the gain level in the system reaches a given threshold. The self-oscillation spasing condition is therefore defined as an unstable point.

This will be reflected by the appearance a non-zero outside field E₂, or in other words, by the appearance of (at least) one non-zero multipolar component . We therefore look for the condition for any given to be non-zero, independently of other modes (p_k = 0, ). For this, we assume steady-state in Eqs (19, 20, 21, 22) by canceling out all time derivatives, as well as zero external field (E₀ = 0), which directly leads to Π₀ = 0 and . Substituting steady-states values for Π₀, and into the boundary condition 23 brings the following equality:

which means that can be non-zero only if

This is the spasing condition for the -th mode. The imaginary part of the spasing condition yields the value for the gain threshold for the -th mode to start self-oscillating: it can be concluded that higher modes require more gain to be triggered. Therefore, the dipolar mode is the one that will be most easily activated:

and it is the dipolar threshold that sets the absolute spasing threshold for the NP.

Note that this is the same condition as discussed in ref. 11 for the singularity of the classical dipolar polarizability formula for a sphere (see also further down). We also have analytically checked that this spasing condition corresponds to the point where losses in the system are exactly balanced by the gain, as holds for the threshold of classical lasers.

We now continue with describing the first stages of the spasing instability, right at the onset of oscillations, when the gain value in the active medium exceeds the above threshold. Because the dipolar mode has the lowest threshold, it will initially have the fastest growing rate among all unstable modes for a given gain value, and hence is expected to dominate this onset.

We therefore consider a dipolar approximation to the set of equations (19, 20, 21, 22, 23), valid as long as all modes with are negligible. In this situation, and still assuming E₀ = 0, the system (19–23) reduces significantly:

Equations (27) and (28) have trivial solutions, and assuming physically meaningful initial conditions, where Π₀|_t=0 = 0 (due to the absence of externally imposed field) and (stable pumping), these variables do not evolve in time: Π₀(t) = 0, . One then only needs to solve the two remaining equations (29) and (30). Using the boundary condition (23) to express p₁, (29) and (30) form a linear system that can be solved analytically:

It appears that the temporal dynamics of the system is dictated by two eigenvalues :

where

and

The parameters are

and the constants c_1,2 are related to initial conditions.

The spasing instability (self-oscillation) will arise if at least one of these eigenvalue has a positive real part. Figure 2 shows that indeed the real part of becomes positive exactly when the gain exceeds the spasing threshold, i.e., when (calculated at the plasmon frequency, ω = ω₀). The real part of the other eigenvalue, , remains always negative.

Figure 2

Spasing threshold and eigenvalues and : becomes positive (instability) exactly when (dipolar spasing threshold).

Parameters: ε_b = 1.85 (ethanol solvent), ħω₂₁ = ħω₀ = 3.21 eV, ħΔ = 0.2 eV.

It is interesting also to consider the span of the unstable region with when the overall gain value varies: this is mapped in Fig. 3 for a 10 nm-size silver nanoparticle. It is at first quite narrow around the plasmon wavelength 2πc/ω₀ = 386 nm (with c the celerity of light) when the gain is close to the threshold, but then reaches spectral widths of a few tens of nanometers, as gain is further increased. This provides a rather large frequency range operation for the nanolaser; note, however, that this spectral range should not be confused with the actual linewidth emission of the lasing system: the latter should be calculated in the non-linear steady-state lasing regime, as the full width at half maximum of the emission intensity, and is expected to be much narrower.

Figure 3

Main figure: map of the stable and linearly unstable region (inside contour), as a function of frequency/wavelength and the overall gain level in the active medium .

Inset: instability region width Δλ as a function of overall gain. Parameters: ε_b = 1.85 (ethanol solvent), ħω₂₁ = ħω₀ = 3.21 eV, ħΔ = 0.2 eV.

Naturally, as the instability develops, the dipolar approximation used in this section will fail: consequences of this failure and later stages are discussed in Section 5 (“Above spasing threshold”) further down. We before consider a situation of practical importance, namely the behaviour of the system when excited by an external field (pump field), both below the spasing threshold and above it, again in the initial stages of the instability.

Response To An External Field - Linear Amplification Regime

Below spasing threshold

We now consider the case where a probe field is shone onto the nanoparticle. We begin by considering situations where the overall gain value remains below the spasing threshold. Since no field is expected to become extremely strong, we consider as before that .

This brings about two physically important consequences: (i) the system (19–22) is linear; (ii) only dipolar terms survive and all multipolar terms are extinct. The latter property can be inferred upon close inspection of the differential system: modes are coupled within pairs of the same order only [Eqs (20, 21, 22)]. Among such pairs, however, only the equations for σ₁ and q₁ exhibit non-zero source terms (proportional to E₀ and Π₀), stemming from the expression of p₁. Hence, after a very short time, all multipolar terms must naturally decay to zero, except dipolar ones which are fed through the above-mentioned source terms. Finally, it is only necessary to solve the equations for Π₀, and the two active modes σ₁ and q₁:

alongside the expression for p₁:

We first look for the (linear) steady states of the system under field; this is best represented by computing the steady-state dipolar polarizability, , of the sphere, which we classically define as the ratio of the dipolar moment p₁ of the NP to the exciting field²⁸: . To this purpose, time derivatives are canceled out in the above system, and with the help of (39), we find:

Unsurprisingly, this no other than the well-known quasi-static formula for the dipolar polarizability of a spherical inclusion²⁸. This simple formula, where the active nature of the host medium is described by letting ε_h have a negative imaginary part [as per Eq. (13)], was indeed studied by some of us in a previous work¹¹ as the simplest step towards loss compensation phenomena: this empirical approach is here revealed as the steady-state, linear approximation of our full time-dynamical model.

The evolution of the steady-state is shown in Fig. 4 for the same 10 nm-silver nanoparticle (see also ref. 11), where the gain level is progressively increased, i.e., goes from zero towards more and more negative values. (Note that we assumed a fixed, large power for the pump: ; The same behaviors can be obtained by fixing the value for and varying from 0 to 1). We observe that first, the plasmon (here observed as a resonance in the polarizability) is gradually amplified with increasing finesse (Fig. 4a–c): this is what we call the “linear amplification regime” where partial loss compensation occurs. This regime by itself, is a very interesting one in practical applications: it requires much less gain to be introduced, is more easily attained than the spasing one, and to a certain extent, is able to compensate the intrinsic losses found in natural plasmonic resonances.

Figure 4

Evolution of the steady-state plasmonic response of a 10-nm silver nanoparticle as gain is increased in the surrounding medium.

Parameters: ε_b = 1.85 (ethanol solvent), ω₂₁ = ω₀ = 3.21 eV, ħΔ = 0.2. (a–c,g–i) Steady-state dipolar polarizability ; (d–f,l–n) Real part of eigenvalues and , corresponding to each of the polarizability curves.

Then, when the gain value reaches the spasing condition [canceling the denominator in 40], is at its sharpest and formally becomes singular (Fig. 4g). This signals the failure of the used linear, dipolar approximation, and new effects, including saturation, will set in as discussed further down. In particular, the singularity shall disappear and all quantities shall remain finite.

Above spasing threshold

We now assume that the gain in the host medium is set above the spasing threshold, and that the system is still excited by an external probe field E₀.

Although it does not represent a fully valid solution anymore, it is worthwhile commenting on the behaviour of the steady-state polarizability above threshold: is finite again, but now exhibiting regions in the spectrum where the imaginary part is negative (Fig. 4g–i), a situation which cannot be observed for purely passive systems: this should here be considered as an indication of the appearance of the spasing regime.

The onset of the spasing regime under field can indeed be studied following the same analysis as in zero-field to demonstrate the onset of an oscillation (instability), by solving the linear system (36–38) analytically. One obtains:

where all parameters are as defined earlier [Eqs (32–35)].

Here again, under field, the spasing instability sets in when becomes positive (the other eigenvalues and remain always stabilizing), as shown in Fig. 4: for low gain values [(d) to (f)], ; then at the spasing threshold ((g)], reaches 0 exactly at the plasmon frequency; finally, from (l) to (n), there is an unstable frequency range around the plasmon frequency, where ; as before, the width of this region shall correspond to the spaser emission width.

A close inspection of Fig. 4(m,n) (zoomed-in detail not provided) reveals an interesting fact: sign changes in always happen extremely near to sign changes of , to a precision better than 0.01 eV. Hence, regions of spasing instability very closely coincide with regions with (“active” polarizability), and therefore, in the present system, an easy rule-of-thumb to predict the extent of the spasing regime is to look for sign changes in the imaginary part of the simple formula (40). Obviously, whenever the instability sets in, it prevents the linear steady-state solution represented by to settle, and the dynamics of the system will become increasingly non-linear as discussed further down. To illustrate the instability, in Fig. 5, we compare the evolution of the dipolar moment p₁(t), as calculated from the exact solution 41 for (as in Fig. 4c) and (as in Fig. 4h), both observed at the same point of the energy spectrum (ħω = 3.23 eV). For the former, (i.e., below spasing threshold) and therefore a steady-state (in the “amplifying regime”). For the latter, so that the spasing instability sets in with a linear growth rate , preventing the linear steady state to appear, and preparing for the coming of a full non-linear temporal dynamics as discussed in the next section.

Figure 5

Time evolution of the dipolar moment of the nanoparticle p₁(t) of a 10-nm silver nanoparticle.

The nanoparticle is under constant excitation from the probe field E₀, while the pump is switched on from t = 0 on. (a) Example of stable regime (amplifying regime, under spasing threshold). (b) Example of onset of the spasing regime (above threshold).

Non-Linear Regime and The Rise of Higher-Order Modes

We now describe qualitatively the later stages of the spasing instability explored in the previous sections.

At the onset of the spasing instability, the amplitude of the plasmonic dipolar mode grows with a typical rise time , both with and without excitation by the field E₀. For times , this regime is linear as described in the previous sections.

For times , since the dipolar mode has been growing exponentially (and so does E₂), the term will no more remain negligible in Eq. 10: the instability enters a regime of non-linear growth. Since both E₂ and Π are spatially non-uniform fields, Eq. 10 indicates that the population inversion N will also take on a non-uniform distribution. In order to account for this, the following expansion for N can be introduced:

where are the mode amplitudes of the population inversion and N₀ is the already presented, spatially uniform (flat) mode. Taking into account this expression transforms the original system (9–11) into a set of non-linear equations displaying intricate cross-coupling terms not present in Eqs (19–22).

The numerical solutions of this system will be explicitly discussed in a forthcoming paper devoted to studying the non-linear amplification regime. Qualitatively, the effect of these new non-linear terms is to couple together and -terms of different -orders via the appearance of the modes. As a consequence, electromagnetic energy will always cascade from the initially unique dipolar mode (which benefit from source terms), towards higher (sourceless) modes , thereby exciting multipoles in the particle.

Therefore we can conclude that, briefly after the spasing instability launches (a few nanoseconds in the example of Fig. 5), the response of the gain-assisted NP always becomes multipolar. Interestingly enough, this is true however small the NP is compared to the probe wavelength, due to local feedback effects with the gain medium. This result means that the response of the spaser, in the geometry studied in this article, is intrinsically a multi-modal one; other geometries for the NP can tailor the response to be exclusively dipolar, for example the “nanoshell” geometry²⁵.

Finally, for times , the population inversion will be sufficiently depleted that it will limit the growth of the instability, and a final state will be reached. The exact nature of this long-term state, as it is well-known from the physics of laser oscillations, can be diverse^24,33,34: steady-state, multistability, stochastic behaviour, burst-like dynamics, etc. Moreover, in the case where the system is subjected to a non-zero external external field E₀, synchronicity issues will rise between the natural resonant frequencies of the NP oscillator (here, the resonance frequencies of the multipoles) and the forcing frequency^25,35. The exact nature of this final state requires a thorough study, depending on the various model parameters, which is out of the scope of this paper.

In case the transient evolution results in a final stationary state, however, the steady-state values of N, Π and P₂ are found easily by canceling time-derivatives in, and solving simultaneously, Eqs (9), (10) and (12). This reveals that the saturated gain medium can be described through a non-linear permittivity^16,25 , as a function of the local field E₂:

with the saturation field .

The non-linear permittivity can then either be used inside numerical simulations¹⁶ to find the field distribution and other properties of the stationary lasing state, or equivalently, the associated non-linear, steady-state value of N can be introduced into our analytical equations (9, 10, 11) to obtain a multipolar mode description of the lasing state.

In conclusion, a full multipolar and dynamical model has been established to describe the optical response of a single metallic nanosphere immersed in an active medium. We have shown that for low gain values in the active medium, the response of the nanoparticle is a steady state and corresponds exactly to the results predicted via the quasi-static classical formulae for dipolar polarizability. This is a regime of amplification, where losses are partially compensated and the plasmonic response is amplified. Then, when gain values exceed a threshold, a spasing instability occurs, where physical quantities such as the dipolar moment of the NP, initially grow exponentially with time. This instability has been demonstrated both in presence, and in absence of an externally-imposed field (i.e., spasing self-oscillation). After a short while, a cascade of non-linear couplings always launches and activates higher-order, multipolar modes, irrespective of how the NP size compares to the exciting (probe) wavelength. Finally, features of the final, non-linear regime were discussed qualitatively.

The importance of geometry for NP-based spasers cannot be overstated, as the nanoparticle’s shape and composition determine the natural “resonator” in the spasing system. In our case of a homogeneous NP in an active medium, our model shows that the optical response of the resonator is unavoidably multi-modal and multipolar, therefore channeling energy into “dark modes” which may or may not suit one needs; then, depending on what type of application is envisioned, various geometries may be sought after, and optimised as a means to control lasing modes of the nano-emitter. Our analysis predicts where in the spectrum and with which specific gain quantity a mode will be activated, allowing in principle to design optical measures aiming at mode-specific amplification regimes, in a very simple geometry.