Semiclassical approach to the description of the basic properties of nanoobjects

Present paper is a review of results, obtained in the framework of semiclassical approach in nanophysics. Semiclassical description, based on Electrostatics and Thomas–Fermi model was applied to calculate dimensions of the electronic shell of a fullerene molecule and a carbon nanotube. This simpli...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2008
Автор: Kornyushin, Y.
Формат: Стаття
Мова:English
Опубліковано: Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України 2008
Назва видання:Физика низких температур
Теми:
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/117564
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Semiclassical approach to the description of the basic properties of nanoobjects / Y. Kornyushin // Физика низких температур. — 2008. — Т. 34, № 10. — С. 1063–1071. — Бібліогр.: 26 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-117564
record_format dspace
spelling irk-123456789-1175642017-05-25T03:03:26Z Semiclassical approach to the description of the basic properties of nanoobjects Kornyushin, Y. Carbon nanotubes, quantum wires and Luttinger liquid Present paper is a review of results, obtained in the framework of semiclassical approach in nanophysics. Semiclassical description, based on Electrostatics and Thomas–Fermi model was applied to calculate dimensions of the electronic shell of a fullerene molecule and a carbon nanotube. This simplified approach yields surprisingly accurate results in some cases. Semiclassical approach provides rather good description of the dimensions of the electronic shell of a fullerene molecule. Two types of dipole oscillations in a fullerene molecule were considered and their frequencies were calculated. Similar calculations were performed for a carbon nanotube also. These results look rather reasonable. Three types of dipole oscillations in carbon nanotube were considered and their frequencies were calculated. Frequencies of the longitudinal collective oscillations of delocalized electrons in carbon peapod were calculated as well. Metallic cluster was modeled as a spherical ball. It was shown that metallic cluster is stable; its bulk modulus and the frequency of the dipole oscillation of the electronic shell relative to the ions were calculated. 2008 Article Semiclassical approach to the description of the basic properties of nanoobjects / Y. Kornyushin // Физика низких температур. — 2008. — Т. 34, № 10. — С. 1063–1071. — Бібліогр.: 26 назв. — англ. 0132-6414 PACS: 73.63.–b http://dspace.nbuv.gov.ua/handle/123456789/117564 en Физика низких температур Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
topic Carbon nanotubes, quantum wires and Luttinger liquid
Carbon nanotubes, quantum wires and Luttinger liquid
spellingShingle Carbon nanotubes, quantum wires and Luttinger liquid
Carbon nanotubes, quantum wires and Luttinger liquid
Kornyushin, Y.
Semiclassical approach to the description of the basic properties of nanoobjects
Физика низких температур
description Present paper is a review of results, obtained in the framework of semiclassical approach in nanophysics. Semiclassical description, based on Electrostatics and Thomas–Fermi model was applied to calculate dimensions of the electronic shell of a fullerene molecule and a carbon nanotube. This simplified approach yields surprisingly accurate results in some cases. Semiclassical approach provides rather good description of the dimensions of the electronic shell of a fullerene molecule. Two types of dipole oscillations in a fullerene molecule were considered and their frequencies were calculated. Similar calculations were performed for a carbon nanotube also. These results look rather reasonable. Three types of dipole oscillations in carbon nanotube were considered and their frequencies were calculated. Frequencies of the longitudinal collective oscillations of delocalized electrons in carbon peapod were calculated as well. Metallic cluster was modeled as a spherical ball. It was shown that metallic cluster is stable; its bulk modulus and the frequency of the dipole oscillation of the electronic shell relative to the ions were calculated.
format Article
author Kornyushin, Y.
author_facet Kornyushin, Y.
author_sort Kornyushin, Y.
title Semiclassical approach to the description of the basic properties of nanoobjects
title_short Semiclassical approach to the description of the basic properties of nanoobjects
title_full Semiclassical approach to the description of the basic properties of nanoobjects
title_fullStr Semiclassical approach to the description of the basic properties of nanoobjects
title_full_unstemmed Semiclassical approach to the description of the basic properties of nanoobjects
title_sort semiclassical approach to the description of the basic properties of nanoobjects
publisher Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
publishDate 2008
topic_facet Carbon nanotubes, quantum wires and Luttinger liquid
url http://dspace.nbuv.gov.ua/handle/123456789/117564
citation_txt Semiclassical approach to the description of the basic properties of nanoobjects / Y. Kornyushin // Физика низких температур. — 2008. — Т. 34, № 10. — С. 1063–1071. — Бібліогр.: 26 назв. — англ.
series Физика низких температур
work_keys_str_mv AT kornyushiny semiclassicalapproachtothedescriptionofthebasicpropertiesofnanoobjects
first_indexed 2025-07-08T12:28:27Z
last_indexed 2025-07-08T12:28:27Z
_version_ 1837081777833771008
fulltext Fizika Nizkikh Temperatur, 2008, v. 34, No. 10, p. 1063–1071 Semiclassical approach to the description of the basic properties of nanoobjects Yuri Kornyushin Maître Jean Brunschvig Research Unit, Chalet Shalva, Randogne CH-3975, Switzerland E-mail: jacqie@bluewin.ch Received November 22, 2007 Present paper is a review of results, obtained in the framework of semiclassical approach in nanophysics. Semiclassical description, based on Electrostatics and Thomas–Fermi model was applied to calculate di- mensions of the electronic shell of a fullerene molecule and a carbon nanotube. This simplified approach yields surprisingly accurate results in some cases. Semiclassical approach provides rather good description of the dimensions of the electronic shell of a fullerene molecule. Two types of dipole oscillations in a fullerene molecule were considered and their frequencies were calculated. Similar calculations were per- formed for a carbon nanotube also. These results look rather reasonable. Three types of dipole oscillations in carbon nanotube were considered and their frequencies were calculated. Frequencies of the longitudinal col- lective oscillations of delocalized electrons in carbon peapod were calculated as well. Metallic cluster was modeled as a spherical ball. It was shown that metallic cluster is stable; its bulk modulus and the frequency of the dipole oscillation of the electronic shell relative to the ions were calculated. PACS: 73.63.–b Electronic transport in nanoscale materials and structures. Keywords: fullerene molecule, carbon nanotube, metallic cluster. 1. Introduction Semiclassical description, based on Electrostatics and Thomas–Fermi model was applied rather successfully to study atomic, molecular and nanoobjects problems [1–5]. Fullerene molecule was studied in [1,3–5]. Carbon nano- tube was studied in [2,3–5]. Carbon peapod was studied in [5]. Fullerene molecule (C60) forms a spherical ball. Let us denote the radius of the sphere, on which the carbon at- oms are situated R f . Chemical bonds determine the value of this radius. The electronic configuration of the constit- uent carbon atom is 1 2 22 2 2s s p . It was assumed [1,3,4] that in a fullerene molecule two 1s electrons of each atom belong to the core (forming the ion itself), two 2s elec- trons form molecular bonds, and two 2 p electrons are delocalized, or free. Hence it was assumed in [1,3,4] that the total number of the delocalized electrons in a fullerene molecule is 120. Another assumption was made, that the delocalized electrons couldn’t penetrate inside the sphere of the ions as they repel each other. The latest experimental data [6] show that R f � 0.354 nm, the total number of the delocalized electrons, N � = 240, the equilibrium internal radius of the electron shell, Rie � 0.279 nm (that is the delocalized electrons do penetrate the sphere of the ions), and the equilibrium ex- ternal radius of the electronic shell, Re � 0.429 nm. In the model considered here it is assumed like in [1,3–5] that the positive charge of the ions, � �eN e( 0 is the electron charge), is distributed homogeneously on the surface of a sphere r R rf� ( is the distance from the cen- ter of a fullerene molecule). The charge of the delocalized electrons, eN , is assumed to be distributed homogene- ously in a spherical layer R r Ri � � (see Fig. 1). Now let us calculate the electrostatic energy of a fullerene molecule [5]. 2. Calculation of the electrostatic energy of a fullerene molecule The electrostatic potential �( )r , arising from electric charge, can be obtained as a solution of the Poisson’s equation [7] �� ��( ) ( )r r� �4 , (1) where �( )r is the charge density. © Yuri Kornyushin, 2008 The uniformly distributed charge density in the spheri- cal layer of the electronic shell of a fullerene molecule is � �( ) ( )r � �3 4 3 3eN/ R Ri . (2) The charge of the delocalized electrons produces the following electrostatic potential inside the spherical layer, where the delocalized electrons are found: � e i ieN/ R R R / r / R /r� � � �[ ( )][( ) ( ) ( )3 3 2 2 33 2 2 . (3) The charge of the ions creates the following potential on the sphere r R f� and outside of it: � i r eN/r( ) � � . (4) At r R f� the positive charge of the sphere creates elec- trostatic potential equal to �eN/R f . According to Eqs. (3), (4), the total electrostatic potential, � � �( ) ( ) ( )r r re i� , is zero at r R� . Hence for r R it is zero also. Inside the neu- tral fullerene molecule, when 0 � �r Ri , the total electro- static potential is as follows: � in i i i feN R R R RR R eN/R� �[ ( ) / ( )]3 2 2 2 . (5) As on the external surface of a neutral fullerene mole- cule, at r R� , the total electrostatic potential is zero, it should be zero on the internal surface, at r Ri� , also. Otherwise the electrons will move to the locations, where their potential energy is smaller. In a state of equilibrium � in � 0. This yields: R R R /i f f� � �0 5 32 2 1 2. { [ ( )] } ,/� � (6) where � � �R Ri is the thickness of the electronic shell. The latest experimental data [6] are: R f � 0 354. nm and � � 0 15. nm. For these values Eq. (6) yields Ri � 0.274 nm. The experimental value is 0.279 nm. The calculated value is in very good agreement with the observed one. The total electrostatic energy of a system, U , consists of the electrostatic energy of the ions, electrostatic energy of the delocalized electrons, and electrostatic energy of the interaction between them. It is described by the fol- lowing expression: U e N / R e N / Rf� �( ) ( )2 2 2 22 2 � � � � [ ( )]( )e N / R R R R R R Rf i f i f 2 2 3 3 2 3 32 3 2 � � �0 1 9 5 53 3 2 5 5 6 3 2. [ ( )] [ ( ) ]eN/ R R R R R /R R Ri i i i . (7) Here the first term describes the electrostatic energy of the ions. The third term describes the interaction energy. The second and the fourth term describe the electrostatic energy of the delocalized electrons. The second term rep- resents the energy of the electrostatic field outside the charged spherical layer. The fourth term represents the energy of the field inside the spherical layer. This term was calculated as an integral over the layer of the square of the gradient of � e r( ), divided by 8�. 3. Calculation of the kinetic and total energies of a fullerene molecule The value of the kinetic energy, T , of a gas of delo- calized electrons, confined in a volume of regarded spher- ical layer 4 33 3�( )R R /i� , and calculated in accord with the ideas of Thomas–Fermi model, is [1] T /m N / R Ri /� �1105 2 5 3 3 3 2 3. ( ) ( )/ � . (8) The total energy, W T U� , is a function of a single variational parameter, �. Indeed, when Eq. (6) is substi- tuted in Eqs. (7), (8) together with R Ri� � �� 0 5. � �{ [ ( )] }/R R /f f� �2 2 1 23 , we see that the total energy is a function of a single variational parameter, �. In a state of equilibrium the total energy reaches a minimum value Wmin � 9.422 keV. The minimum (equilibrium) value of the total energy is reached when � �� �e 0.1374 nm. At that R Ri ie� � 0.2808 nm, and R Re� � 0.4182 nm. The minimum (equilibrium) energy of a positive ion of a fullerene molecule (with N = 239) was calculated to be 10.327 keV. At that in this case Rie = 0.2808 nm and Re � = 0.4182. The difference, appearing in further digits, is very small. It is worthwhile noting that the total energy of a ful- lerene molecule also contains large negative terms, corre- sponding to the energy of the ionic cores and the energy of the chemical bonds. These terms were not calculated here. 1064 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 Yuri Kornyushin Fig. 1. Model, applied to a fullerene molecule and a carbon nanotube [5]. The charge of the ions is distributed homoge- neously on the surface of a sphere of a radius R f , or cylinder of a radius Rn for a carbon nanotube. The delocalized electrons are distributed uniformly in a spherical (or cylindrical for a carbon nanotube) layer R r Ri � � . 4. Collective dipole oscillation of delocalized electrons in a fullerene molecule Linear-response theory was used by G.F. Bertsch et al. to calculate the frequency of a giant collective resonance in a fullerene molecule [8]. The value, calculated in [8] was of about 20 eV. This value, predicted in [8], was mea- sured experimentally and reported in [9]. Two peaks of dipole collective oscillations of delocalized electrons in C60 positive ions, observed ex- perimentally by S.W.J. Scully et al., were reported in [10]. The authors associated the lower peak near 20 eV with a surface plasmon, its frequency being 3 times smaller than the Langmuir frequency [2]. The other peak (about 40 eV) is associated by the authors with some vol- ume plasmon of unknown origin [10]. Mie Oscillation. Mie oscillation is a collective delo- calized electron oscillation in a thin surface layer of a conducting sphere [11,12], where positive and negative charges are not separated in space. So it does not look suitable to regard such an oscillation in the thin surface layer of the electronic shell of a fullerene molecule. Any- way let us estimate possible frequency of such oscillation. Let us consider a thin surface layer of a thickness d, situ- ated on the external surface of the electronic shell of a fullerene molecule. The volume of this layer, 4 2�R de , con- tains electric charge 4 2�Re end (here n is the number of delocalized electrons per unit volume of the electronic shell). Let us shift the layer as a whole by the distance s along the z-axis. This shift creates a dipole moment P R ensde� 4 2� and restoring electric field (in a spherically symmetric object) E P/ R ensd/r e� �4 3 4 32� �( ) corre- spondingly. Restoring force at that is ( )4 32�enR sd/e . The mass of the regarded thin surface layer is 4 2�R nmde (here m denotes the mass of the delocalized electron). Restoring force leads to the dipole oscillation. From Newton equa- tion follows that the frequency of the oscillation consid- ered is M � ( ) /4 32 1 2�e n/ m . This frequency, as was men- tioned by S.W.J. Scully et al. [10], is 3 times smaller than Langmuir frequency. In the model described above n N/ R Re ie� �3 4 3 3�( ). From two last equations written above follows that M e iee N/m R R� �[ ( )] /2 3 3 1 2. This equa- tion yields � M � 21.45 eV when the experimental val- ues of the parameters are used. S.W.J. Scully et al. re- ported the value of one of the two measured peaks, � M , near 20 eV [10]. The agreement is quite reasonable. Using theoretical values of Re = 0.4182 nm and Rie = 0.2808 nm, calculated here, one can obtain the value of the peak 22.72 eV. This result is not too bad either. Simple Dipole Oscillation. Let us consider a dipole type collective quantum oscillation of the electronic shell as a whole (having mass mN ) relative to the ion skeleton. Such an oscillation looks more plausible in a fullerene molecule than Mie oscillation. Regarded system consists of two spherically symmetric objects, the ion skeleton and the electronic shell. Hence this oscillation causes no electric field in the blank interior of a fullerene molecule around its center (cavity). The value of the electrostatic potential inside the electronic shell is � e eeN/ R� �[ ( 3 � � �R R / r / R /rie e ie 3 2 2 33 2 2)][( ) ( ) ( )] [see Eq. (3)]. Let us consider the shift of the electronic shell by a small dis- tance, s, relative to the ion skeleton and calculate the change in the electrostatic energy of a fullerene molecule (in the framework of the accepted model). The first term in Eq. (3) does not contribute to the change in the electro- static energy, as it is constant. The contribution of the third term is also zero because the third term depends on the distance like 1/r, a potential produced by a point charge. Interaction of the ion skeleton with the point charge, situated near the center of the fullerene molecule, does not depend on their relative position, because the electrostatic potential is constant for r R f� . The change in the electrostatic energy does not depend on the direc- tion of the shift (s vector) because of the spherical sym- metry of the problem. This means that the change cannot contain term proportional to s. Let us choose the x-axis along the s vector. Then after the shift we shall have in the third term of Eq. (3) [( ) ]x s y z� 2 2 2 instead of r 2 � � ( )x y z2 2 2 . Taking into account that the term, pro- portional to s in the change of the electrostatic energy is zero, we find that the third term (the only one, which con- tributes) yields the following change in the electrostatic energy of a fullerene molecule due to the shift regarded: U s eNs / R Re ie( ) ( ) ( )� �2 3 32 [5]. This energy, U s( ), is a po- tential of a 3-d harmonic oscillator [13]. For a 3-d har- monic oscillator with mass m NU s mN s /s( ) � 2 2 2 [5,13]. As follows from the last two equations, the frequency, s , is equal to Mie frequency, M . The distance between the quantum levels is � M . Its calculated value is 21.45 eV. Its experimental value is about 20 eV. Quantum excitation to the first excited level manifests itself as 20 eV peak, to the second one as 40 eV peak. Conception of some vol- ume plasmon of unknown origin [10] is not needed. So, the frequency of the dipole oscillation in a ful- lerene molecule is as follows: f e iee N/m R R� �[ ( )] /2 3 3 1 2 . (9) Mie oscillation is also a 3-d (quantum) harmonic oscil- lator. So the results are the same for both types of oscilla- tions. Only Mie oscillation in a fullerene molecule looks less plausible. The frequency of the simple dipole oscillation is equal to Mie frequency in spherically symmetric objects [1]. So it is rather difficult to determine which type of the oscilla- tion occurs in the object under investigation. Semiclassical approach to the description of the basic properties of nanoobjects Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1065 5. Calculation of the electrostatic energy of a carbon nanotube Let us consider a long (comparative to its diameter) carbon nanotube. Let us denote N n the number of delo- calized electrons in a carbon nanotube per unit length. We denote the radius of the cylinder, on which the ions of a carbon nanotube are situated, Rn . The Gauss theorem [7] allows calculating the electro- static field inside the electronic shell of a long carbon nanotube (at R r Ri � � ): E r eN /r r R / R Re n i i( ) ( )( ) ( )� � �2 2 2 2 2 . (10) The electrostatic field, produced by the ions at R r Rn � � can be calculated also [7]: E r eN /ri n( ) � �2 . (11) At r Rn� the electrostatic field, produced by the ions, is zero. The electrostatic potential difference between the internal and external surfaces of the electronic shell in a carbon nanotube should be zero in a state of equilibrium. Otherwise the delocalized electrons will move to the loca- tions, where their potential energy is smaller. The electro- static potential difference is an integral of minus electro- static field between the surfaces mentioned. It is zero when R R /Ri n� 2 0 5exp . . (12) Eq. (12) yields for the thickness of the electronic shell, � � �R Ri , the following result: � � �( ) (exp . )/R R Ri n i 1 2 0 25 . (13) As � � �R Ri , it follows from Eq. (13) that R R R Ri i n� �� ( ) exp ./1 2 0 25 . (14) The total electrostatic energy of a system is equal to the integral of the square of the electrostatic field, divided by 8�. The value of it per unit length of a nanotube is de- scribed by the following expression: U e N / R R R R/R R R /Rn i n i n i� � [ ( ) ][ ln ( ) ln ( )]2 2 2 2 2 4 4 � � �[ ( )]( . . )e N / R R R R Rn i n i 2 2 2 2 2 2 20 75 0 75 (15) Here [see Eq. (14)], R R Ri n /� ( ) exp .1 2 0 25. From this fol- lows that the electrostatic energy of the system is a func- tion of only one variational parameter, Ri . 6. Calculation of the kinetic and total energies of a carbon nanotube The kinetic energy per unit length, T , of a gas of delocalized electrons, confined in a cylindrical layer of a volume per unit length of a nanotube, � ( )R Ri 2 2� , and calculated in accord with the ideas of the Thomas–Fermi model, is as follows [2,3]: T /m N /R R Rn i n i� �1 338 0 52 5 3 2 3 2 3. ( ) [ (exp . ) ]/ / / � . (16) Eq. (14) was used here to express the kinetic energy as a function of Rn and one variational parameter, Ri . Thence, the total energy, W T U� , is a function of a single varia- tional parameter, Ri . In a state of equilibrium the total energy reaches a minimum value Wmin � 71.51 keV/nm at Rn = 0.7 nm and N n = 670 nm–1. Parameter Rn = 0.7 nm was taken from [4]. It was assumed here also that the total number of delocalized electrons in a carbon nanotube is four times larger than the number of carbon atoms. For the same values of Rn = 0.7 nm and N n = 670 nm–1 the mini- mum (equilibrium) value of total energy is reached when Ri = Rnie = 0.577 nm. At that R Rne� � 0.816 nm and � e � = 0.239 nm. It is worthwhile noting that the total energy of a carbon nanotube also contains large negative terms, correspond- ing to the energy of the ionic cores and the energy of the chemical bonds. These terms were not calculated here. 7. Collective dipole oscillations of delocalized electrons in a carbon nanotube Let us consider first a longitudinal surface dipole col- lective oscillation of delocalized electrons in a long carbon nanotube [2]. The electric current can exist only in the layer occupied by collective electrons. In this layer the cur- rent carriers can move in an oscillatory fashion. We assume that the collective electrons move as a whole without de- formation along the axis of a carbon nanotube relative to the ions. We shall call such an oscillation longitudinal. Let us shift the current carriers in the layer by a distance h along the axis of a cylinder. As a result of the shift we have opposite charges on the bases of a cylinder, q eN hn� � . The two bases are situated rather far away from each other, so they should be treated as two separate charged disks. A disk is a limiting case of a flattened ellipsoid of revolution with half-axes a b c� � . So the capacitance of a disk is that of an ellipsoid with a � 0 and b c Rne� � , that is the capa- citance, C R /ne� 2 � [14]. The electrostatic energy of the two disks is U q /Cd � � �2 0 5 2. ( / )� e N R hn ne 2 2 22 . The force, acting on all the electrons as a whole, is a derivative � �dU /dsd �� e N /R hn ne 2 2( ) . Thus we come to the Newton equation for the lN n collective electrons as a whole: mlN d h dr e N /R hn n ne( / ) ( )2 2 2 2 0 �� . (17) Assuming h h tl� 0 sin , and using Eq. (17), we come immediately to the expression: �l n nee N /mlR� ( ) /1 2 . (18) 1066 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 Yuri Kornyushin The angular frequency of the longitudinal oscillation is inversely proportional to the square root of the length of the carbon nanotube [2]. At N n = 670 nm–1, l = 30 nm, Rne = 0.816 nm we have � l = 3.065 eV. Surface dipole collective oscillation of delocalized electrons in a conductive sphere is called, as was men- tioned above, Mie oscillation [11,12]. The frequency of Mie oscillation is square root from 3 times smaller than Langmuir frequency [2], �L e n/m� ( ) /4 2 1 2 (here n is the number of the delocalized electrons per unit volume of the electronic shell). In general case [12] the frequency of Mie oscillation, �M De n/m� ( ) /4 2 1 2 (here D is the depo- larization factor [12]). Since for a sphere [7] D = 1/3, we have for the frequency of Mie oscillation in a sphere �M e n/ m� ( ) /4 32 1 2, as was mentioned above. For a cyl- inder [7] D = 1/2 and n N / R Rn ne nie� ��( )2 2 . Thence we have for the frequency of the transverse surface (Mie) os- cillation in carbon nanotube following result: M n ne niee N /m R R� �[ ( )]2 2 2 2 . (19) For N n = 670 nm–1, Rne = 0.816 nm, and Rnie = 0.577 nm we have � M = 21 eV. Let us consider now a simple dipole oscillation of the electronic shell as a whole relative to the ion skeleton per- pendicular to the carbon nanotube axis. Regarded system consists of two cylindrically symmetric objects, the ion skeleton and the electronic shell. Hence this oscillation causes no electric field in the blank interior of a carbon nanotube around its center (cavity). Like for a fullerene molecule, when the electronic shell of a carbon nanotube is shifted as a whole by a small distance, h, relative to the ion skeleton, the electrostatic energy of a carbon nanotube is changed byW h( ), which is a potential of a 2 � d harmonic oscillator. Let us calculate this change. The electrostatic potential is equal to the in- tegral of minus electric field [7]. Inside the electronic shell of a carbon nanotube, as follows from Eq. (10), the electrostatic potential is as follows: � e en ne nieeN r / R R� � � �[ ( )]2 2 2 � �2 2 2 2 0eN R R R / r/rn nie ne nie[ / ( )] ln( ) , (20) where r0 is some constant. The second term in Eq. (20) co- incides with the potential produced by some linear charge. This term does not contribute to the change in the electrostatic energy of a carbon nanotube, because the in- teraction between the ion skeleton and a linear charge is zero. This is so because the electrostatic potential is con- stant for r Rn� . The change in the electrostatic energy does not contain a term, proportional to h because of the radial symmetry of the problem. So the only contribution comes from the first term of Eq. (20). Let us choose the x-axis along the h vector. Then after the shift we shall have in the first term of Eq. (20) [( ) ]x h y� 2 2 instead of r x y2 2 2� ( ). Taking into account that the term, propor- tional to h in the change of the electrostatic energy is zero, we find that the first term (the only one, which cont- ributes) yields the following value of the change in the electrostatic energy of a carbon nanotube (per unit length) due to the shift regarded: U h eN h R Rn ne nie( ) ( ) / ( )� �2 2 2 [5]. This energy, U h( ), is a potential of a 2 � d harmonic oscillator [5,13]. For a 2 � d harmonic oscillator [13] (mass per unit length of a carbon nanotube mN n )U s( ) � � mN s /n s 2 2 2. As follows from the last two equations, the frequency, s , is equal to Mie frequency, M . The dis- tance between the levels is � M . Its calculated value is 21 eV. So, according to the calculations performed, the quan- tum transition to the first excited level should manifest it- self as 21 eV peak, and to the second one as 42 eV peak. Mie oscillation in a carbon nanotube is also a 2 � d har- monic oscillator. So the results are the same for both types of the oscillations. Only Mie oscillation in a carbon nanotube looks less plausible. Obtained results show that a very simple semiclassical concept of Thomas–Fermi model and Electrostatics often gives rather good agreement between experimental and theoretical results. In particular, this model often gives rather reasonable description of nanoobjects. 8. Carbon peapod As in [15] we assume that the fullerene molecules are encapsulated in a carbon nanotube, and there is an inter- action between collective electrons of a carbon nanotube and fullerene molecules. Let us consider the longitudinal oscillations [4]. We re- mind that Rne denotes equilibrium external radius of the electronic shell of a carbon nanotube and Rnie denotes equi- librium internal radius of a carbon nanotube. Let us calcu- late the part of the electrostatic energy (per unit length), which depends on the shift of the electronic shell of a fullerene molecule s and that of a carbon nanotube h. It con- sists of the energy of the charges arising on the ends of a car- bon nanotube, 2 0 5 22 2 2 2� �. ( / )q /lC e N lR hn ne� , the energy of a fullerene molecule, U s eNs / R Re ie( ) ( ) ( )� �2 3 32 , multi- plied by the number of fullerene molecules in a unit length of a carbon peapod n f , the energy of the interaction of the dipole moments of the fullerene molecules en Nsf with electric field inside a carbon nanotube E q/lCn � �( )2 � �e N /lR hn ne( ) [this interaction energy is equal to �e n NN /lR shf n ne 2( ) ], and the energy of a dipole–dipole in- teraction between the dipoles of fullerene molecules. The last term was calculated for a linear array of fullerene mole- cules with the assumption that the size of the dipoles is smaller than the distance between them. The part of the total electrostatic energy, depending on the shifts s and h is as follows: Semiclassical approach to the description of the basic properties of nanoobjects Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1067 U s h e n N / R Rf e ie( , ) {[ . ( )]� � �2 2 3 30 5 � �2 2 2 2 2 2 2( ) ( )} ( )n /l f n l s e N / lR hf f n ne� � �e n NN /lR shf n ne 2( ) , where f k k /i /i i k i k ( ) ( ) ( )� � � � � � � �1 13 1 1 2 1 1 . (21) The term �2 2 2 3 2e N n /l f n l sf f( ) ( ) in Eq. (21) repre- sents the energy of the dipole–dipole interaction in a peri- odic array of fullerene molecules, calculated in the ap- proximation which neglects the real size of a dipole. This approximation somewhat underestimates the effect of dipole–dipole interaction. The last term in Eq. (21) can be positive or negative depending on the direction of the shifts s and h. It is posi- tive when the shifts are parallel, and negative when they are antiparallel. Let us consider now the transverse oscillations. In ac- cepted here model the shift of the electronic shell (as a whole) relative to the ion skeleton does not produce elec- trostatic field inside the fullerene molecule or carbon nanotube. This is because the ion skeleton and the elec- tronic shell do not lose their symmetry as a result of the shift. So there is no interaction between transverse oscil- lations of a fullerene molecule and carbon nanotube. This means that the two transverse oscillations are not changed in a carbon peapod; they remain the same as they were in a free fullerene molecule and a free carbon nanotube. Let us write the energy of the longitudinal oscillation U s h( , ) in the following form: U s h As Bh Gsh( , ) . .� �0 5 0 52 2 . (22) The specific values of the factors A, B, and G are given in Eq. (21). Now let us write Newton equations for all the collec- tive electrons per unit length of a carbon nanotube, and for all the collective electrons of all the fullerene mole- cules per unit length of a carbon nanotube, for the case of the longitudinal oscillations: mN h/ t Bh Gsn ( )� �2 2 0 � � , mn N s/ t As Ghf f ( )� �2 2 0 � � . (23) Two Eqs. (23) yield the following equations for the amplitudes of the oscillations h0 and s0: ( ) ,B mN h Gsn� � � 2 0 0 0 � � �Gh A mn N sf0 2 0 0( ) . (24) Eqs. (24) yield the frequencies of the oscillations (the same frequencies for both parallel and antiparallel modes): 1 2 2 2 2 2 0 5 2, . {( )� � �f a l � � � [( ) ( ) ] } ,/ � f a l f ne f ln R /lR 2 2 2 2 3 2 2 1 22 4 (25) a f f fe n N f n l /lm 2 2 2 2� ( ) . (26) Here f is the angular frequency of the plasmon in an isolated C60 molecule [Eq. (9)], l is the angular fre- quency of the longitudinal oscillations in an isolated car- bon nanotube [Eq. (18)]; ( ) / f a 2 2 1 22� is the angular frequency of the dipole mode of the plasma longitudinal oscillations of a linear periodic array of fullerene mole- cules. We find that the dipole–dipole interaction between the fullerene molecules reduces the plasmon energy for longitudinal oscillations. We also find in the accepted model that in the limit of infinite length carbon peapod, the coupling between the carbon nanotube and C60 mole- cules vanishes. Eqs. (25), (26) do not take into account spatial disper- sion (dependence on the wave vector q). In [4] Eqs. (25), (26) were generalized, taking into account spatial disper- sion for small wave vectors. Figure 2 shows the coupled longitudinal modes in a finite length carbon peapod (Fig. 1 in [4]), where the dotted line presents l . For small wave vector q, the coupling between the carbon nanotube and C60 molecules has a negligible effect. In the vicinity 1068 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 Yuri Kornyushin 18 16 14 12 10 8 6 4 2 0 1 2 3 4 5 6 q, nm–1 E n er g y, eV Fullerene branch Isolated nanotube Nanotube branch Fig. 2. Plasmon energy as a function of wavevector q along the carbon nanotube direction, for a carbon peapod of finite length (as calculated in [4]). The dotted curve is the plasmon disper- sion of an isolated carbon nanotube, from Fig. 2 of [16]. The lower (higher) solid curve is the carbon nanotube (fullerene) branch, when coupling between the carbon nanotube and C60 molecules is taken into account. Here the following parameters values were used [4]: � f � 20 eV, N = 120, n f � 1.03 nm –1 , Nn � 333 nm –1 , Re � 0.538 nm, R Rf ie� � 0.353 nm, Rne � 0.885 nm, R Rn nie� � 0.7 nm, and l � 29.1 nm (n lf � 30). With these values, we have f n lf( ) � 34.4, � a � 9.07 eV, 4 3�( )n R /lRf e ne � 0.0783, and �( ) / f a 2 2 1 22� = 15.3 eV. of l f a� �( ) /2 2 1 22 , the coupling lowers the carbon nanotube branch while it raises the fullerene branch, forming an anticrossing. 9. Metallic cluster A double-jelly model (delocalized electrons jelly and ions jelly) was applied for the description of metallic clusters [1]. This model does not take into account micro- scopic electrostatic field around the ions. Here we shall take into account this microscopic field. We shall analyze here the degenerate delocalized electrons. A spherical ball models the shape of a cluster. We shall model the conditions as adiabatic ones under given pres- sure. We shall restrict our consideration here by given en- tropy and given pressure condition. In this case the ther- modynamic potential, having a minimum in the state of equilibrium, is the enthalpy (it is a function of thermody- namic variables, entropy S and pressure P: H S P( , ) [17]). Entropy is a function of a radius of a cluster R. So R could be taken as thermodynamic variable instead of entropy. Under different condition the contribution of the entropy term to the Gibbs free energy is negligible. We shall consider here only atmospheric pressure, which can be neglected. So really, we will take into ac- count the energy of a cluster only. 9.1. Electrostatic energy of a separate ion Let us consider a ball of a metallic cluster of a volume V R /� 4 33� (here R is the radius of a cluster), consisting of Nc delocalized electrons. We consider here the ions as point charges, and the delocalized electrons like a nega- tively charged gas. In metallic clusters the delocalized electrons density is so high, that they are degenerate over all the temperature range of the existence of a metallic cluster. Degenerate delocalized electrons screen a long- range electrostatic field of point charges. The screening Thomas–Fermi radius is as follows [18]: 1 6 0 47142 1 2 3 2 1 2/g VE / N e R E /N /eF c F c� � �( ) . ( )/ / /� , (27) where E N /V / mF c� ( ) ( )/3 22 2 3 2� � is the Fermi energy [18]. Fermi energy is proportional to 1 2/R . Thomas–Fermi radius is proportional to R. The electrostatic field around a separate positive ion, submerged into the gas of degenerate delocalized elec- trons, is as follows [18]: � � � �( ) exp( )ze/r gr , (28) where z is the number of delocalized electrons per one atom, an r is the distance from the center of the ion. The electrostatic energy of this field is the integral over the ball volume of its gradient in a second power, di- vided by 8�. The lower limit of the integral on r should be taken as r0 (the radius of the ion), a very small value. Oth- erwise the integral diverges. Calculation yields the fol- lowing expression for the electrostatic energy of a sepa- rate ion: U z e r g gr0 2 2 0 1 00 5 0 5 2� � � . ( . ) exp . (29) For 2 0gr essentially smaller than unity Eq. (29) yields: U z e / r z e g0 2 2 0 2 22 0 75� �( ) . . (30) The first term in the right-hand part of Eq. (30), z e / r2 2 02 , represents the electrostatic energy of the bare ion. It is worthwhile to note that the expansion of a cluster leads to the decrease in the delocalized electron density. From this follows the increase in the screening radius [see Eq. (27)]. The electrostatic energy of a separate ion in- creases concomitantly. One can see it, regarding Eq. (30). 9.2. Total energy of a cluster We regard the ions of a cluster as randomly distrib- uted. It is well known since 1967, that the electrostatic energy of N /zc randomly distributed ions is just U N /z Uc� ( ) 0 [19]. The kinetic energy, T , of a gas of delocalized electrons, confined in a volume ( )4 3 3/ R� , and calculated in accord with the ideas of the Thomas–Fermi model, is as follows [1]: T N /mRc�1105 5 3 2 2. ( ) / � . (31) So the total energy of a cluster, W R( ), is as follows: W R N /mR ze N / r ze N gc c c( ) . ( ) ( ) . . / � �1105 2 0 75 5 3 2 2 2 0 2 � (32) As a function of R,W R( ) has a minimum. So, the cluster is stable. The minimum value of W R( ) is W ze N / r z N me /e c c� �( ) . ( )/2 0 4 3 4 22 0 565 � . (33) This minimum occurs when R Re� : R N /z /mee c� 2 422 2 1 3 2 2. ( ) ( )/ � . (34) The equilibrium volume per one ion is v zR / N /me /ze e c� �4 3 59 532 3 2 2 3� . ( )� . (35) For z ve� � � �1 8 821 10 24. cm3. It is about 3 times smaller than expected. For this value of ve we have average interatomic distance 2.564�10–8 cm and 2 /g � � � �5 27 10 9. cm, which is 2.43 times smaller than average interatomic distance. This means that electrostatic fields of adjacent ions [Eq. (28)] do not overlap. So the electro- static energy of N /zc ions is U N /z Uc� ( ) 0 anyway, are they distributed randomly or not. Overlapping takes place Semiclassical approach to the description of the basic properties of nanoobjects Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1069 when v /me /ze � 0 2871 2 2 3. ( )� only. According to Eq. (35) this is not so at ambient pressure. 9.3. Bulk modulus of a cluster When condensed matter subject is expanded, the in- crease in its elastic energy is as follows [20]: � � �W K V /V KV R/Re e e� �0 5 4 52 2. ( ) . ( ) , (36) where K is the bulk modulus, Ve is the initial equilibrium volume, and Re is the initial equilibrium radius. Here a well-known relation, � �V/V R/Re e� 3 , was used. It is valid when �R is essentially smaller than Re . The change in the total energy of a cluster [Eq. (32)] is � � �W W/ R dRR Re � �0 5 2 2 2. ( ) ( ) . (37) Equations (36), (37) yield: K / R W/ Re R Re � �( ) )( )1 12 2 2� � � . (38) Equations (32), 34), and (38) yield the following ex- pression for the bulk modulus: K z m e //� 0 00105 10 3 4 10 8. ( )� . (39) For z � 1 Eq. (13) yields K � �3 102 1011. erg/cm3. This value is a very reasonable one indeed. 9.4. Collective oscillations Frequency of the dipole oscillation of delocalized electrons, d , in any conductor of a spherical shape is 3 times smaller than Langmuir plasma frequency [12]. For spherical cluster regarded here we have d z me /� 0 2653 1 2 4 3. / � . (40) For z = 1 we have � d = 7.219 eV. This value is close to the classical surface-plasmon frequency in the case of Na [21]. For Na instead of the factor 0.2653 in Eq. (40) there is a factor 0.2497 [21]. So for Na we have � d = = 6.795 eV. This value is only 5.87% less then Eq. (40) yields. Nobody could expect that some simple approach like that used in this paper could compete with modern techniques of Theoretical Physics used in [21,22]. Semi- classical approach works only for simple basic problems. So far it works fine for the problems discussed here. But it cannot describe like in [22] a wide range of photon ener- gies. It was shown in [22] that for lower energies collec- tive effects are essential and for higher energies (higher than 15 eV) a single-particle picture is relevant. We are not discussing here the quadrupole and breath- ing modes, but it is worthwhile mentioning that both of them have an atomic branch of slow oscillations and an electronic one of fast oscillations [1]. 10. Discussion Today many authors use the model of the electronic shell accepted here. This model allows understanding ba- sic experimental results. As a matter of fact this concept is present in papers [6,10]. Two peaks of dipole oscillations of the delocalized electrons in C60 ions, observed experi- mentally by S.W.J. Scully et al., were reported in [10]. The authors ascribe these peaks to the surface (Mie) oscil- lation [11,12] and some bulk dipole mode. The energy values of the quanta of the two peaks were measured to be about 20 and 40 eV. Using experimental data on N , R and Ri and equation for the Mie frequency [2,10], one could calculate the energy of the quantum of the Mie oscilla- tion, 21.45 eV. Simple dipole oscillation of the electronic shell as a whole has the same frequency as Mie oscillation. Oscilla- tion of the electronic shell as a whole relative to the skele- ton of the ions is an oscillation of a 3 � d (for a fullerene molecule) and 2-d (for a carbon nanotube) harmonic os- cillator. The same is for Mie oscillation. As the frequen- cies of the two types of oscillation are the same, cited ex- perimental data could be explained both ways. Concerning carbon nanotubes more sophisticated ob- jects are now under investigation. In [23] collective exci- tations in a linear periodic array of parallel cylindrical carbon nanotubes, consisting of coaxial cylindrical tu- bules, was studied. Multiwall carbon nanostructures were studies in [24]. In [25] dynamic screening of an atom, confined within a finite width fullerene molecule is a subject of research. The authors used a classical dielectric approach. Now let us consider a problem of possibility of the atom confine- ment inside a fullerene molecule. There is no classical reason to find neutral atom confined in the cavity of a fullerene molecule or a neutral carbon nanotube. There is no electric field there and electrostatic potential there is zero like at the infinite distance from a fullerene molecule or a carbon nanotube. It was discussed in Sections 4 and 7 of the present paper. In this problem we have two objects, a fullerene mole- cule and an atom. Let us denote the ionization energy of one object as I and the affinity of the other as �A. Let us transfer one electron with charge e from one object to an- other one. The extra charge �e is situated on the external surface of the electronic shell of a fullerene molecule, r Re� . Atom (ion) is charged with the extra charge of the same value and opposite sign. The electrostatic potential on the external surface of the electronic shell of a ful- lerene molecule is equal to that in the cavity. It was dis- cussed in Sec. 2 of this paper. Both of them are equal to � e/Re . The electrostatic energy of the interaction of the two charges, one in the cavity and the other one on the surface of the electronic shell, is �e /Re 2 . The total energy is changed by I A e /Re� � 2 as a result of a charge transfer. 1070 Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 Yuri Kornyushin When this quantity is negative and minimal, the confine- ment discussed is possible. Either positive or negative ion could be found inside a fullerene molecule. It depends for which ion the total change in the energy is minimal. Using Eq. (3) it is possible to see that in case of a posi- tive ion its energy in the cavity around the center of a fullerene molecule is smaller than that inside the elec- tronic shell. So a positive ion occupies the cavity of a fullerene molecule. Positive ion is not localized in the center of a fullerene molecule, like it was assumed in [25]. It is expected to be found in any place in the cavity. It is opposite for a negative ion. As follows from Eq. (3), the minimal value of its energy is achieved when a negative ion is situated on the surface, where the positive ions of a fullerene molecule are located, r R f� . This en- ergy is smaller than I A e /Re� � 2 , because the positive ion skelton attracts a negative ion. This attraction is not taken into account in the expression I A e /Re� � 2 . In [26] a model of non-homogeneous oscillation of the shell of a fullerene molecule and a carbon nanotube was proposed. This model also can explain experimental re- sults, obtained by S.W.J. Scully et al. [10]. But this model does not look inherent to the topic. That is why it was not discussed here. We modeled the shape of metallic cluster here as a spherical ball. The delocalized electrons in metallic clus- ter are degenerate. Their kinetic energy is calculated in the spirit of Thomas–Fermi model. The delocalized elec- trons screen the electrostatic field of the ions. This field was calculated; it was shown that it contributes essen- tially to the energy of a cluster and its stability. The equilibrium values of the energy of a cluster and its volume were calculated. The bulk modulus of a cluster and simple collective dipole oscillation of the electronic shell as a whole relative to the ion skeleton were calcu- lated also. Obtained results show that a very simple semiclassical concept, based on Electrostatics and Thomas–Fermi mo- del often gives rather good agreement between experi- mental and theoretical results. In particular, this model gives rather reasonable description of nanoobjects. 1. M.Ya. Amusia and Yu. Kornyushin, Contemporary Phys. 41, 219 (2000). 2. Yu. Kornyushin and M.Ya. Amusia, Contemporary Phys. 42, 159 (2001). 3. Yu. Kornyushin, Facta Universitatis. Series: Physics, Che- mistry and Technology 2, 141 (2001). 4. Yu. Kornyushin and W. Que, J. Phys.: Condens. Matter 14, 10203 (2002). 5. Yu. Kornyushin, Facta Universitatis, Series: Physics, Che- mistry and Technology 5, 11 (2007). 6. A. R�del, R. Hentges, U. Becker H. S. Chakraborty, M.E. Madjet, and J.M. Rost, Phys. Rev. Lett. 89, 125503 (2002). 7. L.D. Landau and E.M. Lifshitz, Electrodynamics of Continu- ous Medium, Pergamon Press, Oxford–New York (1997). 8. G.F. Bertsch , A. Bulgac, D. Tom�nek, and Y. Wang, Phys. Rev. Lett. 67, 2690 (1991). 9. I.V. Hertel, H. Steger, J. de Vries, B. Weisser, C. Menzel, B. Kamke, and W. Kamke, Phys. Rev. Lett. 68, 784 (1992). 10. S.W.J. Scully, E.D. Emmons, M.F. Gharaibeh, R.A. Phaneuf, A.L.D. Kilcoyne, A.S. Schlachter, S. Schippers, A. M�ller, H.S. Chakraborty, M.E. Madjet, and J.M. Rost, Phys. Rev. Lett. 94, 065503 (2005). 11. G. Mie, Ann. Phys. 25, 377 (1908). 12. Yu. Kornyushin, Facta Universitatis. Series: Physics, Chem- istry and Technology 3, 35 (2004). 13. L.D. Landau and E.M. Lifshitz, Quantum Mechanics, Per- gamon Press, Oxford–New York (1997). 14. A. Sommerfeld, Electrodynamics, Academic Press, New York (1952). 15. C. Yannouleas, E.N. Bogachek, and U. Landman, Phys. Rev. B50, 7977 (1994). 16. W. Que, J. Phys.: Condens. Matter 14, 5239 (2002). 17. L.D. Landau and E.M. Lifshitz, Statistical Physics, Pergamon Press, Oxford–New York (1997), part 1. 18. C. Kittel, Quantum Theory of Solids, John Wiley & Sons, New York (1963). 19. Yu.V. Kornyushin, Ukr. Phys. J. 12, 1327 (1967). 20. L.D. Landau and E.M. Lifshitz, Theory of Elasticity, Per- gamon Press, Oxford–New York (1997). 21. W. Ekardt, Phys. Rev. Lett. 52, 1925 (1984). 22. O. Frank and J.M. Rost, Phys. Rev. A60, 392 (1999). 23. G. Gumbs and G.R. Aizin, Phys. Rev. B65, 196407 (2002). 24. T. St�ckli et al., Philos. Mag. B79, 1531 (1999). 25. S. Lo, A.V. Korol, and A.V. Solov’yov, J. Phys. B40, 973 (2007). 26. Yu. Kornyushin, Metallofiz. Noveishie Tekhnol. 29, 1313 (2007). Semiclassical approach to the description of the basic properties of nanoobjects Fizika Nizkikh Temperatur, 2008, v. 34, No. 10 1071