Abstract
The changes of dynamical behaviour of a single fullerene molecule inside an armchair carbon nanotube caused by the structural Peierls transition in the nanotube are considered. The structures of the smallest C_{20 }and Fe@C_{20 }fullerenes are computed using the spinpolarized density functional theory. Significant changes of the barriers for motion along the nanotube axis and rotation of these fullerenes inside the (8,8) nanotube are found at the Peierls transition. It is shown that the coefficients of translational and rotational diffusions of these fullerenes inside the nanotube change by several orders of magnitude. The possibility of inverse orientational melting, i.e. with a decrease of temperature, for the systems under consideration is predicted.
Introduction
The structure and elastic properties of carbon nanotubes are studied in connection with the perspectives of their applications in nanoelectronic and nanoelectromechanical devices and composite materials, and are also of fundamental interest, particularly for physics of phase transitions. For example, superconductivity [1], commensurateincommensurate phase transition in doublewalled nanotubes [2], spontaneous symmetry breaking with formation of corrugations along nanotube axis [3] and structural Peierls transition in armchair nanotubes [49] have been considered. In the present Letter, we consider a fundamentally new phenomenon related to phase transitions in nanosystems. In other words, we consider the possibility of inverse orientational melting for molecules encapsulated inside nanotubes caused by structural Peierls transition in the nanotubes.
The possibility of Peierls transition in carbon nanotubes was first considered in [4]. As a result of this transition, armchair nanotubes become semiconducting at low temperature, and Peierls distortions lead to the Kekule structure (see Figure 1) with two essentially different CC bond lengths and a triple translational period (three times more hexagons in the translational unit cell). In previous studies, the Peierls gap [57] and the temperature of the transition to the metallic phase with equal CC bond lengths [4,5,8] were estimated. Recently, the Kekule structure was calculated for the ground state of an infinite armchair (5,5) nanotube by PM3 semiempirical molecular orbital calculations [9]. It was shown that, for the (5,5) nanotube, the difference between CC bond lengths for semiconducting phase is 0.03 Å, whereas the difference between nonequivalent CC bond lengths for metallic phase is only 0.006 Å [9].
Figure 1. Calculated Kekule structure corresponding to the ground state of the (8,8) nanotube. Nanotube axis is shown by the dashed line.
Note that density functional theory (DFT) calculations for the (5,5) nanotube of a finite length [10,11] also gave a 60 atom periodicity of physical properties on the length of nanotube segment which is consistent with Kekule structure for infinite armchair nanotubes. Moreover, Xray crystallographic analysis of chemically synthesized short (5,5) nanotubes [12] shows the Kekule bond length alternation pattern, which was in good agreement with DFT and PM3 calculations also performed in [12]. By the example of the infinite (5,5) nanotube, it was demonstrated that the structural Peierls transition connected with spontaneous symmetry breaking takes place not only with an increase of temperature, but also can be controlled by uniaxial deformation of armchair nanotubes [9].
A dynamical behaviour of molecules encapsulated inside nanotubes can correspond to the following regimes: oscillations about a fixed position and/or a fixed orientation of the molecule (regime A), hindered motion along the nanotube axis and/or rotation of the molecule (regime B) and free motion and/or rotation of the molecule (regime C). In the present Letter, we show the possibility of changes of the dynamical behaviour of molecules encapsulated inside armchair nanotubes as a result of the Peierls transition in the nanotube structure. In other words, these changes can include switching between the regimes A and B, switching between the regimes B and C, and the changes in diffusion coefficients corresponding to the hindered motion and/or rotation of the molecule (regime B). The considered changes are possible in the case where the regime B takes place for at least one phase of the nanotube, i.e. the temperature T_{P }of the Peierls transition should correspond to the temperature range of the regime B (the hindered motion and/or rotation of the molecule). In other words, the temperature T_{P }should be of the same order of magnitude (or a few orders of magnitude less) as energy barriers ΔE for motion and/or rotation of the molecule inside the nanotube at this phase. Note that inverse melting of motion and/or rotation of the molecule is possible, if the Peierls transition from the high to lowtemperature phase of the nanotube occurs with switching from the regime B to the regime C or switching from the regime A to the regime B.
Estimations showed that the Peierls transition temperature is T_{P }≃ 115 K [4,5,8]. According to calculations [13,14], the barriers of the value close to this temperature range were obtained for rotation of the fullerene C_{60 }inside the C_{60}@C_{240 }nanoparticle. It was also found that the changes of bond lengths of the fullerene C_{240}, the outer shell of these nanoparticles, within 0.06 Å lead to an increase of the barriers for rotation by more than an order of magnitude [13,14]. The changes of the nanotube bond lengths caused by the Peierls distortions are of the same order of magnitude (about 0.03 Å for the (5,5) nanotube [9]). Note also that the size of an encapsulated molecule, and therefore, the nanotube radius cannot be too large, since the magnitude of the Peierls distortions decreases with an increase of the armchair nanotube radius [6].
Thus, taking into account the above considerations, we have chosen the smallest fullerene C_{20 }and the magnetic endofullerene Fe@C_{20 }to investigate changes in the dynamical behaviour of molecules inside nanotubes at the Peierls transition. It has been shown that the (8,8) nanotube is the smallest armchair carbon nanotube which can encapsulate the fullerene C_{20 }[15]. A carbon nanotube with the fullerene C_{20 }inside was also used as a model system to simulate a drug delivery via the nanotube [16].
This Letter is organized as follows: "Fullerene and nanotube structures" section presents the DFT calculations of the structure of the C_{20 }and Fe@C_{20 }fullerenes and the PM3 calculations of the structure of the (8,8) nanotube. "Fullerenenanotube interaction" section presents the semiempirical calculations of the barriers for motion and rotation of the fullerenes inside the nanotube. The section that succeeds the latter is devoted to the dynamical behaviour of molecules inside the nanotubes. Our conclusions are summarized in the final section.
Fullerene and nanotube structures
Structures of the C_{20 }and Fe@C_{20 }fullerenes have been calculated using the spinpolarized density functional theory implemented in NWChem 4.5 code [17] with the BeckeLeeYangParr exchangecorrelation functional (B3LYP) [18,19]. Eighteen inner electrons of the iron atom are emulated with the help of the effective core potential  CRENBS ECP [20] (only 8 valence sd electrons are taken into account explicitly). The 631G* basis set is used for describing electrons of the carbon atoms.
The ground state of the fullerene C_{20 }is found to be a singlet state and has D_{2h }symmetry. The calculated energy of the triplet state of the fullerene C_{20 }is found to be 64 meV greater than the energy of the ground state. The ground state of the endofullerene Fe@C_{20 }is found to be a septet state and has C_{2h }symmetry.
The calculated structures of the ground states of the C_{20 }and Fe@C_{20 }fullerenes are shown in Figure 2. The iron atom locates in the centre of the endofullerene. The smallest and the greatest distances between the carbon atoms and the C_{20 }fullerene centre increase by 6 and 1%, respectively, as a result of the iron atom encapsulation.
Figure 2. Calculated structure of the C_{20 }(a) and Fe@C_{20 }(b) fullerenes. The bond lengths are in ångströms. The distances between the carbon atoms and the fullerene centre (in ångströms) are denoted in italics. The atoms which have smaller and greater distances to the fullerene centre are shown by the open and filled circles, respectively. C_{2 }symmetry axis is shown by the dotted line.
The semiempirical method of molecular orbitals modified for onedimensional periodic structures [21] with PM3 parameterization [22] of the Hamiltonian has been used to calculate the structure of the (8,8) nanotube. The method was used previously for calculating the Kekule structure of the (5,5) nanotube ground state and for studying structural transitions controlled by uniaxial deformation of this nanotube [9]. The adequacy of the PM3 parameterization of the Hamiltonian has been demonstrated [23] by the calculation of bond lengths of the C_{60 }fullerene with I_{h }symmetry: the calculated values of the bond lengths agree with the measured ones [24] at the level of experimental accuracy of 10^{3 }Å. The calculated Kekule structure of the (8,8) nanotube ground state is shown in Figure 1. The difference between the lengths of short and long bonds of this Kekule structure of the (8,8) nanotube is close to such a difference of the (5,5) nanotube [9]. The Peierls distortions include also radial distortions of the armchair carbon nanotube with periodicity of half of the translational period of the nanotube (for details see [9]). In the case of the (8,8) nanotube, the longest nanotube radius is 0.547 nm, while the shortest radius is 0.544 nm.
Fullerenenanotube interaction
The structures of the fullerenes C_{20 }and Fe@C_{20 }obtained in "Fullerene and nanotube structures" section have been used for finding the ground state position and for studying motion and rotation of these fullerenes inside the (8,8) nanotube. Two structures of the (8,8) nanotube have been considered: the Kekule structure calculated in "Fullerene and nanotube structures" section, and the structure with all equal bonds 1.423 Å in length (so that to be equal to the average bond length of the calculated the Kekule structure). The calculations of the interaction energy between walls of doublewalled carbon nanotubes showed that the barriers for motion of the short wall relative to the long wall are very sensitive to the length of the long wall [25]. Thus, the size of the system is too large for ab initio calculations. The analogous problem exists also for the considered case of a single fullerene inside a nanotube. Therefore, the interaction between carbon atoms of the fullerenes and the nanotube at the interatomic distance r is described by the LennardJones 126 potential
with the parameters ε = 2.755 meV, σ = 3.452 Å. These parameters of the LennardJones potential for the fullerenenanotube interaction are obtained as the average values of the parameters [26] for fullerenefullerene and fullerenegraphene interactions, in accordance with the procedure described in [26]. Here, the LennardJones potential is used for calculating the potential surface of the interaction energy E_{W }between the fullerene and the infinite nanotube, and we believe that this gives adequate qualitative characteristics of the potential surface shape. The cutoff distance, r = r_{c }of the LennardJones potential is taken equal to r_{c }= 15 Å. For this cutoff distance the errors of calculation of the interaction energy E_{W }between the fullerenes and the (8,8) nanotube and the barriers for relative motion and rotation of the fullerenes inside the nanotube are less than 0.1%. Both the fullerenes and the nanotube are considered to be rigid. An account of structure deformation is not essential for the shape of the potential surface both for the interwall interaction of carbon nanotubes [25,27] and the intershell interaction of carbon nanoparticles [13,14]. For example, the account of the structure deformation of the shells of C_{60}@C_{240 }nanoparticle gives rise to changes of the barriers for relative rotation of the shells which are less than 1% [13,14]. It should also be noted that the symmetry of interaction energy as a function of coordinates describing relative positions of interacting objects is determined unambiguously by symmetries of the isolated objects and does not change if the symmetries of the objects are broken because of their interactions.
The ground state interaction energies between the C_{20 }and Fe@C_{20 }fullerenes, and the (8,8) nanotube with Kekule structure are found to be 1.596 and 1.598 eV, respectively. The angles between the C_{2 }symmetry axes of the C_{20 }and Fe@C_{20 }fullerenes and the nanotube axis at the ground states are 49.6° and 53.1°, respectively. The metastable states with the C_{2 }symmetry axes of the fullerenes perpendicular to the nanotube axis are also found for both C_{20 }and Fe@C_{20}. At the metastable states, the interaction energies are greater by 4.58 and 3.04 meV than the ground state energies, for C_{20 }and Fe@C_{20}, respectively.
The potential surfaces of the interaction energy between the C_{20 }fullerene and the nanotube, E_{W}(φ, z), as functions of the relative displacement of the fullerene along the axis of the nanotube z and the angle of relative rotation of the fullerene about the axis of the nanotube φ are presented in Figure 3 for both the considered structures of the (8,8) nanotube. In the general case, diffusion of a fullerene along the nanotube axis is accompanied by rotation of the fullerene. Our calculations show that the barriers for rotation of both fullerenes about the axes which are perpendicular to the nanotube axis lie between 3 and 23 meV for any orientation of the fullerene and for both the considered structures of the (8,8) nanotube. These barriers are significantly greater than the barriers for rotation of the fullerenes about the axis of the nanotube (shown in Figure 3 for the C_{20 }fullerene). Therefore, diffusion of the fullerenes along the nanotube axis is accompanied only by rotation about the nanotube axis. Thus, the minimal barrier ΔE_{d }for diffusion of the fullerenes along the nanotube axis is the barrier between adjacent minima of the potential surface E_{W}(φ, z). Figure 3 shows that the shapes of the potential surface E_{W}(φ, z), corresponding to the Kekule structure of the (8,8) nanotube and to the structure of metallic phase are essentially different. For a case of the structure corresponding to the metallic phase of the (8,8) nanotube, all the barriers between adjacent minima of the potential surface E_{W}(φ, z) are equivalent. In this case, the same barrier ΔE_{d }= ΔE_{r }should be overcome for diffusion of the fullerenes along the nanotube axis and for rotation of the fullerenes about this axis (see Figure 3b). For a case of the Kekule structure, two different barriers between the adjacent minima exist: the barrier ΔE_{d }for diffusion of the fullerenes along the nanotube axis, and the barrier ΔE_{r }to rotation of the fullerenes about this axis (see Figure 3a). For the Fe@C_{20 }fullerene, the shapes of the potential surfaces E_{W}(φ, z) are qualitatively the same as for C_{20 }for both the considered structures of the nanotube. The calculated values of the barriers ΔE_{d }and ΔE_{r }are listed in Table 1 for both fullerenes. The dependences of the interaction energy between the C_{20 }and Fe@C_{20 }fullerenes and the nanotube on the relative displacement of the fullerene along the axis of the nanotube and the angle of its relative rotation about this axis, corresponding to both the considered structures of the (8,8) nanotube are compared as shown in Figure 4. Figure 4 is a vivid illustration of the significant changes of the barriers ΔE_{d }and ΔE_{r}. The most dramatic change corresponds to the rotation of the Fe@C_{20 }fullerene about the nanotube axis. Let us discuss the reason of the significant changes of the barriers ΔE_{d }and ΔE_{r }at the Peierls transition by the example of the barrier ΔE_{r }for rotation of the Fe@C_{20 }endofullerene about the nanotube axis. Figure 5 presents the dependences E_{Wi}(φ) of the interaction energies between the endofullerene and individual atoms of the nanotube on the angle φ of rotation of the fullerene. Figure 5 shows that maxima of dependences E_{Wi}(φ) for individual atoms of the nanotube occur at different angles φ_{mi }and so the dependence E_{W}(φ) of total energy on the angle of rotation is essentially smoothed. In other words, the barrier ΔE_{r }in the dependence E_{W}(φ) of the total interaction energy between the endofullerene and the nanotube is less by an order of magnitude than the barriers ΔE_{ri }in the dependences of the interaction energy between the endofullerene and only one of the nanotube atoms. Thus, the barrier ΔE_{r }is very sensitive to the values of the barriers ΔE_{ri }and angles φ_{mi}. Therefore, the barrier ΔE_{r }changes considerably at the found changes of the barriers ΔE_{ri }and angles φ_{mi }obtained from the Peierls distortions of the nanotube structure and the change of the nanotube symmetry at the transition. It should also be noted that small barriers to relative motion of nanoobjects resulting from the compensation of contributions of individual atoms to the barriers is a phenomenon well studied by the examples of such systems as doubleshell carbon nanoparticles [13,14], doublewalled carbon nanotubes [25,2732] and a graphene flake in a graphite surface [33]. The considerable changes of barriers for relative rotation of shells at small changes of shell structure were found also for doubleshell carbon nanoparticles [13,14].
Figure 3. The interaction energy E_{W }(in meV) as a function of the fullerene C_{20 }displacement z along the nanotube axis and the angle φ of the fullerene rotation about the nanotube axis. (a) The (8,8) carbon nanotube with the Kekule structure; (b) the (8,8) carbon nanotube with the structure of metallic phase. The energy is given relative to the energy minima. The equipotential lines are drawn at an interval 0.2 meV.
Table 1. Calculated characteristics of the dynamical behaviour of the C_{20 }and Fe@C_{20 }fullerenes inside the (8,8) nanotube of different structure
Figure 4. Interaction energy between the C_{20 }fullerene or Fe@C_{20 }endofullerene and (8,8) nanotube. Dependences of the interaction energy E_{W }between the C_{20 }fullerene (a, b) or Fe@C_{20 }endofullerene (c, d) and the (8,8) nanotube on the fullerene displacement z along the nanotube axis (a, c) and on the angle φ of the fullerene rotation about the nanotube axis (b, d). Solid lines denote nanotube with the Kekule structure, dashed lines denote nanotube with the structure of metallic phase. The energy minimum is positioned at E_{W }= 0, z = 0 and φ = 0.
Figure 5. Interaction energy between the Fe@C_{20 }endofullerene and (8,8) nanotube. Dependences of the interaction energy E_{Wi }between the Fe@C_{20 }endofullerene and individual atoms of the (8,8) nanotube on the angle φ of rotation of the endofullerene about the nanotube axis are denoted by the thin lines. Dependence of the total interaction energy E_{W }between the Fe@C_{20 }endofullerene and the (8,8) nanotube on the angle φ is denoted by the thick line. (a) The (8,8) carbon nanotube with the Kekule structure; (b) the (8,8) carbon nanotube with the structure of metallic phase. All energies are given relative to the energy minima. Only dependences E_{Wi }with high values of the barriers ΔE_{ri }are shown.
The frequencies of small vibrations of the fullerenes along the nanotube axis (ν_{d}), rotational vibrations about the nanotube axis (ν_{z}) and rotational vibrations about two mutually perpendicular lateral axes (ν_{x}, ν_{y}) are also calculated and listed in Table 1. The most remarkable change of frequency as a result of the structural phase transition corresponds to rotational vibrations of the Fe@C_{20 }fullerene about the nanotube axis (this agrees with the changes of the barriers).
Dynamical behaviour of molecules inside nanotube
Let us consider the possible changes of the dynamical behaviour of the C_{20 }and Fe@C_{20 }fullerenes inside the (8,8) nanotube caused by the structural phase transition. The Peierls instability transition temperature T_{P }was estimated for the (5,5) nanotube to correspond to temperature range T_{P }≃ 115 K [4,5,8]. Both barriers ΔE_{d }and ΔE_{r }and the thermal energy k_{B}T_{P }are of the same order of magnitude at the structural Peierls phase transition (see Table 1). Therefore, dramatic changes of the diffusion and drift over these barriers can take place at the Peierls transition for the considered pairs of the encapsulated molecules and the nanotube.
Recently, the diffusion and drift in the periodic potential surface of the interaction energy dependent on the displacement z along the nanotube axis and the angle φ of the rotation about the nanotube axis were considered for the threadlike relative motion of walls of doublewalled carbon nanotubes [28,29]. The expressions for the diffusion coefficient and mobility of a movable wall were obtained [28,29]. In this study, we consider specific cases of diffusion of the molecules along the nanotube axis and rotational diffusion about this axis. In this case, the expressions for diffusion coefficients, D_{d }and D_{r}, mentioned above and corresponding to the diffusion along the nanotube axis and rotational diffusion about this axis, respectively, take the form:
where Ω_{d }and Ω_{r }are the preexponential multipliers in the Arrhenius formula for the frequency of jumps of the molecule between two neighbouring global minima of the potential surface E_{W}(φ, z), δ_{d }is the distance between neighbouring global minima for the motion of the molecule along the nanotube axis, δ_{r }is the angle between neighbouring global minima corresponding to the molecule rotation about the nanotube axis and k_{B }is the Boltzmann constant. The mobility B_{d }for the motion along the axis can be easily obtained from the diffusion coefficient D_{d }using the Einstein ratio D_{d}/B_{d }= k_{B}T. Figure 3 shows that δ_{d }= 0.123 nm and δ_{d }= 0.37 nm for the (8,8) nanotube with the structure of the metallic phase and the Kekule structure, respectively, and δ_{r }= 22.5° for the both structures of this nanotube.
The value of the preexponential multiplier Ω in the Arrhenius formula is usually considered to be related with the frequency ν of corresponding vibrations. We suppose that the ratio Ω/ν remains the same for relative motion of different carbon nanoobjects with graphenelike structure (nanotube walls and fullerenes). For reorientation of the fullerenes of the C_{60}@C_{240 }nanoparticle, the frequency multiplier Ω was estimated by molecular dynamics simulations having the value of 650 ± 350 GHz [13,14]. We expand the potential surface of the intershell interaction energy near the minimum using the same empirical potential as in [13,14], and calculate the frequencies of small relative librations of the shells. The calculated libration frequency has the value ν ≈ 50 GHz, an order of magnitude less than that of the frequency multiplier Ω. In the estimations of this study, we use the values Ω_{d }≈ 10ν_{d }and Ω_{r }≈ 10ν_{z }for the preexponential multipliers.
The temperature dependencies of the diffusion coefficients, D_{d }and D_{r}, estimated using expressions (2) are shown in Figure 6. The dependence of the interaction energy E_{W}(z) on the displacement of the fullerene C_{20 }has two different barriers between the neighbour minima for the case of the nanotube with the Kekule structure (see Figure 4a). Since the frequency of the jumps of the fullerene between the neighbour minima exponentially depends on the barrier, the contribution of jumps over the lower barrier into the total diffusion coefficient is disregarded in these estimations. The temperature range corresponding to the Peierls transition temperature estimates T_{P }≃ 115 K [4,5,8] and the cases ΔE_{d}/k_{B}T < 1 and ΔE_{r}/k_{B}T < 1, where the Arrhenius formula is adequate, are considered. (For rotation of the Fe@C_{20 }inside the nanotube with the Kekule structure the Arrhenius formula is not applicable (ΔE_{r}/k_{B}T > 1) at this temperature range; this case is considered below.) Figure 6 shows that the changes of the diffusion coefficients, D_{d }and D_{r}, at the Peierls transition can be of orders of magnitude. It is of interest that the diffusion coefficient D_{r }for rotational diffusion of the C_{20 }fullerene decreases at the Peierls transition with the increase of temperature.
Figure 6. Dependences of the diffusion coefficients on the reciprocal of temperature 1/T. (a, c) The diffusion coefficient D_{d}; (b) the diffusion coefficient D_{r}. (a, b) The C_{20 }fullerene; (c) the Fe@C_{20 }endofullerene. Solid lines: the (8,8) carbon nanotube with the Kekule structure; dashed lines: the (8,8) carbon nanotube with the structure of metallic phase.
If a molecule is encapsulated inside a nanotube without a structural phase transition, the jump rotational diffusion takes place at low temperatures, ΔE_{r}/k_{B}T > 1, and the free rotation of the molecule occurs at high temperature, ΔE_{r}/k_{B}T < 1. Orientational melting (a loss of the orientational order with an increase of temperature) has a crossover behaviour if the structural phase transition is absent. Firstly, orientational melting was considered for twodimensional clusters with shell structure [3437] and later for doubleshell carbon nanoparticles [13,14], doublewalled carbon nanotubes [25,38,39] and carbon nanotube bundles [40]. In the case where a molecule is encapsulated inside a nanotube with a structural phase transition and the barrier ΔE_{r }for rotation of the molecule is greater for the hightemperature phase than for the lowtemperature phase, an inverse orientational melting (a loss of the orientational order with a decrease of temperature) is possible. In other words, the inverse orientational melting takes place if the Peierls transition temperature lies in the range ΔE_{rl }< k_{B}T_{P }< ΔE_{rh}, where ΔE_{rh }and ΔE_{rl }are the barriers for molecule rotation corresponding to hightemperature and lowtemperature phases of the nanotube, respectively. For the considered molecules inside the (8,8) nanotube, these temperature ranges are estimated to be 3.8 < T_{P }(K) < 10 and 0.58 < T_{P }(K) < 5.1 for the C_{20 }and Fe@C_{20 }fullerenes, respectively (see Table 1). As these temperature ranges are in agreement with the Peierls transition temperature estimates T_{P }≃ 115 K [4,5,8], we predict that the inverse orientational melting is possible for the systems considered. The inverse orientational melting should be more prominent for the case of the Fe@C_{20 }fullerene with the greater ratio of the barriers ΔE_{rh}/ΔE_{rl}.
Let us discuss the possibility of observing the changes of the dynamical behaviour of molecules inside armchair carbon nanotubes at the Peierls transition. We believe that the most promising method is highresolution transmission electron microscopy. This method was used for visualizing dynamics of processes inside nanotubes, such as reactions of fullerene dimerization with monitoring of timedependent changes in the atomic positions [41] and rotation of fullerene chains [42]. The rotational dynamics of C_{60 }fullerenes inside carbon nanotube was studied also by analysing the intermediate frequency mode lattice vibrations using nearinfrared Raman spectroscopy [43]. The orientational melting in a single nanoparticle may be revealed also by IR or Raman study of the temperature dependence of width of spectral lines. A specific heat anomaly in multiwalled carbon nanotubes may be caused by the orientational orderdisorder transition [44]. In the case of encapsulated magnetic molecules (for example, the Fe@C_{20 }endofullerene considered above), the study of the temperature dependence of the electron spin resonance spectra could yield information on the molecule rotational dynamics of these molecules [45].
Conclusive remarks
In this letter, we consider the changes of dynamical behaviour of fullerenes encapsulated in armchair carbon nanotubes caused by the Peierls transition in the nanotube structure by the example of the C_{20 }and Fe@C_{20 }fullerenes inside the (8,8) nanotube. We apply the DFT approach to calculate the structure of the C_{20 }and Fe@C_{20 }fullerenes. The ground state of the (8,8) nanotube is found to be the Kekule structure using the method of molecular orbitals. The LennardJones potential is used for calculating the barriers for motion of the fullerenes along the axis and rotation about the axis of the (8,8) nanotube with the Kekule structure and the structure with all equal bonds corresponding to lowtemperature and hightemperature phases, respectively. We show that the changes in the coefficients of diffusion of the fullerenes along the nanotube axis and their rotational diffusion at the Peierls transition can be as much as several orders of magnitude. The possibility of the inverse orientational melting at the Peierls transition is predicted. The analogous changes of dynamical behaviour are also possible for other large molecules inside armchair nanotubes. We believe that the predicted dynamical phenomena can be observed using highresolution transmission electron microscopy, nearinfrared Raman spectroscopy, specific heat measurements, and by study of electron spin resonance spectra for magnetic molecules.
Abbreviation
DFT: density functional theory.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
NAP and NAV set the problem of dynamical behaviour of small fullerenes inside carbon nanotubes. NAP, SAV, EFK, NNH and ONB performed the calculations of the fullerene and nanotube structures. NAP, NAV and EFK contributed to the analysis of the fullerene and nanotube structures. IVL, AAK, SAV and ONB performed the calculations of the interaction between the fullerenes and the nanotube. AMP and YEL proposed the idea of inverse orientational melting. AMP, YEL, NAP and EFK contributed to the discussion of the dynamical behaviour of the fullerenes inside the nanotube and wrote the original manuscript. All authors commented on the manuscript and approved the final manuscript.
Acknowledgements
This work has been partially supported by the RFBR (Grants 110200604a and 100290021Bel) and BFBR (Grant Nos. F10R062, F11V001). The atomistic calculations are performed on the SKIF MSU Chebyshev supercomputer and on the MVS100K supercomputer at the Joint Supercomputer Center of the Russian Academy of Sciences.
References

