PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nanoreslettspringer open web sitethis articlemanuscript submissionregistrationjournal front pagespringer open web site
 
Nanoscale Res Lett. 2014; 9(1): 174.
Published online 2014 April 10. doi:  10.1186/1556-276X-9-174
PMCID: PMC3996189

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

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 equation M1 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, equation M2 where equation M3 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 p x ,p y , and p z 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: equation M4, 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 r s =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 Figure1.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.

Figure 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 ...

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 left angle bracketn s right angle bracket, variances of the occupation numbers equation M5, and sums of the variances equation M6 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 Figure22 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 Figure2a.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.

Figure 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 Δ ...

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:

equation M7
(1)

where e and m e are the charge and effective mass of the electron, respectively, n e =N/V, ε=p2 / (2 m e ),f(ε) is an equilibrium electron energy distribution function, and τ is the relaxation time. From Equation 1, the classical result equation M8 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:

equation M9
(2)

where left angle bracketn s right angle bracket is the averaged occupation number of the state s. We tested Equation 2 by computing the normalized conductivity defined at constant τ,

equation M10
(3)

The equality equation M11 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 Figure3a3a presents the results obtained at N in the range from 2,000 to 2.5×105,T=300 K. There are pulsations of equation M12vanishing 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 Figure3a3a (at a from 2 to 4.5 nm, i.e., N from 2,000 to 20,000) shows the oscillations of equation M13 with the amplitude increasing with the decrease of a.

Figure 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 ...

The conductivity equation M14at 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 Figure3b,c,3b,c, the value of equation M15is 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 equation M16. The values of Nm and equation M17 are listed in Table Table1.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 Figure11.

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

All the experimental numbers Nm in Figure Figure11 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 equation M19 cluster ions. Other magic numbers of atoms were 137 for equation M20, equation M21, and equation M22 and 139 for equation M23, equation M24, and equation M25. 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 equation M26on 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 equation M27 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 equation M28for 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 equation M29. The normalized conductivity of the neutral sphere is found to be equation M30; equation M31 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

