Paper
24 March 2015 Prediction of stability for carbon nanotori
O. E. Glukhova, A. S. Kolesnikova, M. M. Slepchenkov, G. V. Savostyanov
Author Affiliations +
Abstract
At the present time carbon nanostructures are the main functional material for the development of electronic devices of broad applications. One of the most perspective for practical application forms of carbon nanostructures are nanotori. This paper presents results of prediction the stability of toroidal structures using computer modeling methods. The stability of carbon nanotori are determined by scanning of the local stresses map. Obtained results showed that the highest stability is characterized for carbon nanotori (13, 0) formed by the folding of a zigzag carbon nanotube. These nanotori are characterized by the lowest enthalpy of structure formation (101 eV).

1.

INTRODUCTION

The variety of structural forms and a wide range of physical and chemical properties have caused development of new scientific branch: electronics on carbon nanostructures. Torus-like structure, called carbon nanotori attract special attention among the family of carbon nanostructured materials [1]. For the first time theoretical model of a new class of carbon nanostructures with a torus-like form was proposed in Ref [2].

Geometric characteristics and topology of these forms of carbon were determined unsaturated bonds of carbon nanotubes. The broad prospects for potential applications of nanotori in various scientific and technological fields are caused Interest in the study of nanotori. The use of carbon nanotori as new forms of hydrogen storage is one of these applications. This problem is relevant, because hydrogen is a source of clean energy will play an important role in reducing the global use of hydrocarbon. The results of a theoretical investigation of the adsorption of hydrogen on the carbon nanotorus C120 presented in [3]. Hydrogen was placed inside a nanotorus in two different ways: parallel and perpendicular to the central axis of nanotorus. It was established in consequence of the optimization of nanotorus that all adsorbed hydrogen molecules are evenly distributed along a path that is equidistant path for all points near the inner surface of the nanotorus. The results of calculations of the total energy and the binding energy of the investigated nanotorus C120 showed that hydrogen is advantageous for accumulate carbon nanostructures. The total energy of the system decreases with increasing number of adsorbed hydrogen molecules.

Carbon nanotori have a number of notable magnetic properties due to the peculiarities of the atomic structure. In particular, in 1997 Hadd theoretically predicted that nanotorus C576 has a very large value of the diamagnetic susceptibility, which can be about 130 times higher than that of benzene molecules [4]. Subsequently, a significant paramagnetic moment was found in metal torus-like carbon nanotubes [5]. The existences of a torus-like carbon nanotubes ferromagnetic moment at low temperatures were predicted by the tight-binding method [6]. This effect was caused by the presence of pentagons and heptagons in the structure of the nano torus. Aharonov-Bohm effect is another important magnetic phenomena found in torus-like carbon nanotubes[7].

Evaluation of mechanical properties of nanotori was not carried out despite a number of studies on the stability of the carbon nanotori [8-10]. It is necessary that the material nanotori met all performance requirements for its successful application in the design of electronic nanodevices. The aim of this work is a theoretical study of the strength characteristics of carbon by the tight-binding method and classical molecular dynamics method.

2.

THE MATHEMATICAL APPARATUS FOR THE STUDY OF CARBON NANOTORI

2.1

QUANTUM-CHEMICAL TIGHT-BINDING METHOD

Modification of the tight-binding method, developed in [11] was used in this study. The total energy of the system within the used model is the sum of two terms:

00019_psisdg9339_93390X_page_2_1.jpg

where Erep is repulsive energy, Ebond is filled electronic energy levels. Geometrical and energy parameters of the carbon nano-torus are determined by minimizing its total energy from the characteristic linear parameters of the framework.

Repulsive energy, which takes into account electron-electron and internuclear interaction, is the sum of the repulsive pair potentials

00019_psisdg9339_93390X_page_2_2.jpg

where i, j are number of interacting atoms; ri, rj are Cartesian coordinates. Function Vrep defined by the expression