Takesue I, Haruyama J, Kobayashi N, Chiashi S, Maruyama S, Sugai T, Shinohara H: Superconductivity in entirely endbonded multiwalled carbon nanotubes.
Phys Rev Lett 2006, 96(5):057001. PubMed Abstract  Publisher Full Text

Bichoutskaia E, Heggie MI, Lozovik YuE, Popov AM: Multiwalled nanotubes: Commensurateincommensurate phase transition and NEMS applications.
Fullerenes, Nanotubes and Carbon Nanostructures 2006, 14(2):131140. Publisher Full Text

Connétable D, Rignanese GM, Charlier JC, Blase X: Room temperature Peierls distortion in small diameter nanotubes.
Phys Rev Lett 2005, 94:015503. PubMed Abstract  Publisher Full Text

Mintmire JW, Dunlap BI, White CT: Are fullerene tubules metallic?
Phys Rev Lett 1992, 68(5):631634. PubMed Abstract  Publisher Full Text

Sédéki A, Caron LG, Bourbonnais C: Electronphonon coupling and Peierls transition in metallic carbon nanotubes.

Viet NA, Ajiki H, Ando T: Lattice instability in metallic carbon nanotubes.
J Phys Soc Jpn 1994, 63(8):30363047. Publisher Full Text

Harigaya K, Fujita M: Dimerization structures of metallic and semiconducting fullerene tubules.
Phys Rev B 1993, 47(24):1656316569. Publisher Full Text

