Single plasmon hot carrier generation in metallic nanoparticles

TL;DR

This paper develops a quantum-mechanical approach to describe hot carrier generation from single plasmon decay in metallic nanoparticles, revealing significant differences from classical models and implications for ultrasmall device design.

Contribution

It introduces a novel quantum-specific method for modeling plasmon decay into hot carriers, advancing understanding beyond semiclassical theories.

Findings

Hot carrier generation rates differ from semiclassical predictions.

Decay of non-plasmonic excitations also produces comparable hot carriers.

Provides a foundation for designing ultrasmall plasmonic devices.

Abstract

Hot carriers produced from the decay of localized surface plasmons in metallic nanoparticles are intensely studied because of their optoelectronic, photovoltaic and photocatalytic applications. From a classical perspective, plasmons are coherent oscillations of the electrons in the nanoparticle, but their quantized nature comes to the fore in the novel field of quantum plasmonics. In this work, we introduce a quantum-mechanical material-specific approach for describing the decay of single quantized plasmons into hot electrons and holes. We find that hot carrier generation rates differ significantly from semiclassical predictions. We also investigate the decay of excitations without plasmonic character and show that their hot carrier rates are comparable to those from the decay of plasmonic excitations for small nanoparticles. Our study provides a rigorous and general foundation for further development of plasmonic hot carrier studies in the plasmonic regime required for the design of ultrasmall devices.

Figures

Tables

Table 1. Table 1: Excitation energies, total hot-carrier generation rates and transition dipole moments μ I subscript 𝜇 𝐼 \mu_{I} of the dominant neutral excitations in small sodium nanoparticles. SP refers to electron-hole pair excitations and PL to plasmonic excitations.

System	Excitation	$\hslash {\omega }_I\left(eV\right)$	$N_{\mathrm{tot}}$ ($1/ps$)	${\mu }_I$ (${10}^{-20}$ C nm)
Na40	PL	3.63	27.2	5.6
Na58	PL1	3.67	8.2	6.8
Na58	PL2	3.21	35.2	6.3
Na92	PL	3.07	25.7	9.1
Na92	SP	2.87	30.3	10.3
Na92	Classical PL	3.40	976.8	57.7

Equations

$${\hat{H}}_{\mathrm{eff}}=\sum _i{ϵ}_i{\hat{c}}_i^{†}{\hat{c}}_i+\sum _I\hslash {\omega }_I{\hat{b}}_I^{†}{\hat{b}}_I+\sum _{i,j,I}{g}_{ij}^I{\hat{c}}_i^{†}{\hat{c}}_j({\hat{b}}_I+{\hat{b}}_I^{†}),$$

$$⟨{\varphi }_m\left|{\Sigma }^{\mathrm{GW}}(\omega )\right|{\varphi }_m⟩=\sum _{j,I}\frac{|V_{mj}^I{|}^2}{\hslash \omega -{ϵ}_j-\hslash {\omega }_I{\eta }_j-i\delta },$$

$$V_{mj}^I=\int d\mathbf{r}\int d{\mathbf{r}}^{\prime }{\varphi }_m(\mathbf{r}){\varphi }_j(\mathbf{r})\frac{e^2}{|\mathbf{r}-{\mathbf{r}}^{\prime }|}{\rho }_I({\mathbf{r}}^{\prime })$$

$${\rho }_{\mathrm{I}}(\mathbf{r})=\sum _{vc}{F}_{\mathrm{vc}}^I{\left(\frac{{ϵ}_c-{ϵ}_v}{\hslash {\omega }_I}\right)}^{1/2}{\varphi }_{\mathrm{v}}(\mathbf{r}){\varphi }_{\mathrm{c}}(\mathbf{r})$$

$$⟨{\varphi }_m\left|{\Sigma }^{\mathrm{eff}}(\omega )\right|{\varphi }_m⟩=\sum _{j,I}\frac{|g_{mj}^I{|}^2}{\hslash \omega -{ϵ}_j-\hslash {\omega }_I{\eta }_j-i\delta }.$$

$$g_{mj}^I=V_{mj}^I.$$

$${\Gamma }_I=\frac{2\pi }{\hslash }\sum _{vc}|g_{vc}^I{|}^2\delta ({ϵ}_c-{ϵ}_v-\hslash {\omega }_I).$$

