Open Access Nano Express

Statistical and electrical properties of the conduction electrons of a metal nanosphere in the region of metal-insulator transition

Vitaly V Datsyuk* and Iryna V Ivanytska

Author Affiliations

Department of Physics, Taras Shevchenko National University of Kyiv, 64, Volodymyrska Str., Kyiv 01601, Ukraine

For all author emails, please log on.

Nanoscale Research Letters 2014, 9:174  doi:10.1186/1556-276X-9-174


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/9/1/174


Received:10 December 2013
Accepted:25 March 2014
Published:10 April 2014

© 2014 Datsyuk and Ivanytska; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Abstract

Statistical and electrical properties of the conduction electrons of a silver or gold sphere with a radius from 1 to 2 nm are shown to differ drastically from the properties of electrons in a bulk metal sample. If the radius of a noble metal sphere decreases from 10 to 1 nm, its conductivity oscillates around the bulk metal value with increasing amplitude and drops at the 'magic’ numbers of electrons. These numbers are equal to 186, 198, 254, 338, 440, 556, 676, 832, 912, 1,284, 1,502, and 1,760, in agreement with various experimental data. We show that the conductivity and capacitance of a metal nanosphere can be changed by several orders of magnitude by adding or removing just a few electrons.

Keywords:
Nanoparticle; Metal-insulator transition; Conductivity; Capacitance

Background

In the past, a measurement of optical absorption by silver nanoparticles embedded in glass showed that the particles had normal metallic properties when their diameters were decreased down to 2.2 nm [1]. Contrary to this finding, metal particles with sizes below 2 nm cannot be conducting according to more recent papers [2,3]. Very recently, it was understood that the metal-insulator transition (MIT) is gradual so that particles with certain 'magic’ numbers of electrons become insulating while others remain conducting [4]. If electrons move inside a sphere, then the numbers 186, 198, 254, 338, 440, 556, 676, 832, 912, 1,284, 1,502, and 1,760 are known to be 'magic’. It was experimentally found that the above numbers are indeed magic for clusters of many metals [5-16]. This allows one to consider the motion of electrons in a spherical jellium [8,12,17,18].

We recently studied statistical properties of 500 to 2,000 free electrons confined in a spherical potential well with a radius from 1.2 to 2 nm. The averaged occupation numbers of the electron energy levels and the variances of the occupation numbers were computed for both isolated metal nanoparticles and those in equilibrium with an electron bath. The sum of the variances of all occupation numbers was found to depend on the number of electrons nonmonotonically dropping by a few orders of magnitude at 'magic numbers’ of electrons. Here, we show how the statistical properties of the conduction electrons are related with the electrical properties of metal nanoparticles. Calculations of the DC conductivity and capacitance of single nanometer-sized noble metal spheres are reported. We predict a transistor-like behavior of a single nanoparticle when an additional charge of the particle drastically changes its conductivity and capacitance.

Methods

Statistical and transport models

The electron statistics and capacitance of metal nanoparticles are investigated by the Gibbs ensemble method. The transport properties of electrons are studied using a quantum-mechanic generalization of a formula that can be obtained by solving the Boltzmann transport equation.

The Boltzmann transport equation (BTE) is widely used in modern physical kinetics [19,20] even at nanoscales [21,22]. In our previous work, we solved BTE and determined the response of the conduction electrons on an electromagnetic wave allowing for dependence of the distribution function on the wavenumber. Then, we obtained the spatially dispersive permittivity of metal and applied a generalized Mie theory [23] to calculate spectra of light extinction by silver nanospheres. Because of excitation of the rotationless (longitudinal) plasmon-polariton waves, the frequency of the Fröhlich resonance ω=3.5 eV [24] was found to increase with decreasing the radius a, so ω=3.635 eV for particles with a=1.5 nm and ω=3.73 eV at a=1 nm (see Figure two of [25]; a similar dependence of ω on a as well as the formulas of [23] is found in a recent paper [26]). The blueshift by 0.15 eV and the width of the plasmon resonance calculated for a particle beam created by Hilger, Tenfelde, and Kreibig [27] were in excellent agreement with the experimental data. We concluded that silver clusters with a<1 nm do not contribute into the extinction spectra of the beams. However, it was not possible to establish whether the contribution of these ultrathin particles into the integral extinction spectrum vanished due to MIT.

Energies of the conduction electrons