Huang Y, Okada M, Tanaka K, Yamabe T: Estimation of Peierlstransition temperature in metallic carbon nanotube.
Solid State Commun 1996, 97(4):303307. Publisher Full Text

Poklonski NA, Kislyakov EF, Hieu NN, Bubel' ON, Vyrko SA, Popov AM, Lozovik YuE: Uniaxially deformed (5,5) carbon nanotube: Structural transitions.
Chem Phys Lett 2008, 464(46):187191. Publisher Full Text

Zhou Z, Steigerwald M, Hybertsen M, Brus L, Friesner RA: Electronic structure of tubular aromatic molecules derived from the metallic (5,5) armchair single wall carbon nanotube.
J Am Chem Soc 2004, 126(11):35973607. PubMed Abstract  Publisher Full Text

Matsuo Y, Tahara K, Nakamura E: Theoretical studies on structures and aromaticity of finitelength armchair carbon nanotubes.
Org Lett 2003, 5(18):31813184. PubMed Abstract  Publisher Full Text

Nakamura E, Tahara K, Matsuo Y, Sawamura M: Synthesis, structure, and aromaticity of a hoopshaped cyclic benzenoid [10]cyclophenacene.
J Am Chem Soc 2003, 125(10):28342835. PubMed Abstract  Publisher Full Text