$$N_e^I(E)=\frac{2\pi }{\hslash }\sum _{vc}|g_{vc}^I{|}^2\delta ({ϵ}_c-{ϵ}_v-\hslash {\omega }_I)\delta (E-{ϵ}_c).$$

$$N_{\mathrm{tot}}^I={\int }_{E_F+\delta E}^{\infty }dE{N}_e^I(E),$$

$$\Phi_{sc}(\mathbf{r},\omega)=E_{0}\frac{\epsilon(\omega)-1}{\epsilon(\omega)+2}\left\{\begin{tabular}[]{c}$r\cos\theta$ for $r\leq R$,\\ $R^{3}\frac{\cos\theta}{r^{2}}$ for $r>R$.\end{tabular}\right.$$

$$ϵ(\omega )=1-\frac{{\omega }_0^2}{{\omega }^2+i\omega {\gamma }_P/\hslash },$$

$$g_{ij}^{SC}=⟨{\varphi }_i|{\Phi }_{sc}|{\varphi }_j⟩.$$

$$E_0^{1PL}=\frac{{\gamma }_P}{{\mu }_P},$$

$$\left[-\frac{{\hslash }^2}{2m}{\nabla }^2+V_{KS}(\mathbf{r})\right]{\varphi }_i(\mathbf{r})={ϵ}_{\mathrm{i}}{\varphi }_i(\mathbf{r}),$$

$$\sum _{v^{\prime }{c}^{\prime }}{\Omega }_{vc,v^{\prime }{c}^{\prime }}{F}_{v^{\prime }{c}^{\prime }}^{\mathrm{I}}={\hslash }^2{\omega }_{\mathrm{I}}^2{F}_{vc}^{\mathrm{I}},$$

$${\Omega }_{vc,v^{\prime }{c}^{\prime }}={\delta }_{\mathrm{c}{\mathrm{c}}^{\prime }}{\delta }_{\mathrm{v}{\mathrm{v}}^{\prime }}({ϵ}_{\mathrm{c}}-{ϵ}_{\mathrm{v}}{)}^2+2\sqrt{f_{vc}({ϵ}_{\mathrm{c}}-{ϵ}_{\mathrm{v}})}{K}_{vc,v^{\prime }{c}^{\prime }}\sqrt{f_{v^{\prime }{c}^{\prime }}({ϵ}_{{\mathrm{c}}^{\prime }}-{ϵ}_{{\mathrm{v}}^{\prime }})},$$

$$K_{vc,v^{\prime }{c}^{\prime }}=\int d\mathbf{r}\int d{\mathbf{r}}^{\prime }{\varphi }_{\mathrm{v}}(\mathbf{r}){\varphi }_{\mathrm{c}}(\mathbf{r})\frac{e^2}{|\mathbf{r}-{\mathbf{r}}^{\prime }|}{\varphi }_{{\mathrm{v}}^{\prime }}({\mathbf{r}}^{\prime }){\varphi }_{{\mathrm{c}}^{\prime }}({\mathbf{r}}^{\prime }).$$

$$f_{\mathrm{I}}=\frac{2m}{{\hslash }^2}{\left|\sum _{vc}\sqrt{f_{vc}({ϵ}_c-{ϵ}_v}){\mu }_{vc}{F}_{vc}^I\right|}^2\equiv \frac{2m}{\hslash }{\omega }_I{\mu }_I^2,$$

Full text

{affiliations}

Department of Physics, Imperial College London, London SW7 2AZ

Departments of Physics and Materials, Imperial College London, London SW7 2AZ

The Thomas Young Centre for Theory and Simulation of Materials

1 Introduction

Localized surface plasmons (LSPs) in metallic nanoparticles facilitate drastic electric field enhancements and large light absorption cross sections that can be harnessed in nanophotonic applications, such as plasmon-enhanced biosensing [1], surface-enhanced Raman scattering [2], data storage [3] or nanoheaters [4, 5]. Recently, there has been significant interest in the decay of the LSPs into electrons and holes. The resulting carriers are energetic or ”hot” and can be used in solar energy conversion applications, including solar cells [6, 7] or photocatalysts [8, 9, 10]. For example, Mukherjee and coworkers demonstrated that hot electrons can induce challenging chemical reactions, such as the dissociation of hydrogen molecules on gold surfaces [11]. Moreover, the fast decay of LSPs can be used for new quantum information devices and in nanocircuitry [12, 13].

To provide insight and guidance in this rapidly evolving field, a detailed theoretical understanding of hot electron processes, including plasmon decay, hot carrier thermalization and recombination dynamics, is needed. Using semiclassical approaches, which combine a classical description of the LSP with a quantum-mechanical description of hot carriers, several groups analyzed the distribution of hot carriers resulting from the plasmon decay and studied its dependence on the nanoparticle size, material and environment [14, 15, 16, 17]. Providing general insight, the semiclassical approach is frequently based on a bulk dielectric function (such as a Drude model) and therefore cannot be used to describe small nanoparticles where quantum confinement effects play an important role [18, 19]. In addition, it has been observed that non-plasmonic excitations, such as electron-hole pairs, can take place at similar energies as the plasmon resonance, but such excitations are not described accurately on the basis of a the semiclassical approach [20, 21, 22]. Finally, a classical description of plasmons is less appropriate in the limit of low plasmon densities, where their quantized character must be captured [23, 24, 25]. Other groups have employed first-principles real-time time-dependent density-functional theory (TDDFT) to study plasmon [26, 27] and hot carrier properties in small metallic nanoparticles [28, 29, 21, 30]. While this method allows the study of nonlinear properties and scales favorably with the system size, it does not include a quantized treatment of the plasmon. Few attempts have been made to describe the effect of electron-plasmon interactions using quantized plasmons in metallic nanoparticles. Notably, Gerchikov and coworkers [31] and Weick et al. [32] used a separation of centre-of-mass motion and relative motion of the electrons to derive a quantized electron-plasmon Hamiltonian. However, their approach cannot be used to study the decay of other neutral excitations, such as electron-hole pairs.

In this paper, we thus present a fully quantum-mechanical approach to calculating the properties of hot carriers resulting from the decay of neutral excitations, such as LSPs or electron-hole pairs. In particular, we employ an effective Hamiltonian which describes the interaction of fermionic quasiparticles with bosonic neutral excitations and determine its parameters, including quasiparticle and plasmon energies and electron-plasmon coupling constants, using quantum-mechanical calculations. Most importantly, the electron-plasmon coupling strength is derived by comparing the electronic self energy of the effective Hamiltonian with the first-principles self energy within the GW approximation (where the electron self energy is approximated as the product of the electron Green’s function G and the screened Coulomb interaction W). After identifying the dominant plasmonic and non-plasmonic neutral excitations, we first calculate the hot carrier generation rates for spherical nanoparticles with different radii and study the effect of quantum confinement on the hot carrier distributions. We find that a larger fraction of the plasmon energy is distributed to the electrons rather than the holes. Secondly, we compare hot carrier rates from the decay of plasmonic and non-plasmonic excitations. We also compare our quantum-mechanical results with semiclassical calculations and show that there is a significant discrepancy in the hot carriers rates for small nanoparticles.

2 Results

2.1 Electron-plasmon coupling

To describe the interaction between charged fermionic quasiparticles and neutral bosonic excitations, such as plasmons and electron-hole pairs, in metallic nanostructures, we employ the following effective Hamiltonian

[TABLE]

where $\hat{b}_{I}^{{\dagger}}$ ($\hat{b}_{I}$) creates (annihilates) a neutral excitation in state $I$ with energy $\hbar\omega_{I}$ and $\hat{c}_{i}^{\dagger}$ ($\hat{c}_{i}$) creates (destroys) a quasiparticle in state $i$ with energy $\epsilon_{i}$. Also, $g_{ij}^{I}$ is the electron-plasmon coupling, i.e. the matrix element that describes the scattering of quasiparticles from state $i$ into state $j$ via the emission or absorption of a collective excitation in state $I$. See Supplementary note 1 for the definition of the single plasmon operator within the second quantisation formalism.

A key challenge in using $H_{\rm eff}$ to describe metallic nanoparticles is the accurate determination of the various parameters, including energies of neutral and quasiparticle excitations and their coupling. Quasiparticle energies are formally defined as the poles of the one-electron Green’s function and can be measured in photoemission and inverse photoemission experiments[33]. Calculating the Green’s function, for example via the GW method, is very challenging for metallic nanoparticles [34] and therefore quasiparticle energies are often approximated using Kohn-Sham energies obtained from density-functional theory (DFT)[35]. Similarly, energies of neutral excitations are defined as poles of a two-particle Green’s function and can be obtained by solving the Bethe-Salpeter equation or from time-dependent density-functional theory (TDDFT). Plasmons arise because of the long-ranged nature of the Coulomb interaction. This is captured by the Hartree contribution to the total energy [36], while exchange-correlation effects often play a minor role for plasmon properties. Neglecting exchange-correlation effects results in the well-known random-phase approximation (RPA) for neutral excitations[37].

To determine the electron-plasmon coupling, we follow Lundqvist’s approach [38] for the homogeneous electron gas and compare the second-order electron self-energy of $H_{\rm eff}$ with the correlation contribution of the ab initio GW self energy, see Figure 1a. The latter is given by [39]

[TABLE]

where the coefficient $\eta_{j}$ has the value $+1$ for unoccupied orbitals and $-1$ for occupied orbitals, $\delta$ represents a positive infinitesimal and $V^{I}_{mj}$ denotes the fluctuation potential given by

[TABLE]

with the transition density

[TABLE]

characterizing the $I$-th neutral excitation. Note that $\phi_{v}(\mathbf{r})$ ($\phi_{c}(\mathbf{r})$) denote single-particle wavefunctions of occupied (empty) states and the coefficients $F^{I}_{vc}$ are obtained by solving the Casida equation, see Methods section.

For the second-order electron self energy of $\hat{H}_{\rm eff}$, we find [40]

[TABLE]

Comparing Eq. (2) and Eq. (5) yields the final expression for the electron-plasmon coupling which is given by

[TABLE]

This result has the intuitive interpretation that the coupling between the neutral excitation $I$ and the quasiparticles states $m$ and $j$ is due to the electric potential induced by the transition density of the neutral excitation, see Eq. (3). Note that the expression is completely general and does not depend on the dimensionality, material or size of the system under consideration.

2.2 Hot carrier generation

To model the decay of a neutral excitations, such as a plasmon, into hot electrons and holes, we calculate the self energy of the neutral excitation and evaluate the Feynman diagram shown in Figure 1(b). The resulting decay rate $\Gamma_{I}$ is given by [15]

[TABLE]

In our numerical calculations, the delta function was replaced by a Gaussian with a standard deviation of 0.12 eV.

The generation rate $N^{I}_{e}(E)$ of hot electrons with energy $E$ from the decay of the $I$-th neutral excitation is given by

[TABLE]

A similar formula describes the generation rate of hot holes $N^{I}_{h}(E)$. In our numerical calculations, the second delta function was replaced by a Gaussian with a standard deviation of 0.05 eV.

Finally, the total generation rate $N^{I}_{\rm tot}$ of hot electrons resulting from the decay of the $I$-th collective excitation is given by[14]

[TABLE]

where $E_{F}$ denotes the Fermi energy which for metallic nanoparticles is defined as the middle of the gap between the highest occupied state and the lowest unoccupied one and $\delta E$ is a threshold energy which is typically chosen larger than the available thermal energy. As we are dealing with small nanoparticles, the energy gaps between occupied and unoccupied states are always larger than the thermal energy at room temperature and we therefore set $\delta E$ to zero. Note that the total rate of excited electrons is equal to the total rate of excited holes.

2.3 Semiclassical approach

To model plasmon decay and hot carrier generation in nanoplasmonic systems, a semiclassical approach is usually employed [15, 14]. In this approach, the metallic nanoparticle of radius $R$ is assumed to be exposed to an incident electric field (along the z-direction) with strength $E_{0}$. The total potential due to the perturbing field and the induced polarization in the nanoparticle is calculated using the quasistatic approximation according to

[TABLE]

Here, $\epsilon(\omega)$ denotes the dielectric function of the bulk material which is often described using a Drude model

[TABLE]

where $\gamma_{P}$ denotes the plasmon width and $\omega_{0}$ is the bulk plasmon frequency $\omega_{0}=\sqrt{4\pi ne^{2}/m}$ where $e$, $m$ and $n$ denote the electron charge, mass and density, respectively. Improved results can be obtained by using non-local approximations to the dielectric function [41, 42].

The hot-carrier generation rates are then obtained by evaluating Fermi’s Golden Rule in Eq. (7) with the electron-plasmon coupling strength

[TABLE]

Importantly, the semiclassical expression for the electron-plasmon coupling depends on the strength of the incident light field which excites the plasmon. In contrast, the quantum definition, Eq. (6), does not depend on $E_{\rm 0}$. To resolve this discrepancy, we note that the semiclassical result contains the interaction of the electrons with both the induced charge density and the external light field, while the quantum result only describes the interaction with the (light induced) neutral excitations. Moreover, the semiclassical result captures the effect of exciting multiple plasmons (in fact, the expression becomes exact in the limit of a large number of plasmon quanta), while the quantum result describes the decay of a single plasmon excitation.

In order to meaningfully compare the results of the two approaches, we derive an expression for the electric field strength $E^{1PL}_{0}$ which is required to excite a single plasmon, see Supplementary note 2. We find that

[TABLE]

where $\mu_{P}$ denotes the dipole moment of the plasmonic state. The corresponding semiclassical transition dipole moment is given by $\mu_{SC}(\omega)=R^{3}\left(\frac{\epsilon(\omega)-1}{\epsilon(\omega)+2}\right)E^{1PL}_{0}$. When evaluated at the classical Mie frequency $\omega_{cl}=\omega_{0}/\sqrt{3}$, we obtain $\mu_{SC}(\omega_{cl})=R^{3}\hbar\omega_{cl}/\mu_{P}$, which is independent of the plasmon linewidth $\gamma_{P}$.

3 Numerical results

We present results for three sodium nanoparticles: Na40, Na58 and Na92, which consist of 40, 58 and 92 Na atoms, respectively. These systems are modelled as jellium spheres with diameters of 1.4 nm, 1.6 nm and 1.9 nm, respectively. We first analyze the neutral excitations of these systems and identify those with plasmonic character and then calculate the hot carrier distributions resulting from their decay. We also compare our results to semiclassical calculations and obtain hot carrier rates for nanoparticles in different dielectric environments.

3.1 Identifying plasmonic excitations

Distinguishing plasmon-like excitations, which are a typically thought of as collective oscillations of all electrons in the nanoparticle, from other neutral excitations, such as bound or unbound electron-hole pairs, is difficult. Recently, several approaches have been proposed to address this problem, for example by scaling the electron-electron interactions[43, 44], by studying the time-dependent occupations of the Kohn-Sham states[45] or by investigating the optical response of the nanostructure[46, 47] that led to the introduction of the generalized plasmonicity index. Below, we present an alternative method for identifying the plasmon which is based on a graphical analysis of various physically transparent quantities (such as the excitation energy and oscillator strength of a neutral excitations) which are directly obtained from a solution of Casida’s equation. Importantly, this approach does not require multiple solutions of Casida’s equation at different strengths of the electron-electron interaction and shares many advantages of the generalized plasmonicity index, which is proportional to the total induced dipole moment of a transition.

Importantly, our quantum method for calculating hot-carrier rates can be used for both plasmon-like and electron-hole pair-like excitations. However, to enable a comparison with semiclassical calculations, we have employed the following approach to identify excitations with a plasmonic character. In particular, we analyze three properties of neutral excitations: (i) the energy $\hbar\omega_{I}$ of the excitation, (ii) its oscillator strength $f_{I}$ and (iii) its collectivity $C_{I}$, see Methods section. For plasmon-like excitations, we expect both high oscillator strengths and collectivities as well as energies not too far from the classical plasmon energy $\hbar\omega_{cl}$. For sodium nanoparticles, the classical plasmon energy is given by $\hbar\omega_{cl}=3.40$ eV.

Figure 2 a)-c) summarizes the properties of neutral excitations for the three sodium nanoparticles under consideration. In these graphs, each circle represents an excited state with its radius being proportional to $C_{I}$. For Na40, we find one state (denoted PL for plasmon in Figure2a) with a high oscillator strength, large collectivity and energy close to $\hbar\omega_{cl}$. This state gives rise to a strong peak in the absorption spectrum (see inset) and its transition density has a dipolar shape, see Figure 3. We therefore assign this state a plasmonic character. For Na58, we find two states with large oscillator strengths and high collectivities. The energies of both states are within 0.3 eV from $\hbar\omega_{cl}$ and their transition densities are both plasmon-like. We therefore assign both excitations plasmonic character (and denote them by PL1 and PL2 in in Figure2b)). For Na92, we find two states with high oscillator strengths. The state with energy $\hbar\omega_{SP}=2.87$ eV has a low collectivity and therefore electron-hole pair character (denoted SP for single pair in in Figure 2c)), while the state with energy $\hbar\omega_{P}=3.07$ eV has a high collectivity and plasmonic character (denoted PL in the figure).