It is a common practice to assume that the conduction electron energy ε in metal is a function of the continuous quasi-momentum <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M1">View MathML</a> This approach can be used to model the properties of bulk metals in which an electron state s can be described by a set of four quantum numbers, <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M2">View MathML</a> where <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M3">View MathML</a> is a projection of the electron spin. However, if an electron moves inside a microscopic sphere, the vector p can no longer describe an electron state. The set px,py, and pz has to be replaced by a set of the radial q, angular momentum l, and magnetic m quantum numbers. Then, integrals over p should be replaced by sums over the electron states s, according to the following rule: <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M4">View MathML</a>, where V is the volume of the body and h is the Plank constant. In this study, we used discrete sets of the electron energies εs.

Shape of metal nanoparticles

Every set of the electron energies εs was obtained at a certain number N of electrons confined in a potential well of an ideal spherical shape. We used the parameters of bulk silver and gold [28], namely, the depth of the potential well U=9.8 eV and the Wigner-Seitz radius rs=0.16 nm. Because of the spherical symmetry of the electronic wave functions, different physical processes have peculiar features at magic numbers of electrons. The magic numbers, which we determined for perfect noble metal spheres, are presented in the 'Background’ Section of this paper. Surprisingly, almost the same values were determined in various experiments working with particles of different metals, even though no special efforts were made to produce perfect spheres. In addition, different theoretical papers also reported similar magic numbers, according to Figure 1. This means that effects associated with the peculiarities of the spacing of εs in spherical nanoparticles are sensitive neither to surface distortions nor the values of the parameters U and rs.

thumbnailFigure 1. Experimental (centered boxes with error bars) and theoretical (crosses) 'magic’ numbers of electrons in metal clusters. Solid grid lines indicate Nm= 186, 198, 254, 338, 440, 556, 676, 760, 832, 912, 1,012, 1,100, 1,284, 1,502, and 1,760. Dashed grid lines indicate Nm= 268, 542, 1,074, and 1,206.

Results and discussion

Variances of the occupation numbers

In our previous work [29], we reported statistical properties of the conduction electrons in isolated metal nanospheres. To study the systems with a fixed number of electrons, the method of the canonical ensemble was applied. The averaged occupation numbers 〈ns〉, variances of the occupation numbers <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M5">View MathML</a>, and sums of the variances <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M6">View MathML</a> were computed and discussed. In [29], we also examined the properties of the conduction electrons in grand canonical ensembles where the chemical potential μ0 was fixed. Figure 2 represents the values of Δ calculated at fixed N (canonical ensembles) and μ0 (grand canonical ensembles). The sum of the variances depends on the number of electrons nonmonotonically dropping by several orders of magnitude at magic numbers of electrons. The decrease in Δ can occur if (i) the distance between the Fermi level and the neighboring higher energy level, εf+1-εf, is large compared to the thermal energy and (ii) the Fermi level is fully occupied at absolute zero temperature. Addition of one atom to a particle with Nm conduction electrons results in a substantial increase in the Fermi energy, as is evident from Figure 2a. If a particle has a magic number of electrons, the chemical potential lies in the gap between the distant energy levels, so the number of the current carriers is greatly reduced. The influence of this effect on the electrical properties of the metal nanoparticles is studied below.

thumbnailFigure 2. Fermi energies and variances of the occupation numbers of electronic states of single Ag or Au spheres. (a) Fermi energy as a function of the number N of conduction electrons. (b) Sums of the variances Δ normalized to the bulk metal value Δb in canonical ensembles (points) and grand canonical ones (crosses). The grid lines are the same as in Figure 1.

Conductivity

The response of the conduction electrons of metals to an infrared and far infrared radiation is well described by a Drude dielectric function [30]. In the corresponding limit of small emission wavenumbers, this function can be derived by using either a quantum theory by Lindhard [31] or the classical Boltzmann transport equation [32] (see derivations in [20]). The permittivity of metal is expressed through the conductivity σ that, in the DC case, can be calculated from the following equation:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M7">View MathML</a>

(1)

where e and me are the charge and effective mass of the electron, respectively, ne=N/V, ε=p2 / (2 me),f(ε) is an equilibrium electron energy distribution function, and τ is the relaxation time. From Equation 1, the classical result <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M8">View MathML</a> is obtained at τ=const and any f(ε) finite at ε=0 and vanishing at ε. The formula for σb can also be derived by substituting a zero-temperature Fermi-Dirac distribution function into Equation 1.

A generalization of Equation 1 for discrete energy levels gives the following formula:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M9">View MathML</a>

(2)

