Abstract

Using either quasistatic approximation or exact Mie expansion, we characterize the localized surface plasmons supported by a metallic spherical nanoparticle. We estimate the quality factor and define the effective volume of the th mode in such a way that coupling strength with a neighbouring dipolar emitter is proportional to the ratio (Purcell factor). The role of Joule losses, far-field scattering, and mode confinement in the coupling mechanism is introduced and discussed with simple physical understanding, with particular attention paid to energy conservation.

1. Introduction

Metallic nanoparticles support localized surface plasmonpolaritons (SPP) strongly confined at the metal surface ensuring efficient electromagnetic coupling with neighbouring materials, offering a variety of applications such as surface-enhanced spectroscopies [1, 2], photochemistry [3], or optical nanoantennas [4]. This also opens the way towards control of light emission at the nanoscale [510]. The ratio is generally used to quantify the coupling strength between a dipolar emitter and a cavity mode (polariton). and refers to the mode quality factor and effective volume, respectively. High ratio may lead to strong coupling regime with characteristics Rabi oscillations revealing cycles of energy exchange between the emitter and the cavity. In the weak coupling regime, the emitter energy dissipation into the cavity mode is nonreversible and follows the Purcell factor: where is the emitter spontaneous decay rate into the cavity compared to its free-space value , is the optical index inside the cavity, and is the emission wavelength. Spontaneous emission rate in complex systems follows the more general Fermi’s golden rule that expresses as a function of the density of electromagnetic mode [11]. However, for describing the coupling to a cavity mode, the Purcell factor is usually preferred since it clearly introduces the cavity resonance quality factor and the mode extension on which coupling remains efficient, bringing therefore a clear physical understanding of the coupling process.

In this context, we propose to determine the quality factors and effective volumes of localized SPPs supported by a metallic nanosphere. Indeed, these quantities are useful parameters to understand and evaluate the coupling mechanisms between dipolar emitters and plasmonic nanostructures [12, 13]. This would help for achieving strong coupling regime [14, 15] or for designing plasmonic nanolasers [1618]. We have chosen a spherical particle, a highly symmetrical system, since it is fully analytical and simple expressions can be derived with clear physical meaning [1923]. Moreover, our results could be extended to more complex structures [24, 25].

In Section 2, we focus on the dipolar mode and present in details the derivation of its quality factor and effective volume. We extend our approach to each mode of the particle in Section 3. Finally, we discuss the coupling efficiency to one of the particle modes in the last section. For the sake of simplicity, all the analytical expressions are derived assuming a particle in air and dipolar emitter perpendicular to the particle surface. The generalisation to arbitrary emitter orientation and a background medium of optical index is provided in the appendix. Exact calculations presented in the main text using Mie expansion correctly include both the background medium and dipolar orientation.

2. Dipolar Mode

We first characterize the dipolar mode of a spherical particle. For sphere radius small compared to the excitation wavelength , the electric field is considered uniform over the metallic particle. The metallic particle is polarized by the incident electric field and behaves as a dipole where is the nanoparticle quasistatic (dipolar) polarisability and is the metal dielectric constant. The dipole plasmon resonance appears at such that . In case of Drude metal, the dipolar resonance is with the bulk metal plasma angular frequency. However, expression (3) does not satisfy the optical theorem (energy conservation). It is wellknown that this apparent paradox is easily overcome by taking into account the finite size of the particle and leads to define the effective polarisability [26]: The corrective term () is the so-called radiative reaction correction and microscopically originates from the radiation emitted by the charges oscillations induced inside the nanoparticle by the excitation field [27].

The dipolar polarisability presents a simple shape near the resonance if the metallic dielectric constant follows Drude model [20] ( and refers to metal plasma frequency an Ohmic loss rate, resp.) (for the sake of clarity, we use lowercase notation for the emitter decay rates and uppercase for the quantities associated with the metallic particle); where is the decay rate of the particle dipolar mode and includes both the Joule () and radiative [] losses rates.

2.1. Quality Factor

The dipolar response can be described by either the extinction efficiency , proportional to , or scattering efficiency , proportional to , and therefore follows a Lorentzian profile centered at and with a full width at half maximum (FWHM) : The quality factor of this resonance is therefore

As an example, we consider a silver sphere in air. The Drude model parameters are and that lead to a resonance peak at () and radiative energy . The quality factor is then . Note that the metal optical properties are better described when including the bound electrons in the Drude model: , ( for silver, see also the appendix). In that case, we obtain (), , and .

Although Drude model qualitatively explains the shape of the resonance, a more representative value of the quality factor can only be determined using tabulated data for the dielectric constant of the metal [28]. The extinction efficiencies associated with the dipolar resonance are represented in Figure 1. For silver particle (Figure 1(a)), the extinction efficiency closely follows a Lorentzian shape, as expected, with a quality factor , in good agreement with the value obtained using Drude model (with the contribution of the bound electrons included). In case of gold, the resonance profile is not well defined due to interband transitions for but a quality factor of 7 can be estimated (Figure 1(b)). Figure 1(c) represents the extinction efficiency for a gold particle embedded in polymethylmethacrylate (PMMA) into which metallic particles are routinely dispersed. This leads to a small redshift of a resonance, avoiding therefore the resonance disturbance by interband absorption, and we recover partly the Lorentzian profile [29].