equation M32
(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

equation M33
(5)

where Δ is the sum of the variances of the occupation numbers shown in Figure Figure22 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 [similar, equals] 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

  • Kreibig U. Electronic properties of small silver particles: the optical constants and their temperature dependence. J Phys F. 1974;9:999–1014. doi: 10.1088/0305-4608/4/7/007. [Cross Ref]
  • Roldughin VI. Quantum-size colloid metal systems. Russ Chem Rev. 2000;9(10):821–844. doi: 10.1070/RC2000v069n10ABEH000605. [Cross Ref]
  • Chen M, Cai Y, Yan Z, Goodman DW. On the origin of the unique properties of supported Au nanoparticles. J Am Chem Soc. 2006;9(19):6341–6346. doi: 10.1021/ja0557536. [PubMed] [Cross Ref]
  • Mikkellä M-H. Experimental study of nanoscale metal clusters using synchrotron radiation excited photoelectron spectroscopy. Academic dissertation, University of Oulu; 2013.
  • 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;9:229–236. doi: 10.1016/0168-1176(85)80021-5. [Cross Ref]
  • 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;9:33–41. doi: 10.1016/0168-1176(86)85021-2. [Cross Ref]
  • Göhlich H, Lange T, Bergmann T, Martin TP. Electronic shell structure in large metallic clusters. Phys Rev Lett. 1990;9:748–751. doi: 10.1103/PhysRevLett.65.748. [PubMed] [Cross Ref]
  • 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;9:1627–1630. doi: 10.1103/PhysRevLett.65.1627. [PubMed] [Cross Ref]
  • 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;9:209–213. doi: 10.1016/0009-2614(90)85389-T. [Cross Ref]
  • Pedersen J, Bjørnholm S, Borggreen J, Hansen K, Martin TP, Rasmussen HD. Observation of quantum supershells in clusters of sodium atoms. Nature. 1991;9:733–735. doi: 10.1038/353733a0. [Cross Ref]
  • 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;9:53–57. doi: 10.1016/0009-2614(91)80191-Y. [Cross Ref]
  • 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;9:215–222. doi: 10.1016/S0009-2614(91)85132-G. [Cross Ref]
  • Alameddin G, Hunter J, Cameron D, Kappes MM. Electronic and geometric structure in silver clusters. Chem Phys Lett. 1992;9:122–129. doi: 10.1016/0009-2614(92)85439-H. [Cross Ref]
  • 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;9:4741–4745. doi: 10.1103/PhysRevB.48.4741. [PubMed] [Cross Ref]
  • 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;9:2271–2277. doi: 10.1103/PhysRevB.47.2271. [PubMed] [Cross Ref]
  • 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;9:013203.
  • Genzken O, Brack M. Temperature dependence of supershells in large sodium clusters. Phys Rev Lett. 1991;9:3286–3289. doi: 10.1103/PhysRevLett.67.3286. [PubMed] [Cross Ref]
  • Koch E, Gunnarsson O. Density dependence of the electronic supershells in the homogeneous jellium model. Phys Rev B. 1996;9:5168–5177. doi: 10.1103/PhysRevB.54.5168. [PubMed] [Cross Ref]
  • Lundstrom M. Fundamentals of Carrier Transport. Cambridge: Cambridge University Press; 2000.
  • Dressel M, Grüner G. Electrodynamics of Solids. Optical Properties of Electrons in Matter. New York: Cambridge University Press; 2002.
  • 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;9(3):727–736.
  • Narumanchi SVJ, Murthy JY, Amon CH. Boltzmann transport equation-based thermal modeling approaches for hotspots in microelectronics. Heat and Mass Transfer. 2006;9(6):478–491. doi: 10.1007/s00231-005-0645-6. [Cross Ref]
  • Datsyuk VV. A generalization of the Mie theory for a sphere with spatially dispersive permittivity. Ukr J Phys. 2011;9(2):122–129.
  • Bohren CF, Huffman DR. Absorption and Scattering of Light by Small Particles. New York: Wiley; 1983.
  • Datsyuk VV, Tovkach OM. Optical properties of a metal nanosphere with spatially dispersive permittivity. J Opt Soc Am B. 2011;9(5):1224–1230. doi: 10.1364/JOSAB.28.001224. [Cross Ref]
  • Raza S, Yan W, Stenger N, Wubs M, Mortensen NA. Blueshift of the surface plasmon resonance in silver nanoparticles: substrate effects. Opt Expr. 2013;9:27344–27355. doi: 10.1364/OE.21.027344. [PubMed] [Cross Ref]
  • Hilger A, Tenfelde M, Kreibig U. Silver nanoparticles deposited on dielectric surfaces. Appl Phys B. 2001;9(4):361–372. doi: 10.1007/s003400100712. [Cross Ref]
  • Ashcroft NW, Mermin ND. Solid State Physics. New York: Holt, Rinehart and Winston; 1979.
  • Datsyuk VV, Ivanytska IV. Statistical properties of conduction electrons in an isolated metal nanosphere. J Stat Phys. 2013;9:969–978. doi: 10.1007/s10955-013-0798-5. [Cross Ref]
  • 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;9:4493–4499. doi: 10.1364/AO.24.004493. [PubMed] [Cross Ref]
  • Mermin ND. Lindhard dielectric function in the relaxation-time approximation. Phys Rev B. 1970;9:2362–2363. doi: 10.1103/PhysRevB.1.2362. [Cross Ref]
  • Kliewer KL, Fuchs R. Lindhard dielectric functions with a finite electron lifetime. Phys Rev. 1969;9:552–558. doi: 10.1103/PhysRev.181.552. [Cross Ref]
  • Wood DM. Classical size dependence of the work function of small metallic spheres. Phys Rev Lett. 1981;9:749–749. doi: 10.1103/PhysRevLett.46.749. [Cross Ref]

Articles from Nanoscale Research Letters are provided here courtesy of Springer