where 〈ns〉 is the averaged occupation number of the state s. We tested Equation 2 by computing the normalized conductivity defined at constant τ,

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M10">View MathML</a>

(3)

The equality <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M11">View MathML</a> should hold for 'large’ particles since properties of a macroscopic body are independent of the boundary conditions for the electron wave function. The calculations were performed by using sets of εs for N free electrons confined in a spherical potential well with the radius a=rsN1/3, where rs=0.16 nm. Figure 3a presents the results obtained at N in the range from 2,000 to 2.5×105,T=300 K. There are pulsations of <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M12">View MathML</a>vanishing as sphere radii increase above 9 nm that corresponds to N>2×105. Therefore, Equation 2 works well, and particles with a≥10 nm can be regarded as macroscopic. The left-hand side of the curve in Figure 3a (at a from 2 to 4.5 nm, i.e., N from 2,000 to 20,000) shows the oscillations of <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M13">View MathML</a> with the amplitude increasing with the decrease of a.

thumbnailFigure 3. Normalized DC conductivity.(a) Normalized DC conductivity vs rigid-wall sphere radius a=rsN1/3 at N from 2,000 to 2.5×105. Normalized DC conductivity of a neutral silver or gold sphere at (b)N= 180 to 382 and (c)N= 382 to 2,000. The grid lines are the same as in Figure 1.

The conductivity <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M14">View MathML</a>at N= 200 to 2,000 was calculated by using more realistic values of εs found for a spherical potential well with the parameters of silver and gold. According to Figure 3b,c, the value of <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M15">View MathML</a>is not a monotonic function of N and drops sharply when N is equal to one of the magic numbers Nm. The appearance of magic numbers is a general property of fermionic systems. In this paper, the magic numbers of the conduction electrons are identified by the dips in the conductivity <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M16">View MathML</a>. The values of Nm and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M17">View MathML</a> are listed in Table 1. The found values of Nm are in excellent agreement with the experimental and theoretical magic numbers of clusters of many metals according to Figure 1.

Table 1. Normalized conductivity (%) calculated for an Ag or Au particle with a magic number of atoms

All the experimental numbers Nm in Figure 1 were obtained by using the mass spectroscopy from dips in the mass spectra. For example, Katakuse and co-workers [6] found magic numbers of atoms equal to 197 for negative cluster ions of silver (Ag)n- and 199 for positive <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M19">View MathML</a> cluster ions. Other magic numbers of atoms were 137 for <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M20">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M21">View MathML</a>, and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M22">View MathML</a> and 139 for <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M23">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M24">View MathML</a>, and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M25">View MathML</a>. In these cases, the negatively and positively charged and neutral particles had the same magic numbers of electrons, Nm= 198 and 138. Thus, the anomalous properties of the metal nanoparticles in the experiments [5-15] are determined by electron motion [29] but not their atomic structure. Moreover, the model of single electrons trapped in a spherical potential well was shown to be adequate [6] though the shape of the clusters obtained by the bombardment of metal sheets with Xe ions was not controlled.

A nonlinear dependence of <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M26">View MathML</a>on N can occur even in a single sphere if N varies around Nm. To examine electric properties of a single charged nanoparticle, let us consider a sphere in thermal equilibrium with a reservoir of electrons, so the electrochemical potential μ=μ0+eϕ is constant inside the sphere; here μ0 is the chemical potential of the neutral sphere and ϕ is the electric potential. For a fixed μ, we determined <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M27">View MathML</a> by using the Fermi-Dirac occupation numbers and computed the charge of the sphere Q=e (N-N0), where N0 is the number of electrons in the neutral particle. We calculated the quantities Q and <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M28">View MathML</a>for a charged 336-atom Ag or Au nanoparticle. We found that the 336-atom particle holds two extra electrons when the value ϕ changes in a wide range of about 0.6 V. If the mean number of electrons in the particle is equal to 338, then <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M29">View MathML</a>. The normalized conductivity of the neutral sphere is found to be <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M30">View MathML</a>; <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M31">View MathML</a> In the considered example, the neutral sphere is conductive, but the charged one with two extra electrons turns out to be an insulator.

Capacitance

A parameter that describes the dependence of Q on ϕ is the electric capacitance

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M32">View MathML</a>

(4)

A straightforward calculation of the derivative of Q gives the capacitance of the charged particle with 338 electrons C=6.1×10-22 F that is much lower than C=1.1×10-17 F of the neutral 336-atom sphere. The change in the capacitance C(338)/C(336)=5.3×10-5 is similar to the the correspondent change in the conductivity.

