Advances in Physical Chemistry
Volume 2016 (2016), Article ID 1862959, 4 pages
http://dx.doi.org/10.1155/2016/1862959

Research Article

Theoretical Studies of the Stone-Wales Defect in C₃₆ Fullerene Embedded inside Zigzag Carbon Nanotube

Konstantin S. Grishakov,1,2 Konstantin P. Katin,1,2 and Mikhail M. Maslov1,2

¹Department of Condensed Matter Physics, National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Kashirskoe Shosse 31, Moscow 115409, Russia
²Laboratory of Computational Design of Nanostructures, Nanodevices and Nanotechnologies, Research Institute for the Development of Scientific and Educational Potential of Youth, Aviatorov Street 14/55, Moscow 119620, Russia

Received 21 July 2016; Accepted 20 October 2016

Academic Editor: Sylvio Canuto

Copyright © 2016 Konstantin S. Grishakov et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We apply density functional theory at PBE/6-311G(d) level as well as nonorthogonal tight-binding model to study the Stone-Wales transformation in C₃₆ fullerene embedded inside the (14,0) zigzag carbon nanotube. We optimize geometries of two different isomers with the and the symmetries and the transition state dividing them. The mechanism of Stone-Wales transformation from to symmetry for the encapsulated C₃₆ is calculated to be the same as for the isolated one. It is found that the outer carbon wall significantly stabilizes the isomer. However, carbon nanotube reduces the activation barrier of Stone-Wales rearrangement by 0.4 eV compared with the corresponding value for the isolated C₃₆.

1. Introduction

Carbon peapods are the host-guest compounds (endohedral complexes) that consist of fullerenes encapsulated inside the single-walled carbon nanotubes. Original peapods containing only C₆₀ fullerenes were synthesized in 1998 through the catalytic pulsed laser vaporization of graphite [1]. At present peapods are of special interest as they are forming the basis for future nanoelectronic devices, since the fullerenes encapsulated have strong effect on the electronic properties of the outer nanotube [2]. Raman spectroscopy and electronic transport measurements on individual peapod confirm this significant doping impact produced by the encapsulated C₆₀s [3]. In addition, ab initio calculations clearly predict the dependence of charge transfer in the carbon peapod on the mutual orientation of the fullerenes [4] and their distribution [5].

The Stone-Wales (SW) defects [6] in the cages of encapsulated fullerenes are the subject of special interest, because they can lead to the significant changes in the electronic structure of peapod [3]. The SW defects occur during the peapods synthesis or as a result of their thermal and ultraviolet treatment [7]. This type of defects can be regarded as a rotation of single C–C bond around its center through an angle of 90° [6]. It should be noted that SW defects are characterized by the lowest activation energy barrier compared with the other possible atomic rearrangements, but they considerably influence on the physicochemical characteristics of fullerenes and other sp ²-hybridized carbon systems [8–12]. Furthermore, they play a crucial role in the processes of fullerenes isomerization [13] and coalescence [7]. For example, it is possible to transform peapod to the double-walled carbon nanotube by means of the partial or complete coalescence of the encapsulated fullerenes into the nanotube with smaller diameter under the high temperatures (~3600°C) through the SW mechanism [14, 15]. Recent experiments confirmed that such transformation resulted in the drastic reduction of the conductivity of the whole system [3]. Thus, the introduction of isolated SW defect as well as the initiation of series of SW mechanisms can be regarded as an effective method for the tuning of peapod electronic characteristics.