Table 1 shows an overall redshift of the plasmonic excitation energy as the nanoparticle size increases. This is caused by quantum confinement effects and the spill-out of the electron wave functions beyond the geometrical radius of the nanoparticle [32]. However, the redshift is not monotonic which has also been observed experimentally for nanoparticles with a diameter smaller than 3.5 nm [48, 49] and was predicted theoretically for small nanoparticles [32]. Moreover, the presence of more than one plasmonic resonance for a given nanoparticle could explain the experimentally observed scattered distribution of energies in this size regime [48].

Figure 3 shows the transition densities of the plasmonic excitations in Na40, Na58 and Na92. The transition densities exhibit a clear dipolar shape, but are significantly more complex than predicted by classical electrodynamics. In particular, the transition densities are not only localized at the surface of the nanoparticles and show multiple regions of positive and negative charge. While the transition density of the smallest system (Na40) is extended throughout the nanoparticle, it becomes more localized in the surface regions for the larger nanoparticles in agreement with the classical result. In contrast, the transition densities of single-pair excitations can exhibit different features. For example, the transition density of the dominant single-pair state of Na92, see inset of Figure 4, is localized in the centre of the nanoparticle.

3.2 Hot carrier distributions

Figure 5 shows the energy distribution of electrons and holes resulting from the decay of a plasmon in the three sodium nanoparticles under consideration. Energy conservation imposes that the distribution of electrons (blue curve) has the same shape as the distribution of holes (red curve), but shifted by the plasmon energy. For the smaller nanoparticles, Na40 and Na58, the distributions exhibit only a few peaks. This is a consequence of the large energy level spacing (due to quantum confinement) which makes it difficult to find energy-conserving transitions, see insets of Figures. 5 a)-d). In contrast, a large number of peaks are observed in the hot carrier distributions of Na92, where the energy level spacing is significantly reduced. Interestingly, the larger fraction of the plasmon energy is transferred to the hot electrons in these systems.