00019_psisdg9339_93390X_page_2_3.jpg

where 00019_psisdg9339_93390X_page_2_10.jpg = 10.92 eV. Energy filled levels determined by the formula

00019_psisdg9339_93390X_page_2_4.jpg

where εn is energy electron state with the number n (eigenvalues of the Hamiltonian). «2» takes into account the electron spin. Interatomic matrix elements of the Hamiltonian are calculated by the following formula

00019_psisdg9339_93390X_page_2_5.jpg

where r is the distance between the atoms; i,j are orbital momenta of the wave functions; a is an index indicating a type of bond (σ and π). The values of the atomic terms εs, εp and equilibrium overlap integrals 00019_psisdg9339_93390X_page_2_6.jpg, 00019_psisdg9339_93390X_page_2_7.jpg, 00019_psisdg9339_93390X_page_2_8.jpg, 00019_psisdg9339_93390X_page_2_9.jpg are shown in Table 1 [11].

Table 1.

Parameter values for tight-inding method

εs, eVεp, eV, eV, eV, eV, eVp1p2,nmp3,nmp4p5,eVp6
-10.932-5.991-4.3443.9695.457-1.9382.7690.2320.1542210.924.455

2.2.

MOLECULAR DYNAMICS METHOD

Dynamic compression of carbon nano tori we investigated by molecular dynamics method. Newton’s equations of motion for a system of interacting particles are solved within of this method.

00019_psisdg9339_93390X_page_3_1.jpg

where ri is radius vector i atom, mi is atom mass, Fi is the total force acting on the i-th atom from the other particles:

00019_psisdg9339_93390X_page_3_2.jpg

where U (r) is the potential energy, which depends on the mutual arrangement of all the atoms; n - number of atoms.

Solution of the equations of motion was carried by a predictor-corrector scheme of the third order [12]. Modeling the compression process was carried out in steps of 1 fs. The forces acting on the atoms of the studied system, were calculated by quantum-chemical tight-binding method.

3.

TOPOLOGICAL MODEL NANOTORI

Several models of nanotori corresponding equilibrium configuration are discussed in this paper. The geometry of topological models is described by the point group of symmetry D5d.

Nanotorus C120

C120 is the smallest torus of the group. Basic segment of this torus is presented in Fig. 1a. Atomic coordinates of framework of nanotorus were generated by the symmetry operations of C5, which acted on the atomic coordinates of the torus segment as a result of the four transformations over atoms segment. Atomic coordinates of segment are given in the initial approximation. Atomic framework of molecule C120 are shown in Fig. 1b. Figure 1 shows that ten pentagons (light gray shading in Fig. 1b) form the outer circle, and ten heptagons (dark gray shading) form the inner circle.

Fig. 1.

Nanotorus C120: a) reference framework segment; b) the atomic framework.

00019_psisdg9339_93390X_page_4_1.jpg

Atomic coordinates nanotorus corresponding to ground state are determined by minimizing the total energy of the coordinates of the atoms within the tight-binding method. Various lengths of bonds are given numbers on the Figure 1: 1 – 1,38 Å, 2 – 1,46 Å, 3 – 1,46 Å, 4 – 1,47 Å, 5 – 1,46 Å, 6 – 1,47 Å, 7 – 1,46 Å, 8 – 1,43 Å, 9 – 1,46 Å, 10 – 1,40 Å, 12 – 1,46 Å. It can be seen, the shortest bond length is in the heptagon on the inner circle, where the framework structure in greatest deformation. The radius of the inner circle is equal to 2,03 Å. The largest bond lengths are located in the Pentagon of outside circle whose radius is equal to 5,86 Å.

Nanotorus C340