Here we present a theoretical study of the SW defect formation in the C₃₆ fullerene cage embedded inside the zigzag carbon nanotube. It should be noted that C₃₆ possesses a number of distinctive features compared with C₆₀. It has more strained cage and higher reactivity [16]. In addition, its smaller effective diameter allows one to introduce it inside the thinner carbon nanotubes for peapod formation. The SW activation energy barrier in the isolated C₃₆ is equal to 6.23 eV [17], which is lower than the corresponding value 7.16 eV for C₆₀ [18]. Therefore, the C₃₆ isomerization process can be initiated more easily. Various C₃₆ isomers possess different physicochemical properties; for example, they may vary in electronic characteristics (their HOMO-LUMO gaps lie in the range of 0.20–0.83 eV [19]) and in thermodynamic stability. One of the most stable isomers of C₃₆ has symmetry; thus, this fullerene is stretched along its symmetry axis. Therefore, it can be expected that the encapsulated -C₃₆s are arranged inside the carbon nanotube so that their symmetry axes are parallel to the main axis of the nanotube. Thus, the inner side of the carbon nanotube wall is able to stabilize the -C₃₆ fullerene and to prevent its isomerization. Charge transfer between the tube and the fullerene can also effect on the SW defect formation in C₃₆. Note that for the neutral and positively charged C₃₆ the most energetically favorable isomer has symmetry, whereas for the negatively charged C₃₆ the isomer is more thermodynamically stable [17]. On the other hand, it was demonstrated that both positive and negative charges changed weakly the SW activation barrier of isolated C₃₆ [17]. So, the main purpose of this work is to study the influence of the host zigzag carbon nanotube on the formation of SW defect in the encapsulated C₃₆ guest fullerene.

2. Materials and Methods

In our study we consider C₃₆ fullerenes with the and the symmetries both isolated and encapsulated inside the center of the (14,0) zigzag carbon nanotube. Previous ab initio calculations confirmed that such type of nanotubes with the encapsulated C₃₆s demonstrated the highest binding energy among all zigzag [20] and armchair [21] tubes. In our simulations carbon nanotube is constructed from 224 carbon atoms. To avoid the dangling bonds hydrogen passivation of the tube edges is performed. So, carbon nanotube considered is 11 Å in diameter with the length of 16 Å. Note that the same representation was used in [20]. Finally, our atomistic model of peapod has the chemical view C₃₆@C₂₂₄H₂₈ as it is shown in Figure 1.

(a)

(b)

Figure 1: Atomistic model of carbon peapod C₃₆@C₂₂₄H₂₈: side view (a) and front view (b).

The geometries of minima and transition states of C₃₆ and C₃₆@C₂₂₄H₂₈ molecular systems are obtained within the nonorthogonal tight-binding model (NTBM) [22]. Previously, we successfully used this model to simulate the C₆₀s coalescence in peapods [3], to determine the SW mechanism in isolated C₃₆ [23, 24] and other fullerenes [23–25], and so forth. Initial molecular structures are relaxed without any symmetry constrains until the residual forces acting on atoms are less than 10−3 eV/Å. For each configuration obtained the Hessian matrix and normal frequencies are also calculated in the framework of NTBM model to confirm the energy minima (imaginary frequencies are absent) or transition state (the frequency spectrum has only one imaginary frequency). Next, using the geometries found at NTBM level the total energies of various C₃₆ configurations (both minima and transition state) are computed by means of the density functional theory (DFT) approach via the single-point energy calculation procedure. All the quantum-chemical DFT calculations are performed in the frame of GAMESS software package [26]. The hybrid Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional [27] is used along with the 6-311G(d) basis set [28]. It was shown that this basis set was good enough to provide an adequate description of SW defect formation in isolated C₃₆ [17]. Finally, the peapod total energy is calculated according to the following formula [29]:where , , and are the total energies of the peapod sample and encapsulated C₃₆ obtained at the corresponding level of theory (NTBM or DFT, resp.). Grimme empirical dispersion corrections [30] are taken into account for all energies. GAMESS program package is used for their calculation. Note that formula (1) is applicable both to reactant/product configurations and to transition state.

3. Results and Discussion