Table 1 shows the total number of hot carriers generated by the decay of a single excitation in the three Na nanoparticles. The total generation rate for the dominant plasmonic excitation (PL in Na40 and Na92 and PL2 in Na58) does not show a strong dependence on the nanoparticle radius. It is important to recall, however, that these rates are calculated for a single excitation quantum. As the transition dipole moment of these excitations increases with system size, see Table 1, more excitation quanta are excited per incoming photon and this gives rise to a larger number of hot carriers in larger nanoparticles.

As discussed in the Methods section, our approach can also be used to study the hot-carrier distribution resulting from the decay of neutral excitations without plasmonic character. Figure 4 shows the hot-carrier distributions resulting from decay of the prominent single-pair excitations in Na92 (denoted SP in Figure 2 c). Similar to the plasmonic excitation in this system, the larger fraction of the excitation energy is transferred to the hot electrons. Interestingly, the total rate of hot carriers from the decay of the electron-hole pair excitation is larger than from the decay of the plasmonic excitation in this system, see Table 1.

The energy of the localized surface plasmon can be easily modified by placing the nanoparticle in different dielectric environments. Specifically, the plasmon energy is reduced as the dielectric constant of the environment increases. To study the effect of the dielectric environment on the hot-carrier generation rates resulting from the decay of plasmonic excitations, we evaluated Eq. (7) with a reduced plasmon energy $\tilde{\omega}_{P}$ (but keeping the orbital energies and coupling strengths unchanged). Figure 6 shows the resulting total hot-carrier generation rates for three Na nanoparticles as function of the environment-screened plasmon energy. We observe that the reduction of the plasmon energy can lead to significant enhancements in the total hot-carrier rates when the reduced plasmon energy matches the energy of multiple transitions from occupied to empty states in the quasiparticle spectrum. A particularly strong increase is found for plasmon energies smaller than half of the unscreened value $\omega_{P}$. This is caused by a strong increase of the electron-plasmon coupling strength for low-energy transitions, see inset of Figure 6.