C340 is the next class of symmetry D5d of nanotori. Its difference consists not only in increasing the radii of the inner and outer circles, but also in positioning heptagons and pentagons. Basic segment of nanotorus C340 is shown in Fig. 2a, and the atomic cell generally shown in Fig. 2b. The difference framework of nanotori C120 and C340 is clearly manifested in the process of formation of the cells of the inner circle heptagon. If heptagon of nanotorus C120 are positioned adjacent to each, the heptagon of nanotorus C340 are separated by one hexagon. Naturally, the tension of the inner circle in this case framework is markedly reduced.

Fig. 2.

Nanotorus C340: a) reference framework segment; b) the atomic framework.

00019_psisdg9339_93390X_page_5_1.jpg

Technique generate coordinate framework remains the same (the coordinate transformation of the base segment) (Fig. 2a). Let us consider the atomic structure nanotorus C340. Lengths of bond to the base segment are equal: 1 – 1,38 Å, 2 – 1,40 Å, 3 – 1,40 Å, 4 –1,40 Å, 5 – 1,40 Å, 6 – 1,40 Å, 7 – 1,41 Å, 8– 1,44 Å, 9– 1,41 Å, 10 – 1,40 Å, 11 – 1,44 Å, 12 – 1,46 Å, 13 – 1,43 Å, 14 – 1,43 Å, 15 – 1,42 Å, 16 – 1,47 Å, 17 – 1,45 Å. Numbers of the various of bond lengths are given numbers in Fig. 2a. The shortest bond length of the torus C340 is a heptagon on the inner circle as in the nanotorus C120. The radius of the inner circle is equal to 6,09 Å. The largest bond lengths are located in the pentagon outside circle whose radius is equal to 10,68 Å.

Nanotorus C460

Feature of atomic carbon framework nanotorus C460 consists in that the radius of the inner circle 4.97 nm and inside part of a framework formed by heptagon. Heptagon is separated hexagons (as in nanotora C340). The radius of the outer circle of the nanotorus C460 is much higher than on the radius at the nanotorus C340: 14,88 Å. This increase in size of the outer circle of torus is provided significant removal of from each other pentagons. Basic segment of nanotorus C460 are shown in Fig. 3a. Structure of the entire framework is shown in Fig. 3b. The number of different of bond lengths also was increased by increasing the number of atoms of the base segment. We present just a few of them: lowest bond length 1,40 Å is observed in hexagons, forming the surface of the torus, the highest 1,55 Å observed in the Pentagon of outer circle.

Fig. 3.

Nanotorus C460: a) reference framework segment; b) the atomic framework.

00019_psisdg9339_93390X_page_5_2.jpg

We calculated the reaction enthalpy of formation of torus-like structures for a quantitative estimate of stability formed torus-like structures. The calculation results are in Table 2. It can be seen from the Table 2 that the most stable structure is nanotorus C340. Enthalpy and energy per atom in these nanotori are practically identical with the same parameters C60.

Table 2.

Some energy characteristics of nano tori

The number of atoms in the structureEnthalpy, kcal • mol/atom
12013,38
34010,22
46034,63

4.

STRENGTH CHARACTERISTICS OF CARBON NANO TORI

Strength characteristics of carbon nanotori were studied based on the behavior of the objects under deformation in real time. We have studied the process of axial compression. Deformation investigated nanotori was carried out at a speed of 20 m / s along the axis Z (Fig. 4). The pressure in the work was seen as applied load. Pressure was created by the graphene plate approaching to torus with a velocity of 20 m / sec. Carbon nanotori during the study were exposed to longitudinal compression by 1-5%. The value of the compressive force applied to the objects is fixed at each stage of compression. The calculation results for the case of deformation of nanotorus C120, C340, C460 were shown in Fig. 5-7.

Fig.4.

The compression scheme of the carbon nanotorus.

00019_psisdg9339_93390X_page_6_1.jpg

Fig.5.

Dependence deforming force on the magnitude of axial compression nanotorus C120

00019_psisdg9339_93390X_page_7_1.jpg

Fig. 6.

Dependence deforming force on the magnitude of axial compression nanotorus C340

00019_psisdg9339_93390X_page_7_2.jpg

Fig. 7.