2.2. Effective Volume

Coupling rate of a dipolar emitter to a dipolar particle expresses for very short emitter-particle distances ( is the distance to the particle center) as [1921] for a dipole emitter orientation perpendicular to the nanoparticle surface. Using (6), we obtain, in case of a dipolar emitter emission tuned to the dipolar particle resonance In order to determine the dipolar mode effective volume, we now identify the coupling rate to the Purcell factor ((1), assuming ), so that we obtain in full agreement with expression recently derived by Greffet et al. by defining the optical impedance of the nanoparticle antenna [30]. For a nm radius sphere in air, we estimate the mode effective volume for an emitter 10 nm away from the particle surface ().

Unlike here, usual definition of the mode effective volume does not include the emitter position. For instance, in case of cavity quantum electrodynamics (cQED) applications, it can be expressed as , so that it directly characterizes the mode extension. However, in that case, Purcell factor expression (1) is only valid for an emitter located at the cavity center where the mode-emitter coupling is maximum (mode antinode). In close analogy with the definition used in cQED, we derive later another expression for the SPP effective volume, independent of the emitter position (see (22), and the discussion below).

3. Multipolar Modes

If the excitation field is generated by a dipolar emitter (fluorescent molecules, quantum dots,…), it cannot be considered uniform anymore and the dipolar approximation fails. One needs therefore to consider the coupling strength to high-order modes (Figure 2). Here again, we discuss the mode quality factors and volume first using quasistatic approximation and then discuss their quantitative behaviour using exact Mie theory.

The th multipole tensor moment of the metallic particle is given by with and is the vector differential operator. As discussed above, the dipole moment is the unique mode excited in an uniform field. However, the dipolar emitter near field behaves as and strongly varies spatially so that higher modes can be excited in the particle.

For a Drude metal, the th resonance appears at . Therefore, higher-order modes accumulate near . Moreover, as discussed for the dipolar case, quasistatic expression of the th mode polarisability (14) does not obey energy conservation. Applying the optical theorem, we recently extend the radiative correction to all the spherical particles modes [22]. This leads to

3.1. Quality Factors

It is now a simple matter to generalize the dipolar mode analysis reported in the previous section to all the particle modes. The th mode polarisability can be approximated near resonance. A simple expression for the polarisability is achieved considering a Drude metal: where is the total decay rate of the th mode, that includes both ohmic losses and radiative scattering. As expected, for a given mode , the radiative scattering rate increases with the particle size since it couples more efficiently to the far field. For instance, we obtain ( with ) for the quadrupolar mode of a silver particle in air. As expected, the radiative rate of the quadrupolar is strongly reduced compared to the dipolar mode.

The quality factor associated with the th mode is therefore

Figure 3 details the quality factor of the two first modes of a silver sphere in PMMA using tabulated value for . The quadrupolar mode presents a quality factor almost 5 times higher than the dipolar mode since it has limited radiative losses. Indeed, quadrupolar mode poorly couples to the far field. Finally, similar values () are obtained for all the higher modes (). Here again, assuming a Drude metal and in the quasistatic approximation, we qualitatively explains this result. Actually, Q factors of high-order modes are absorption loss limited (, so called dark modes) and tend to . However this overestimates the mode quality factor () as compared to the value deduced using tabulated value for and exact Mie expansion.

3.2. Effective Volumes

Last, we express the total coupling strength of a dipolar emitter to the spherical metallic particle for very short separation distances [21, 22]: so that the coupling strength to the th mode is easily deduced as where we have used approximated expression (16) for the th polarisability. The mode effective volume is then straightforwardly deduced from comparison to the Purcell factor  (1) We obtain for instance quadrupole effective volume for an emitter 10 nm away from a 25 nm radius metallic particle in air. The calculated effective volumes of the three first modes are shown on Figure 4. At long distance, the smaller effective volume (most efficient coupling to an emitter) is associated with the dipole mode since it presents the largest extension. When the emitter-particle distance decreases the coupling strength to the quadrupolar (), then hexapolar () mode becomes stronger as revealed by their lower effective volume. Most efficient coupling between an emitter and one of the particle modes obviously occurs at contact since plasmons modes are confined near the particle surface.

Finally, it is worthwhile to note that mode effective volume is generally defined independently on the particle-emitter distance so that it gives an estimation of the mode extension. This is done by evaluating the maximum effective volume available and for a random emitter orientation. In case of localized SPP, this is achieved for contact (). Since the decay rate of a randomly oriented emitter expresses , it comes where is the metallic sphere volume and refers to a dipole parallel to the sphere surface (see (B.4) in Appendix B, with ). Dipolar and quadrupolar mode volumes are and , respectively. The expression (22), derived for a Drude metal, quantifies an extremely important property of localized SPPs; their effective volume does not depend on the wavelength and is of the order of the particle volume. Metallic nanoparticles therefore support localized modes of strongly subwavelength extension. Highest-order modes have negligible extension ( for , see Figure 5). As a comparison, photonics cavity modes are generally limited by the diffraction limit so that their effective volume is at best of the order of (see [31]). Nevertheless, the extremely reduced SPP volume is achieved at the expense of the mode quality factor. It should be noticed that a low-quality factor indicates a large cavity resonance FWHM so that emitter-SPP coupling can occur on a large spectrum range.