3.3 Comparison to the semiclassical approach

To compare the hot-carrier rates from our quantum approach to the semiclassical approximation, we carry out semiclassical calculations using the electric field strength $E_{\rm 0}^{1PL}$ which generates a single plasmon in the nanoparticle. Table 1 shows the resulting total hot-carrier rate for Na92 for the classical plasmon energy. Surprisingly, we find that the semiclassical rate is more than an order of magnitude larger than the quantum result.

To understand this discrepancy, we compare the transition dipole moment of the LSP obtained from the quantum calculation with the results from the semiclassical approach. Table 1 shows that the semiclassical transition dipole moment of the plasmon is six times larger than the quantum mechanical one. This is consequence of two factors: (i) as discussed above, the transition density of the plasmons in small nanoparticles deviates significantly from the perfect dipolar shape predicted by classical electrodynamics (see Figure 3) and (ii) in small nanoparticles, the plasmon is one of many excitations which share the total available oscillator strength (according to the f-sum rule), while in the semiclassical approach the total oscillator strength is concentrated in the plasmon, see Supplementary note 3. Because of their smaller transition dipole moments, the quantum plasmons couple less strongly to the transition dipole moments of the hot electron-hole pairs resulting in significantly reduced generation rates compared to the semiclassical calculations.

4 Discussion