First of all, we optimize the geometries of isolated C₃₆ isomers with the and the symmetries and the transition state (TS) dividing them (see Figure 2). The choice of these isomers is caused by their high thermodynamic stability [16] and by their interconversion from one to another through the SW mechanism [17]. We obtain that the length of “rotated” C–C bond (see Figure 2) is equal to 1.455, 1.296, and 1.459 Å for -C₃₆, TS-C₃₆, and -C₃₆ atomic configurations, respectively. These data coincide well with the corresponding values of 1.443, 1.292, and 1.445 Å, calculated in [17] at DFT/B3LYP/6-311G(d) level of theory. Moreover, NTBM model predicts that the energy difference between the -C₃₆ and -C₃₆ neutral isomers is equal to 0.14 eV [23]. This result is in good agreement with the value of 0.11 eV also reported in [17]. Note that the simplified single-point energy calculations also give the small absolute value of energy difference of 0.07 eV, which indicates that -C₃₆ and -C₃₆ states lie very close in energy. Thus, the activation energy barrier obtained for the --C₃₆ SW transformation of 6.7 eV agrees with the more accurate value of 6.3 eV [17]. So, NTBM model provides adequate geometries of the C₃₆ local minima and transition state.

Figure 2: Potential energies Δ of endohedral complex C₃₆@C₂₂₄H₂₈ (solid line) and isolated C₃₆ fullerene (dashed line) during the SW transformation --C₃₆. The insets show the C₃₆ isomers with the and the symmetries and the transition state (TS). Rotation of C–C bond during the SW transformation is indicated in black. The atoms of carbon nanotube are not displayed for clarity.

We also calculate the energy and structure of the analogous -C₃₆ and -C₃₆ isomers as well as TS-C₃₆ embedded inside the C₂₂₄H₂₈ sample nanotube. It is confirmed that there are no intermediate states along the isomerization reaction path as it was earlier established for the isolated neutral and charged C₃₆s [17]. Nevertheless, negligible changes in bond lengths of the encapsulated C₃₆ are identified compared with the isolated one. The “rotated” C–C bond length decreases by only 0.003, 0.001, and 0.002 Å for the -C₃₆, TS-C₃₆, and -C₃₆ atomic configurations, respectively. This slight reduction in bond lengths confirms that the (14,0) carbon nanotube possesses too large diameter for significant compression of the C₃₆ cage. On the other hand, we obtain the carbon nanotube effects on the comparative thermodynamic stability of two C₃₆ isomers considered. For the encapsulated C₃₆ atomic configuration, axially stretched isomer lies lower in energy than the more spherical one by 0.57 eV. Thus, the embedded -C₃₆ is more thermodynamically stable than the -C₃₆.

The activation energy barrier for the --C₃₆ SW transformation inside the carbon nanotube is found to be 6.3 eV. This value is less than the activation energy barrier for the corresponding isomerization reaction for the isolated C₃₆ by Δ = 0.4 eV. For estimation of the sensitivity of Δ to the DFT functional used, we recalculate the values of transition barriers for the --C₃₆ SW transformation (for isolated C₃₆ and embedded one) in the frame of wB97X-D functional, which contains Grimme empirical dispersion corrections by default [31, 32]. The same corrections are applied to NTBM energies as well. As a result, we obtain that the energy barrier for SW transformation in isolated fullerene is equal to 7.8 eV, and the analogous barrier for the encapsulated C₃₆ is equal to 7.5 eV. Although these values are higher than the corresponding ones obtained within the PBE functional, one can see that Δ depends weakly on the functional choice.

The results obtained demonstrate that the activation of SW mechanism in the encapsulated C₃₆ becomes ~exp(Δ) times faster compared with the isolated fullerene (here is the temperature and is the Boltzmann constant). For example, the factor exp(Δ) is equal to 3.6 at = 3600 K (this value is typical for the fullerenes coalescence [14]) and achieves ~5·10⁶ at room temperature (300 K).

4. Conclusions

