Abstract
With the synthesis of nanocar structures the idea of transporting energy and payloads on the surface became closer to reality. To eliminate the concern of diffusive surface motion of nanocars, in this study, we evaluate the motion of C60 and C60-based nanovehicles on graphene and hexagonal boron-nitride (BN) surfaces using molecular dynamics simulations and potential energy analysis. Utilizing the graphene-hBN hybrid substrate, it has been indicated that C60 is more stable on boron-nitride impurity regions in the hybrid substrate and an energy barrier restricts the motion to the boron-nitride impurity. Increasing the temperature causes the molecule to overcome the energy barrier frequently. A nanoroad of boron-nitride with graphene sideways is designed to confine the surface motion of C60 and nanovehicles at 300 K. As expected, the motion of all surface molecules is limited to the boron-nitride nanoroads. Although the motion is restricted to the boron-nitride nanoroad, the diffusive motion is still noticeable in lateral directions. To obtain the unidirectional motion for C60 and nanocars on the surface, a temperature gradient is applied to the surface. The unidirectional transport to the nanoroad regions with a lower temperature occurs in a short period of time due to the lower energies of molecules on the colder parts.
Similar content being viewed by others
Introduction
Understanding the motion of molecules on surfaces is of great importance in developing nano-transportation systems on the surface1. Using experimental and computational studies, many researchers have tried to understand the surface motion of various molecules in the past decade2,3,4,5,6,7,8,9. The experimental tools such as scanning tunneling microscopy (STM) require highly specialized equipment for investigating the surface phenomena at nanometer scale10,11. On the other hand, computational methods permit us to evaluate the effect of different parameters on the surface motion of nano-materials with lower costs.
In the primary studies, it was shown that the nanocars experience a diffusive motion on the substrate at different temperatures12,13. Recent computational attempts are mostly focused on proposing a controllable nanocar as a nanocarrier in the delivery of molecules14,15. Rectilinearity of motion has been achieved in some investigations using an external agent such as temperature gradient16 or electrical field17,18. In another group of studies, the directional surface motion of molecules has been obtained using surface modifications19,20. Using a hybrid substrate of gold and silver, Nemati et al.21 demonstrated that, the motion of molecules could be restricted to gold regions, bringing forth the idea of nanoroads. In other studies, it has been illustrated that C60 and fullerene-based nanocars experience directed movements using a thermal gradient on the substrate22,23,24. Kianezhad et al.25 have also studied the unidirectionality of the motion of carbon nanotubes on a golden substrate by evaluating the effects of chirality, diameter, temperature, and length.
The boron-nitride and graphene monolayers are capable of forming surface ripples (up to 1 nm out of plane, surface deformation in z-direction26). As a result, boron-nitride and graphene sheets demonstrate a wavy structure which can affect the surface motions27,28. Considering the structural and electrical characteristics29,30, graphene sheets can be considered as a potential candidate for nanoscale transportation applications. Mofidi et al.26 have studied the effect of thermal vibrations of graphene and surface ripples on the dynamical behavior of C60. They realized that surface ripples could result in decreasing desorption temperature due to the increase of the fluctuations of the molecule.
Similar to the atomic structure of graphene, a boron-nitride sheet is composed of a hexagonal structure of three boron and three nitrogen atoms. The boron-nitrogen bond length is 1.44 Å31, which is quite similar to the carbon–carbon bond length of 1.42 Å32 in graphene. Based on their structural similarity, graphene and boron-nitride sheets can be combined to form a hybrid substrate for transportation purposes at the nanoscale33. Vaezi et al.34 studied the mechanism of the motion of nanovehicles on hexagonal boron-nitride. They indicated that nanovehicles obtain higher mobility on the boron-nitride compared to the metallic surfaces. In another study, they investigated the motion of a C60 molecule on a hexagonal boron-nitride monolayer using molecular dynamics and potential energy approach to evaluate the effect of temperature, where they showed that C60 has lower potential energy on boron-nitride33.
Controlling the motion of molecules on the surface plays important role in the surface phenomena such as chemical reactions and catalysis35,36. The directional motions of nanomaterial are beneficial in the fabrication of nanostructures on the surface through the bottom-up approach37. On the other hand, the nanocars have been synthesis to carry payloads on the surface38. These nanomachines have shown the potential of transporting cargoes on the surface39. As a result, the directed movement of these nano-vehicles are also useful for designing nano-scale transportation and delivery systems on the surface.
In the present study, we evaluated the motion of C60 and C60-based nanovehicles on the graphene/hBN hybrid substrates to steer the surface motion, using the molecular dynamics method and the potential energy analysis. An initial assessment of the surface motion of molecules is presented based on the variation of the C60 potential energy on the BN and graphene surfaces. Using molecular dynamics simulations, the effect of impurities is investigated, suggesting C60 inclines to remain on the boron-nitride as it is more stable. Afterwards, using a nanoroad of boron-nitride inside a graphene substrate, we try to control the motion of objects on the surface. Finally, a combination of boron-nitride nanoroad and the temperature gradient is used to further steer the motion.
Methods
In order to reach a comprehensive understanding of the dynamic behavior of nanomachines on the surface, we used both molecular dynamics simulations and potential energy analysis. In this study, the nanocar and the nanotruck results are compared with C60 since both are composed of four fullerene wheels. C60 is a rigid spherical molecule with Ih symmetry (Fig. 1c). It weighs \(720.6\) g/mol, made of only carbon, and it is stable at high temperatures40,41. The primary difference between C60-based nanocar and nanotruck is the chassis22. The C60-based nanocar (Fig. 1a) is made of carbon and hydrogen, weighs \(3757.587 \frac{\mathrm{g}}{\mathrm{mol}}\), and has C2 symmetry42. On the other hand, the chassis structure in a C60-based nanotruck (Fig. 1b) is different; besides carbon and hydrogen, there are four nitrogen atoms, making the chassis more rigid compared to the nanocar. The C60-based nanotruck weighs \(3602.202 \frac{\mathrm{g}}{\mathrm{mol}}\) and it has a D2 symmetry which makes it more symmetrical than nanocar. Also, C60-based nanotruck is noticeably smaller than a nanocar38. Neither C60 nor nanocar and nanotruck has electric dipoles21. A schematic view of C60 and C60-based nanocar and nanotruck is indicated in Fig. 1, showing the atomic structure, molecular size, and chemical formulation of each molecule.
In this study, Tersoff potential43 has been employed to model the inner interactions of the graphene and boron-nitride monolayer. Also, substrate ripples are taken into account since the wavy surface can change the dynamical behavior of the nanocar on the graphene and boron-nitride monolayer due to the variation of contact level29,34. The internal interaction of fullerene is modeled using the Tersoff potential43 as well. Molecular mechanics (MM3) force field44,45 has been used to study the internal interactions of the C60-based nanocar and nanotruck. The bond and angle terms are considered by the harmonic style, while the dihedral style of OPLS force field was employed to calculate the torsional terms46,47. It is worth mentioning that improper terms are not included in the calculations. Two different types of carbon atoms are considered in the simulations as "\(CA\)" and "\(C2\)" based on their hybridization. "\(CA\)" atoms have sp2 or sp3 hybridization while "\(C2\)" carbon atoms have sp34. Different atom types are observable in Fig. 1. Supplementary material includes different parameters employed in the simulation of nanovehicles and fullerene.
Utilizing the Restricted Hartree–Fock (RHF) method, Vaezi et al.34 determined the partial atomic charge of each particle in the C60-based nanocar and nanotruck. Akimov et al.17 illustrated that with the absence of an external electric field, the charge transfer between non-bonded atoms is negligible. This conclusion has also been used in various studies17,21,22,23,48,49. As a result, we neglected the electrostatic forces to decrease the computational cost of the investigation since the charge transfer is not large enough to noticeably affect the motion. In order to account for the interactions between molecules (C60 and C60-based nanocar and nanotruck) with the graphene and boron-nitride monolayers, the 6–12 Lennard–Jones (LJ) potential with a cut-off radius of 13 Å is used. To be more specific, the LJ potential function employed to calculate the non-bonded interactions is as follow50:
where \(\varepsilon\) and \(\sigma\) are non-bonding parameters of LJ potential. Table S4 indicates the parameters used to model the non-bonding interactions. As shown in Table S4, identical LJ non-bonding parameters are used to model the interactions between the substrate and both carbon types in the nanovehicles34,51,52,53. Also, the non-bonding interactions between C60 and the substrate are considered the same as the “C2” (and “CA”) parameters, indicated in Table S4.
Noose–Hover (NVT ensemble) thermostat is employed to keep the desired temperature of the system constant. The temperature damping parameter of the thermostat (Tdamp), which determines the heat exchange rate between the thermostat and the system, is set to 100 fs. Then, before the primary simulation, the simulation system is allowed to relax freely at the desired temperatures for 1 ns to let the system reach an equilibration condition. The equations of motion were integrated using the velocity Verlet algorithm, and the time-step was set to 1 fs to achieve precise simulations. LAMMPS54 package was used for all the simulations.
The efficient application of nanomachines such as nanocars and nanotrucks requires a directional movement20. In this study, mean square displacement (MSD) is used to describe the diffusive motion of the surface-moving molecules. The diffusive coefficient is another parameter that can be derived from the MSD in order to characterize the diffusive motion. MSD is defined as55:
where r is the center of the mass position vector. In statistical mechanics, < > Expresses the ensemble averaging on the phase space, which is equal to the time integration of the system in ergodic terms56,57. In other words, < > defines as averaging over all atoms. The slope of the MSD is defined as diffusion coefficient55.
Since we are investigating the horizontal motion on the substrate in the XY plane, disregarding the vertical vibration, b is considered 4. However, for three-dimensional simulations, b will be 6.
For an accurate estimation of MSD, many simulations must be performed and then averaging should be done regarding the trajectory of all of them. Since this method is time-consuming, another way of calculating the MSD is to perform a relatively long simulation and then divide the simulation into smaller sub-trajectories. By assuming that each short-time trajectory is independent of others, averaging could be done over all of them58. In the present research, for MSD calculations, the data in each simulation is divided into 30 smaller simulations with 500 ps duration.
Results and discussion
The analysis of the potential energy
In this section, the potential energy between the molecules (i.e., nanocar, nanotruck, and C60) and the substrates have been evaluated. Considering the molecules and the substrate as rigid objects, the potential energy is calculated using the 6–12 Lennard–Jones potential. Additionally, because the results of this section are compared to those of molecular dynamics simulations, the initial configuration and the parameters of the Lennard–Jones potential are considered similar to those employed in the molecular dynamics analysis.
First, the intermolecular interactions between molecules and a pure boron-nitride and graphene substrate are measured separately, and their variation in Z-direction (Perpendicular to the substrate) is demonstrated in Fig. 2a–c. As depicted in Fig. 2a–c, the nanocars and fullerene experience more attraction on the pure boron-nitride substrate than on the pure graphene substrate. As a result, we can employ a substrate consisting of a mixture of both in some particular patterns to limit the molecular movement on the substrate. Moreover, the interaction between the nanocar and the substrate is stronger than that of the nanotruck and the substrate. The interaction between the nanotruck and substrate is stronger than the interaction between the C60 and the substrate, which arise from the fact that the nanocar is composed of more atoms than that of the nanotruck, and C60.
In the next step, to further survey the variation of potential energy, we investigated the potential energy in the case when objects are transferred from the pure boron-nitride substrate to the pure graphene substrate. The schematic representation of the setup considered for this evaluation is depicted in Fig. 2d. Each object is located on the boron-nitride side at first, and then by 0.5 Å steps, it shifts in the X-direction, and its minimum potential energy in Z-direction is computed. It is worth mentioning that, the vertical distance corresponding to the minimum potential energy is computed using the line search algorithm, and the minimum potential energy at the equilibrium point is recorded. The result of the preceding explained process for all objects is illustrated in Fig. 2e.
In Fig. 2e, the potential energy decay at the middle of the diagrams shows the transition of nanocars and C60 over the boundary, which depends on the structure of the molecules. In case of smallest molecule (i.e., C60), we observe smaller variation of potential energy. On the other hand, the nanocar has the largest descent of the potential energy, and since its wheels pass the border separately, the variation of potential energy does not have a constant slope. The difference between the potential energies of molecules, when they are placed on the boron-nitride and graphene substrate is the energy barrier against the transition. The highest and lowest energy barriers are corresponding to the transitions of nanocar and C60, respectively.
In the next step, eight types of local impurity are considered to evaluate the potential energy of C60 when it moves on these regions. Toward this aim, four circle-like impurities are considered in the middle of the graphene substrate, as shown in Fig. 3a. The impurity regions include 1, 3, 7, and 19 cells (hexagons) of boron-nitride, including 6, 13, 24, and 54 atoms of boron and nitrogen, respectively. Contrary to the former circumstance, the boron-nitride substrate with graphene impurities is also studied, whose schematic figures are depicted in Fig. 3b. The C60 was displaced horizontally in the X–Y plane with 0.1 Å steps, and the minimum potential energy in every point was calculated with respect to the vertical position of C60. It should be mentioned that, at each point, the hexagon face of the fullerene molecule is parallel to the surface. The final result is presented in Fig. 3c,d to clarify the effect of local impurity on the variation of potential energy.
The variation of the potential energy of C60 reveals that, increasing the size of boron-nitride impurity in the graphene substrate increases the depth of the potential energy well on the boron-nitride regions. As a consequence, the C60 requires higher thermal energy to overcome the energy barrier and move out of the impurity. On the other hand, as it can be seen in Fig. 3d, increasing the size of graphene impurity in the boron-nitride layer leads to a higher energy barrier in the graphene regions, which reduces the possibility of C60 entering this region. Also, it is presumable that C60 is unstable on the larger graphene impurity. In addition, the energy level in the center of 19 cell impurity is close to the state where the C60 is on the pure substrate. In other words, when C60 is located at the center of 19 cell impurity of both boron-nitride and graphene impurity, it almost does not interact with the atoms outside the impurity.
Molecular dynamics of C60 on the local impurities
Using molecular dynamics simulations, the trajectories of the motions of C60 on the substrates with the impurities are obtained (Fig. 3b). Each simulation is conducted in the temperature range of 100–700 K, with 100 K intervals. A 50 Å \(\times\) 50 Å layer was considered as substrates, which includes the impurity at the center. After relaxation process, the atoms which are located at the edges of the substrate were considered fixed to prevent excessive distortion of the substrate. The trajectories of the motion of C60 on graphene substrate with 4 types of boron-nitride impurity are presented in Fig. 4, during the 15 ns simulations. Here, only the trajectories of the simulations at the temperatures of 100, 300, 500, and 700 K are plotted in Fig. 4a, while complete trajectories are available in the Supplementary material. It is noteworthy to mention that, the following trajectories are obtained from the position of the C60 center of mass (COM) relative to the substrate center of mass.
According to Fig. 4, C60 motion is completely restricted in all four impurities at 100 K, except for a short period of time in 1 cell impurity. At 300 K, the larger impurity limited the C60 motion on the substrate. In 500 K, the motion is not as limited as in lower temperatures except for 19 cells of boron-nitride impurity, where a noticeable limitation can be seen. However, at 700 K, the C60 receives sufficient thermal energy to overcome the energy barrier, and its trajectories do not show any restriction imposed by the impurity. The constraint movement of C60 induced by boron-nitride impurity is intensified by the size of the impurity and diminishes with the simulation temperature.
For a quantitative evaluation of the C60 motion on the graphene substrate, the simulation time corresponding to the motion of C60 on the graphene substrate was estimated, and by dividing it by the total simulation time, we obtained the portion of the simulation time in which the C60 was on the graphene substrate. Figure 4b approves that by increasing the temperature and decreasing the impurity size, the C60 is more likely to overcome the energy barrier of the impurity.
In the opposite state of the previous situation, Fig. S3 (see Supplementary material) exhibits trajectories of a C60 on the boron-nitride substrate with four types of graphene impurities at the temperature range of 100–700 K. As Fig. S3 illustrates, all four kinds of graphene impurity prevent the C60 to enter the impurity area, at the temperature of 100 K. At the temperature of 500 K, only large impurities noticeably prevent the C60 to enter the impurity regions (i.e., 7 and 19 cells). Eventually, the C60 can overcome the energy barrier and move on the graphene parts more frequently at the temperature of 700 K.
The intermolecular interactions between the C60 and substrate are studied during the molecular dynamics simulations. According to the potential energy between C60 and substrates (Fig. 5), the vdW interactions fluctuate around a constant level of energy when C60 remains on the boron-nitride impurity. In case of 1 cell impurity, the C60 is able to overcome the energy barrier for a short period of time, at 100 K. As a result, the potential energy level increases when the fullerene leave the boron-nitride impurity, and the intermolecular interactions decreased. It can be seen in Fig. 5 that, by enlarging the impurity at a constant temperature, the average potential energy of the C60 decreases, which shows the stronger attraction of C60 and the larger BN impurities. This observation was mentioned earlier in the potential energy analysis. By raising the temperature in a particular impurity size, the range of the energy fluctuation increases based on the higher thermal fluctuations of the fullerene at higher temperatures.
Another point that can be deduced from Fig. 5 is that, desorption occurs at the temperature of 700 K. In other words, the thermal fluctuations of C60 increase to 0 eV, which means that C60 has a long distance from the substrate at 700 K. Since the distance of fullerene to the substrate becomes more than the cut-off radius, the interaction between the C60 and the substrate becomes trivial. According to the result of the simulations at 700 K, it seems that the impurity size has no effect on the desorption of C60 on the impurity, which can explain the irregular change of percentages observed in the last row of Fig. 4b.
The desorption temperature of C60 on the graphene substrate with 1, 3, 7, and 19 cells of boron-nitride impurity is also studied. Based on Fig. 2c, when the distance of C60 COM from the substrate becomes more than 14 Å, the LJ potential energy becomes zero. Therefore, the 14 Å vertical distance can be considered as a threshold at which the fullerene and the substrate no longer interact with each other. Accordingly, desorption occurs when the C60 distance to the substrate reaches this threshold and never interacts again during the simulation. To report the desorption temperature in this investigation, we increased the temperature of the simulations by 25 K intervals to obtain the lowest temperature that the C60 desorbs. The results of the mentioned process for the C60 on the boron-nitride impurities are presented in Table 1. Although there is no regular relationship between the impurity size and the desorption temperature, the desorption temperature corresponding to the 19-cell impurity has a slightly higher value indicating that higher thermal energy is required for the desorption of fullerene due to the higher attraction of 19-cell impurity.
Steering the surface motion by designing a nanoroad
In this section, we aim to restrict the motion of nanocars and C60 in a unidirectional path by using the difference of the potential energy of the molecules on boron-nitride and graphene substrates. Therefore, both sides of boron-nitride substrate are bonded to the strips of graphene to limit the lateral movement of the molecules. The schematic representation of the initial conditions of the C60, nanocar and nanotruck on this particular type of substrate (i.e., nanoroads) are depicted in Fig. 6. The width of the boron-nitride part of the substrate is chosen proportional to the width of nanovehicles and fullerene. More precisely, the width of the boron-nitride part of the substrate is about 8.5 Å larger than the width of nanovehicles and C60.
As demonstrated in the simulation setup, the length of the nanoroad is 280 Å, and contrary to the Y-direction, the boundary condition of the nanoroad is periodic in the X-direction. The temperature of the simulations was adjusted to 300 K. The simulations were performed for 15 ns, and the nanovehicles and fullerene were allowed to move freely on the substrate. The trajectory of nanocar, nanotruck, and fullerene during 15 ns simulations are shown in Fig. 6d–f. The blue strips of the substrate in Fig. 6 indicate the graphene substrate, and the pink part in the middle indicates the boron-nitride nanoroads.
As expected, C60 and nanocars demonstrate limited motion in the Y-direction, and0 none of the molecules entered the graphene region. Although fullerene and nanovehicles oscillate in Y-direction, they experience directional movements in X-direction. As a result, using a nanoroad, a directional motion is available for molecules on the surface. Nonetheless, the molecules still have diffusive motion in the Y-direction, especially C60, which has less interaction with the substrate due to the fewer atoms of molecule.
Unidirectional motion by applying temperature gradient on the nanoroads
In this section, to decrease the fluctuations of molecules and reach a more unidirectional motion, the effect of temperature gradient on the surface is surveyed. Since the free energy of the system decreases in the physical phenomena59, the nanocars and fullerene are expected to move from the beginning of the substrate with higher temperature to the colder end of the substrate.
In order to create models with an appropriate uniform temperature gradient, the following setup is implemented. First of all, a similar substrate with a length of 500 Å is considered, as employed in the last section. The temperature of the left side and the right side of the substrate are considered 600 and 300 K, respectively. Then, the substrate is allowed to relax with an NVT ensemble at 450 K for 100 ps. Afterwards, as shown in Fig. 7, a 10 Å strips of the beginning and the end of the substrate is fixed to prevent large displacements of the substrate. Next, the NVT ensembles are employed to apply the temperature gradient, as demonstrated in Fig. 7. To obtain a regular thermal distribution, five NVT ensembles are applied to five groups of substrate atoms at equal distances. Also, the NVE ensembles are applied to other atoms between NVT groups. Finally, the system was relaxed for 1 ns to find the desired temperature gradient on the surface.
The C60 and nanovehicles are placed at the left part of the substrate. The simulation continues until the COM of molecules reaches the right edge of the substrate with a temperature of 300 K. The trajectory of the motion of molecules is exhibited in Fig. 8 during the main simulations.
As we observe in Fig. 8, the fluctuation in the Y-direction is dramatically reduced for C60, nanocar and nanotruck compared to the simulations of previous section, in which the temperature was constant throughout the substrate. Because the nanocars and C60 tend to move toward the cold region on the right side of the substrate, most of their movements are in the X-direction.
In the following, a quantitative comparison is provided to better comprehend the influence of the temperature gradient on the surface motion of molecules. The average speed of the nanocar, nanotruck, and fullerene in the X-direction on the nanoroad with uniform temperature are 0.0314 Å/ps, 0.0263 Å/ps, and 0.0794 Å/ps, respectively. On the other hand, the average speed of these molecules in the X-direction on the nanoroad with the temperature gradient are 1.6734 Å/ps, 1.4982 Å/ps, and 0.8237 Å/ps, respectively. Accordingly, it can be seen that the average speed of the molecules in the X-direction has increased significantly in the presence of the temperature gradient compared to the motion at constant temperature. This is because the molecules spend less time oscillating in the Y-direction in the presence of temperature gradient. Correspondingly, nanocar, nanotruck and C60 tend to move in one direction, which increases the distance traveled in the same period of time (speed) in comparison with the motion of molecules at constant temperature. Although C60 has a higher average speed in the X-direction when subjected to a constant temperature compared to nanovehicles, its average speed in the X-direction is lower than the nanovehicles when it is subjected to the thermal gradient. This is due to the fluctuations of the C60 in the Y-direction; thus, the fullerene requires more time to pass the road.
To better investigate the directionality of the motion in both conditions of constant temperature and temperature gradient, the horizontal velocity of the molecules COM is calculated in each time step. The positive direction of the X-axis is considered as a reference, and the angle between the horizontal velocity vector and reference is obtained in each time step. For a particular angle of \({\theta }_{i}\), the magnitude of the velocity where the angle between the horizontal velocity vector and reference is equal to \({\theta }_{i}\) is accumulated. It should be noted that the angles between the velocity vector and the reference vector are between 0 and 180 degrees. Using the abovementioned process, the distribution of the horizontal velocities of molecules are obtained (Fig. 8b,c). According to Fig. 8b,c, unidirectional motion is not perceptible when the substrate is subjected to a constant temperature, and the velocity is distributed at different angles. Nevertheless, when a thermal gradient is applied to the substrate, the velocity vector is clearly directional, showing that the motion is toward the right side of the substrate where the temperature diminishes.
In the next step, we defined the unidirectional motion index (UMI) to check the unidirectionality of the motion quantitatively. UMI is defined as follows:
where \({V}_{acc}\) is the accumulated velocity vector at different angles. This parameter finds values between 0 and 1. When the horizontal velocity is uniformly distributed in different angles, the index is equal to 0, and when the velocity is limited to a smaller range of angles, the UMI index will be closer to 1. The obtained UMIs of C60 and nanovehicles are given in Table S5 of Supplementary material for both substrates with different thermal conditions. According to the UMI values provided in Table S5, the molecules have a higher UMI when the substrate is subjected to a temperature gradient, and their motions are more unidirectional compared to the motions on the nanoroads with a uniform temperature.
Conclusion
To provide the controllable motion of molecules on the surface, we evaluated the motion of C60 and C60-based nanocar and nanotruck on graphene/boron-nitride hybrid surfaces, using MD simulations and the potential energy analysis. It has been shown that C60 and nanovehicles are more stable on the boron-nitride impurity in a graphene monolayer unless enough thermal energy is provided to prevail the energy barrier. On the other hand, the molecules are not stable on a graphene impurity in a boron-nitride substrate. Considering different sizes of impurity reveals that at a constant temperature, larger boron-nitride impurities are more capable of restricting the motion.
Considering the energy barrier between two monolayers, we have tried to restrict the motion using a nanoroad of boron-nitride with graphene sideways. As a result of the higher attraction of molecules with boron-nitride regions, the motion of C60 and nanovehicles was limited to the boron-nitride pathways; however, a diffusive motion was still dominant, especially in the lateral direction.
To further control the motion, a thermal gradient was applied to the nanoroads. The surface molecules showed unidirectional motion towards the lower temperature part of the substrate to reach lower free energy. The rectilinear motion of C60 and nanovehicles induced by hybrid surfaces exposed to the temperature gradient benefits targeted nano-transportation on the surface.
Data availability
The computational protocols used in this work, the input and output files, are available upon request from corresponding author.
References
Otero, R. et al. Lock-and-key effect in the surface diffusion of large organic molecules probed by STM. Nat. Mater. 3, 779–782 (2004).
Hou, Q.-W., Cao, B.-Y. & Guo, Z.-Y. Thermal gradient induced actuation in double-walled carbon nanotubes. Nanotechnology 20, 495503 (2009).
van Venrooy, A., García-López, V., Li, J. T., Tour, J. M. & Dubrovskiy, A. V. Nanocars with permanent dipoles: preparing for the second international nanocar race. J. Org. Chem. 85, 13644–13654 (2020).
Mofidi, S. M., Nejat Pishkenari, H., Ejtehadi, M. R. & Akimov, A. V. Locomotion of the C60-based nanomachines on Graphene surfaces. Sci. Rep. 11, 2576 (2021).
Jin, T., García-López, V., Chen, F., Tour, J. & Wang, G. Imaging single molecular machines attached with two BODIPY dyes at the air–solid interface: High probability of single-step-like photobleaching and nonscaling intensity. J. Phys. Chem. C 120, 26522–26531 (2016).
Lin, X., Han, Q. & Huang, J. Effect of defects on the motion of carbon nanotube thermal actuator. Nanotechnol. Rev. 8, 79–89 (2019).
Akimov, A. & Kolomeisky, A. B. Dynamics of single-molecule rotations on surfaces that depend on symmetry, interactions, and molecular sizes. J. Phys. Chem. C 115, 125–131 (2011).
Jin, T. et al. Diffusion of nanocars on an air-glass interface. J. Phys. Chem. C 122, 19025–19036 (2018).
Vaezi, M., Nejat Pishkenari, H. & Ejtehadi, M. R. Collective movement and thermal stability of fullerene clusters on the graphene layer. Phys. Chem. Chem. Phys. 24, 11770–11781 (2022).
Khoei, A. R. & Kianezhad, M. A machine learning-based atomistic-continuum multiscale technique for modeling the mechanical behavior of Ni3Al. Int. J. Mech. Sci. 239, 107858 (2023).
Khoei, A. R., Youzi, M. & Eshlaghi, G. T. Mechanical properties and γ/γ’ interfacial misfit network evolution: A study towards the creep behavior of Ni-based single crystal superalloys. Mech. Mater. 171, 104368 (2022).
Lavasani, S. M. H., Pishkenari, H. N. & Meghdari, A. A closer look at the motion of p-carborane on gold surface. In 2016 International Conference on Manipulation, Automation and Robotics at Small Scales (MARSS) 1–6 (2016). https://doi.org/10.1109/MARSS.2016.7561702.
Pishkenari, H. N., Nemati, A., Meghdari, A. & Sohrabpour, S. A close look at the motion of C60 on gold. Curr. Appl. Phys. 15, 1402–1411 (2015).
Nemati, A. & Meghdari, A. Directing the diffusive motion of fullerene-based nanocars using nonplanar gold surfaces †. Phys. Chem. Chem. Phys. 20, 332–344 (2017).
Rad, A. S. & Aghaei, S. M. Potential of metal–fullerene hybrids as strong nanocarriers for cytosine and guanine nucleobases: A detailed DFT study. Curr. Appl. Phys. 18, 133–140 (2018).
Becton, M. & Wang, X. Thermal gradients on graphene to drive nanoflake motion. J. Chem. Theory Comput. 10, 722–730 (2014).
Akimov, A. V. & Kolomeisky, A. B. Unidirectional rolling motion of nanocars induced by electric field. J. Phys. Chem. C 116, 22595–22601 (2012).
Grill, L. et al. Rolling a single molecular wheel at the atomic scale. Nat. Nanotechnol. 2, 95–98 (2007).
Nemati, A., Pishkenari, H. N., Meghdari, A. & Ge, S. S. Influence of vacancies and grain boundaries on diffusive motion of surface rolling molecules. J. Phys. Chem. C https://doi.org/10.1021/acs.jpcc.0c03697 (2020).
Nemati, A., Nejat Pishkenari, H., Meghdari, A. & Sohrabpour, S. Directing the diffusive motion of fullerene-based nanocars using nonplanar gold surfaces. Phys. Chem. Chem. Phys. 20, 332–344 (2017).
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, 26018–26030 (2019).
Nemati, A., Nejat Pishkenari, H., Meghdari, A. & Ge, S. S. Directional control of surface rolling molecules exploiting non-uniform heat-induced substrates. Phys. Chem. Chem. Phys. 22, 26887–26900 (2020).
Akimov, A. V., Nemukhin, A. V., Moskovsky, A. A., Kolomeisky, A. B. & Tour, J. M. Molecular dynamics of surface-moving thermally driven nanocars. J. Chem. Theory Comput. 4, 652–656 (2008).
Lohrasebi, A., Neek-Amal, M. & Ejtehadi, M. R. Directed motion of C60 on a graphene sheet subjected to a temperature gradient. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 83, 3–6 (2011).
Kianezhad, M., Youzi, M., Vaezi, M. & Nejat, H. International journal of mechanical sciences rectilinear motion of carbon nanotube on gold surface. Int. J. Mech. Sci. 217, 107026 (2022).
Mofidi, S. M., Nejat Pishkenari, H., Ejtehadi, M. R. & Akimov, A. V. Role of graphene surface ripples and thermal vibrations in molecular dynamics of C60. J. Phys. Chem. C 123, 20026–20036 (2019).
Li, Y. et al. Differences in self-assembly of spherical C60 and planar PTCDA on rippled graphene surfaces. Carbon N. Y. 145, 549–555 (2019).
Bao, W. et al. Controlled ripple texturing of suspended graphene and ultrathin graphite membranes. Nat. Nanotechnol. 4, 562–566 (2009).
Li, Z., Liu, Z., Sun, H. & Gao, C. Superstructured assembly of nanocarbons: Fullerenes, nanotubes, and graphene. Chem. Rev. 115, 7046–7117 (2015).
Ozmaian, M., Fathizadeh, A., Jalalvand, M., Ejtehadi, M. R. & Allaei, S. M. V. Diffusion and self-assembly of C60 molecules on monolayer graphyne sheets. Sci. Rep. 6, 21910 (2016).
Li, L. H. et al. Large-scale mechanical peeling of boron nitride nanosheets by low-energy ball milling. J. Mater. Chem. 21, 11862–11866 (2011).
Chen, D. M., Shenai, P. M. & Zhao, Y. Tight binding description on the band gap opening of pyrene-dispersed graphene. Phys. Chem. Chem. Phys. 13, 1515–1520 (2011).
Vaezi, M., Nejat Pishkenari, H. & Nemati, A. Mechanism of C60 rotation and translation on hexagonal boron-nitride monolayer. J. Chem. Phys. 153, 234702 (2020).
Vaezi, M., Nejat Pishkenari, H. & Nemati, A. Mechanism of the motion of nanovehicles on hexagonal boron-nitride: A molecular dynamics study. Comput. Mater. Sci. 207, 111317 (2022).
Simpson, G. J., García-López, V., Daniel Boese, A., Tour, J. M. & Grill, L. How to control single-molecule rotation. Nat. Commun. 10, 4631 (2019).
Besenbacher, F., Lægsgaard, E. & Stensgaard, I. Fast-scanning STM studies. Mater. Today 8, 26–30 (2005).
Vaezi, M., Nejat Pishkenari, H. & Ejtehadi, M. R. Nanocar swarm movement on graphene surfaces. Phys. Chem. Chem. Phys. 24, 27759–27771 (2022).
Sasaki, T., Morin, J. F., Lu, M. & Tour, J. M. Synthesis of a single-molecule nanotruck. Tetrahedron Lett. 48, 5817–5820 (2007).
Vives, G. & Tour, J. M. Synthesis of single-molecule nanocars. Acc. Chem. Res. 42, 473–487 (2009).
Cuberes, M. T., Schlittler, R. R. & Gimzewski, J. K. Room-temperature repositioning of individual C60 molecules at Cu steps: Operation of a molecular counting device. Appl. Phys. Lett. 69, 3016–3018 (1996).
Gaisch, R. et al. Internal structure of C60 on Au(110) as observed by low-temperature scanning tunneling microscopy. J. Vac. Sci. Technol. B Microelectron. Nanom. Struct. Process. Meas. Phenom. 12, 2153–2155 (1994).
Shirai, Y. et al. Surface-rolling molecules. J. Am. Chem. Soc. 128, 4854–4864 (2006).
Tersoff, J. Modeling solid-state chemistry: Interatomic potentials for multicomponent systems. Phys. Rev. B 39, 5566–5568 (1989).
Allinger, N. L., Yuh, Y. H. & Lii, J. H. Molecular mechanics. The MM3 force field for hydrocarbons. J. Am. Chem. Soc. 111, 8551–8566 (1989).
Rappé, A. K. & Casewit, C. J. Molecular Mechanics Across Chemistry (University Science Books, 1997).
Watkins, E. K. & Jorgensen, W. L. Perfluoroalkanes: Conformational analysis and liquid-state properties from ab initio and Monte Carlo calculations. J. Phys. Chem. A 105, 4118–4125 (2001).
Allinger, N. L. Conformational analysis. 130. MM2. A hydrocarbon force field utilizing V1 and V2 torsional terms. J. Am. Chem. Soc. 99, 8127–8134 (1977).
Mofidi, S., Torkashvand, A. M., Moeini, A., Pazira, E. & Ahmadi, H. Relationship between fractal dimensions of stream and morphometric characteristics of basin for the soil conservation from water erosion. Int. J. River Basin Manage. 5, 1–11 (2021).
Nemati, A., Meghdari, A., Nejat Pishkenari, H. & Sohrabpour, S. Investigation into thermally activated migration of fullerene-based nanocars. Sci. Iran. 25, 1835–1848 (2018).
Wang, X., Ramírez-Hinestrosa, S., Dobnikar, J. & Frenkel, D. The Lennard–Jones potential: When (not) to use it. Phys. Chem. Chem. Phys. 22, 10624–10633 (2020).
Shayeganfar, F. & Shahsavari, R. Oxygen- and lithium-doped hybrid boron-nitride/carbon networks for hydrogen storage. Langmuir 32, 13313–13321 (2016).
Leven, I., Azuri, I., Kronik, L. & Hod, O. Inter-layer potential for hexagonal boron nitride. J. Chem. Phys. 140, 104106 (2014).
Lee, J. H. A study on a boron-nitride nanotube as a gigahertz oscillator. J. Korean Phys. Soc. 49, 172–176 (2006).
Plimpton, S. & Hendrickson, B. Parallel molecular dynamics algorithms for simulation of molecular systems. ACS Symp. Ser. 592, 114–132 (1995).
Qian, H., Sheetz, M. P. & Elson, E. L. Single particle tracking. Analysis of diffusion and flow in two-dimensional systems. Biophys. J. 60, 910–921 (1991).
Greiner, W., Neise, L., Stocker, H., Stöcker, H. & Rischke, D. Thermodynamics and Statistical Mechanics (Springer, 1995).
Attard, P. Thermodynamics and Statistical Mechanics: Equilibrium by Entropy Maximisation (Elsevier Science, 2002).
Ernst, D. & Köhler, J. Measuring a diffusion coefficient by single-particle tracking: Statistical analysis of experimental mean squared displacement curves. Phys. Chem. Chem. Phys. 15, 845–849 (2013).
Lohrasebi, A., Neek-Amal, M. & Ejtehadi, M. R. Directed motion of C60 on a graphene sheet subjected to a temperature gradient. Phys. Rev. E- Stat. Nonlinear Soft Matter Phys. 83, 1–10 (2011).
Author information
Authors and Affiliations
Contributions
M.K.: performing simulations, writing manuscript, proposing the idea, preparing figures, and providing discussion. M.Y.: writing manuscript, preparing figures, and performing simulations. M.V.: editing the manuscript, data analysis, providing discussions. H.N.P.: proposing the idea of the project, editing the manuscript and provided the discussions.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/.
About this article
Cite this article
Kianezhad, M., Youzi, M., Vaezi, M. et al. Unidirectional motion of C60-based nanovehicles using hybrid substrates with temperature gradient. Sci Rep 13, 1100 (2023). https://doi.org/10.1038/s41598-023-28245-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-28245-4
This article is cited by
-
Toward steering the motion of surface rolling molecular machines by straining graphene substrate
Scientific Reports (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.