We have presented a new approach for studying hot carrier generation in metallic nanoparticles which takes quantum plasmonic effects, such as the bosonic nature of the plasmon, fully into account. We employ an effective fermion-boson Hamiltonian and determine its parameters for specific nanoparticles using (time-dependent) density-functional theory and many-body perturbation theory within the GW approximation. In particular, an expression for the coupling strength of quasiparticles to neutral excitations is obtained by comparing the self energy of the effective Hamiltonian with the first-principles GW self energy. We have used this approach to study the decay of single plasmons into hot carriers in small sodium nanoparticles with different radii. We find that some systems exhibit multiple plasmonic excitations, while others have electron-hole pair excitations with large oscillator strengths. The hot carrier distributions from the decay of these excitations exhibit a molecular character with discrete peaks that are more closely spaced as the nanoparticle size increases. Interestingly, we find that in all systems a large fraction of the plasmon energy is transferred to the hot electrons. We also compare our results to semiclassical calculations and find that the semiclassical results provide qualitative insights but overestimate hot carrier rates for the small nanoparticles under consideration. Our approach opens the possibility to study hot carrier generation in the quantum plasmonic regime with potential application to quantum-controlled devices, including single-photon sources, transistors and ultra-compact circuitry at the nanoscale.

4.1 Methods