Dependence deforming force on the magnitude of axial compression nanotorus C460

00019_psisdg9339_93390X_page_8_1.jpg

It is seen from Figure 8 that with compression C120 deforming force acting on the torus increases monotonically, approaching the saturation near values 120 nN. A similar pattern is observed during deformation carbon nanotori C340 and C460. Deforming force with increasing compression nanotorus C340 increases under the linear law. The small drop in value of deforming force is observed when compression of the structure 3%. Deforming force at a higher compression is almost unchanged. Deforming force of nanotorus C460, is growing rapidly in compression of object by 1-2%, and during subsequent deformation (3.5%) changes insignificant in the range of (160-180 nN).

Young’s modulus of compression was calculated for quantitative assessment strength properties of the investigated nanotori. Algorithm for computing pseudo Young’s modulus is as follows.

  • 1) At the initial stage, we optimize the structure by minimizing the total energy of the structure of with respect to coordinates of the atoms by quantum-chemical tight-binding method.

  • 2) Optimized structure was shrinking by 1%. After that, the ends of structure are firmly fixed and the structure was optimized again.

  • 3) Young’s modulus of compression of carbon nanotorus calculated by the formula

00019_psisdg9339_93390X_page_8_2.jpg

where S is the cross section nanotorus; the force required to compress the torus defined by the formula

00019_psisdg9339_93390X_page_9_1.jpg

where ΔE is the energy of elastic tension, ΔL is reducing the value of the structure.

Young’s modulus of compression for all studied in this paper nanotori was calculated by the algorithm described above. Results are presented in Table 3.

Table 3.

Calculated values of the Young’s modulus of compression for nanotori

Structure of the nanotorusYoung’s modulus of compression, TPa
C1200,54
C3400,54
C4600,65

It can be concluded from table data that nano torus C460 has the highest modulus of elasticity. The order of the calculated modulus tells us that nanotori possess high strength properties comparable with the strength properties of carbon nanotubes [13]. Consequently, this carbon material can be considered as basic material for the base element of modern electronics.

5.

CONCLUSION

We have installed the most stable carbon nanotori during the theoretical study their electron energy and strength characteristics. Nanotorus C340 has the most stable structure among the investigated objects. This fact is confirmed by the results of calculations of the enthalpy of reaction of structure formation (10.22 kcal • mol / atom) and the energy per atom (-43.02 eV). The values obtained are practically identical with the same parameters of the C60 fullerene, which is one of the most stable members of the family of carbon nanoclusters.

It is shown that the process of dynamic compression nanotori is accompanied by constant increasing deforming force. Numerical estimation of Young’s modulus in compression of carbon nanotori possible to establish that the carbon nanotorus C460 characterized by the greatest modulus of elasticity (0.65 TPa). The growth of the elastic modulus is observed with an increase in the geometric dimensions of the nanotorus.

The results obtained in this study indicate that carbon torus-like structures are characterized by stability and high strength characteristics. These structures can be used as the material of element base of modern nanoelectronics. In particular, their use as a material for the creation of cold cathodes is promising use in nanoelectronics.

ACKNOWLEDGEMENTS

The work was supported by the grant of the President of the Russian Federation for young scientists and graduate students engaged in advanced research and development in priority areas of modernization of the Russian economy (grant number SP-2302.2013.1), supported by the Ministry of Education and Science of the Russian Federation within the framework of the project of the state task in the field of scientific work (project №3.1155.2014/K) and supported by grants of RFBR (projects №14-01-31429, 14-01-31508).

REFERENCES

[1]. 

Sano, M., Kamino, A., Okamura, J., Shinkai, S., “Ring Closure of Carbon Nanotubes,” Science, 293 (5533), 1299 –1301 (2001). https://doi.org/10.1126/science.1061050 Google Scholar

[2]. 

Avron, J. E., Berger, J., “Tiling rules for toroidal molecules,” Phys. Rev. A, 51 (2), 1146 –1159 (1995). https://doi.org/10.1103/PhysRevA.51.1146 Google Scholar