Lozovik YuE, Popov AM: Orientational melting of twoshell carbon nanoparticles: molecular dynamics study.
Chem Phys Lett 2000, 328(46):355362. Publisher Full Text

Lozovik YuE, Popov AM: Molecular dynamics study of orientational melting and thermodynamic properties of C_{60}@C_{240 }nanoparticles.
Phys Solid State 2002, 44:186194. Publisher Full Text

Zhou L, Pan ZY, Wang YX, Zhu J, Liu TJ, Jiang XM: Stable configurations of C_{20 }and C_{28 }encapsulated in single wall carbon nanotubes.
Nanotechnology 2006, 17(8):18911894. Publisher Full Text

Chen M, Zang J, Xiao D, Zhang C, Liu F: Nanopumping molecules via a carbon nanotube.
Nano Res 2009, 2(12):938944. Publisher Full Text

Straatsma TP, Aprá E, Windus TL, Dupuis M, Bylaska EJ, de Jong W, Hirata S, Smith DMA, Hackler MT, Pollack L, Harrison RJ, Nieplocha J, Tipparaju V, Krishnan M, Brown E, Cisneros G, Fann GI, Früchtl H, Garza J, Hirao K, Kendall R, Nichols JA, Tsemekhman K, Valiev M, Wolinski K, Anchell J, Bernholdt D, Borowski P, Clark T, Clerc D, Dachsel H, Deegan M, Dyall K, Elwood D, Glendening E, Gutowski M, Hess A, Jaffe J, Johnson B, Ju J, Kobayashi R, Kutteh R, Lin Z, Littlefield R, Long X, Meng B, Nakajima T, Niu S, Rosing M, Sandrone G, Stave M, Taylor H, Thomas G, van Lenthe J, Wong A, Zhang Z: NWChem, A Computational Chemistry Package for Parallel Computers, Version 4.5. Richland, WA: Pacific Northwest National Laboratory;

Becke AD: Densityfunctional thermochemistry. III. The role of exact exchange.
J Chem Phys 1993, 98(7):56485652. Publisher Full Text

Lee C, Yang W, Parr RG: Development of the ColleSalvetti correlationenergy formula into a functional of the electron density.
Phys Rev B 1988, 37(2):785789. Publisher Full Text

Ross RB, Powers JM, Atashroo T, Ermler WC, LaJohn LA, Christiansen PA: Ab initio relativistic effective potentials with spinorbit operators. IV. Cs through Rn.
J Chem Phys 1990, 93(9):66546670. Publisher Full Text

Stewart JJP: MNDO cluster model calculations on organic polymers.

Stewart JJP: Optimization of parameters for semiempirical methods.
J Comp Chem 1989, 10(2):209264. Publisher Full Text

Bubel' ON, Vyrko SA, Kislyakov EF, Poklonski NA: Totally symmetric vibrational modes of fullerene C_{60}.

Leclercq F, Damay P, Foukani M, Chieux P, BellissentFunel MC, Rassat A, Fabre C: Precise determination of the molecular geometry in fullerene C_{60 }powder: A study of the structure factor by neutron scattering in a large momentumtransfer range.
Phys Rev B 1993, 48(4):27482758. Publisher Full Text

Belikov AV, Lozovik YuE, Nikolaev AG, Popov AM: Doublewall nanotubes: classification and barriers to walls relative rotation, sliding and screwlike motion.
Chem Phys Lett 2004, 385(12):7278. Publisher Full Text

Girifalco LA, Hodak M, Lee RS: Carbon nanotubes, buckyballs, ropes, and a universal graphitic potential.
Phys Rev B 2000, 62(19):1310413110. Publisher Full Text

Kolmogorov AN, Crespi VH: Smoothest bearings: Interlayer sliding in multiwalled carbon nanotubes.
Phys Rev Lett 2000, 85(22):47274730. PubMed Abstract  Publisher Full Text

Lozovik YuE, Minogin AV, Popov AM: Nanomachines based on carbon nanotubes.
Phys Lett A 2003, 313(12):112121. Publisher Full Text

Lozovik YuE, Minogin AV, Popov AM: Possible nanomachines: Nanotube walls as movable elements.
JETP Lett 2003, 77(11):631635. Publisher Full Text

Damnjanović M, Dobardžić E, Milošević I, Vuković T, Nikolić B: Lattice dynamics and symmetry of double wall carbon nanotubes.

Damnjanović M, Vuković T, Milošević I: Superslippery carbon nanotubes.

Vuković T, Damnjanović M, Milošević I: Interaction between layers of the multiwall carbon nanotubes.

Matsushita K, Matsukawa H, Sasaki N: Atomic scale friction between clean graphite surfaces.
Solid State Commun 2005, 136:5155. Publisher Full Text

Lozovik YuE: Ion and electron clusters.
Sov Phys Usp 1987, 30(10):912913. Publisher Full Text

Lozovik YuE, Mandelshtam VA: Coulomb clusters in a trap.
Phys Lett A 1990, 145(5):269271. Publisher Full Text

Bedanov VM, Peeters FM: Ordering and phase transitions of charged particles in a classical finite twodimensional system.
Phys Rev B 1994, 49(4):26672676. Publisher Full Text

Lozovik YuE, Rakoch EA: Energy barriers, structure, and twostage melting of microclusters of vortices.
Phys Rev B 1998, 57(2):12141225. Publisher Full Text

Kwon YK, Tománek D: Electronic and structural properties of multiwall carbon nanotubes.
Phys Rev B 1998, 58(24):R16001R16004. Publisher Full Text

Bichoutskaia E, Heggie MI, Popov AM, Lozovik YuE: Interwall interaction and elastic properties of carbon nanotubes.
Phys Rev B 2006, 73(4):045435. Publisher Full Text

Kwon YK, Tománek D: Orientational melting in carbon nanotube ropes.
Phys Rev Lett 2000, 84(7):14831486. PubMed Abstract  Publisher Full Text

Koshino M, Niimi Y, Nakamura E, Kataura H, Okazaki T, Suenaga K, Iijima S: Analysis of the reactivity and selectivity of fullerene dimerization reactions at the atomic level.
Nat Chem 2010, 2(2):117124. PubMed Abstract  Publisher Full Text

Warner JH, Ito Y, Zaka M, Ge L, Akachi T, Okimoto H, Porfyrakis K, Watt AAR, Shinohara H, Briggs GAD: Rotating fullerene chains in carbon nanopeapods.
Nano Lett 2008, 8(8):23282335. PubMed Abstract  Publisher Full Text

Zou Y, Liu B, Wang L, Liu D, Yu S, Wang P, Wang T, Yao M, Li Q, Zou B, Cui T, Zou G, Wagberg T, Sundqvist B, Mao HK: Rotational dynamics of confined C_{60 }from nearinfrared Raman studies under high pressure.
Proc Natl Acad Sci 2009, 106(52):2213522138. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Jorge G, Bekeris V, Escobar M, Goyanes S, Zilli D, Cukierman A, Candal R: A specific heat anomaly in multiwall carbon nanotubes as a possible sign of orientational orderdisorder transition.
Carbon 2010, 48(2):525530. Publisher Full Text

Simon F, Kuzmany H, Náfrádi B, Fehér T, Forró L, Fülöp F, Jánossy A, Korecz L, Rockenbauer A, Hauke F, Hirsch A: Magnetic fullerenes inside singlewall carbon nanotubes.
Phys Rev Lett 2006, 97(13):136801. PubMed Abstract  Publisher Full Text