By calculating the derivative of Q in Equation 4 at N defined through the Fermi-Dirac occupation numbers, we get

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/174/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/174/mathml/M33">View MathML</a>

(5)

where Δ is the sum of the variances of the occupation numbers shown in Figure 2 by crosses. Equation 5 expresses the relation between the reaction of the conduction electrons to the electric field and the fluctuations of the occupation numbers of the electron states. Thus, the peculiarities of spacing and degeneracy of the electronic energy levels have similar effects on the statistical and electrical properties of a nanometer-sized particle.

During the calculations we neglected Coulomb effects. These effects are as follows. When an electron leaves a neutral metal sphere, it overcomes the attraction of the positive charge remaining on the sphere. Consequently, the work function increases by the value ΔU = 0.54/a(nm) eV [33]. For example, ΔU ≃ 0.5 eV for a 338-atom noble-metal sphere. In addition, Coulomb forces should change the binding energy of the electron in a charged sphere compared to that in the neutral one. Though the change in U can be large, it should not be critical to the effects studied in this paper. Indeed, they depend on the value of εf+1-εf, but this difference is a weak function of U. For example, for a noble metal sphere with 338 conduction electrons, we get εf+1-εf=0.69 eV at U=9.8 eV, and εf+1-εf=0.74 eV if U.

Conclusion

In conclusion, the statistical properties, conductivity, and capacitance of a single nanometer-sized metal sphere depends very strongly on the number of conduction electrons N in the range from 200 to 2,000. In particular, the DC conductivity drops by several orders of magnitude if N is equal to one of the magic numbers. For instance, addition of two electrons to a 336-atom noble metal sphere should reduce both the conductivity and capacitance of the particle by four orders of magnitude.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

VVD and IVI together carried computations, analyzed results, and prepared the manuscript. Both authors read and approved the final manuscript.

