Next Article in Journal
Mechanical Performance of Multidirectional Buckling-Based Negative Stiffness Metamaterials: An Analytical and Numerical Study
Next Article in Special Issue
Dynamics of Frenkel Excitons in Pentacene
Previous Article in Journal
Spray-Drying of Electrode Materials for Lithium- and Sodium-Ion Batteries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microscopic Electron Dynamics in Metal Nanoparticles for Photovoltaic Systems

1
Department of Quantum Technologies, Faculty of Fundamental Problems of Technology, Wrocław University of Science and Technology, 50-370 Wrocław, Poland
2
Madrid Institute for Advanced Studies in Nanoscience (IMDEA Nanoscience), C/Faraday 9, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Materials 2018, 11(7), 1077; https://doi.org/10.3390/ma11071077
Submission received: 15 May 2018 / Revised: 19 June 2018 / Accepted: 20 June 2018 / Published: 25 June 2018
(This article belongs to the Special Issue Advanced nanostructures for Photonics and Photovoltaics)

Abstract

:
Nanoparticles—regularly patterned or randomly dispersed—are a key ingredient for emerging technologies in photonics. Of particular interest are scattering and field enhancement effects of metal nanoparticles for energy harvesting and converting systems. An often neglected aspect in the modeling of nanoparticles are light interaction effects at the ultimate nanoscale beyond classical electrodynamics. Those arise from microscopic electron dynamics in confined systems, the accelerated motion in the plasmon oscillation and the quantum nature of the free electron gas in metals, such as Coulomb repulsion and electron diffusion. We give a detailed account on free electron phenomena in metal nanoparticles and discuss analytic expressions stemming from microscopic (Random Phase Approximation—RPA) and semi-classical (hydrodynamic) theories. These can be incorporated into standard computational schemes to produce more reliable results on the optical properties of metal nanoparticles. We combine these solutions into a single framework and study systematically their joint impact on isolated Au, Ag, and Al nanoparticles as well as dimer structures. The spectral position of the plasmon resonance and its broadening as well as local field enhancement show an intriguing dependence on the particle size due to the relevance of additional damping channels.

Graphical Abstract

1. Introduction

An accurate description of microscopic properties of metal nanoparticles (metal NPs—MNPs) is important to predict the optical response of e.g., molecules in close proximity to metal surfaces and resulting field enhancement and quenching effects. Nanoparticles as part of functionalized layers in sensing, spectroscopy [1] and light harvesting applications, photovoltaics [2,3,4,5,6,7] and photocatalysis [8,9,10,11,12], can improve the performance of such devices. They are efficient subwavelength scatterers improving the light trapping effect and MNPs provide, in particular, large local fields enhancing charge carrier generation, absorption, and light-induced effects from other nanostructures such as spectral conversion [13] or photoluminescence [14].
For over a hundred years, modeling of the optical properties of MNPs relies on classical electrodynamics. In highly symmetric cases (spherical and cylindrical NPs) analytic solutions are obtained within Mie scattering theory [15] using corresponding basis functions. The electric part E of the electromagnetic field creates a polarization field P = α ( ϵ 0 , ϵ ) E in solid matter, expressed in terms of the permittivities ϵ 0 ( ω ) and ϵ ( ω ) of the environment and the bulk material, respectively. This polarizability α , depending only on the optical response at a frequency ω , neglects microscopic electron interaction effects at the ultimate nanoscale arising not only from the quantum nature of the free electron gas in metals, but also from accelerated motion in the plasmon oscillation.
Light-matter interaction involves processes within the electron subsystem in solids, crystals and molecules. Inhomogeneities on the length scale of the de Broglie wavelength λ e = h 2 m E produce scattering and interference effects of electrons which mutually interact with incoming light, see Figure 1a. Hereby, h is Planck’s constant, m is the (effective) electron mass which depends on the bulk material, and E is the energy of the electron wave. Typically, this wavelength is about 7.5 nm in solids at room temperature T = 300 K, where E = k B T with the Boltzmann constant k B . For MNPs, the main source of electron scattering is the particle surface, see Figure 1b, where the surface-to-volume ratio indicates the relevance of such scattering events.
Microscopic interaction effects of electrons in metals are accurately described using first-principle methods, e.g., Density Functional Theory (DFT) [16,17,18]. These solve Schrödinger’s equation for a large, but finite number of electron wave functions from all atoms in the considered system. Unfortunately, even with strong approximations such as the Time Dependent Local Density Approximation (TDLDA), time-consuming algorithms limit their applicability to particles of a few nanometers in size [19,20,21]. Moreover, advances in fabrication of nanostructures along with experimental access to particle sizes and interparticle spacings below 10 nm led to the possibility of direct or indirect observation of such effects [22,23,24,25,26,27,28,29]. The situation described above resulted in increased interest in semi-classical approaches towards the incorporation of damping and interaction effects stemming from the quantum nature of charge carriers, illustrated in Figure 1. In this article, we present two such semi-classical approaches, the Random Phase Approximation (RPA) and Generalized Nonlocal Optical Response (GNOR), and ultimately combine them into a single framework to study their joint impact on MNPs of different materials, sizes and in different environments.
The original formulation of light scattering by a sphere by Gustav Mie [15] excludes microscopic dynamics of the conduction band electrons in bulk and surface effects. However, efforts to extend have been made since the 1970s [30,31,32,33,34,35,36,37,38,39]. Advanced semi-classical material models can be derived from perturbative theories [40,41], by separating the free electron dynamics from the core electron polarization via the hydrodynamic equation for an electron plasma [41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59], and from microscopic theories [60,61,62,63,64]. It should be noted that a major advantage of ab initio methods lies in their capability to account for the electron spill-out (evanescent tail of the electron wave functions) of the electron density into the surrounding dielectric medium. It was shown within the hydrodynamic framework that the electron spill-out can be adequately incorporated [57,65] and a current-dependent potential can be accounted for [66], which is, however, out of scope of the present study.
In this article, we combine two semi-classical approaches towards microscopic electron dynamics into a single feasible framework to address quantum corrections in MNPs allowing the description of isolated particles, clusters and large-scale (two- or three-dimensional) devices via the integration of analytical expressions into standard procedures. We hereby focus on results on damping in MNPs derived from the microscopic Random Phase Approximation (RPA), stemming from Lorentz friction, and spatial dispersion (nonlocal) effects obtained with the hydrodynamic approach. We discuss briefly the separate ingredients of these approaches in the next sections and give more details in the methods section. Moreover, we compare and combine the different processes of mesoscale electron dynamics stemming from scattering, Figure 1a,b, irradiation (Lorentz friction), Figure 1c, and nonlocal interaction, Figure 1d, and study their impact on the optical response of isolated MNPs and dimers. An emphasis is put on the size regimes where these effects are dominant for the materials silver, Figure 2, as well as for aluminum and gold, Figure 3.

2. Results