In this paper we presented the study of SW defect formation in C₃₆ fullerene placed inside the (14,0) zigzag carbon nanotube. We obtain that the isomerization reaction passes through the in-plane C–C bond rotation without any intermediate steps. The presence of carbon nanotube leads to the negligible changes in the geometry of reactants, products, and transition state. It is found that the most thermodynamically stable C₃₆ isomer in the case of host-guest compound has symmetry whereas the most thermodynamically stable unrelated C₃₆ isomer possesses symmetry. The activation energy barrier for the --C₃₆ SW transformation inside the nanotube is less than the corresponding value for the isolated C₃₆ isomerization reaction by 0.4 eV. However, the carbon nanotube significantly stabilizes the -C₃₆ isomer. We believe that the data reported in this article will be useful for the further theoretical and experimental studies on the C₃₆s isomerization or their coalescence inside the carbon peapods.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The reported study was funded by RFBR according to the research Project no. 16-32-00588 mol_a.

References

B. W. Smith, M. Monthioux, and D. E. Luzzi, “Encapsulated C₆₀ in carbon nanotubes,” Nature, vol. 396, no. 6709, pp. 323–324, 1998.

M. Monthioux, “Filling single-wall carbon nanotubes,” Carbon, vol. 40, no. 10, pp. 1809–1823, 2002.

V. Prudkovskiy, M. Berd, E. Pavlenko et al., “Electronic coupling in fullerene-doped semiconducting carbon nanotubes probed by Raman spectroscopy and electronic transport,” Carbon, vol. 57, pp. 498–506, 2013.

O. Dubay and G. Kresse, “Density functional calculations for C₆₀ peapods,” Physical Review B—Condensed Matter and Materials Physics, vol. 70, no. 16, Article ID 165424, pp. 1–10, 2004.

J. Chen and J. Dong, “Electronic properties of peapods: effects of fullerene rotation and different types of tube,” Journal of Physics: Condensed Matter, vol. 16, no. 8, pp. 1401–1408, 2004.

A. J. Stone and D. J. Wales, “Theoretical studies of icosahedral C₆₀ and some related species,” Chemical Physics Letters, vol. 128, no. 5-6, pp. 501–503, 1986.

Y. Zhao, Y. Lin, and B. I. Yakobson, “Fullerene shape transformations via Stone-Wales bond rotations,” Physical Review B, vol. 68, no. 23, Article ID 233403, 2003.

J. Wei, H. Hu, H. Zeng, Z. Wang, L. Wang, and P. Peng, “Effects of nitrogen in Stone-Wales defect on the electronic transport of carbon nanotube,” Applied Physics Letters, vol. 91, no. 9, Article ID 092121, 2007.

M. Kabir, S. Mukherjee, and T. Saha-Dasgupta, “Substantial reduction of Stone-Wales activation barrier in fullerene,” Physical Review B, vol. 84, no. 20, Article ID 205404, 2011.

M. Kabir and K. J. Van Vliet, “Kinetics of topological stone-wales defect formation in single-walled carbon nanotubes,” The Journal of Physical Chemistry C, vol. 120, no. 3, pp. 1989–1993, 2016.

A. I. Podlivaev and L. A. Openov, “Dynamics of the Stone-Wales defect in graphene,” Physics of the Solid State, vol. 57, no. 4, pp. 820–824, 2015.

L. A. Openov and A. I. Podlivaev, “Interaction of the stone–wales defects in graphene,” Physics of the Solid State, vol. 57, no. 7, pp. 1477–1481, 2015.

Y. E. Lozovik and A. M. Popov, “Formation and growth of carbon nanostructures: fullerenes, nanoparticles, nanotubes and cones,” Physics-Uspekhi, vol. 40, no. 7, pp. 717–737, 1997.

E. Hernández, V. Meunier, B. W. Smith et al., “Fullerene coalescence in nanopeapods: a path to novel tubular carbon,” Nano Letters, vol. 3, no. 8, pp. 1037–1042, 2003.

F. Ding, Z. Xu, B. I. Yakobson et al., “Formation mechanism of peapod-derived double-walled carbon nanotubes,” Physical Review B—Condensed Matter and Materials Physics, vol. 82, no. 4, Article ID 041403, 2010.