References

  1. Kreibig U: Electronic properties of small silver particles: the optical constants and their temperature dependence.

    J Phys F 1974, 4:999-1014. Publisher Full Text OpenURL

  2. Roldughin VI: Quantum-size colloid metal systems.

    Russ Chem Rev 2000, 69(10):821-844. Publisher Full Text OpenURL

  3. Chen M, Cai Y, Yan Z, Goodman DW: On the origin of the unique properties of supported Au nanoparticles.

    J Am Chem Soc 2006, 128(19):6341-6346. PubMed Abstract | Publisher Full Text OpenURL

  4. Mikkellä M-H: Experimental study of nanoscale metal clusters using synchrotron radiation excited photoelectron spectroscopy. Academic dissertation, University of Oulu; 2013. OpenURL

  5. Katakuse I, Ichihara T, Fujita Y, Matsuo T, Sakurai T, Matsuda H: Mass distributions of copper, silver and gold clusters and electronic shell structure.

    Int J Mass Spectrom Ion Process 1985, 67:229-236. Publisher Full Text OpenURL

  6. Katakuse I, Ichihara T, Fujita Y, Matsuo T, Sakurai T, Matsuda H: Mass distributions of negative cluster ions of copper, silver, and gold.

    Int J Mass Spectrom Ion Process 1986, 74:33-41. Publisher Full Text OpenURL

  7. Göhlich H, Lange T, Bergmann T, Martin TP: Electronic shell structure in large metallic clusters.

    Phys Rev Lett 1990, 65:748-751. PubMed Abstract | Publisher Full Text OpenURL

  8. Bjørnholm S, Borggreen J, Echt O, Hansen K, Pedersen J, Rasmussen HD: Mean-field quantization of several hundred electrons in sodium metal clusters.

    Phys Rev Lett 1990, 65:1627-1630. PubMed Abstract | Publisher Full Text OpenURL

  9. Martin TP, Bergmann T, Göhlich H, Lange T: Observation of electronic shells and shells of atoms in large Na clusters.

    Chem Phys Lett 1990, 172:209-213. Publisher Full Text OpenURL

  10. Pedersen J, Bjørnholm S, Borggreen J, Hansen K, Martin TP, Rasmussen HD: Observation of quantum supershells in clusters of sodium atoms.

    Nature 1991, 353:733-735. Publisher Full Text OpenURL

  11. Martin TP, Bjørnholm S, Borggreen J, Bréchignac C, Cahuzac P, Hansen K, Pedersen J: Electronic shell structure of laser-warmed Na clusters.

    Chem Phys Lett 1991, 186:53-57. Publisher Full Text OpenURL

  12. Persson JL, Whetten RL, Cheng H-P, Berry RS: Evidence for quantized electronic level structure for 100-1300 electrons in metal-atomic clusters.

    Chem Phys Lett 1991, 186:215-222. Publisher Full Text OpenURL

  13. Alameddin G, Hunter J, Cameron D, Kappes MM: Electronic and geometric structure in silver clusters.

    Chem Phys Lett 1992, 192:122-129. Publisher Full Text OpenURL

  14. Ito H, Sakurai T, Matsuo T, Ichihara T, Katakuse I: Detection of electronic-shell structure in divalent-metal clusters (Hg)n.

    Phys Rev B 1993, 48:4741-4745. Publisher Full Text OpenURL

  15. Bréchignac C, Cahuzac P, Carlier F, de Frutos M, Roux JP: Temperature effects in the electronic shells and supershells of lithium clusters.

    Phys Rev B 1993, 47:2271-2277. Publisher Full Text OpenURL

  16. Bonatsos D, Karoussos N, Lenis D, Raychev PP, Roussev RP, Terziev PA: Unified description of magic numbers of metal clusters in terms of the three-dimensional q-deformed harmonic oscillator.

    Phys Rev A 2000, 62:013203. OpenURL

  17. Genzken O, Brack M: Temperature dependence of supershells in large sodium clusters.

    Phys Rev Lett 1991, 67:3286-3289. PubMed Abstract | Publisher Full Text OpenURL

  18. Koch E, Gunnarsson O: Density dependence of the electronic supershells in the homogeneous jellium model.

    Phys Rev B 1996, 54:5168-5177. Publisher Full Text OpenURL

  19. Lundstrom M: Fundamentals of Carrier Transport. Cambridge: Cambridge University Press; 2000. OpenURL

  20. Dressel M, Grüner G: Electrodynamics of Solids. Optical Properties of Electrons in Matter. New York: Cambridge University Press; 2002. OpenURL

  21. Jin S, Tang T-W, Fischetti MV: Simulation of silicon nanowire transistors using Boltzmann transport equation under relaxation time approximation.

    Electron Devices IEEE Trans 2008, 55(3):727-736. OpenURL

  22. Narumanchi SVJ, Murthy JY, Amon CH: Boltzmann transport equation-based thermal modeling approaches for hotspots in microelectronics.

    Heat and Mass Transfer 2006, 42(6):478-491. Publisher Full Text OpenURL

  23. Datsyuk VV: A generalization of the Mie theory for a sphere with spatially dispersive permittivity.

    Ukr J Phys 2011, 56(2):122-129. OpenURL

  24. Bohren CF, Huffman DR: Absorption and Scattering of Light by Small Particles. New York: Wiley; 1983. OpenURL

  25. Datsyuk VV, Tovkach OM: Optical properties of a metal nanosphere with spatially dispersive permittivity.

    J Opt Soc Am B 2011, 28(5):1224-1230. Publisher Full Text OpenURL

  26. Raza S, Yan W, Stenger N, Wubs M, Mortensen NA: Blueshift of the surface plasmon resonance in silver nanoparticles: substrate effects.

    Opt Expr 2013, 21:27344-27355. Publisher Full Text OpenURL

  27. Hilger A, Tenfelde M, Kreibig U: Silver nanoparticles deposited on dielectric surfaces.

    Appl Phys B 2001, 73(4):361-372. Publisher Full Text OpenURL

  28. Ashcroft NW, Mermin ND: Solid State Physics. New York: Holt, Rinehart and Winston; 1979. OpenURL

  29. Datsyuk VV, Ivanytska IV: Statistical properties of conduction electrons in an isolated metal nanosphere.

    J Stat Phys 2013, 152:969-978. Publisher Full Text OpenURL

  30. Ordal MA, Bell RJ, Long LL, Querry MR, Alexander RW Jr: Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W.

    Appl Opt 1985, 24:4493-4499. PubMed Abstract | Publisher Full Text OpenURL

  31. Mermin ND: Lindhard dielectric function in the relaxation-time approximation.

    Phys Rev B 1970, 1:2362-2363. Publisher Full Text OpenURL

  32. Kliewer KL, Fuchs R: Lindhard dielectric functions with a finite electron lifetime.

    Phys Rev 1969, 181:552-558. Publisher Full Text OpenURL

  33. Wood DM: Classical size dependence of the work function of small metallic spheres.

    Phys Rev Lett 1981, 46:749-749. Publisher Full Text OpenURL