In this section, we describe how the various parameters of the effective electron-plasmon Hamiltonian $H_{\rm eff}$, Eq. (1), are obtained.

4.2 Quasiparticle energies

To model the electronic structure of spherical metallic nanoparticles of radius $R=r_{\rm s}N^{1/3}$ (where $r_{\rm s}$ denotes the Wigner-Seitz radius of the bulk material and $N$ is the number of electrons in the nanoparticle), we employ the jellium approach which has been widely used to describe metallic clusters [50, 51, 52, 53]. In this parameter-free method, the positive charge of the atomic nuclei is smeared out homogeneously throughout the volume of the nanoparticle.The jellium approach often yields accurate electronic properties for metals with s- or p-electron bands, but at a significantly reduced computational cost compared to full atomistic descriptions. Other authors [14, 54] have used phenomenological spherical well models to describe the electronic structure of metallic nanoparticles, but we have found that such models do not reproduce the experimentally observed ordering of energy levels [55]when electron-electron interactions are included. For example, photoelectron spectroscopy of sodium clusters reveals the following ordering of energy levels (in order of increasing energy): 1s, 1p, 1d, 2s, 1f and 2p [55]. This ordering is reproduced by the jellium approach, but not by the interacting spherical well approach.

Ground state properties, such as the total energy or the electron density, are obtained by solving the Kohn-Sham equations