Note that mode effective volume is generally defined by its energy confinement [12]. Khurgin and Sun used this definition and obtained [32] that is in agreement with our expression for the dipolar mode volume () but leads to slightly different values for other modes (e.g., instead of ). Recently, Koenderink showed that defining the mode volume on the basis of energy density could lead to underestimating the Purcell factor near plasmonic nanostructures [33] (however, he defined the coupling rate to the whole system rather than considering the coupling into a single mode). Oppositely, we adopt here a phenomenological approach where the mode volume is defined so that Purcell factor remains valid. Nevertheless, both methods lead to very similar results for the effective volumes of localized SPPs supported by a nanosphere. Since mode volume is generally a simple way to qualitatively characterize the mode extension, expressions derived here or by Sun et al. could be used.

4. -Factor

Purcell factor quantifies the coupling strength between a quantum emitter and a (plasmon) mode but lacks information on the coupling efficiency as compared to all the other emitter relaxation channels. For the sake of clarity, it has to be mentioned that coupling strength to one SPP mode corresponds to the total emission decay rate induced by this coupling. It does not permit to distinguish radiative () and nonradiative () coupling into a single mode. This could be done by numerically cancelling the imaginary part of the metal dielectric function (see also [34] for similar discussion in case of coupling to delocalized SPPs). However, one can determine the coupling efficiency into a single mode as compared to all other modes. This coupling efficiency, or the so-called -factor, is easily estimated in case of a spherical metallic particle since all the available channels are taken into account in the Mie expansion. Coupling efficiency into th mode writes where is the decay rate of a randomly oriented molecule. -factor is represented on Figure 6 for the first three modes. A maximum efficiency of can be achieved in the dipolar mode () and into quadrupolar mode (). The coupling efficiency into the hexapolar mode is lower ( around ) since it has a very low extension, as indicated in Figure 5. For very short distances, all the coupling efficiencies drop down to zero since all the higher-order modes accumulate in this region, opening numerous alternative decay channels. Moreover, it is possible to efficiently couple the emitter to either the dipolar or quadrupolar mode, by matching the emitter and mode wavelengths. This is of strong importance when designing a SPASER (or plasmon laser) [5] so that the active mode can be tuned on the dipolar or quadrupolar mode. This last spasing mode would consist of an extremely localized and ultrafast nanosource [18].

5. Conclusion

We explicitly determined the effective volumes for all the SPP modes supported by a metallic nanosphere. Their quality factor is also approximated in the quasistatic case or calculated using exact Mie expansion. Rather low-quality factors ranging from 10 to 100 can be achieved, associated with extremely confined effective volume of nanometric dimensions. This results in high ratios indicating efficient coupling strength with a quantum emitter. At the opposite to cavity quantum electrodynamics where the coupling strength is obtained on diffraction limited volumes thanks to ultrahigh factors, plasmonics structures allow efficient coupling on nanometric scales with reasonable factors, defining therefore high-bandwidth interaction. Finally, high coupling efficiencies (80%–90%) to dipolar or quadrupolar mode can be achieved and are of great interest for nanolasers realization.

Appendices

In this appendix, we derive the expression of mode quality factor and effective volume for a spherical particle embedded in homogeneous background of optical index . For a better description of the metal optical properties, we include the contributions of the bound electrons into the metal dielectric constant .

A. Effective Polarisability Associated to th SPP Mode

The th multipole tensor moment of the metallic particle is given by The resonance angular frequency is then (). Finally, the effective polarisability, including finite size effects writes [22] with the wavenumber in the background medium.

Considering a Drude metal, we achieve a simple approximated expression for near a resonance so that the quality factor expression remains valid but with the corrected expressions for resonance frequency and total dissipation rate .

B. Effective Volumes

The total coupling strength of a dipolar emitter to the spherical metallic particle, embedded in medium, expresses [21, 22] The coupling strength to the th mode is and the mode effective volume is deduced from comparison to the Purcell factor  (1) Finally, we define the mode volume as All the expressions obtained in this appendix reduce to the simple analytical case discussed in the main text for (bound electrons contribution neglected) and (background medium is air).

Acknowledgments

This work is supported by the Agence Nationale de la Recherche (ANR) under Grants Plastips (ANR-09-BLAN-0049-01), Fenoptiχs (ANR-09-NANO-23), and HYNNA (ANR-10-BLAN-1016). This paper was also motivated by the preparation of a lecture for the 2nd Summer School On Plasmonics organized by Nicolas Bonod in October 2011.