K. M. Kadish and R. S. Ruoff, Fullerenes: Chemistry, Physics, and Technology, John Wiley & Sons, New York, NY, USA, 2000.

Y.-F. Jin and C. Hao, “Computational study of the Stone-Wales transformation in C36,” The Journal of Physical Chemistry A, vol. 109, no. 12, pp. 2875–2877, 2005.

W. I. Choi, G. Kim, S. Han, and J. Ihm, “Reduction of activation energy barrier of Stone-Wales transformation in endohedral metallofullerenes,” Physical Review B—Condensed Matter and Materials Physics, vol. 73, no. 11, Article ID 113406, 2006.

K. H. Kim, Y.-K. Han, and J. Jung, “Basis set effects on relative energies and HOMO-LUMO energy gaps of fullerene C₃₆,” Theoretical Chemistry Accounts, vol. 113, no. 4, pp. 233–237, 2005.

Y. Bao-Hua, W. Yang, and H. Yuan-He, “Structures and electronic properties of C₃₆ encapsulated in single-wall carbon nanotubes,” Chinese Journal of Chemistry, vol. 23, no. 4, pp. 370–376, 2005.

Y. Baohua, W. Yang, and H. Yuanhe, “Theoretical studies of C₃₆ encapsulated in zigzag single-wall carbon nanotubes,” Chinese Science Bulletin, vol. 51, no. 1, pp. 25–30, 2006.

M. M. Maslov, A. I. Podlivaev, and K. P. Katin, “Nonorthogonal tight-binding model with H–C–N–O parameterisation,” Molecular Simulation, vol. 42, no. 4, pp. 305–311, 2016.

A. I. Podlivaev, K. P. Katin, D. A. Lobanov, and L. A. Openov, “Specific features of the stone-wales transformation in the C₂₀ and C₃₆ fullerenes,” Physics of the Solid State, vol. 53, no. 1, pp. 215–220, 2011.

A. I. Podlivaev and K. P. Katin, “On the dependence of the lifetime of an atomic cluster on the intensity of its heat exchange with the environment,” JETP Letters, vol. 92, no. 1, pp. 52–56, 2010.

K. P. Katin and A. I. Podlivaev, “Dynamic characteristics of the low-temperature decomposition of the C₂₀ fullerene,” Physics of the Solid State, vol. 52, no. 2, pp. 436–438, 2010.

M. W. Schmidt, K. K. Baldridge, J. A. Boatz et al., “General atomic and molecular electronic structure system,” Journal of Computational Chemistry, vol. 14, no. 11, pp. 1347–1363, 1993.

J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 77, no. 18, pp. 3865–3868, 1996, Erratum to Physical Review Letters, vol. 78, no. 7, p. 1396, 1997

R. Krishnan, J. S. Binkley, R. Seeger, and J. A. Pople, “Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions,” The Journal of Chemical Physics, vol. 72, no. 1, pp. 650–654, 1980.

M. Svensson, S. Humbel, R. D. J. Froese, T. Matsubara, S. Sieber, and K. Morokuma, “ONIOM: A multilayered integrated MO + MM method for geometry optimizations and single point energy predictions. A test for Diels-Alder reactions and Pt(P(t-Bu)₃)₂ + H₂ oxidative addition,” The Journal of Physical Chemistry, vol. 100, no. 50, pp. 19357–19363, 1996.

S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, “A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu,” Journal of Chemical Physics, vol. 132, no. 15, Article ID 154104, pp. 1–9, 2010.

J.-D. Chai and M. Head-Gordon, “Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections,” Physical Chemistry Chemical Physics, vol. 10, no. 44, pp. 6615–6620, 2008.

R. Peverati and K. K. Baldridge, “Implementation and performance of DFT-D with respect to basis set and functional for study of dispersion interactions in nanoscale aromatic hydrocarbons,” Journal of Chemical Theory and Computations, vol. 4, no. 12, pp. 2030–2048, 2008.