[3]. 

Cruz-Torres, A., Castillo-Alvarado, F. DE L., Ortiz-Lopez, J., Arellano, J. S., “Arellano Hydrogen Storage Inside a Toroidal Carbon Nanostructure C120: Density Functional Theory Computer Simulation,” Intern. J. Quantum Chemistry, 110 2495 –2508 (2010). Google Scholar

[4]. 

Haddon R. C., “Electronic properties of carbon toroids,” Nature, 388 (6637), 31 –32 (1997). https://doi.org/10.1038/40302 Google Scholar

[5]. 

Liu, L., Guo, G. Y.,Jayanthi, C. S., Wu, S. Y., “Colossal Paramagnetic Moments in Metallic Carbon Nanotori,” Phys. Rev. Let, 88 (21), 217206 (2002). https://doi.org/10.1103/PhysRevLett.88.217206 Google Scholar

[6]. 

Rodriguez-Manzo, J. A., Lopez-Urias, F., Terrones, M., Terrones, H., “Magnetism in Corrugated Carbon Nanotori: The Importance of Symmetry, Defects, and Negative Curvature,” Nano Let, 4 (11), 2179 –2183 (2004). https://doi.org/10.1021/nl0486968 Google Scholar

[7]. 

Lin, M. F., Chuu, D. S., “Persistent currents in toroidal carbon nanotubes,” Phys. Rev. B, 57 (11), 6731 –6737 (1998). https://doi.org/10.1103/PhysRevB.57.6731 Google Scholar

[8]. 

Liu, C. P., Chen, H. B., Ding, J. W., “Magnetic response of carbon nanotori: the importance of curvature and disorder,” J. Physics: Condensed Matter, 20 (1), 015206 (2008). Google Scholar

[9]. 

Liu, C., P., Xu, N., “Magnetic response of chiral carbon nanotori: The dependence of torus radius,” Physica B: Condensed Matter, 403 (17), 2884 –2887 (2008). https://doi.org/10.1016/j.physb.2008.02.030 Google Scholar

[10]. 

Wang X., Wang Z., Liu Yq., Wang C., Bai C., Zhu D., “Ring formation and fracture of a carbon nanotube,” Chemical Physics Letters, 339 (1), 36 –40 (2001). https://doi.org/10.1016/S0009-2614(01)00291-3 Google Scholar

[11]. 

Glukhova, O. E., Terentev, O.A., “Theoretical study of the dependence of the Young’s modulus and torsion of thin single-walled carbon nanotubes zigzag and armchair on the geometrical parameters,” Solid State Physics.. V, 48 (7), 1329 –1335 (2006). Google Scholar

[12]. 

Glukhova, O. E., “Dimerization of miniature C20 and C28 fullerenes in nanoautoclave,” J. Molecular Modeling, 17 (3), 573 –576 (2011). https://doi.org/10.1007/s00894-010-0763-z Google Scholar

[13]. 

Wang, Y., Tomanek, D., Bertsh, G. F., “Stiffness of a solid composed of C60 clusters,” Phys. Rev. B, 44 (12), 6562 –5665 (1991). https://doi.org/10.1103/PhysRevB.44.6562 Google Scholar
© (2015) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
O. E. Glukhova, A. S. Kolesnikova, M. M. Slepchenkov, and G. V. Savostyanov "Prediction of stability for carbon nanotori", Proc. SPIE 9339, Reporters, Markers, Dyes, Nanoparticles, and Molecular Probes for Biomedical Applications VII, 93390X (24 March 2015); https://doi.org/10.1117/12.2080102
Lens.org Logo
CITATIONS
Cited by 2 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Carbon

Chemical species

Nanostructures

Hydrogen

Carbon nanotubes

Molecules

Motion models

RELATED CONTENT


Access provided by
Saratov State Univ.
Back to Top