[TABLE]

where $\phi_{i}(\mathbf{r})$, $V_{KS}(\mathbf{r})$ and $\epsilon_{i}$ denote the Kohn-Sham orbitals, Kohn-Sham potential and Kohn-Sham energies, respectively. The Kohn-Sham potential consists of a nuclear, a Hartree and an exchange-correlation contribution, for which we employ the Perdew-Zunger parametrization of the local density approximation (LDA) [56].

We only carry out calculations for nanoparticles with closed electronic shells. For such systems, the Kohn-Sham potential exhibits spherical symmetry and the Kohn-Sham orbitals can be expressed as the product of a spherical harmonic and a radial function which is obtained by direct integration on a real-space grid. To approximate quasiparticle energies which correspond to electron addition or removal energies, we have calculated the ionization potential of the nanoparticles using the $\Delta$-SCF approach (i.e., by calculating the total energy difference of the neutral and ionized nanoparticles) and shifted all Kohn-Sham energies such that the energy of the highest occupied orbital agrees with the calculated ionization potential. Alternatively, the quasiparticle energies can be obtained from GW calculations, but at a significantly larger computational cost. Figure 7 shows the resulting quasiparticle energy levels and Kohn-Sham potential for a sodium nanoparticle ($r_{s}=4.00$ Bohr) consisting of 92 atoms.

4.3 Neutral excitations

Neutral excitations are obtained using the random phase approximation (RPA). In particular, we solve Casida’s equation [57]

[TABLE]

where $\hbar\omega_{I}$ denotes the energy of the neutral excitations and $F^{I}_{vc}$ is the corresponding eigenvector, which determines the transition density, see Eq. (4). The Casida matrix is given by

[TABLE]

where $K_{vc,v^{\prime}c^{\prime}}$ denote Coulomb matrix elements

[TABLE]

The Coulomb integrals were computed using the LIBERI library [58] which we modified to perform integrals using real spherical harmonics. Note that we have chosen to work within a linear-response framework because it allows exploitation of spherical symmetry which is broken by the perturbing electric field in a real-time framework.

Besides their energy, other important properties of neutral excitations are their oscillator strength $f_{I}$ and their collectivity $C_{I}$. The oscillator strength is obtained from the eigenvectors of Casida’s equation via

[TABLE]

where $\mu_{vc}=e\int d\mathbf{r}z\phi_{v}(\mathbf{r})\phi_{c}(\mathbf{r})$ denotes the dipole moment matrix element for the vc-transition in the z-direction, which is parallel to the perturbing electric field, $\mu_{I}$ is the transition dipole moment of the excitation $I$ and $f_{vc}$ is the occupation number difference.

To define the collectivity of a neutral excitation, we first discuss the two extreme cases. When an excitation is perfectly collective, we expect that all components of $|F^{I}_{vc}|^{2}$ are equal to $1/N_{\rm pair}$, where $N_{\rm pair}$ denotes the total number of electron-hole pair states (note that this is finite as we only consider bound electron states). For an ideal electron-hole pair excitation, on the other hand, all components of $|F^{I}_{vc}|^{2}$ are zero except for one whose value is unity. For a general excitation, we first test if any component of $|F^{I}_{vc}|^{2}$ is larger than $\lambda/N_{\rm pair}$, where $\lambda$ is set to $500$ (we have tested that our results do not depend strongly on the choice of $\lambda$). If such a component is found, the state is identified as an electron-hole pair state and $C_{I}$ is set equal to the number of such components. Otherwise, the state is identified as collective and the collectivity is calculated as the total number of non-zero components of $|F^{I}_{vc}|^{2}$. In practice, we find that excitations with a plasmonic character have collectivity values corresponding to $\sim 10$ percent of non-zero components, while for electron-hole pairs typical $C_{I}$ is less than one percent. All DFT and TDDFT calculations were carried out using in-house computer codes.