We briefly discuss classical electrodynamics and mesoscopic electron dynamics obtained from the RPA and GNOR theories. In summary, we compare quantum correction models stemming from microscopic RPA derivations with the following, semi-classical damping expressions
( 1 a ) γ = γ p , ( Mie ) ( 1 b ) γ = γ p + C v F a , ( Kreibig ) ( 1 c ) γ = γ p + C v F a + ω 1 3 ω 1 a c 3 , ( perturbative ) ( 1 d ) γ = Im ( Ω 2 ) = 1 3 l + 1 + 6 l q 2 2 3 3 l A + A 2 1 3 6 l , ( Lorentz ) ,
and nonlocal interaction effects. Both approaches are described in more detail in the next sections and the methods section. The advantage in the analytic formulation is the straightforward integration with existing computational tools for nanospheres using modified Mie simulations and multiple scattering techniques [67] for clusters thereof or commercial software such as COMSOL (http://www.comsol.com).

2.1. Classical and Phenomenological Approaches

Typically, the optical response of a metal is described with the Drude model via the frequency-dependent permittivity
ϵ ( ω ) = ϵ b ω p 2 ω ( ω + i γ p ) ,
where ϵ b is the background permittivity given by bound (valence band) electrons, ω p 2 = 4 π n 0 e 2 / m is the plasmon frequency, determined by the material dependent electron density n 0 and mass m, and γ p is the inherent (bulk) damping rate. This widely used Drude model applies only to bulk material and should be modified for nanostructures to include effects due to the finite size of the system. One of the corrections considered by Kreibig and von Fragstein [68] is the inclusion of an additional damping due to the scattering on the physical particle boundaries, depicted in Figure 1b. This is in particular important in particles of sizes equal or smaller than the mean free path λ b of electrons in bulk metal. In such a case, the electrons will experience (in the classical picture) additional scattering from the boundary of the system. Mathematically, it is described as γ K = v F / L e f f , where v F is the Fermi velocity of the electron gas and L e f f is the effective mean free path of electrons resulting from collisions with the particle surface [38,68,69]. The common feature is that L e f f reflects the volume (proportional to the number of electrons inside the nanoparticle) to surface ratio of the particle. According to this, we get the γ K ( a ) = C v F / a , where a is the radius of the nanoparticle and C is a constant of the order of unity which depends on the scattering type and particle radius. Similarly, collision effects in the bulk, depicted in Figure 1a, can be described via the damping term γ p v F / 2 λ b .

2.2. Random Phase Approximation

Nevertheless, this phenomenological approach neglects the microscopic dynamics of electrons inside the MNP. Their accelerated movement (plasmon oscillation) leads to energy loss via irradiation of the electromagnetic field, see Figure 1c. In case of nanoparticles much smaller than the incident wavelength, this effect can be expressed by the Lorentz friction, an effective field stemming from the plasmon induced dipole field D ( t ) as E L = 2 / 3 c 3 3 D ( t ) / t 3 , with c being the speed of light [70]. The dynamics of the electron density can be described using a driven, damped oscillator, with the incident electromagnetic wave being the driving force and the damping arising form electron scattering (bulk γ p and Kreibig damping γ K ) and electromagnetic field irradiation (Lorentz friction).
An analytical form of the exact solution for the damping γ and self-frequency ω L (the exponents Ω i of solution ∼ e i Ω i t for self-modes i) including Lorentz friction exists [61], which is discussed in more detail in the methods section. They can be summarized as follows
Ω 1 = i 3 l i 2 1 / 3 ( 1 + 6 l q ) 3 l A i A 2 1 / 3 3 l Im , Ω 2 = i 3 l + i ( 1 + i 3 ) ( 1 + 6 l q ) 2 2 / 3 3 l A + i ( 1 i 3 ) A 2 1 / 3 6 l = ω L + i γ , Ω 3 = ω L + i γ = Ω 2 * ,
where A = B + 4 ( 1 6 l q ) 3 + B 2 1 / 3 , B = 2 + 27 l 2 + 18 l q , q = 1 τ 0 ω 1 , l = 2 3 ϵ 0 a ω p c 3 3 and 1 / τ 0 = γ p . Exact inclusion of the Lorentz friction indicates that the radiative losses and the self-frequencies are a complicated function of particle radius as given by Equation (3), see the methods section for a detailed discussion.
Direct comparison to experimental work for this framework is available within Refs. [61,62,63,64] and good agreement has been found.

2.3. Nonlocal Optical Response

Aside from electron irradiation due to Lorentz friction, we discuss spatial dispersion (nonlocality) which denominates the effects of electron coupling over a short distance, see Figure 1d [40]. Such interactions are inherent to the solution for the displacement field D of the Coulomb equation
D ( ω , r ) = 0 D ( ω , r ) = d r ϵ ( ω , r , r ) E ( ω , r ) .
In homogeneous media, we can assume a dependence on the distance r r rather than on the specific position of electrons, which allows solving Maxwell’s equations in Fourier space D ( ω , k ) = ϵ ( ω , k ) E ( ω , k ) .
The dependence on the wave vector k enables us to describe nonlocal electron-electron interaction (Coulombic force) and electron diffusion effects. It is important to note that the large-k response that originates in the subwavelength oscillations of plasmonic excitations is not only an inherent prerequisite for many intriguing wave phenomena, but also particularly sensitive to nonlocality. However, the common Mie result has no upper wavelength cut-off and does suppress short-range electron interactions which can strongly dampen the response beyond ω / v F . We show in the corresponding section below that accounting for nonlocal response leads to longitudinal pressure waves as additional solutions to the combined system of differential equations of the electromagnetic wave equation and (linearized) Navier-Stokes equation. This is in contrast to the damping expressions derived by Kreibig and for Lorentz friction. Such additional waves offer further damping channels, however, they can also support resonant enhancement effects [12,51,59,71].
Experimental work focusing on the blueshift found for nanoparticles decreasing in size, as well as the influence of the electron-spill out has been studied in Refs. [22,23,24,25,26,27,28,29], including comparisons with the hydrodynamic model.

2.4. Remarks on Retardation, Multipolar Response and Computational Feasibility

Both of the presented semi-classical approaches towards microscopic corrections in the mesoscale electron dynamics in metal nanoparticles have the advantage of analytic expressions fully compatible with existing computational procedures. For the quantum confinement picture of Kreibig and the mesoscopic RPA result for the Lorentz friction, modified damping terms were derived, see Equations (1a)–(1d), which can be used to directly replace the damping in the Drude expression for the permittivity given in Equation (2) and subsequently be used in standard Mie calculations and procedures to calculate optical properties of complex structures, e.g., with a multiple scattering approach [67] or within commercial software such as COMSOL.
It is important to note that although all electrons participate in plasmon oscillations, part of their irradiation is absorbed by other electrons in the system. This is in analogy with the skin-effect [72] in metals and introduces an effective radiation active electron layer of the depth h 1 / σ ω ( σ is the conductivity) underneath the particle surface. Therefore, the effective energy transfer outside of the nanoparticle will be reduced by the factor 4 π 3 a 3 a h 3 / 4 π a 3 3 . According to this, we expect a decrease of radiative damping, especially for larger particles.
The nonlocal theory introduces a novel type of electron motion, longitudinal pressure waves, in addition to the transversal modes stemming from the classical electromagnetic wave equation. This additional electronic excitation offers further damping channels due to the energy lost in dampened motion. Here, the Mie coefficients are derived from the coupled system of optical and electronic excitation yielding modified scattering matrices that can again be implemented in existing methods. The properties of the longitudinal wave are given by analytic expressions such as their wave vector and their importance with respect to the common Mie solution is entirely captured in a single additional term, see the methods section for details.
Retardation is important when either the particle radius or the overall system size becomes large, i.e., for particle dimers, clusters and arrays. Although the presented microscopic effects are highly localized, they can have a strong impact on a larger particle or system in the interplay with long-range retardation effects. In addition, particle layer modes can couple to nonlocal modes within particle arrays and thus increase their impact on a larger scale [59,71]. It is thus noteworthy that the hydrodynamic theory and the damping terms stemming from microscopic analysis within the RPA allow fully retarded calculations; equally for planar geometries (nonlocal Fresnel coefficients) [51] and regular, two-dimensional particle arrays [41,59,71] and even charge carriers in electrolytes (Nonlocal Soft Plasmonics) [12].

3. Discussion

3.1. Single Metal Nanoparticles

We compare the quantum correction models introduced in the previous section, see Figure 1, as well as the combined effect of Kreibig damping Equation (1b), Lorentz friction Equation (1d) and spatial dispersion to classical Mie calculations for the materials gold, aluminum and silver in Figure 2 and Figure 3. Hereby, we show the effect on the Localized Surface Plasmon Resonance (LSPR) for all materials in Figure 2a, confirming that the modified damping rates do not alter the resonance position predicted by the classical calculations, whereas nonlocal response—and in combination with any damping model—does predict an increasing blueshift of the nanoparticle resonance with decreasing particle size. Looking at the extinction cross section as a function of particle radius in Figure 2b for silver and Figure 3 for gold and aluminum, we find that all correction models result in a reduction of the optical response in dependence of both the material and particle size, typically yielding a different optimized particle size. Hereby, Kreibig damping with a ∼ 1 / a dependence drastically attenuates the optical response for the smaller size regime below the maxima (15 nm for Ag, 20 nm for Au, and 10 nm for Al), while the complex size dependence of Lorentz friction results in a greater effect above this particle size. The diffusion coefficient in the hydrodynamic (GNOR) model (imaginary part of the nonlocal parameter β GNOR ) is chosen thus that its dampening effect captures the Kreibig result [56]. This is best seen in Figure 3a for Au. The hydrodynamic pressure (real part of the nonlocal parameter β GNOR ) describes Coulomb interaction between electrons and results in the blueshift observed in Figure 2a at very small particle sizes below 5 nm. We can further incorporate the analytical expressions for Lorentz friction. This combined result shows the strongest attenuation since all different damping channels are included. At a larger particle size (60 nm for Ag, 80 nm for Au, and 40 nm for Al) all material models converge with classical Mie theory where the mesoscale electron dynamics cease to have an impact.
The damping associated with the Lorentz friction can be approximated to the simpler perturbative expression Equation (1c) in a narrow size window, see the methods section for a detailed discussion. Since the exact solution can be obtained with analytical expressions which can be incorporated into standard calculation schemes, we discuss exclusively exact Lorentz friction results.
We study the (maximum) field enhancement factor EF = | E | 2 / | E 0 | 2 just outside of the NP ( r a + ) for the different damping models in Figure 4 for gold nanospheres. Hereby, Figure 4a shows the spectral position of the field maximum. The local field enhancement reveals the size dependence of the field resonance with the damping rates. It should be emphasized that Kreibig damping shows a strong redshift for small particle sizes of the spectral position of local field enhancement maxima in contrast to experimental findings [25,26,27] and approaches the Mie result for larger sizes. Nonlocal optical response agrees with the blueshift of the plasmon resonance found experimentally for noble metals, as already seen in the extinction cross section, Figure 2a. However, in order to correctly describe simple metals, the inclusion the electron spill-out region [52,57,65] is crucial. Furthermore, advances towards the spatial dispersion found in (doped) semiconductors were made recently [73,74], which is of further interest when using dielectric nanoparticles to enhance the performance of photovoltaic devices.
Lorentz friction is closest to the classical calculation for smaller sizes and deviates stronger at larger sizes. This is in agreement with the findings of Figure 2 and Figure 3. The corresponding field enhancement, shown in Figure 4b for gold MNPs in water, is strongly suppressed for the considered particle size range when including the damping models while spatial dispersion by itself reduces the predicted field enhancement mostly for smaller particle sizes and converges with the classical Mie result rapidly with increasing particle size. This behavior is corrected by incorporating Lorentz friction into the GNOR result.
Figure 4c,d shows the (maximum) field enhancement of gold nanoparticles in dependence of the refractive index (RI) of the surrounding medium (from air n = 1 to Si n = 3.4 ) for two particle sizes. This is accompanied with a linear (in case of the nonlocal theory approximately linear) shift in the resonance wavelength towards longer wavelengths (not shown). With increasing RI of the host medium, the enhancement factor reaches a saturation value which for increasing particle size converges for all material models discussed. The discrepancy between the local field enhancement values predicted remains similar for small MNPs in different host media spanning several orders of magnitude.
The complexity of the Lorentz friction makes it necessary to restrict ourselves to the dipolar response of the plasmon oscillation. It is therefore important to consider the material, particle size and wavelength regime in order to assess whether the dipolar response model is adequate for the system under study. We show in Figure 5 for Au NPs the dipolar and the converged result of local field enhancement obtained from classical Mie calculations at a fixed frequency close to the respective plasmon resonance. Here, the dipolar approximation is valid up to ca. 100 nm in particle radius which in general covers the discussed microscopic effects well. The inset in Figure 5 compares this for the combined theories showing small differences already for particles above 25 nm radius.

3.2. Dimers

For particle dimers, in addition to their size, the particle distance becomes important and retardation effects cannot be neglected for larger particles in close proximity. This can transfer the impact of localized microscopic electron dynamics onto a larger structure. Figure 6 shows the (maximum) field enhancement at the center of a gold dimer in water as a function of both particle size and distance for the different theories considered. The impact of nonlocal response, Figure 6b, on the classical Mie theory, Figure 6a, is visible as strong quenching of the local fields. It is worth remembering that one main effect is a blueshift in the position of the maximum enhancement factor, see again Figure 4a and Ref. [41]. In addition, the maximum field enhancement within the parametric area of particle and gap size is EF 9000 for the Mie calculations and EF 3000 for the nonlocal theory, showing that indeed there is an impact of the longitudinal waves found. The damping observed within Kreibig theory, Figure 6c, is dramatic for the dimer setup and the dominant contribution in the combined theory as seen in Figure 6d. This is also evidenced by comparing the Lorentz friction with and without nonlocal damping, see Figure 6e,f, respectively. The Lorentz friction has a strong impact on the optical response for larger particle sizes, but also dampens the dimer setup for increasing gap size, which points towards retardation and the increasing structural size as the main source for this damping effect. This leads to slightly stronger damping when combined with the additional plasmon quenching within GNOR in Figure 6f.
The strong field quenching poses limitations to the photovoltaic effect in solar cells. However, considering different materials for MNPs and their environment, the size regimes where local field quenching is dominant can be avoided with the presented theory of combined damping.

3.3. Summary

In conclusion, we have presented a number of semi-classical corrections to incorporate electron dynamics and non-classical interaction effects into optical response calculations for nanoparticles. Hereby, pure damping models, such as the Kreibig damping and Lorentz friction, derived from microscopic RPA theory, show an intriguing dependence on the particle size, where the material influences relevant size regimes. On the other hand, semi-classical nonlocal theories allow evoking additional modes in the system by explicitly considering mesoscopic dynamics of free electrons. This results in a correction of the spectral position of resonant phenomena and introduces additional, implicit damping channels. The phenomenological Kreibig damping does yield a plasmon broadening that agrees with experiments [38], however, it also introduces a redshift of the resonance with respect to the classical Mie result contrary to measurements on nanoparticles [25,26,27,29]. This is addressed by using the hydrodynamic GNOR (generalized nonlocal optical response) approach, i.e., by introducing a diffusion parameter, able to reproduce the Kreibig damping while fully capturing the observed plasmon broadening.
An important aspect is that the resulting analytical expressions can be implemented into existing computational procedures in a straightforward manner, as isolated theories or combined, allowing the comparison to experiments with little added numerical effort. We have studied the combined effect of these mesoscopic electron interaction effects for single nanospheres and gold dimers and have evidenced the importance of retardation as a way to communicate localized quantum effects and impact a larger structure.
The straightforward inclusion of electro-optical effects at the nanoscale into (metal) nanoparticle systems is of importance in nanostructures employed for photovoltaics and catalysis as well as in spectroscopy and sensing applications.

4. Methods

4.1. Electron Dynamics within the RPA

The model of electron dynamics inside MNPs [60,61,62] presented here is an extension to the RPA theory developed by Pines and Bohm [75] for bulk metals. In our model, a finite, rigid jellium defines the shape of a nanoparticle. The plasmon oscillations are described as local electron density fluctuations ρ ^ ( r , t ) obtained from the Heisenberg equation
d 2 ρ ^ ( r , t ) d t 2 = 1 ( i ) 2 [ [ ρ ^ ( r , t ) , H ^ e ] H ^ e ]
with a corresponding Hamiltonian H ^ e for electrons inside the MNP in the jellium model taking the following form
H ^ e = j = 1 N e 2 j 2 2 m e 2 n e ( r ) d 3 r | r j r | + 1 2 j j e 2 | r j r j | .
The operator of the local electron density is defined as
ρ ( r , t ) = Ψ e ( t ) | j δ ( r r j ) | Ψ e
where Ψ e is the electron wave function, N e is the number of collective electrons, r j and m are their positions and mass. The ion field is approximated as averaged background charge density and described as n e ( r ) | e | = n e Θ ( a r ) | e | , where Θ is the Heaviside step function, a is the radius of the MNP and n e = N e / V .
The first term in the Hamiltonian stands for the kinetic energy of electrons, the second for interaction between electrons and positive background charges (approximating the ion lattice potential) and the last for electron-electron Coulomb interaction.
Taking into account the sharp form of the positive charge density n e ( r ) , one can decompose Equation (5) into two parts corresponding to the inside and outside of the NP, which leads to two separate solutions describing the surface and bulk plasmons. This description is valid for NPs larger than ca. 5 nm for which the surface is well defined and the spill-out effect is negligible.
δ ρ ˜ ( r , t ) = δ ρ 1 ˜ ( r , t ) ,   for   r < a , δ ρ 2 ˜ ( r , t ) ,   for   r a , ( r a + ) .
The electron density fluctuations are then described with the formulas
2 δ ρ 1 ˜ ( r , t ) t 2 = 2 3 ϵ F m 2 δ ρ 1 ˜ ( r , t ) ω p 2 δ ρ 1 ˜ ( r , t )
and
2 δ ρ 2 ˜ ( r , t ) t 2 = [ 2 3 ϵ F m r r δ ρ 2 ˜ ( r , t ) + ω p 2 4 π r r d 3 r 1 1 r r 1 ( δ ρ 1 ˜ ( r 1 , t ) Θ ( a r 1 ) + δ ρ 2 ˜ ( r 1 , t ) Θ ( r 1 a ) ) ] δ ( a + ε r ) 2 3 m 3 5 ϵ F n e + ϵ F δ ρ 2 ˜ ( r , t ) r r δ ( a + ε r ) .
where ϵ F is the Fermi energy.
The structure of the above equations is of an harmonic oscillator, which allows including a damping term in phenomenological manner by adding to the right hand side 2 / τ 0 ρ ˜ 1 ( 2 ) ( r , t ) / t . The damping 2 / τ 0 = γ p + γ K includes collision effects and Kreibig damping due to the particle boundary.
Assuming homogeneity of the external electric field E ( t ) inside the NP (dipole approximation), the solution for surface modes reduces to a single dipole mode
δ ρ ˜ ( r , t ) = m = 1 1 Q 1 m Y 1 m ( Ω ) , for r a , ( r a + )
and for bulk modes δ ρ ˜ ( r , t ) = 0 where r < a .
The function Q 1 m ( t ) ( m = 1 , 0 , 1 ) represents dipole modes, Y l m ( Ω ) is the spherical function. The former can be related to the vector q ( t ) via Q 11 = 8 π / 3 q x ( t ) , Q 10 = 4 π / 3 q x ( t ) , Q 1 1 = 8 π / 3 q y ( t ) satisfying the equation
2 t 2 + 2 τ 0 t + ω 1 2 q ( t ) = e n e m E ( t ) .
Then the plasmon dipole can be defined as
D ( t ) = e d 3 r r δ ρ ( r , t ) = 4 π 3 e q ( t ) a 3 .
Knowing this, the damping caused by electric field irradiation can be simply added to the right hand side of Equation (12) as additional field E L = 2 / 3 c 3 3 D ( t ) / t 3 hampering charge oscillations and can be rewritten in the form
2 t 2 + ω 1 2 D ( t ) = t 2 τ 0 D ( t ) + 2 3 ω 1 ϵ 0 ω p a c 3 3 2 t 2 D ( t ) .
The above equation is a third order linear differential equation and the exponents ∼ e i Ω i t of its solutions are given in Equation (3). A perturbation approach can be applied to Equation (14) for small particles using 2 D ( t ) / t 2 = ω 1 2 D ( t ) . Then the resulting damping term takes the form γ = 2 / τ 0 + ω 1 / 3 ϵ 0 ω a / c 3 3 . The comparison of both damping terms is shown in Figure 7 justifying the usage of the perturbation formulation for (gold) particles with radii up to ca. 30 nm, where the second term proportional to ∼ a 3 still fulfills the perturbation constrain.
For larger radii, the discrepancy between both solutions grows rapidly since the irradiation losses within the perturbation approach scale as a 3 . Therefore, the radiative losses dominate plasmon damping for large nanospheres. On the other hand, scattering is more important for smaller nanospheres scaling as 1 a . One can observe thus the size-dependent crossover in Figure 3a of the damping at ca. 12 nm for gold.

4.2. Electron Dynamics with the Hydrodynamic Model

In recent years, a great effort to theoretically [41,43,44,46,47,48,49,50,51,53,54,55,56,59] describe and subsequently to experimentally [25,26,27,29] verify the effect of spatial dispersion in metals was made. In the hydrodynamic approach, coupling the electromagnetic wave equation
× × E k 2 ϵ b E = 4 π i k 2 ω j ind
to the (linearized) Navier-Stokes equation
j ind = i ω + i γ p ω p 2 4 π E β 2 + D ( γ p i ω ) ρ ind
allows treating the conduction band electrons as a plasma subject to short-ranged interaction such as the Coulomb force included in the pressure term p = β 2 ρ ind and electron diffusion via the diffusion coefficient D. It is convenient to abbreviate β GNOR 2 = β 2 + D ( γ p i ω ) (where GNOR refers to the Generalized Nonlocal Optical Response model [55,56]). With this, we can write the wave equation in a compact form
2 E + k 2 ϵ E = η ρ ind ,
where η = 4 π 1 ϵ b k 2 β GNOR 2 ω ( ω + i γ ) and ϵ = ϵ b ω p 2 / ω ( ω + i γ p ) . Together with the continuity equation j ind = i ω ρ ind , we readily obtain a separate wave equation for the induced charges
β GNOR 2 2 ρ ind = ϵ ϵ b ω ( ω + i γ p ) ρ ind ,
where E = 4 π / ϵ b ρ ind was used. This yields the wave vector of the longitudinal field and motion of electrons
q = 1 β GNOR ϵ ϵ b ω ( ω + i γ p ) .
Nonlocal theories predict finite distributions of induced charges at an illuminated metal surface—in contrast to classical electrodynamics—with a characteristic penetration depth Im ( 1 / q ) comparable to the electron spill-out [41,76].
Thus, this system of coupled equations yields an additional wave solution, longitudinal in character, and can be solved for different geometries leading to nonlocal extensions of Mie [41,48] and Fresnel coefficients [51], including for charge carriers in electrolytes [12]. Typically, hard-wall boundary conditions are assumed for the additional boundary condition n ^ j ind 0 prohibiting electrons to trespass through the particle surface into the dielectric surrounding, using a uniform electron density n 0 = ω p 2 m / 4 π e 2 inside the material and neglecting the electron spill-out. However, it was shown that a smooth surface distribution of electrons can be taken into account accurately [57,65] and that the hydrodynamic model is capable of dealing with the spill-out by solving the above equations with position-dependent material parameters ω p ( z ) 2 = 4 π n 0 ( z ) e 2 / m .
The main observations of nonlocal theories are a blueshift of the plasmon resonance with respect to the common local approximation and plasmon broadening, in particular tied to the diffusion coefficient which can be set to fully capture the broadening found with Kreibig damping [55,56]. In the present work, we have adopted the diffusion coefficients as deduced in Ref. [56] for the different materials, reflected, for instance, in the correspondence between the Kreibig and GNOR result for gold in Figure 3a. Moreover, we add the Lorentz friction result from the RPA technique summarized in Equation (1d) to our GNOR calculations.
Next, we present the derivation of nonlocal Mie scattering coefficients of individual spheres and nanoshells described with the hydrodynamic model [41] starting from Equation (17) which describes the evolution of the electric field, together with Equation (18) which is the wave equation for the induced charge. The resulting scattering matrices can be used to investigate interacting spheres with a multiple scattering method [67]. The hydrodynamic model has no free parameters which makes the resultant nonlocal response for the short distances involved in the interaction (Coulomb force, diffusion) between the charges of MNPs the sole source of these effects, in contrast to the quantum-confinement picture for plasmon broadening presented by Kreibig.
It is convenient to use an expansion of the electric field into scalar functions [77] as
E = ( 1 / k ) ψ L + L ψ M + × L k i ψ E ,
where L = i r × is the angular momentum operator, and the superscripts E, M, and L indicate electric, magnetic, and longitudinal components, respectively. The additional boundary condition, Equation (16), becomes with r ^ j = 0
β GNOR 2 r ρ ind = e 2 n 0 m k r ψ L + 1 r l ( l + 1 ) ψ E
in terms of the scalar functions and the angular momentum number l using the identity r · ( × L ) = ( i r × ) · L = L 2 = l ( l + 1 ) . The boundary conditions for the electric and magnetic field components result in the continuity of ψ M , ( 1 + r r ) ψ M , ψ L + ( 1 + r r ) ψ E , and ϵ ψ E for the scalar functions.
The magnetic and electric scalar functions ψ ν ( ν = { E , M } ) obey a Helmholtz equation of the form ( 2 + k 2 ϵ ) ψ ν = 0 and can therefore be expanded in terms of spherical Bessel functions ψ ν = L ψ L ν j L ( k r ) . Similarly, the electron density is expanded into ρ ind ( r , ω ) = L ρ L j L ( q r ) , with the longitudinal wave vector q given by Equation (19). The longitudinal scalar function satisfies a different wave equation, namely 2 ψ L = 4 π k / ϵ b , which we find from the Coulomb law ϵ b E = 4 π ρ ind .
Note that the above analysis is needed for the metal region, where the electric ( ν = E ) and magnetic ( ν = M ) field are given by A l ν j L , with j L = j l m ( k r ) . Outside the particle, the longitudinal scalar function vanishes since there are no induced charges in the dielectric surrounding. Therefore, the electric scalar field is given by j l m ( k 0 r ) + t l ν h l m + ( k 0 r ) with unknown parameters A l ν and scattering matrix t l ν . Exploiting the boundary conditions stated above, we find a set of linear equations for the magnetic and electric scattering matrices. Interestingly, the magnetic scattering matrix is unchanged with respect to the local theory, indicating that magnetic modes are not sensitive to the induced longitudinal modes. The scattering matrix for the electric scalar function is more complicated than in the local approximation due to the appearance of ψ L in the metal region that contains information on the nonlocal response. The additional boundary condition yields a prescription to calculate ρ L .
The local scattering matrix can then be extended by a single parameter describing nonlocal behavior of the electron motion in the conduction band
g l = l ( l + 1 ) j l ( θ ) j l ( q a ) q a j l ( q a ) ϵ ϵ b 1
and becomes with θ 0 = k a ϵ 0 and θ = k a ϵ .
t l E = ϵ j l ( θ ) [ θ 0 j l ( θ 0 ) ] + ϵ 0 j l ( θ 0 ) ( [ θ j l ( θ ) ] + g l ) ϵ j l ( θ ) [ θ 0 h l + ( θ 0 ) ] ϵ 0 h l + ( θ 0 ) ( [ θ j l ( θ ) ] + g l ) ,
where the primes indicate differentiation with respect to the θ variables. The scattering coefficients t l ν fully contain the optical response of the particle for an external observer.
Note that the nonlocal parameter g vanishes under the assumption of local response ( β GNOR 0 g l 0 ) fully recovering the original Mie coefficients [15,78]. This allows us to study the electro-optical properties of NPs with only a small correction in available numerical procedures, see for instance Figure 5.
Likewise, for a nonlocal metal nanoshell, the magnetic response is insensitive to the nonlocal properties of the material. The electric part, however, mixes with the longitudinal components from the two interfaces of the metal intermediate layer. For the electric scalar functions, we obtain a linear system of six equations and analytical solutions exist for the metal nanoshell [41,79].

4.3. Simulations

The modeling presented in this article was obtained by both using the commercial software COMSOL Multiphysics (http://www.comsol.com) and in-house numerical code to evaluate Mie coefficients, from Equation (23).
To make predictions that can be compared to experiments, the expressions obtained are used to calculate e.g., the extinction cross section of an individual sphere via
σ ext = 2 π k 2 ϵ 0 l ( 2 l + 1 ) Im ( t l E + t l M ) .
Note that only the electric scattering matrix is sensitive to nonlocal contributions.
The scalar electric field is obtained from j l m ( k 0 r ) + t l E h l m + ( k 0 r ) outside the particle, with the corresponding spherical Bessel and Hankel functions and the related vector field from Equation (20).
The analytic damping expressions Equations (1a)–(1d) are directly introduced as damping terms in the permittivity of the different material permittivities, Equation (2).
For dimers, we use a multiple elastic scattering approach [67].

Author Contributions

Investigation, K.K. and C.D.; Methodology, W.J.; Validation, K.K. and C.D.; Visualization, K.K. and C.D.; Writing—original draft, C.D.; Writing—review & editing, L.J. and W.J.

Funding

This research was funded by the European Cooperation in Science and Technology (COST) under the COST Action MP1406 MultiscaleSolar.

Acknowledgments

C.D. thanks the Comunidad de Madrid (Ref. 2017-T2/IND-6092) and acknowledges financial support by the Spanish Ministry of Economy, Industry and Competition (MINECO) via funding of the Centers of Excellence Severo Ochoa (Ref. SEV-2016-0686).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Weidinger, I.; Ly, K.H.; Öner, H.I.; Querebillo, C.; David, C.; Gernert, U.; Walter, C.; Driess, M.; Leimkühler, S. High electromagnetic field enhancement of TiO2 nanotubes electrodes. Angew. Chem. 2018, 57, 7225–7229. [Google Scholar] [CrossRef]
  2. Atwater, H.A.; Polman, A. Plasmonics for improved photovoltaic devices. Nat. Mater. 2010, 9, 205–213. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, C.; Zhang, G.; Li, H.M.; Yoo, W.J. Localized Surface Plasmon Resonances Caused by Ag Nanoparticles on SiN for Solar Cell Applications. J. Korean Phys. Soc. 2010, 56, 1488–1491. [Google Scholar]
  4. David, C.; Connolly, J.P.; Chaverri Ramos, C.; García de Abajo, F.J.; Sánchez Plaza, G. Theory of random nanoparticle layers in photovoltaic devices applied to self-aggregated metal samples. Sol. Energy Mater. Sol. Cells 2013, 109, 294–299. [Google Scholar] [CrossRef]
  5. Cortés-Juan, F.; Chaverri Ramos, C.; Connolly, J.P.; David, C.; García de Abajo, F.J.; Hurtado, J.; Mihailetchi, V.D.; Ponce-Alcántara, S.; Sánchez, G. Effect of Ag nanoparticles integrated within antireflection coatings for solar cells. J. Renew. Sustain. Energy 2013, 5, 033116. [Google Scholar] [CrossRef] [Green Version]
  6. David, C. Multi-type particle layer improved light trapping for photovoltaic applications. Appl. Opt. 2016, 55, 7980–7986. [Google Scholar] [CrossRef] [PubMed]
  7. Jang, Y.H.; Jang, Y.J.; Kim, S.; Quan, L.N.; Chung, K.; Kim, D.H. Plasmonic Solar Cells: From Rational Design to Mechanism Overview. Chem. Rev. 2016, 116, 14982–15034. [Google Scholar] [CrossRef] [PubMed]
  8. Cushing, S.K.; Li, J.; Meng, F.; Senty, T.R.; Suri, S.; Zhi, M.; Li, M.; Bristow, A.D.; Wu, N. Photocatalytic Activity Enhanced by Plasmonic Resonant Energy Transfer from Metal to Semiconductor. J. Am. Chem. Soc. 2012, 134, 15033–15041. [Google Scholar] [CrossRef] [PubMed]
  9. Akimov, A.V.; Neukirch, A.J.; Prezhdo, O.V. Theoretical Insights into Photoinduced Charge Transfer and Catalysis at Oxide Interfaces. Chem. Rev. 2013, 113, 4496–4565. [Google Scholar] [CrossRef] [PubMed]
  10. Zhang, X.; Chen, Y.L.; Liu, R.S.; Tsai, D.P. Plasmonic photocatalysis. Rep. Prog. Phys. 2013, 76, 046401. [Google Scholar] [CrossRef] [PubMed]
  11. Ma, X.C.; Dai, Y.; Yu, L.; Huang, B.B. Energy transfer in plasmonic photocatalytic composites. Light Sci. Appl. 2016, 5, e16017. [Google Scholar] [CrossRef]
  12. David, C. Two-fluid, hydrodynamic model for spherical electrolyte systems. Sci. Rep. 2018, 8, 7544. [Google Scholar] [CrossRef] [PubMed]
  13. Lakhotiya, H.; Nazir, A.; Madsen, S.P.; Christiansen, J.; Eriksen, E.; Vester-Petersen, J.; Johannsen, S.R.; Jeppesen, B.R.; Balling, P.; Larsen, A.N.; et al. Plasmonically enhanced upconversion of 1500nm light via trivalent Er in a TiO2 matrix. Appl. Phys. Lett. 2016, 109, 263102. [Google Scholar] [CrossRef]
  14. Yuan, Z.; Pucker, G.; Marconi, A.; Sgrignuoli, F.; Anopchenko, A.; Jestin, Y.; Ferrario, L.; Bellutti, P.; Pavesi, L. Silicon nanocrystals as a photoluminescence down shifter for solar cells. Sol. Energy Mater. Sol. Cells 2011, 95, 1224–1227. [Google Scholar] [CrossRef]
  15. Mie, G. Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen. Ann. Phys. 1908, 330, 377–445. [Google Scholar] [CrossRef] [Green Version]
  16. Liebsch, A. Surface-plasmon dispersion and size dependence of Mie resonance: Silver versus simple metals. Phys. Rev. B 1993, 48, 11317–11328. [Google Scholar] [CrossRef]
  17. Larkin, I.A.; Stockman, M.I. Imperfect Perfect Lens. Nano Lett. 2005, 5, 339–343. [Google Scholar] [CrossRef] [PubMed]
  18. Zuloaga, J.; Prodan, E.; Nordlander, P. Quantum Description of the Plasmon Resonances of a Nanoparticle Dimer. Nano Lett. 2009, 9, 887–891. [Google Scholar] [CrossRef] [PubMed]
  19. Lermé, J.; Palpant, B.; Prével, B.; Cottancin, E.; Pellarin, M.; Treilleux, M.; Vialle, J.L.; Perez, A.; Broyer, M. Optical Properties of gold metal clusters: A time-dependent local-density-approximation investigation. Eur. Phys. J. D 1998, 4, 95–108. [Google Scholar] [CrossRef]
  20. Savage, K.J.; Hawkeye, M.M.; Esteban, R.; Borisov, A.G.; Aizpurua, J.; Baumberg, J.J. Revealing the quantum regime in tunnelling plasmonics. Nature 2012, 491, 574–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Esteban, R.; Borisov, A.G.; Nordlander, P.; Aizpurua, J. Bridging quantum and classical plasmonics with a quantum-corrected model. Nat. Commun. 2012, 3, 825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Prodan, E.; Nordlander, P.; Halas, N.J. Electronic Structure and Optical Properties of Gold Nanoshells. Nano Lett. 2003, 3, 1411–1415. [Google Scholar] [CrossRef]
  23. Danckwerts, M.; Novotny, L. Optical Frequency Mixing at Coupled Gold Nanoparticles. Phys. Rev. Lett. 2007, 98, 026104. [Google Scholar] [CrossRef] [PubMed]
  24. Kumar, P.S.; Pastoriza-Santos, I.; Rodríguez-González, B.; de Abajo, F.J.G.; Liz-Marzán, L.M. High-yield synthesis and optical response of gold nanostars. Nanotechnology 2008, 19, 015606. [Google Scholar] [CrossRef] [PubMed]
  25. Scholl, J.A.; Koh, A.L.; Dionne, J.A. Quantum plasmon resonances of individual metallic nanoparticles. Nature 2012, 483, 421–427. [Google Scholar] [CrossRef] [PubMed]
  26. Ciracì, C.; Hill, R.T.; Mock, J.J.; Urzhumov, Y.; Fernández-Domínguez, A.I.; Maier, S.A.; Pendry, J.B.; Chilkoti, A.; Smith, D.R. Probing the Ultimate Limits of Plasmonic Enhancement. Science 2012, 337, 1072–1074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Raza, S.; Stenger, N.; Kadkhodazadeh, S.; Fischer, S.V.; Kostesha, N.; Jauho, A.P.; Burrows, A.; Wubs, M.; Mortensen, N.A. Blueshift of the surface plasmon resonance in silver nanoparticles studied with EELS. Nanophotonics 2013, 2, 131–138. [Google Scholar] [CrossRef]
  28. Haberland, H. Looking from both sides. Nature 2013, 494, E1–E2. [Google Scholar] [CrossRef] [PubMed]
  29. Raza, S.; Kadkhodazadeh, S.; Christensen, T.; Di Vece, M.; Wubs, M.; Mortensen, N.A.; Stenger, N. Multipole plasmons and their disappearance in few-nanometer silver nanoparticles. Nat. Commun. 2015, 6, 8788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Fuchs, R.; Kliewer, K.L. Surface Plasmon in a Semi-Infinite Free-Electron Gas. Phys. Rev. B 1971, 3, 2270–2278. [Google Scholar] [CrossRef]
  31. Ruppin, R. Optical Properties of a Plasma Sphere. Phys. Rev. Lett. 1973, 31, 1434–1437. [Google Scholar] [CrossRef]
  32. Ruppin, R. Optical properties of spatially dispersive dielectric spheres. J. Opt. Soc. Am. 1981, 71, 755–758. [Google Scholar] [CrossRef]
  33. Dasgupta, B.B.; Fuchs, R. Polarizability of a small sphere including nonlocal effects. Phys. Rev. B 1981, 24, 554–561. [Google Scholar] [CrossRef]
  34. Fuchs, R.; Claro, F. Multipolar response of small metallic spheres: Nonlocal theory. Phys. Rev. B 1987, 35, 3722–3727. [Google Scholar] [CrossRef]
  35. Rojas, R.; Claro, F.; Fuchs, R. Nonlocal Response of a small coated sphere. Phys. Rev. B 1988, 37, 6799–6808. [Google Scholar] [CrossRef]
  36. Ruppin, R. Optical absorption by excitons in microcrystals. J. Phys. Chem. Solids 1989, 50, 877–882. [Google Scholar] [CrossRef]
  37. Ruppin, R. Optical absorption by a small sphere above a substrate with inclusion of nonlocal effects. Phys. Rev. B 1992, 45, 11209–11215. [Google Scholar] [CrossRef]
  38. Kreibig, U.; Vollmer, M. Optical Properties of Metal Clusters; Springer: Berlin, Germany, 1995. [Google Scholar]
  39. Ruppin, R. Extinction properties of thin metallic nanowires. Opt. Commun. 2001, 190, 205–209. [Google Scholar] [CrossRef]
  40. García de Abajo, F.J. Nonlocal Effects in the Plasmons of Strongly Interacting Nanoparticles, Dimers, and Waveguides. J. Phys. Chem. C 2008, 112, 17983–17987. [Google Scholar] [CrossRef] [Green Version]
  41. David, C.; García de Abajo, F.J. Spatial Nonlocality in the Optical Response of Metal Nanoparticles. J. Phys. Chem. C 2011, 115, 19470–19475. [Google Scholar] [CrossRef] [Green Version]
  42. Bloch, F. Zur Bremsung rasch bewegter Teilchen beim Durchgang durch Materie. Z. Phys. 1933, 81, 363–376. [Google Scholar] [CrossRef]
  43. Eguiluz, A.; Ying, S.C.; Quinn, J.J. Influence of the electron density profile on surface plasmons in a hydrodynamic model. Phys. Rev. B 1975, 11, 2118–2121. [Google Scholar] [CrossRef]
  44. Sipe, J.E.; So, V.C.Y.; Fukui, M.; Stegeman, G.I. Analysis of second-harmonic generation at metal surfaces. Phys. Rev. B 1980, 21, 4389–4402. [Google Scholar] [CrossRef]
  45. Leung, P.T. Decay of molecules at spherical surfaces: Nonlocal effects. Phys. Rev. B 1990, 42, 7622–7625. [Google Scholar] [CrossRef]
  46. Aizpurua, J.; Rivacoba, A. Nonlocal effects in the plasmons of nanowires and nanocavities excited by fast electron beams. Phys. Rev. B 2008, 78, 035404. [Google Scholar] [CrossRef]
  47. McMahon, J.M.; Gray, S.K.; Schatz, G.C. Optical Properties of Nanowire Dimers with a Spatially Nonlocal Dielectric Function. Nano Lett. 2010, 10, 3473–3481. [Google Scholar] [CrossRef] [PubMed]
  48. Raza, S.; Toscano, G.; Jauho, A.P.; Wubs, M.; Mortensen, N.A. Unusual resonances in nanoplasmonic structures due to nonlocal response. Phys. Rev. B 2011, 84, 121412. [Google Scholar] [CrossRef]
  49. Wiener, A.; Fernández-Domínguez, A.I.; Horsfield, A.P.; Pendry, J.B.; Maier, S.A. Nonlocal Effects in the Nanofocusing Performance of Plasmonic Tips. Nano Lett. 2012, 12, 3308–3314. [Google Scholar] [CrossRef] [PubMed]
  50. De Ceglia, D.; Campione, S.; Vincenti, M.A.; Capolino, F.; Scalora, M. Low-damping epsilon-near-zero slabs: Nonlinear and nonlocal optical properties. Phys. Rev. B 2013, 87, 155140. [Google Scholar] [CrossRef]
  51. David, C.; Mortensen, N.A.; Christensen, J. Perfect imaging, epsilon-near zero phenomena and waveguiding in the scope of nonlocal effects. Sci. Rep. 2013, 3, 2526. [Google Scholar] [CrossRef] [PubMed]
  52. Teperik, T.V.; Nordlander, P.; Aizpurua, J.; Borisov, A.G. Robust subnanometric plasmon ruler by rescaling of the nonlocal optical response. Phys. Rev. Lett. 2013, 110, 263901. [Google Scholar] [CrossRef] [PubMed]
  53. Luo, Y.; Fernandez-Dominguez, A.I.; Wiener, A.; Maier, S.A.; Pendry, J.B. Surface Plasmons and Nonlocality: A Simple Model. Phys. Rev. Lett. 2013, 111, 093901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Christensen, T.; Yan, W.; Raza, S.; Jauho, A.P.; Mortensen, N.A.; Wubs, M. Nonlocal Response of Metallic Nanospheres Probed by Light, Electrons, and Atoms. ACS Nano 2014, 8, 1745–1758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Mortensen, N.A.; Raza, S.; Wubs, M.; Søndergaard, T.; Bozhevolnyi, S.I. A generalized non-local optical response theory for plasmonic nanostructures. Nat. Commun. 2014, 5, 3809. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Raza, S.; Bozhevolnyi, S.I.; Wubs, M.; Mortensen, N.A. Nonlocal optical response in metallic nanostructures. J. Phys. Condens. Matter 2015, 27, 183204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Toscano, G.; Straubel, J.; Kwiatkowski, A.; Rockstuhl, C.; Evers, F.; Xu, H.; Mortensen, N.A.; Wubs, M. Resonance shifts and spill-out effects in self-consistent hydrodynamic nanoplasmonics. Nat. Commun. 2015, 6, 7132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Yan, W.; Wubs, M.; Mortensen, N.A. Projected Dipole Model for Quantum Plasmonics. Phys. Rev. Lett. 2015, 115, 137403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. David, C.; Christensen, J.; Mortensen, N.A. Spatial dispersion in two-dimensional plasmonic crystals: Large blueshifts promoted by diffraction anomalies. Phys. Rev. B 2016, 94, 165410. [Google Scholar] [CrossRef]
  60. Jacak, J.; Krasnyj, J.; Jacak, W.; Gonczarek, R.; Chepok, A.; Jacak, L. Surface and volume plasmons in metallic nanospheres in semiclassical RPA-type approach; Near-field coupling of surface plasmons with semiconductor substrate. Phys. Rev. B 2010, 82, 035418. [Google Scholar] [CrossRef]
  61. Jacak, W.A. Lorentz Friction for Surface Plasmons in Metallic Nanospheres. J. Phys. Chem. C 2015, 119, 6749–6759. [Google Scholar] [CrossRef]
  62. Kluczyk, K.; Jacak, W. Damping-induced size effect in surface plasmon resonance in metallic nano-particles: Comparison of RPA microscopic model with numerical finite element simulation (COMSOL) and Mie approach. J. Quant. Spectr. Radiat. Transf. 2016, 168, 78–88. [Google Scholar] [CrossRef]
  63. Jacak, W.; Popko, E.; Henrykowski, A.; Zielony, E.; Gwozdz, K.; Luka, G.; Pietruszka, R.; Witkowski, B.; Wachnicki, L.; Godlewski, M.; et al. On the size dependence and spatial range for the plasmon effect in photovoltaic efficiency enhancement. Sol. Energy Mater. Sol. Cells 2016, 147, 1–16. [Google Scholar] [CrossRef]
  64. Kluczyk, K.; David, C.; Jacak, W.A. On quantum approach to modeling of plasmon photovoltaic effect. JOSA B 2017, 34, 2115–2127. [Google Scholar] [CrossRef]
  65. David, C.; García de Abajo, F.J. Surface Plasmon Dependence on the Electron Density Profile at Metal Surfaces. ACS Nano 2014, 8, 9558–9566. [Google Scholar] [CrossRef] [PubMed]
  66. Ciracì, C. Current-dependent potential for nonlocal absorption in quantum hydrodynamic theory. Phys. Rev. B 2017, 95. [Google Scholar] [CrossRef]
  67. De Abajo, F.J.G. Multiple scattering of radiation in clusters of dielectrics. Phys. Rev. B 1999, 60, 6086–6102. [Google Scholar] [CrossRef]
  68. Kreibig, U.; Fragstein, C.V. The Limitation of Electron Mean Free Path in Small Silver Particles. Z. Phys. 1969, 224, 307–323. [Google Scholar] [CrossRef]
  69. Coronado, E.A.; Schatz, G.C. Surface plasmon broadening for arbitrary shape nanoparticles: A geometrical probability approach. J. Chem. Phys. 2003, 119, 3926–3934. [Google Scholar] [CrossRef]
  70. Landau, L.D.; Lifshits, E.M. Field Theory; Nauka: Moscow, Russia, 1973. [Google Scholar]
  71. David, C.; Christensen, J. Extraordinary optical transmission through nonlocal holey metal films. Appl. Phys. Lett. 2017, 110, 261110. [Google Scholar] [CrossRef]
  72. Abrikosov, A.A. Fundamentals of the Theory of Metals; North-Holland: Amsterdam, The Netherlands, 1988. [Google Scholar]
  73. Maack, J.R.; Mortensen, N.A.; Wubs, M. Size-dependent nonlocal effects in plasmonic semiconductor particles. EPL (Europhys. Lett.) 2017, 119, 17003. [Google Scholar] [CrossRef] [Green Version]
  74. Maack, J.R.; Mortensen, N.A.; Wubs, M. Two-fluid hydrodynamic model for semiconductors. Phys. Rev. B 2018, 97, 115415. [Google Scholar] [CrossRef] [Green Version]
  75. Pines, D. Elementary Excitations in Solids; Advanced Book Program; Perseus Books: Reading, MA, USA, 1999. [Google Scholar]
  76. Lang, N.D.; Kohn, W. Theory of metal Surfaces: Charge Density and Surface Energy. Phys. Rev. B 1970, 1, 4555–4568. [Google Scholar] [CrossRef]
  77. Low, F.E. Classical Field Theory: Electromagnetism and Gravitation; John Wiley & Sons: New York, NY, USA, 1997. [Google Scholar]
  78. Jackson, J.D. Classical Electrodynamics, 2nd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 1975. [Google Scholar]
  79. Tserkezis, C.; Stefanou, N.; Wubs, M.; Mortensen, N.A. Molecular fluorescence enhancement in plasmonic environments: Exploring the role of nonlocal effects. Nanoscale 2016, 8, 17532–17541. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Illustration of sources of plasmon damping and electron interaction phenomena. (a) Electron-electron collisions in the bulk material; (b) Electron-surface collisions due to confinement; (c) Electron irradiation due to acceleration during plasmon oscillation; (d) Short-ranged electron-electron interactions, such as Coulomb force and electron diffusion.
Figure 1. Illustration of sources of plasmon damping and electron interaction phenomena. (a) Electron-electron collisions in the bulk material; (b) Electron-surface collisions due to confinement; (c) Electron irradiation due to acceleration during plasmon oscillation; (d) Short-ranged electron-electron interactions, such as Coulomb force and electron diffusion.
Materials 11 01077 g001
Figure 2. Impact of quantum corrections on single nanoparticles. (a) Spectral position of the localized surface plasmon resonance (LSPR) for gold, silver and aluminum; (b) Extinction cross section normalized to the surface of a hemisphere for silver evaluated at the respective LSPR wavelengths from (a).
Figure 2. Impact of quantum corrections on single nanoparticles. (a) Spectral position of the localized surface plasmon resonance (LSPR) for gold, silver and aluminum; (b) Extinction cross section normalized to the surface of a hemisphere for silver evaluated at the respective LSPR wavelengths from (a).
Materials 11 01077 g002
Figure 3. Extinction cross section normalized to the surface of a hemisphere for isolated (a) gold and (b) aluminum nanoparticles evaluated at the respective LSPR wavelengths from Figure 2a.
Figure 3. Extinction cross section normalized to the surface of a hemisphere for isolated (a) gold and (b) aluminum nanoparticles evaluated at the respective LSPR wavelengths from Figure 2a.
Materials 11 01077 g003
Figure 4. Maximum enhancement factor EF = | E | 2 / | E 0 | 2 at the particle surface for gold. Dependence of (a) the maximum EF and (b) its wavelength position for the different quantum corrections on the particle radius in water. (c), (d) The same as a function of the permittivity ϵ 0 of the surrounding medium for nanospheres of (c) R = 10 nm and (d) R = 50 nm.
Figure 4. Maximum enhancement factor EF = | E | 2 / | E 0 | 2 at the particle surface for gold. Dependence of (a) the maximum EF and (b) its wavelength position for the different quantum corrections on the particle radius in water. (c), (d) The same as a function of the permittivity ϵ 0 of the surrounding medium for nanospheres of (c) R = 10 nm and (d) R = 50 nm.
Materials 11 01077 g004
Figure 5. Size regime for multipolar response in metal nanoparticles. Enhancement factor EF = | E | 2 / | E 0 | 2 at the particle surface, where the EF is maximized, for λ = 500 nm close to the corresponding Mie resonance of gold with classical Mie coefficients and as inset with combined microscopic corrections. The calculations are based on the dipolar response (black), the first three multipoles (red) and the converged result (blue).
Figure 5. Size regime for multipolar response in metal nanoparticles. Enhancement factor EF = | E | 2 / | E 0 | 2 at the particle surface, where the EF is maximized, for λ = 500 nm close to the corresponding Mie resonance of gold with classical Mie coefficients and as inset with combined microscopic corrections. The calculations are based on the dipolar response (black), the first three multipoles (red) and the converged result (blue).
Materials 11 01077 g005
Figure 6. Impact of microscopic electron dynamics on gold dimers in water. We show the maximum field enhancement at the gap center of gold dimers dispersed in water in dependence of their radius ( a > 0.5 nm) and separation (>0.1 nm) for (a) classical Mie calculations, (b) spatial dispersion with GNOR, (c) Kreibig damping, (d) all RPA corrections combined, (e) Lorentz friction and (f) GNOR with Lorentz friction. The incident field is polarized along the dimer axis and the maximum EF is evaluated at the respective resonance frequency calculated for each case.
Figure 6. Impact of microscopic electron dynamics on gold dimers in water. We show the maximum field enhancement at the gap center of gold dimers dispersed in water in dependence of their radius ( a > 0.5 nm) and separation (>0.1 nm) for (a) classical Mie calculations, (b) spatial dispersion with GNOR, (c) Kreibig damping, (d) all RPA corrections combined, (e) Lorentz friction and (f) GNOR with Lorentz friction. The incident field is polarized along the dimer axis and the maximum EF is evaluated at the respective resonance frequency calculated for each case.
Materials 11 01077 g006
Figure 7. Comparison of RPA damping rates. The perturbative solution (red) and exact Lorentz friction (blue) for a Au nanoparticle in water.
Figure 7. Comparison of RPA damping rates. The perturbative solution (red) and exact Lorentz friction (blue) for a Au nanoparticle in water.
Materials 11 01077 g007

Share and Cite

MDPI and ACS Style

Kluczyk, K.; Jacak, L.; Jacak, W.; David, C. Microscopic Electron Dynamics in Metal Nanoparticles for Photovoltaic Systems. Materials 2018, 11, 1077. https://doi.org/10.3390/ma11071077

AMA Style

Kluczyk K, Jacak L, Jacak W, David C. Microscopic Electron Dynamics in Metal Nanoparticles for Photovoltaic Systems. Materials. 2018; 11(7):1077. https://doi.org/10.3390/ma11071077

Chicago/Turabian Style

Kluczyk, Katarzyna, Lucjan Jacak, Witold Jacak, and Christin David. 2018. "Microscopic Electron Dynamics in Metal Nanoparticles for Photovoltaic Systems" Materials 11, no. 7: 1077. https://doi.org/10.3390/ma11071077

APA Style

Kluczyk, K., Jacak, L., Jacak, W., & David, C. (2018). Microscopic Electron Dynamics in Metal Nanoparticles for Photovoltaic Systems. Materials, 11(7), 1077. https://doi.org/10.3390/ma11071077

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop