
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

72873
10.1038/s41598-024-72873-3
Article
Self-propelled directed transport of C60 fullerene on the surface of the cone-shaped carbon nanotubes
Vaezi Mehran 1
Nejat Pishkenari Hossein 2
1 https://ror.org/024c2fq17 grid.412553.4 0000 0001 0740 9747 Center for Nanoscience and Nanotechnology, Institute for Convergence Science and Technology, Sharif University of Technology, Tehran, 14588-89694 Iran
2 https://ror.org/024c2fq17 grid.412553.4 0000 0001 0740 9747 Nano Robotics Laboratory, Mechanical Engineering Department, Sharif University of Technology, Tehran, Iran
16 9 2024
16 9 2024
2024
14 2163025 6 2024
11 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Directed transportation of materials at molecular scale is important due to its crucial role in the development of nanoelectromechanical devices, particularly the directional movements along the carbon nanotubes (CNTs), due to the applications of CNTs as nano-manipulators, confined reactors, and drug or other materials delivery systems. In the present investigation, we evaluate the movements of C60 fullerenes on the surface of the cone-shaped CNTs. The fullerene molecules indicate directed motion toward the narrower end of CNTs, which is due to the potential energy gradient along the nanotube length. A continuum model is proposed to evaluate the mechanism of the directed motion and the results of the theoretical model are compared with numerical simulations. Directed movements have been examined at various opening angles of CNTs, considering the trajectories of motions, variation of potential energy, and diffusion coefficients. At smaller opening angles, the driving force on the C60 increases and the molecule experiences more directed transport along the nanotube. The motion of fullerene has also been simulated inside the cone-shaped CNTs, with similar opening angle, and different average radius. At lower average radius of the cone-like nanotubes, the motion of C60 is comparatively more rectilinear. Directional transport of fullerene has been observed in the opposite direction, when the molecule moves on the external surface of the cone-like CNTs, which is due to the stronger interaction of C60 with the parts of the external surface with larger radius. The effect of temperature has been evaluated by performing the simulations at the temperature range of 100 to 400 K. The direction of the velocity reveals that the thermal fluctuations at higher temperatures hinder the directed motion of molecule along the cone-shaped CNTs. The results of the present study propose a new method to obtain directed transport of molecules which can be helpful in different applications such as drug delivery systems.

Subject terms

Carbon nanotubes and fullerenes
Nanoscience and technology
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Achieving directed motion at nano-scale has a crucial role for several applications in the field of nanotechnology1,2. Manufacturing functional nano-devices through the directed assembly of smaller nanoparticles3, is among the applications of controlling the movements of low-dimensional materials. Steered locomotion of nanostructures can also be exploited to perform precise tasks at molecular dimensions4. As an instance, Regan et al.5 reported the fabrication of a mass conveyor based on carbon nanotubes, or Barreiro and her colleagues6 could transport the gold particle cargoes with sub-nanometer accuracy along the nanotube length.

Controlling the molecular motions along the nanotube is of particular importance, due to the special applications of this nanostructure. Due to the hollow structure of nanotubes, they are employed as nano-reactors7–9, in which several products can be achieved. Kuzmany et al.10 synthesized sub-nanometer graphene nano-ribbons by the reactions of ferrocene molecules inside the carbon nanotubes. The reaction of fullerene molecules has been studied inside the CNT, as well. At elevated temperatures, C60 fullerenes transform into different products in carbon nanotube such as inner nanotube11, and C60 oligomers and polymers12. Consequently, controllable transport of materials inside the CNT can help to control the kinetics of the reactions confined within the nanotubes. Moreover, carbon nanotubes have shown considerable potential as the atomic force microscope tips13–15. As a result, controlling the maneuverability of particles on the CNT can provide the opportunity to design precise nano-manipulators16,17.

Directional motion at molecular scale has been reported in previous investigations18–21, using various techniques. Rectilinear transport of materials has been obtained on the surfaces subjected to strain gradient22–25, on which adsorbed molecules move from regions with higher strain to less strained parts26. Constructing vacancies on the surface is another method to direct the motion of adsorbed materials27. Youzi et al.28 demonstrated the one-dimensional motion of molecular machines, surrounded by the series of vacancies, on the silicene surface. Applying temperature gradient on the substrate is another technique to steer the molecular motions29. The adsorbed materials migrate to the lower temperature area, due to the lower free energy level in these regions30. The nanoscale directed movement has also been obtained by angustotaxis effect31, where a trapped capped nanotube is derived toward the narrower end of a channel. The underlying physics is attributed to the higher contact area between the nanotube and the channel, at narrower side.

Several attempts were made so far to control the motion of nanoparticles along the carbon nanotubes. Previous studies reported the fabrication of artificial nanomotor consists of a short CNT which can translate relative to a coaxial inner longer nanotube6. The electric voltage applied to the ends of longer nanotube leads to a temperature gradient along the CNT. The shorter nanotube prefers to move toward the ends of longer CNT which have a lower temperature. The electromigration phenomenon has also been employed for achieving reversible mass transport along the carbon nanotubes5. By establishing electric current through the carbon nanotube, the indium nanoparticles experience directed movements on the external surface of the nanotube. The migration of the electrons can move indium particles through the collisions and momentum transfer. Using molecular dynamics simulations, Legoas et al.32 showed the presence of a retraction force on the extruded core of a double-walled carbon nanotube. After the release of extruded core, it is completely retracted inside the outer shell. The retraction force is due to the increase of contact area between the core and the outer shell. Machado and the colleagues33 reported the curvature driven motion of the encapsulated CNT inside a spiral-shaped carbon nanotube. Since the bending energy of the inner nanotube is proportional to the square of curvature, the encapsulated CNT experiences a driving force toward the regions of lower curvature. Despite proposing different techniques for nano-manipulations along the CNTs, using cone-shaped nanotubes did not receive attention in the previous investigations. The conical carbon nanotubes are experimentally synthesized through different methods such as arc discharge34, vapor deposition35, pyrolysis36 and bottom-up organic synthesis37,38. Exploiting conical nanotubes for controlling the movements seems more practical than the previous methods such as using a spiral-shaped CNT or an extruded core of the double-walled nanotube.

In the present study, we investigate the directed transport of fullerene molecule along the cone-like carbon nanotubes. Directed motions of C60 fullerenes are studied inside the CNTs with different opening angles. For this purpose, we consider the trajectories of motion, variation of the potential energy, potential energy surface (PES) analysis, and the diffusion coefficients. A continuum model has been proposed to evaluate the mechanism of the directed motion. At a similar opening angle of the cone-like CNTs, we evaluate the effect of average radius, by calculating the distribution of the fullerene velocity, and the trajectories of motion. The transportation of molecules has also been studied outside the nanotubes. In case of moving outside the cone-shaped nanotube, the fullerene molecule shows directed motion in opposite direction compared with the motion inside the CNT. Ultimately, we consider the effect of temperature on the directionality of the motion.

Computational methods

The motion of C60 fullerene has been simulated on the surface of the cone-like CNTs, using an all-atom molecular dynamics method. The fullerene molecule moves inside/outside the nanotubes, with different opening angles, as illustrated in Fig. 1. The radius of the CNTs decreases linearly along their length. Different cone angles (hereafter called as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\alpha\:$$\end{document} angle) have been chosen including 112.9, 60.0, and 19.2 degrees, while the slant length of CNTs is 80Å. The atomic structures of the carbon nanotube were obtained from Nanotube Modeler39. The fullerene nanostructure contains 60 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:s{p}^{2}$$\end{document} carbon atoms, arranged in 20 hexagonal and 12 pentagonal faces40.

Fig. 1 One of the cone-shaped nanotubes used in the simulations. The slant length of the nanotube is 80 Å, and the opening angle is 60.0 degrees. The color bar indicates the variations of the CNT radius. The axes of nanotubes are parallel to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z$$\end{document}-axis of the coordinate system.

Different simulations have been performed in the canonical NVT ensemble, which was implemented by the Nose-Hoover thermostat41,42. The temperature of the simulations was set to 100 K, and the damping parameter of the thermostats was 100 fs. It should be mentioned that, in Sect. 3.4, where we evaluate the effect of temperature, the temperature of the system changes from 100 K to 400 K. The inter-molecular interactions of C60 and nanotubes were described by Lennard-Jones (LJ) potential, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\varphi\:=4\epsilon\:\left({\left(\frac{\sigma\:}{l}\right)}^{12}-{\left(\frac{\sigma\:}{l}\right)}^{6}\right)$$\end{document}. When the distance between a pair of atoms (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:l$$\end{document}) is equal to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{2}^{1/6}\sigma\:$$\end{document}, the atoms find the minimum potential energy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\epsilon\:$$\end{document}. In our simulations, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\epsilon\:$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\sigma\:$$\end{document} parameters were assumed, 2.4 meV and 3.4Å, respectively43; which is suitable for describing the interactions between carbon atoms44. The cut-off radius of the pair potential was considered 12 Å. The inner-interactions of C60 atoms were captured by Tersoff potential45. To understand the effect of CNTs characteristics on the transportation more precisely, the nanotubes had fixed structures during the simulations.

The fullerene molecule was initially placed on the top part of CNTs with the largest radius (red region in Fig. 1), and the potential energy was minimized to find its local stable position. The initial velocity of C60 fullerene corresponds to the temperature of the simulation systems. The reproducibility of the results was examined by repeating the simulations, at different initial conditions of C60 (i.e., different distribution of the initial velocity). All of the simulations were conducted by Large-scale Atomic/Molecular Massively Parallel Simulator (i.e., LAMMPS package)46. The equations of the motions were integrated by velocity Verlet method47, and the time step of the integrations was adjusted to 1 fs. Each simulation continued until the C60 reached the narrower part of the CNT (blue parts in Fig. 1).

Results and discussion

Mechanism of the motion: Continuum approach and numerical simulations

To understand the mechanism of the directed motion of C60, the analysis of the potential energy has been performed. In this analysis, the fullerene molecule is translated to different points in horizontal plane (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:x$$\end{document}-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:y$$\end{document} plane). At each horizontal position (e.g., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:({x}_{0},{y}_{0})$$\end{document}), the molecule is displaced to different heights along the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z$$\end{document}-axis. The potential energy between C60 and carbon nanotube is calculated at each height, and the minimum potential energy of the molecule is captured. It should be noted that, the translation of C60 occurs with the step-length of 0.5 Å in the directions of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:x$$\end{document}-, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:y$$\end{document}-, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z$$\end{document}-axis, and the nanotube has the opening angle of 38.9 degrees. Figure 2 shows the results of the analysis of potential energy. Based on this analysis, as we move from the wider end of the CNT (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{R}_{\text{m}\text{a}\text{x}}$$\end{document}) to the narrower end (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{R}_{\text{m}\text{i}\text{n}}$$\end{document}), the reduction of the potential energy is observed. The potential energy gradient along the length of CNT provides the driving force on the fullerene. Hence, the C60 molecule demonstrates directed motion inside the nanotube. At the wider end of the CNT, the potential energy of fullerene is almost − 0.6 eV, which is similar to the case where the C60 is on a graphene surface48. The potential energy level decreases to -1.1 eV, at the narrower end of the CNT, which is approximately equal to the potential energy of a C60, sandwiched between two graphene layers. Since the fullerene is more stable at the narrower end of the CNT, it experiences a directed movement along the cone-like nanotubes.

Fig. 2 Analysis of the potential energy of C60. At each horizontal position such as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:({x}_{0},{y}_{0})$$\end{document}, we displaced the C60 in the direction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z$$\end{document}-axis, and the minimum energy has been captured which is corresponding to a specific height of molecule (e.g., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{z}_{0}$$\end{document}).

The variation of the potential energy of C60 has been analyzed, using a continuum model. The continuum model is indeed a valuable tool for predicting the interactions between the materials. The continuum model assumes that matter is continuous and completely fills the space it occupies, while the materials consist of discrete atoms. In the present study, the surface of the nanotubes is considered continuous. This assumption is true if the materials include a large number of atoms. In case of nanomaterials with few number of atoms the continuous assumption is not right. Here, since the conical nanotubes contain thousands of atoms we are allowed to use the assumption. Considering a discrete structure for the surface of nanotubes, we should calculate the summation of the vdW interactions between the adsorbed atom (an atom of fullerene) and all of the atoms of the cone-shaped CNT. On the other hand, when we use continuum model, the interaction of adsorbed atom and nanotube is obtained by calculating an integral over the surface of the cone-shaped CNT. As a result, the computational costs decrease as the continuum model is employed. However, the continuum model may lead to the complex integrals which should be solved by the numerical methods.

As mentioned in the model, it is assumed that the carbon atoms are uniformly distributed on the surface of the cone-shaped CNT. As a result, the potential energy between an adsorbed atom (an atom of C60) and CNT is calculated by integrating the Lennard-Jones potential as follows,1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{adatom}=\iint \varphi \left(r,\theta\:,z\right)\:{\rho\:}_{A}\frac{1}{\text{sin}\alpha\:/2}rdrd\theta\:.$$\end{document}

In the recent equation, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\varphi\:$$\end{document} is the LJ potential energy which is a function of the coordinates of CNT’s atoms (i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:(r,\theta\:,z)$$\end{document}), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{\rho\:}_{A}$$\end{document} is the areal density of the atoms over the CNT surface, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\frac{1}{\text{sin}\alpha\:/2}rdrd\theta\:$$\end{document} is surface element of a cone-like CNT with the opening angle of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\alpha\:$$\end{document}. The 6–12 Lennard-Jones potential function is given by2 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\varphi\:=4\epsilon\:\left({\left(\frac{\sigma\:}{l}\right)}^{12}-{\left(\frac{\sigma\:}{l}\right)}^{6}\right)$$\end{document}

where, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\epsilon\:$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\sigma\:$$\end{document} are the LJ parameters and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:l$$\end{document} is the distance between the adatom and the surface element of CNT. Considering the position of adatom as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:({r}_{0},{\theta\:}_{0},{z}_{0})$$\end{document} and the position of the surface element as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:(r,\theta\:,z)$$\end{document}, the separation distance can be written as,3 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{l}^{2}={r}^{2}+{r}_{0}^{2}-2r{r}_{0}\text{cos}\left(\theta\:-{\theta\:}_{0}\right)+{\left(z-{z}_{0}\right)}^{2}.$$\end{document}

By substituting Eq. 7 into Eq. 1, the potential energy between adatom and con-like surface is rewritten as,4 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{adatom}={U}_{R}-{U}_{A}$$\end{document}

where, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{R}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{A}$$\end{document} are expressed as,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{R}=\frac{4\epsilon\:{\rho\:}_{A}{\sigma\:}^{12}}{\text{sin}\left(\alpha\:/2\right)}{\int\:}_{0}^{2\pi\:}{\int\:}_{{R}_{min}\:}^{{R}_{max}}\frac{rdrd\theta\:}{{\left[{r}^{2}+{r}_{0}^{2}-2r{r}_{0}\text{cos}\left(\theta\:-{\theta\:}_{0}\right)+{\left(rcot(\alpha\:/2)-{z}_{0}\right)}^{2}\right]}^{6}}$$\end{document}

5 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{A}=\frac{4\epsilon\:{\rho\:}_{A}{\sigma\:}^{6}}{\text{sin}\left(\alpha\:/2\right)}{\int\:}_{0}^{2\pi\:}{\int\:}_{{R}_{min}\:}^{{R}_{max}}\frac{rdrd\theta\:}{{\left[{r}^{2}+{r}_{0}^{2}-2r{r}_{0}\text{cos}\left(\theta\:-{\theta\:}_{0}\right)+{\left(rcot(\alpha\:/2)-{z}_{0}\right)}^{2}\right]}^{3}}.$$\end{document}

The cone equation in cylindrical coordinate is considered as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z=rcot\left(\alpha\:/2\right)$$\end{document}, and the apex of the cone-like CNT is located at the origin of the coordinate system as shown in Fig. 3d. The total potential energy between C60 and cone-shaped CNT is calculated by adding up the potential energies of all atoms in the C60 molecule.6 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{total}=\sum\limits_{i=1}^{60}{U}_{R}\left(i\right)-{U}_{A}\left(i\right)$$\end{document}

In the next step, the C60 was displaced along the nanotube length, while the molecule is located in equilibrium distance to the CNT wall. The trajectory of the translation of C60 inside the nanotube is demonstrated in Fig. 3d. At each position of C60 inside the cone-like CNT, we calculate the potential energy of molecule (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{U}_{total}$$\end{document}) by numerical integration of the Eq. 5. Figure 3 demonstrates the potential energy of C60 along the nanotubes length from the theoretical continuum model as well as MD simulations. In the theoretical model, C60 is displaced to different locations along the nanotube, and the potential energy of molecule is calculated at each position from the integrals shown in Eq. 5. Since the tube lies between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:15<z<55$$\end{document}, there is no surface atoms in the left vicinity of fullerene when its position is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z<15$$\end{document} and hence the potential energy between fullerene and tube increases. Similar increase in potential energy is observed when fullerene position is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z>55$$\end{document}. These sudden changes in potential energy act like energy barrier and prevent fullerene to exit from the conical tube surface except at sufficiently higher temperatures.

The variation of the potential energy of C60 has been shown in Fig. 3 as the molecule travels the length of cone-shaped CNT. The potential energy level decreases as C60 moves toward the narrower end of CNT, which confirms the presence of a driving force on the molecule. By decreasing the opening angle from 112.9∘ (Fig. 3a) to 60.0∘ (Fig. 3b), the potential energy level decays to lower values, which is attributed to the increase of contact area between C60 and CNT, at smaller opening angles. It is found that, the fullerene experiences larger energy gradient by decreasing the opening angles (e.g., 19.2∘ in Fig. 3c), which is equal to 0.28 eV. As a result, the decrease of opening angle leads to the faster transport of C60 along the nanotube length. Due to the increase of potential energy gradient (driving force) on the fullerene at smaller opening angles, we expect to observe more directed movements. The simulation results are also illustrated in Fig. 3a and c by hollow circle. The simulation results are completely compatible with the results obtained from the continuum model (dashed lines). The sudden increase of the potential energy in the results of the continuum model refers to energy level when the molecule is out of the cone-shaped nanotube.

Fig. 3 The potential energy of C60 as a function molecule’s height. C60 travels the length of cone-shaped CNT, and it is located in equilibrium distance to the CNT wall. The opening angles of nanotubes are (a) 112.9∘, (b) 60.0∘, and (c) 19.2∘. The variation of the potential energy has been obtained by presenting a continuum model and numerical simulations. (d) The trajectory of the displacements of C60 has been shown inside the nanotube.

Cone-shaped CNTs with different opening angles

In this section, we aim to study the effect of opening angles of CNTs, on the motion of fullerene molecules inside the nanotubes. The opening angles of imaginary cones are considered 112.9, 60.0, and 19.2 degrees. The radius of these CNTs decreases linearly, along their axes such as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:76.6\:\text{\AA\:}\to\:10.0\:\text{\AA\:}$$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:50.0\:\text{\AA\:}\to\:10.0\:\text{\AA\:}$$\end{document}, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:23.3\:\text{\AA\:}\to\:10.0\:\text{\AA\:}$$\end{document}. The slant length in different carbon nanotubes is equal to 80 Å. In Fig. 4, the snapshots of the motion of C60 fullerenes are demonstrated inside the cone-like CNTs from the top view. As we observe in this figure, the molecules experience directed motions to the narrower end of the nanotubes. However, the rectilinearity of the movements depends on the opening angles of the CNTs. According to Fig. 4, at larger opening angles (e.g., 112.9∘), although the molecule moves toward the narrower end of CNT, it does not walk on a straight path. As we decrease the cone angle to lower values (e.g., 60.0∘), the C60 shows more directed motion to the end of CNT. The translocation time also lessens by decreasing the opening angle of CNTs. At the opening angle of 112.9∘, the fullerene travels the length of CNT in 560 picoseconds; while, in case of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\alpha\:=60.0^\circ\:$$\end{document}, C60 passes the same length in 73 picoseconds. It seems that, by decreasing the opening angle, the driving force on the fullerene increases. The validity of the last claim would be examined by the analysis of potential energy.

Fig. 4 Snapshots of the motion of C60 inside the cone-like nanotubes with the opening angles of 112.9, 60.0 and 19.2 degrees. The color bars indicate the variation of the nanotubes radius.

The reproducibility of the directed motion has been examined, by performing the simulations with different initial conditions of C60. The simulations have been repeated five times, while the initial coordinates of fullerenes are similar, but the distribution of the initial velocity is different in each simulation. Figure 5 illustrates the trajectories of the motion of fullerene molecules inside the cone-like CNTs, from the top view. Based on the results of the simulations, at each opening angle of CNT, the fullerene prefers to move toward the narrower end of the nanotube. When the opening angle of the nanotube is 112.9∘ (Fig. 5a), the motion of C60 is almost diffusive; however, it reaches the narrower end of nanotube during the simulation. As the opening angle decreases to 60.0∘ and 19.2∘ (Fig. 5b and c), we observe more directed movements, and the fullerene travels the length of CNT in more straight paths. In case of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\alpha\:=19.2^\circ\:$$\end{document}, in most of the simulations, the molecule shows spiral motions as it approaches the end of nanotube (Fig. 5c). According to the analysis of potential energy (Fig. 3), the driving force on the C60 molecule increases by decreasing the opening angle of the cone-shaped nanotubes. Since a larger driving force is acting on the C60 at smaller opening angles, the molecule is less affected by the random forces of the thermal fluctuations. As a result, the fullerene molecule maintains the direction of the initial velocity, particularly in the azimuth direction. The preserved direction of the azimuth velocity leads to the observation of spiral motions at smaller opening angles.

The average trajectory of the motions was obtained for each opening angle (black curves in Fig. 5). At different opening angles of cone-shaped CNTs, the average trajectories almost show the directional transport of C60 molecules to the narrower end of nanotubes. However, the standard deviations are larger at different points for CNTs with larger opening angles (e.g., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\alpha\:=112.9^\circ\:$$\end{document}). On the other hand, the standard deviation decreases at smaller opening angles such as\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:19.2^\circ\:$$\end{document}, especially in the direction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:y$$\end{document}-axis. As a result, the fullerene experiences more directional movements inside the cone-shaped CNTs with smaller opening angles.

Fig. 5 Trajectories of the motion of fullerene on the CNTs with the opening angles of (a) 112.9, (b) 60.0 and (c) 19.2 degrees. At each opening angle of CNT, the distribution of the initial velocity is different on the fullerene atoms.

Using the following equation, mean square displacements has been measured for the directed motions, at different opening angles of CNTs.7 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:MSD=\:\langle{\left(z\left(t\right)-z\left(0\right)\right)}^{2}\rangle=2D{t}^{\beta\:}\:$$\end{document}

In Eq. 7, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z\left(t\right)$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:z\left(0\right)$$\end{document} are the coordinates of the C60 in the direction of nanotube length at time \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:t$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:0$$\end{document}, respectively. The angle brackets in this equation also refer to the ensemble averaging, which is conducted by averaging the square displacements over five simulations with different seed numbers. Figure 6 illustrates the MSDs of the directed movements of fullerene as the functions of simulation time. As we decrease the opening angles of CNTs, the MSD of the directed motion grows faster with time.

Fig. 6 Mean square displacements of fullerene have been obtained during the motion inside the CNTs with different opening angles.

The mean square displacements is related to the diffusion coefficient (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:D$$\end{document}) through the power-law relation presented in Eq. 7. In this equation, the diffusion regime is determined by the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\beta\:$$\end{document} parameter, which is named as anomaly parameter. Table 1 indicates that the diffusion coefficient of the directed motion increases by decreasing the opening angle of cone-like nanotubes. The larger diffusion coefficient of C60 is related to the growth of driving force at smaller opening angles (previously concluded from Fig. 3). According to Table 1, since the anomaly parameter is larger than 1, the fullerene motion follows a super-diffusive regime, at different angles of cone-like nanotubes. However, the anomaly parameter relatively finds lower values at larger opening angles of the CNTs.

Table 1 Diffusion coefficients and anomaly parameters of the directed motion C60 inside the nanotubes with different opening angles.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\alpha\:\:$$\end{document}(degree)	D( Å2/ps)	\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\beta\:$$\end{document}	
19.2	0.055	2.721	
60.0	0.016	2.567	
112.9	0.007	2.295	

Effect of average radius on the directed transportation

In this section, we evaluate the effect of average radius on the directed movements of fullerene. For this purpose, three types of nanotubes have been chosen, which are similar in opening angle, but vary in the average radius. The opening angle of three nanotubes is 38.9∘, while the radius of CNTs linearly changes from 7.9 Å to 34.5 Å, from 34.5 Å to 61.2 Å and from 61.2 Å to 87.9 Å. In these nanotubes, the average values of the radius are 21.2 Å, 47.85 Å, and 74.55 Å, respectively. Figure 7 demonstrates the trajectories of the motion fullerene inside these nanotubes, which have the same opening angle of 38.9∘ and the same slant length of 80 Å, but they are different in average radius. At the average radius of 21.2 Å (Fig. 7a), the fullerene travels the length of CNT in a more directed path. By increasing the average radius to 47.85 and 74.55 Å (Fig. 7b and c), we observe more random movements in the trajectories of the fullerene. As a result, the increase of average radius in cone-like nanotubes hinders the directed motion of molecules inside the tube. This result seems reasonable, because in the case of a very large average radius, the fullerene interacts with a lower number of CNTs atoms, and it hardly senses the curvature of the nanotube. On the other hand, when the radius of CNT is small enough (compared with molecule size), the directed transport is better observed along the cone-like nanotube.

Fig. 7 Trajectories of the motion of fullerene in cone-like nanotubes with similar opening angle of 38.9 degrees. The radius of nanotubes linearly changes as (a) 7.9 Å → 34.5 Å, (b) 34.5 Å → 61.2 Å and (c) 61.2 Å → 87.9 Å.

To better quantify the directionality of the motion, we calculated the distribution of the fullerene horizontal velocity relative to a unit vector that perpendicularly pints to the nanotube axis (named as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{n}$$\end{document}). At each point such as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:({x}_{0},{y}_{0},{z}_{0})$$\end{document}, the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{n}$$\end{document} vector is defined as,8 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{n}=-\widehat{r}=-{(x}_{0},{y}_{0},0)/\| \left({x}_{0},{y}_{0},0\right)\|$$\end{document}

where, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{r}$$\end{document} is the normal radial vector in the cylindrical coordinate system, and the brackets calculate the length of the vector. In case of a rectilinear motion to the apex of a cone-like CNT, the horizontal velocity of fullerene (i.e., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{(v}_{x},{v}_{y},0)$$\end{document}) is completely in the direction of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{n}$$\end{document}, at each step of the simulation. The distribution of the C60 horizontal velocity around the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{n}$$\end{document} has been obtained for different average radii of the CNT (Fig. 8). It should be mentioned that, analogous to the previous analysis in Fig. 7, the opening angle is similar in all of the nanotubes (38.9 degrees), and they are repeated five times to ensure the obtained results.

At the average radius of 21.2 Å (Fig. 8a), the horizontal velocity is mostly distributed in the directions around 0 degree, which means that C60 moves in the direction of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{n}$$\end{document} most of the time. Consequently, the molecule almost experiences a directed movement toward the apex of the cone-like CNT, when the average radius is 21.2 Å. At larger average radius of nanotube such as 47.85 Å, and 74.55 Å (Fig. 8b and c), we observe the growth of velocity distribution in the directions other than the direction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:\widehat{n}$$\end{document}. However, a larger portion of the velocity is still distributed around 0 degree, which implies the directed movement of the molecule. As a result, in case of larger average radius of CNT, the fullerene spends more time moving in the random directions, but it still prefers to reach the narrower end of the nanotube. The last results are in agreement with the trajectories of the motion, which have already been shown in Fig. 7.

Fig. 8 Distribution of the horizontal velocity of C60, in different angles relative to the direction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:-\widehat{r}$$\end{document}, when the radius of nanotubes changes as (a) 7.9 Å → 34.5 Å, (b) 34.5 Å → 61.2 Å and (c) 61.2 Å → 87.9 Å. The results are extracted from five simulations with different initial conditions of C60.

Moving outside the cone-shaped nanotubes

Using carbon nanotubes as the tips of the atomic force microscope49 leads to the higher imaging resolution50 as well as, the precise manipulation of molecules at nano-scale dimensions51. As a result, controlling the motion of adsorbed molecules outside the CNT is of interest in the field of nano-manipulation.

The variation of the potential energy of C60 molecule has been evaluated using the continuum model present in Sect. 3.1, as well as the numerical simulations. According to Fig. 9, the potential energy of the molecule decreases as it moves toward the wider end of the cone-shaped CNT. Compared to the case where C60 moves inside the nanotube (Fig. 3), the direction of the energy gradient is different. The molecule prefers to move from the narrower end of CNT to wider regions, because the molecule has lower potential energy on wider regions. This conclusion seems reasonable, because when the fullerene is on the wider end of the nanotube there is larger area of CNT under the C60 and the fullerene interacts with more CNTs atoms. The fullerene experiences a larger energy gradient when it moves outside the cone-shaped nanotubes with smaller opening angles. For example, at the opening angle of 19.2∘ (Fig. 9c) the potential energy of C60 changes by 0.1 eV from narrower end to wider side of the CNT. The energy gradients in Fig. 9 are smaller than the ones observed when C60 was inside the nanotube (Fig. 3). It should be mentioned that, as we observe in Fig. 9, the results of the theoretical model are in agreement with the variations of the energy recorded from the simulations.

Fig. 9 Potential energy of C60 when it moves outside the cone-shaped nanotubes with the opening angles of (a) 112.9∘, (b)60.0∘ and (c) 19.2∘. The variation of the potential energy was obtained from the continuum model as well as the numerical simulations.

The trajectories of the motion of C60 has been captured during the motion outside the cone-shaped nanotubes with different opening angles. As we observe in Fig. 10, the trajectories of the motion confirm the migration of fullerene molecule from the parts with smaller radius to the regions with larger radius. To examine the repeatability of the result, five simulations have been performed by changing seed number. At the opening angles of 112.9∘ and 60.0∘ (Fig. 10a and b), the average path of the motion is more straight. While at the opening angle of 19.2∘, since the molecule rotates around the CNT, the average trajectory indicates an almost spiral path.

Fig. 10 Trajectories of the motion of C60 outside the cone-like carbon nanotubes with the opening angles of (a) 112.9∘, (b)60.0∘ and (c) 19.2∘. The averages and standard deviations of the trajectories are demonstrated in the figure.

Directed transportation at different temperatures

To understand the effect of temperature on the directionality of motion inside the cone-like nanotubes, the simulations have been conducted at the temperature range of 100 to 400 K. At each temperature, we calculate different parameters such \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{v}_{r}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{v}_{\theta\:}$$\end{document}, which are the average values of radial and angular speeds of fullerene, in the cylindrical coordinate system. In the rectilinear transport of C60 toward the narrower end of CNT, the speed of molecule is in radial direction, while the angular speed has a negligible share. As a result, the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{v}_{r}/{v}_{\theta\:}$$\end{document} ratio has also been calculated, which is another criterion to understand the directionality of the fullerene motion.

According to Fig. 11a, the radial speed is more affected by the opening angle of CNTs, rather than the temperature of the environment. At the opening angle of 112.9∘, the radial speed increases by raising the temperature, which is related to the growth of thermal energy at higher temperatures. However, the increase of radial speed with the temperature is not observable clearly at smaller opening angles of CNTs. Figure 11b demonstrates the variation of angular speed at different temperatures, and opening angles of nanotubes. As the temperature rises, the angular speed of fullerene also increases, especially at smaller opening angles such as 38.9∘. The angular speed growth refers to the non-directional movements of the C60 molecule at higher temperatures. The ratio of radial speed to angular one (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{v}_{r}/{v}_{\theta\:}$$\end{document}) indicates the proper conditions for achieving directional motion inside the cone-like CNTs. According to Fig. 11c, the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{v}_{r}/{v}_{\theta\:}$$\end{document} ratio considerably increases as the opening angle of nanotube decreases, and when the temperature of the system is lower. The last conclusion is consistent with the trajectories of motion and potential energy gradients in Figs. 3 and 5, respectively. At smaller opening angle of CNT, larger driving force is acting on the C60, which leads to a directional movement inside the cone-like nanotube. Moreover, as we increase the temperature, the thermal fluctuations hinder the directed transport of fullerene and the molecule experiences non-directional movements in the random directions.

Fig. 11 The average values of (a) radial and (b) angular speeds of C60 at different temperatures and opening angles of cone-like nanotubes. (c) The ratio of radial to angular speed has been obtained as a criterion for the directional transport of molecule inside the CNTs.

Previous investigations on the one-dimensional directed transport of nanomaterials can be categorized to different groups. In some of the studies the one-dimensional movements occur on the surface. It is possible to steer the motion of materials on the surface, by applying some modifications on the substrates. Nemati et al.52 demonstrated the restricted movements of fullerene molecule along the edge of the atomic steps on the gold surface. Using hybrid substrates is another method to restrict the lateral motion of molecules on the surface53. According to the previous investigations, on the gold-silver hybrid substrate, the adsorbed molecules prefers to move on gold regions due to the lower potential energy of molecule on this area53. Using the atomic vacancies on the substrate is another technique which has been proposed to constrain the surface motion in one-dimension54. Although these modifications of the substrates help to restrict the surface motion in one-dimension, the direction of the movements is not under control in these methods. In some of the previous works, an external agent has been employed for achieving directional motions. The electric field of the scanning electron microscope is able to manipulate the particles on the surface in a desired direction55. This type of manipulation is commonly implemented in the experimental investigations56. By applying transverse vibration on the graphene nanoribbon, Chen et al.57 demonstrated that the adsorbed nanoparticle moves away from the fluctuating regions, a phenomenon named as fluctuotaxis. The difference in atomic fluctuations of the substrate behind and ahead of the nanoparticle gives rise to a driving force acting on the adsorbed particle.

Another group of investigations utilize carbon nanotubes as a track for the one-dimensional directed transportation. The nanotubes subjected to temperature gradient can direct the encapsulated nanoparticles toward the regions with lower temperature58, because the particle finds lower free energy on this area59. The directed transportation of nanomaterials on the CNTs is also available by passing the electric current through the carbon nanotubes60. According to the electromigration phenomenon, the direction of motion is controllable by changing the current direction. The directed motions caused by electromigration and temperature gradient are utilized to fabricate nanomachines based on carbon nanotubes5,6. Directed transport of materials along the carbon nanotubes is also achievable by performing some modifications on the structure of nanotubes. As previously mentioned, the spiral-shaped nanotubes are able to direct the motion of encapsulated nano-objects33. The curvature gradient in these nanotubes leads to a driving force toward the areas with lower curvatures. In this study, we showed that the conical structure of the nanotubes steers the motion of molecules located on the internal and external surfaces of the CNTs. The synthesis of cone-shaped nanotubes is experimentally feasible61. Therefore, it is conceivable to put this novel approach for directed transport of nanomaterials into practice.

Conclusions

Directed transport of C60 fullerene has been studied along the cone-like nanotubes, using molecular dynamics simulations. The radius of CNTs decreases linearly as we move along the nanotube length. According to the results of the simulations, the fullerene experiences directed movements toward the narrower end of the CNTs. The potential energy of the molecule is lower at the narrower end of the nanotube, which is due to the increase of interactions between the C60 molecule and the CNT surface.

Different opening angles were chosen for the nanotubes including 112.9, 60.0, and 19.2 degrees. The trajectories of the motion of fullerene revealed that, the C60 molecule has more directional movement inside the CNTs with smaller opening angles (e.g., 19.2∘). To understand the reason for this observation, variations of the potential energy were obtained at different opening angles of CNTs. Based on the results of the simulations, by decreasing the opening angle of nanotube, the driving force on the fullerene becomes larger. At smaller opening angles, the molecule finds a larger potential energy gradient, which leads to the rectilinear motion of C60. Based on the mean square displacement analysis, the fullerene showed larger diffusion coefficients as the opening angle of CNT decreases.

At a constant opening angle, the directionality of the motion was examined at different average radius of CNT. For this purpose, the opening angle of nanotubes was chosen 38.9∘, while the radius of CNTs linearly changed as 7.9 Å → 34.5 Å, 34.5 Å → 61.2 Å and 61.2 Å → 87.9 Å. The trajectories of the fullerene motion demonstrated that, the C60 has more diffusive motions inside the nanotubes with larger average radius (e.g. 61.2 Å → 87.9 Å). On the other hand, at smaller average radius of nanotubes (e.g. 7.9 Å → 34.5 Å), the fullerene showed a directed motion toward the narrower end of CNT. The distribution of the horizontal velocity of fullerene also revealed that, when the average radius of CNT is smaller, the directed transport is better achieved along the cone-like nanotube.

The motion of fullerene molecules has been studied outside the cone-like nanotube. The directed transport of C60 was observed in the opposite direction, when the molecule moves outside the CNT. It has been shown that, when the fullerene is located outside the CNT, it has stronger interactions with the parts with larger radius. The fullerene motion inside the cone-like nanotube has been evaluated at the temperature range of 100 to 400 K. The ratio of radial to angular speed (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\:{v}_{r}/{v}_{\theta\:}$$\end{document} ratio) indicated that the C60 molecule has more directed movement, at smaller opening angle of nanotube, and lower temperatures. As we increase the temperature, the thermal fluctuations hinder the directed transport of fullerene and the molecule experiences movements in random directions.

Author contributions

M.V.: performing simulations, writing manuscript, preparing figures, and providing discussion. H.N.P.: supervision, editing the manuscript, data analysis, providing discussions, proposing the idea of the project.

Data availability

The data that support the findings of this study are available on request from the corresponding author.

Declarations

Competing interests

The authors declare no competing interests.

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Sheikhpour M An effective nano drug delivery and combination therapy for the treatment of tuberculosis Sci. Rep. 2022 12 1 9591 10.1038/s41598-022-13682-4 35688860
Sheikhpour, M. et al. An effective nano drug delivery and combination therapy for the treatment of tuberculosis. Sci. Rep. 12(1), 9591 (2022).35688860 10.1038/s41598-022-13682-4
2. Kamazani FM Sotoodehnejad nematalahi F Siadat SD Pornour M Sheikhpour M A success targeted nano delivery to lung cancer cells with multi-walled carbon nanotubes conjugated to bromocriptine Sci. Rep. 2021 11 1 24419 10.1038/s41598-021-03031-2 34952904
Kamazani, F. M., Sotoodehnejad nematalahi, F., Siadat, S. D., Pornour, M. & Sheikhpour, M. A success targeted nano delivery to lung cancer cells with multi-walled carbon nanotubes conjugated to bromocriptine. Sci. Rep. 11(1), 24419 (2021).34952904 10.1038/s41598-021-03031-2
3. Chai Z Childress A Busnaina AA Directed assembly of nanomaterials for making nanoscale devices and structures: mechanisms and applications ACS nano 2022 16 11 17641 17686 10.1021/acsnano.2c07910 36269234
Chai, Z., Childress, A. & Busnaina, A. A. Directed assembly of nanomaterials for making nanoscale devices and structures: mechanisms and applications. ACS nano. 16(11), 17641–17686 (2022).36269234 10.1021/acsnano.2c07910
4. Wang, J. & Nakano, T. (eds) Collective rotational motion of bio-nanomachines via chemical and physical interactions. GLOBECOM 2022–2022 IEEE Global Communications Conference; : IEEE. (2022).
5. Regan B Aloni S Ritchie R Dahmen U Zettl A Carbon nanotubes as nanoscale mass conveyors Nature 2004 428 6986 924 927 10.1038/nature02496 15118721
Regan, B., Aloni, S., Ritchie, R., Dahmen, U. & Zettl, A. Carbon nanotubes as nanoscale mass conveyors. Nature. 428(6986), 924–927 (2004).15118721 10.1038/nature02496
6. Barreiro A Subnanometer motion of cargoes driven by thermal gradients along carbon nanotubes Science 2008 320 5877 775 778 10.1126/science.1155559 18403675
Barreiro, A. et al. Subnanometer motion of cargoes driven by thermal gradients along carbon nanotubes. Science. 320(5877), 775–778 (2008).18403675 10.1126/science.1155559
7. Miners SA Rance GA Khlobystov AN Chemical reactions confined within carbon nanotubes Chem. Soc. Rev. 2016 45 17 4727 4746 10.1039/C6CS00090H 27301444
Miners, S. A., Rance, G. A. & Khlobystov, A. N. Chemical reactions confined within carbon nanotubes. Chem. Soc. Rev. 45(17), 4727–4746 (2016).27301444 10.1039/C6CS00090H
8. Eskandari, S., Koltai, J., László, I., Vaezi, M. & Kürti, J. Formation of nanoribbons by carbon atoms confined in a single-walled carbon nanotube—A molecular dynamics study. J. Chem. Phys. 158(22). (2023).
9. Nakanishi R Thin single-wall BN-nanotubes formed inside carbon nanotubes Sci. Rep. 2013 3 1 1 6 10.1038/srep01385
Nakanishi, R. et al. Thin single-wall BN-nanotubes formed inside carbon nanotubes. Sci. Rep. 3(1), 1–6 (2013).10.1038/srep01385
10. Kuzmany H Well-defined sub-nanometer graphene ribbons synthesized inside carbon nanotubes Carbon 2021 171 221 229 10.1016/j.carbon.2020.08.065
Kuzmany, H. et al. Well-defined sub-nanometer graphene ribbons synthesized inside carbon nanotubes. Carbon. 171, 221–229 (2021).10.1016/j.carbon.2020.08.065
11. Khlobystov AN Carbon nanotubes: from nano test tube to nano-reactor ACS nano 2011 5 12 9306 9312 10.1021/nn204596p 22200191
Khlobystov, A. N. Carbon nanotubes: from nano test tube to nano-reactor. ACS nano. 5(12), 9306–9312 (2011).22200191 10.1021/nn204596p
12. Terrones M Visualizing fullerene chemistry Nat. Chem. 2010 2 2 82 83 10.1038/nchem.526 21124394
Terrones, M. Visualizing fullerene chemistry. Nat. Chem. 2(2), 82–83 (2010).21124394 10.1038/nchem.526
13. Wilson NR Macpherson JV Carbon nanotube tips for atomic force microscopy Nat. Nanotechnol. 2009 4 8 483 491 10.1038/nnano.2009.154 19662008
Wilson, N. R. & Macpherson, J. V. Carbon nanotube tips for atomic force microscopy. Nat. Nanotechnol. 4(8), 483–491 (2009).19662008 10.1038/nnano.2009.154
14. Siria A Niguès A Electron Beam detection of a nanotube scanning force microscope Sci. Rep. 2017 7 1 11595 10.1038/s41598-017-11749-1 28912433
Siria, A. & Niguès, A. Electron Beam detection of a nanotube scanning force microscope. Sci. Rep. 7(1), 11595 (2017).28912433 10.1038/s41598-017-11749-1
15. Gao Z Optical detection of individual ultra-short carbon nanotubes enables their length characterization down to 10 nm Sci. Rep. 2015 5 1 17093 10.1038/srep17093 26603487
Gao, Z. et al. Optical detection of individual ultra-short carbon nanotubes enables their length characterization down to 10 nm. Sci. Rep. 5(1), 17093 (2015).26603487 10.1038/srep17093
16. Wei H Control of length and spatial functionality of single-wall carbon nanotube AFM nanoprobes Chem. Mater. 2008 20 8 2793 2801 10.1021/cm7031465
Wei, H. et al. Control of length and spatial functionality of single-wall carbon nanotube AFM nanoprobes. Chem. Mater. 20(8), 2793–2801 (2008).10.1021/cm7031465
17. Vaezi, M. Nanomachines Based on Carbon Nanotubes. (2024).
18. Vaezi M Pishkenari HN Nemati A Mechanism of the motion of nanovehicles on hexagonal boron-nitride: a molecular dynamics study Comput. Mater. Sci. 2022 207 111317 10.1016/j.commatsci.2022.111317
Vaezi, M., Pishkenari, H. N. & Nemati, A. Mechanism of the motion of nanovehicles on hexagonal boron-nitride: a molecular dynamics study. Comput. Mater. Sci. 207, 111317 (2022).10.1016/j.commatsci.2022.111317
19. Vaezi, M., Nejat Pishkenari, H. & Nemati, A. Mechanism of C60 rotation and translation on hexagonal boron-nitride monolayer. J. Chem. Phys. 153(23). (2020).
20. Wang J Manesh KM Motion control at the nanoscale Small 2010 6 3 338 345 10.1002/smll.200901746 20013944
Wang, J. & Manesh, K. M. Motion control at the nanoscale. Small. 6(3), 338–345 (2010).20013944 10.1002/smll.200901746
21. Kianezhad M Youzi M Vaezi M Pishkenari HN Rectilinear motion of carbon nanotube on gold surface Int. J. Mech. Sci. 2022 217 107026 10.1016/j.ijmecsci.2021.107026
Kianezhad, M., Youzi, M., Vaezi, M. & Pishkenari, H. N. Rectilinear motion of carbon nanotube on gold surface. Int. J. Mech. Sci. 217, 107026 (2022).10.1016/j.ijmecsci.2021.107026
22. Zhang B Rapid programmable nanodroplet motion on a strain-gradient surface Langmuir 2019 35 7 2865 2870 10.1021/acs.langmuir.8b03774 30600996
Zhang, B. et al. Rapid programmable nanodroplet motion on a strain-gradient surface. Langmuir. 35(7), 2865–2870 (2019).30600996 10.1021/acs.langmuir.8b03774
23. Vaezi M Nejat Pishkenari H Ejtehadi MR Programmable transport of C60 by straining graphene substrate Langmuir 2023 39 12 4483 4494 10.1021/acs.langmuir.3c00180 36926912
Vaezi, M., Nejat Pishkenari, H. & Ejtehadi, M. R. Programmable transport of C60 by straining graphene substrate. Langmuir. 39(12), 4483–4494 (2023).36926912 10.1021/acs.langmuir.3c00180
24. Vaezi M Nejat Pishkenari H Toward steering the motion of surface rolling molecular machines by straining graphene substrate Sci. Rep. 2023 13 1 20816 10.1038/s41598-023-48214-1 38012233
Vaezi, M. & Nejat Pishkenari, H. Toward steering the motion of surface rolling molecular machines by straining graphene substrate. Sci. Rep. 13(1), 20816 (2023).38012233 10.1038/s41598-023-48214-1
25. Huang Y Zhu S Li T Directional transport of molecular mass on graphene by straining Extreme Mech. Lett. 2014 1 83 89 10.1016/j.eml.2014.12.006
Huang, Y., Zhu, S. & Li, T. Directional transport of molecular mass on graphene by straining. Extreme Mech. Lett. 1, 83–89 (2014).10.1016/j.eml.2014.12.006
26. Vaezi, M. Programmable oscillation of C60 inside carbon nanotubes subjected to strain gradient. J. Appl. Phys. 134(23). (2023).
27. Nirmalraj P Daly R Martin N Thompson D Motion of Fullerenes around Topological defects on metals: implications for the Progress of Molecular Scale devices ACS Appl. Mater. Interfaces 2017 9 9 7897 7902 10.1021/acsami.7b00408 28233982
Nirmalraj, P., Daly, R., Martin, N. & Thompson, D. Motion of Fullerenes around Topological defects on metals: implications for the Progress of Molecular Scale devices. ACS Appl. Mater. Interfaces. 9(9), 7897–7902 (2017).28233982 10.1021/acsami.7b00408
28. Youzi M Kianezhad M Vaezi M Pishkenari HN Motion of nanovehicles on pristine and vacancy-defected silicene: implications for controlled surface motion Phys. Chem. Chem. Phys. 2023 25 42 28895 28910 10.1039/D3CP02835F 37855185
Youzi, M., Kianezhad, M., Vaezi, M. & Pishkenari, H. N. Motion of nanovehicles on pristine and vacancy-defected silicene: implications for controlled surface motion. Phys. Chem. Chem. Phys. 25(42), 28895–28910 (2023).37855185 10.1039/D3CP02835F
29. Kianezhad M Youzi M Vaezi M Nejat Pishkenari H Unidirectional motion of C60-based nanovehicles using hybrid substrates with temperature gradient Sci. Rep. 2023 13 1 1100 10.1038/s41598-023-28245-4 36670148
Kianezhad, M., Youzi, M., Vaezi, M. & Nejat Pishkenari, H. Unidirectional motion of C60-based nanovehicles using hybrid substrates with temperature gradient. Sci. Rep. 13(1), 1100 (2023).36670148 10.1038/s41598-023-28245-4
30. Lohrasebi A Neek-Amal M Ejtehadi M Directed motion of C 60 on a graphene sheet subjected to a temperature gradient Phys. Rev. E 2011 83 4 042601 10.1103/PhysRevE.83.042601
Lohrasebi, A., Neek-Amal, M. & Ejtehadi, M. Directed motion of C 60 on a graphene sheet subjected to a temperature gradient. Phys. Rev. E. 83(4), 042601 (2011).10.1103/PhysRevE.83.042601
31. Leng J Hu Y Chang T Nanoscale directional motion by Angustotaxis Nanoscale 2020 12 9 5308 5312 10.1039/C9NR10108J 31872851
Leng, J., Hu, Y. & Chang, T. Nanoscale directional motion by Angustotaxis. Nanoscale. 12(9), 5308–5312 (2020).31872851 10.1039/C9NR10108J
32. Legoas S Molecular-dynamics simulations of carbon nanotubes as gigahertz oscillators Phys. Rev. Lett. 2003 90 5 055504 10.1103/PhysRevLett.90.055504 12633370
Legoas, S. et al. Molecular-dynamics simulations of carbon nanotubes as gigahertz oscillators. Phys. Rev. Lett. 90(5), 055504 (2003).12633370 10.1103/PhysRevLett.90.055504
33. Machado LD Bizao RA Pugno NM Galvão DS Controlling movement at nanoscale: curvature driven mechanotaxis Small 2021 17 35 2100909 10.1002/smll.202100909
Machado, L. D., Bizao, R. A., Pugno, N. M. & Galvão, D. S. Controlling movement at nanoscale: curvature driven mechanotaxis. Small. 17(35), 2100909 (2021).10.1002/smll.202100909
34. Berkmans J Jagannatham M Reddy R Haridoss P Synthesis of thin bundled single walled carbon nanotubes and nanohorn hybrids by arc discharge technique in open air atmosphere Diam. Relat. Mater. 2015 55 12 15 10.1016/j.diamond.2015.02.004
Berkmans, J., Jagannatham, M., Reddy, R. & Haridoss, P. Synthesis of thin bundled single walled carbon nanotubes and nanohorn hybrids by arc discharge technique in open air atmosphere. Diam. Relat. Mater. 55, 12–15 (2015).10.1016/j.diamond.2015.02.004
35. Chen I-C Fabrication of high-aspect-ratio carbon nanocone probes by electron beam induced deposition patterning Nanotechnology 2006 17 17 4322 10.1088/0957-4484/17/17/007
Chen, I-C. et al. Fabrication of high-aspect-ratio carbon nanocone probes by electron beam induced deposition patterning. Nanotechnology. 17(17), 4322 (2006).10.1088/0957-4484/17/17/007
36. Saito Y Arima T Growth of cone-shaped carbon material inside the cell lumen by heat treatment of wood charcoal J. Wood Sci. 2002 48 451 454 10.1007/BF00770709
Saito, Y. & Arima, T. Growth of cone-shaped carbon material inside the cell lumen by heat treatment of wood charcoal. J. Wood Sci. 48, 451–454 (2002).10.1007/BF00770709
37. Shoyama K Würthner F Synthesis of a carbon nanocone by cascade annulation J. Am. Chem. Soc. 2019 141 33 13008 13012 10.1021/jacs.9b06617 31380635
Shoyama, K. & Würthner, F. Synthesis of a carbon nanocone by cascade annulation. J. Am. Chem. Soc. 141(33), 13008–13012 (2019).31380635 10.1021/jacs.9b06617
38. Zhu Z-Z Rational synthesis of an atomically precise carboncone under mild conditions Sci. Adv. 2019 5 8 eaaw0982 10.1126/sciadv.aaw0982 31467971
Zhu, Z-Z. et al. Rational synthesis of an atomically precise carboncone under mild conditions. Sci. Adv. 5(8), eaaw0982 (2019).31467971 10.1126/sciadv.aaw0982
39. Melchor S Dobado JA CoNTub An algorithm for connecting two arbitrary carbon nanotubes J. Chem. Inf. Comput. Sci. 2004 44 5 1639 1646 10.1021/ci049857w 15446821
Melchor, S., Dobado, J. A. & CoNTub An algorithm for connecting two arbitrary carbon nanotubes. J. Chem. Inf. Comput. Sci. 44(5), 1639–1646 (2004).15446821 10.1021/ci049857w
40. Baskar AV Self-assembled fullerene nanostructures: synthesis and applications Adv. Funct. Mater. 2022 32 6 2106924 10.1002/adfm.202106924
Baskar, A. V. et al. Self-assembled fullerene nanostructures: synthesis and applications. Adv. Funct. Mater. 32(6), 2106924 (2022).10.1002/adfm.202106924
41. Nosé S A unified formulation of the constant temperature molecular dynamics methods J. Chem. Phys. 1984 81 1 511 519 10.1063/1.447334
Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 81(1), 511–519 (1984).10.1063/1.447334
42. Hoover WG Canonical dynamics: equilibrium phase-space distributions Phys. Rev. A 1985 31 3 1695 10.1103/PhysRevA.31.1695
Hoover, W. G. Canonical dynamics: equilibrium phase-space distributions. Phys. Rev. A. 31(3), 1695 (1985).10.1103/PhysRevA.31.1695
43. Vaezi M Pishkenari HN Ejtehadi MR Nanocar swarm movement on graphene surfaces Phys. Chem. Chem. Phys. 2022 24 45 27759 27771 10.1039/D2CP03856K 36349663
Vaezi, M., Pishkenari, H. N. & Ejtehadi, M. R. Nanocar swarm movement on graphene surfaces. Phys. Chem. Chem. Phys. 24(45), 27759–27771 (2022).36349663 10.1039/D2CP03856K
44. Rafii-Tabar H Computational modelling of thermo-mechanical and transport properties of carbon nanotubes Phys. Rep. 2004 390 4–5 235 452 10.1016/j.physrep.2003.10.012
Rafii-Tabar, H. Computational modelling of thermo-mechanical and transport properties of carbon nanotubes. Phys. Rep. 390(4–5), 235–452 (2004).10.1016/j.physrep.2003.10.012
45. Tersoff J Modeling solid-state chemistry: interatomic potentials for multicomponent systems Phys. Rev. B 1989 39 8 5566 10.1103/PhysRevB.39.5566
Tersoff, J. Modeling solid-state chemistry: interatomic potentials for multicomponent systems. Phys. Rev. B. 39(8), 5566 (1989).10.1103/PhysRevB.39.5566
46. Thompson AP LAMMPS-a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales Comput. Phys. Commun. 2022 271 108171 10.1016/j.cpc.2021.108171
Thompson, A. P. et al. LAMMPS-a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comput. Phys. Commun. 271, 108171 (2022).10.1016/j.cpc.2021.108171
47. Lippert, R. A. et al. A common, avoidable source of error in molecular dynamics integrators. J. Chem. Phys. 126(4). (2007).
48. Vaezi M Pishkenari HN Ejtehadi MR Collective movement and thermal stability of fullerene clusters on the graphene layer Phys. Chem. Chem. Phys. 2022 24 19 11770 11781 10.1039/D2CP00667G 35506871
Vaezi, M., Pishkenari, H. N. & Ejtehadi, M. R. Collective movement and thermal stability of fullerene clusters on the graphene layer. Phys. Chem. Chem. Phys. 24(19), 11770–11781 (2022).35506871 10.1039/D2CP00667G
49. Cheng B Controlled growth of a single carbon nanotube on an AFM probe Microsystems Nanoengineering 2021 7 1 80 10.1038/s41378-021-00310-w 34721888
Cheng, B. et al. Controlled growth of a single carbon nanotube on an AFM probe. Microsystems Nanoengineering. 7(1), 80 (2021).34721888 10.1038/s41378-021-00310-w
50. Cheung, C. L., Hafner, J. H. & Lieber, C. M. Carbon nanotube atomic force microscopy tips: Direct growth by chemical vapor deposition and application to high-resolution imaging. Proceedings of the National Academy of Sciences. ;97(8):3809-13. (2000).
51. Firouzi MM Pishkenari HN Mahboobi S Meghdari A Manipulation of biomolecules: a molecular dynamics study Curr. Appl. Phys. 2014 14 9 1216 1227 10.1016/j.cap.2014.06.014
Firouzi, M. M., Pishkenari, H. N., Mahboobi, S. & Meghdari, A. Manipulation of biomolecules: a molecular dynamics study. Curr. Appl. Phys. 14(9), 1216–1227 (2014).10.1016/j.cap.2014.06.014
52. Nemati A Pishkenari HN Meghdari A Sohrabpour S Directing the diffusive motion of fullerene-based nanocars using nonplanar gold surfaces Phys. Chem. Chem. Phys. 2018 20 1 332 344 10.1039/C7CP07217A
Nemati, A., Pishkenari, H. N., Meghdari, A. & Sohrabpour, S. Directing the diffusive motion of fullerene-based nanocars using nonplanar gold surfaces. Phys. Chem. Chem. Phys. 20(1), 332–344 (2018).10.1039/C7CP07217A
53. Nemati A Nejat Pishkenari H Meghdari A Ge SS Controlling the diffusive motion of fullerene-wheeled nanocars utilizing a hybrid substrate J. Phys. Chem. C 2019 123 42 26018 26030 10.1021/acs.jpcc.9b08335
Nemati, A., Nejat Pishkenari, H., Meghdari, A. & Ge, S. S. Controlling the diffusive motion of fullerene-wheeled nanocars utilizing a hybrid substrate. J. Phys. Chem. C. 123(42), 26018–26030 (2019).10.1021/acs.jpcc.9b08335
54. Nemati A Nejat Pishkenari H Meghdari A Ge SS Influence of vacancies and grain boundaries on the diffusive motion of surface rolling molecules J. Phys. Chem. C 2020 124 30 16629 16643 10.1021/acs.jpcc.0c03697
Nemati, A., Nejat Pishkenari, H., Meghdari, A. & Ge, S. S. Influence of vacancies and grain boundaries on the diffusive motion of surface rolling molecules. J. Phys. Chem. C. 124(30), 16629–16643 (2020).10.1021/acs.jpcc.0c03697
55. Akimov AV Kolomeisky AB Unidirectional rolling motion of nanocars induced by electric field J. Phys. Chem. C 2012 116 42 22595 22601 10.1021/jp306938b
Akimov, A. V. & Kolomeisky, A. B. Unidirectional rolling motion of nanocars induced by electric field. J. Phys. Chem. C. 116(42), 22595–22601 (2012).10.1021/jp306938b
56. Shirai Y Osgood AJ Zhao Y Kelly KF Tour JM Directional control in thermally driven single-molecule nanocars Nano Lett. 2005 5 11 2330 2334 10.1021/nl051915k 16277478
Shirai, Y., Osgood, A. J., Zhao, Y., Kelly, K. F. & Tour, J. M. Directional control in thermally driven single-molecule nanocars. Nano Lett. 5(11), 2330–2334 (2005).16277478 10.1021/nl051915k
57. Chen Y Fluctuotaxis: Nanoscale directional motion away from regions of fluctuation Proc. Natl. Acad. Sci. 2023 120 31 e2220500120 10.1073/pnas.2220500120 37487105
Chen, Y. et al. Fluctuotaxis: Nanoscale directional motion away from regions of fluctuation. Proc. Natl. Acad. Sci. 120(31), e2220500120 (2023).37487105 10.1073/pnas.2220500120
58. Wei N Wang H-Q Zheng J-C Nanoparticle manipulation by thermal gradient Nanoscale Res. Lett. 2012 7 1 9 10.1186/1556-276X-7-154 22214494
Wei, N., Wang, H-Q. & Zheng, J-C. Nanoparticle manipulation by thermal gradient. Nanoscale Res. Lett. 7, 1–9 (2012).22214494 10.1186/1556-276X-7-154
59. Lohrasebi A Neek-Amal M Ejtehadi M Directed motion of C 60 on a graphene sheet subjected to a temperature gradient Phys. Rev. E—Statistical Nonlinear Soft Matter Phys. 2011 83 4 042601 10.1103/PhysRevE.83.042601
Lohrasebi, A., Neek-Amal, M. & Ejtehadi, M. Directed motion of C 60 on a graphene sheet subjected to a temperature gradient. Phys. Rev. E—Statistical Nonlinear Soft Matter Phys. 83(4), 042601 (2011).10.1103/PhysRevE.83.042601
60. Kulshrestha N Misra A Misra D Electrical transport and electromigration studies on nickel encapsulated carbon nanotubes: possible future interconnects Nanotechnology 2013 24 18 185201 10.1088/0957-4484/24/18/185201 23575106
Kulshrestha, N., Misra, A. & Misra, D. Electrical transport and electromigration studies on nickel encapsulated carbon nanotubes: possible future interconnects. Nanotechnology. 24(18), 185201 (2013).23575106 10.1088/0957-4484/24/18/185201
61. Zhang Q The synthesis of conical carbon Small Methods 2021 5 3 2001086 10.1002/smtd.202001086
Zhang, Q. et al. The synthesis of conical carbon. Small Methods. 5(3), 2001086 (2021).10.1002/smtd.202001086
