وبلاگ صنایع شیمیایی 
Dynamic friction force in a carbon peapod oscillator Haibin Su^{1,2, William A. Goddard III1}, Yang Zhao3 1)Materials and Process Simulation Center, California Institute of Technology, Pasadena, CA 91125, U.S.A. 2)Division of Materials Science, Nanyang Technological University, Singapore 3)Department of Chemistry, University of Hong Kong, Pokfulam Road, Hong Kong, China Abstract We investigate a new generation of fullerene nanooscillators: a singlewalled carbon nanotube with one buckyball inside with an operating frequency in the tensofgigahertz range. A quantitative characterization of energy dissipation channels in the peapod pair has been performed via molecular dynamics simulation. Edge effects are found to the dominant cause of dynamic friction in the carbonpeapod oscillators. A comparative study on energy dissipation also reveals significant impact of temperature and impulse velocity on the frictional force. 1 Nanoscale fabrication technologies have made pervasive impact in the past 20 years^{1. One of the manifestations has been in the area of nanoelectromechanicalsystems (NEMS)2, which broadly refers to the application of nanofabrication technologies to construct sensors, actuators, and nanoscale integrated systems for a variety of applications. Very recently, Zettl’s group reported that frictional forces are very small, c.a. in the magnitude of 1014 N per إ2, during the controlled and reversible telescopic extension of multiwall carbon nanotubes3. Furthermore, it has been proposed that the transit time for complete nanotube core retraction (on the order of 1 to 10 ns) implies the possibility of exceptionally fast electromechanical switches4. In fact, oscillating crystals have a long history, dating back to 1880 when the piezoelectric effect was discovered by the Curie brothers5. Quartz crystal oscillators are widely used to provide regular pulses to synchronize various parts of an electronic system. But a typical crystal is millimeters in size, which could not be directly integrated into a computer chip. Motivated by the observation of Zettl’s group’s, Zheng and Jiang4 proposed a new type of nanooscillators operating completely differently from conventional quartz oscillators. Since then, designing this type of nanooscillator has been carried out actively. Legoas and collaborators6 first simulated an 38GHz nanooscillator consisting of a (9,0) carbon nanotube (CNT) inside of an (18,0) CNT. Zhao }et al 7 found that offaxial rocking motion of the inner nanotube and wavy deformation of the outer nanotube are responsible for energy dissipation in a doublewalled nanotube oscillator. So far, no successful experimental realization of the bitube oscillators has been reported. This is probably due to the considerable amount of energy dissipation, and the difficulty of preparing bitube type oscillator unit from multiwall carbon nanotubes with high quality. Fortunately, we have effective ways to place buckyballs inside nanotubes. For instance, singlewalled carbon nanotube (SWNT) can be synthesized by pulsedlaser vaporization route, whereby the sublimation of solid C_{60 in the presence of open SWNT causes the fullerenes to enter SWNT and selfassemble into 1D chains8. Therefore, it is feasible for single C60 to enter nanotube by van der Walls interactions. The peapod formation process has been widely studied912. Many studies are also focused on the effect of the nanotube diameter on binding properties9, 10. Filling SWNT with C60 is } 2 exothermic or endothermic depending on the size of the nanotube. C_{60@(10,10) is found to be stable (exothermic) while other peapods with smaller radius such as the (9,9) and (8,8) tubes are endoethermic10. Among many proposed interesting applications of peapod structures, one recent nanomechanical resonance study has been performed on C60filled carbon nanotube bundles towards developing next generation resonating systems13. The topic of this article is another feasible application of carbonpeapod system. Comparing with bitube structure models, here we propose that it is much more practical to replace the inner tube with a buckyball since the dynamic friction force between two objects is indeed proportional to the area of the overlapping sections from the perspective of modern tribology. Similar carbonpeapod osillators were investigated previously14, 15, focusing on elastic properties14 and lengthdependent oscillating frequencies15. It was also found that the frictional behavior of the carbon peapod depends on the diameter and chirality of the nanotube15. However, little attention has been paid so far to detailed energy dissipation mechanisms in the carbonpeapod oscillators. It is the aim of this article to investigate energy dissipation channels and effects of temperature and impulse velocity in the peapod oscillators via molecular dynamics simulation. } The structure model consists of one (10,10) SWNT of length 50.05إ with one C60 molecule inside. The edges of SWNT are passivated with hydrogen atoms. The diameter of C60 is 6.83إ, which fits nicely into (10,10) tube with a diameter 13.56إ (see Fig. 1). The force field used here for sp2 carbon centers was developed by fitting experimental lattice parameters, elastic constants and phonon frequencies for graphite16. This force field uses LennardJones 126 van der Waals interactions (Rv = 3.8050, Dv = 0.0692), Morse bond stretches (Rb = 1.4114, kb = 720, Db = 133.0), cosine angle bends (θa = 120, kθθ = 196.13, krθ =72.41, krr = 68), and a twofold torsion (Vt = 21.28), where all distances are in إ, angles in degrees, energies and force constants in kcal/mol. This force field correctly predicts that the energetically favorable packing for C60 is facecentered cubic and for C70 is hexagonal closepacking (h.c.p.). A more interesting result relevant to this work is the good agreement of the sublimation energy between calculation (40.9 kcal/mol at 739 K) and measurement (40.1 + 1.3 kcal/mol at 739 K)16. This force field is employed here to study the binding energy Eb 3 60)10,10()10,10@(60cEEcEbE−−= (Eq. 1). It is expected that buckyball molecules are more attracted to a SWNT than to each other due to a larger contact area with SWNTs, and therefore, more carboncarbon van der Waals interactions. Once a fullerene enters a nanotube, the van der Waals attraction keeps it inside. From our forcefield parameters, the binding energy is – 81.4 kcal/mol for putting one C_{60 inside a (10,10) SWNT, which is consistent with previous results reported by the Girifalco group17 and Ulbricht et al11. However, the reported binding energies by Okada et al10 (c.a. –11.8 kcal/mol) and the Louie group12 (c.a. –23.1 kcal/mol) are different from the above results, which is not surprising as it is wellknown that longrange attractive (London dispersion) interactions are not adequately described by the DFT methods18 based on the local density approximation and the generalized gradient approximation.. } As a preparation, our system is first equilibrated at 120K, 180K, 240K, 300K, and 360K with the NVT dynamics for 30 ps. After an impulse velocity is added to buckyball, dynamics simulation is carried out with the thermostat only attached to the tube. It is clear that the total energy of the C60@(10,10) system is not conserved. It is desirable to study the energy dissipation channels which allow energy flow from buckyball to tube, and then to the thermostat. The tube serves as a potential well to confine the motion of the ball. It is easy to estimate the maximallyallowed impulse velocity, which is 960 m/s, using the relation602maxcbmEv=. In our simulation, the C60 molecule is given an initial translational velocity of 480 m/s. Various impulse velocities between 100 m/s and 700 m/s also have been applied for T = 300K in order to study the relation between the dynamic friction force and the initial velocity. Kinetic energy data of the buckyball from the molecular dynamics simulation are smoothened by averaging over every five time periods (see Fig. 2a). As the buckyball moves to the opening end of the tube, it is slowed down until the translational velocity vanishes thanks to the van der Waals attraction 4 between tube and ball. At this moment, the potential energy reaches its local maximum, after which the ball returns back to the tube (see Fig. 2b and 2c). The period is 20 ps, corresponding to a frequency of 50GHz, which is quite encouraging for its potential applications to the nanofabrication field. In this system, there are two important channels for energy dissipation: offaxial wavy motion and the edge effect. The former means that the ball moves offaxially inside the tube; the latter refers to the case when the ball moves to the edges of the tube. Macroscopic models of friction between solids dictate that friction is proportional to normal force, independent of contact area. This is socalled Amontons’s law. For dynamic friction force, Coulomb’s law states that it is independent of velocity. It is interesting to note that recent experiments^{1921 by atomic force microscope (AFM)22 suggest that microscopic friction does not always behave according to traditional Amontons’s law and Coulomb’s law, thereby suggesting the need for new laws that account for atomic scale phenomena23, 24. So far, both area and velocity dependences have been demonstrated. It should be emphasized that edge effects should be paid more attention at nanoscale. This will be illustrated further in this Letter. } The dynamic friction force can be readily computed from the energy data lEfΔΔ−=, (Eq. 2) where Δl is the distance the buckyball travel while its kinetic energy is decreased by an amount ΔE. Comparing with previous studies on doublewalled oscillators, here only the tube is attached to the thermostat, a setting arguably more convenient to study energy dissipation in nanotube oscillators. Under this circumstance, the simulation yields the upper bound of dynamic friction force. It is interesting to compare the dynamic friction force at 300K evaluated in this study with those from the literature. In Zettl’s paper, the dynamic friction force per area is estimated to be less than 4.3 x 1015 N/Å2. In our study, this dynamic friction force is 0.17 pN, which is far less than that in bitubelike system (usually in the magnitude of nN). Note that the contact area between ball and tube is 97.9 5 إ^{2. This means the dynamic friction force per unit area of our system is 1.8 x 1015 N/إ2 which contains contributions from both the offaxial motion and the edge effects. } We now address the temperature effect on dynamic friction. The results, plotted in Fig. 3a and 3b, show that the dynamic friction force increases as the raising temperature. The source of this dynamic friction is the energy flow from the ball to the tube. This is caused by both the offaxial motion between the ball and the tube, and the edge effects. We define the offaxial angle, θ, of buckyball with respect to the tube as ))((cos21211rrrrabsvvvv×•=−θ, (Eq. 3). where 1rv is the vector from origin to the center of ball, 2rvis the axial vector from origin to the right edge of tube, abs means taking absolute value so that the value of θ is in the first quadrant (see Fig. 3c). No constraint is applied on the shape of the SWNT during the simulation. Indeed, there exists an energy transfer channel via the coupling between the offaxial wavy motion of the buckyball and the radial breathing mode of the SWNT. At the atomic scale, the shape of SWNT cannot be perfectly rigid due to thermal vibrations. Especially, when the buckyball moves near the edges of the SWNT, the axial vector of SWNT tilts away by about 0.4o from xaxis (assuming the SWNT initially is aligned along the xaxis). Thus, this variation should be taken into account. In practice, we use the mass center of the right edge of tube as the reference point to compute the axial vector, and then calculate the offaxial angle based on Eq. 3 for each ps during the simulation. The average offaxial angles tabulated in Table 1 are analyzed following this procedure. In Table (1.a), the average offaxial angles for different temperatures with a fixed impulse velocity (v = 480m/s) are collected. There is a clear trend that θ increases monotonically with the temperature, indicating that the ball takes on larger offaxial motion at higher temperatures which leads to higher energy dissipation. In addition, the energy transfer between the ball and the tube at the tube edges is also enhanced. Therefore, the dynamic friction force increases (see Fig. 3a). In this case, the contributions to dynamic friction from two dissipation channels are mixed together. However, as we stated before, when the buckyball moves near the edges of the SWNT (see Fig. 3d), there is significant energy transfer between them. Therefore, we must study 6 quantitatively energy dissipation due to the edge effect. To do this, we have examined the relation between the dynamic friction force and the impulse velocity in this system by varying impulse velocities (between 100 m/s and 700 m/s) at T = 300K. The results are plotted in Fig. 3b. Indeed, the velocity dependence of dynamic friction has been reported by Gnecco et al^{21 for a silicon tip sliding on a NaCl(100) surface. In our simulation, the velocity is in the order of 102 m/s, which is ten orders of magnitude higher than those in Gnecco et al’s paper21. Thus, we are in two completely different velocity regimes. In addition, the friction force calculated here arises from smooth sliding between a bukcyball and the inner walls of a SWNT without applied normal forces, while in Gnecco et al’s paper, strong forces are applied between the tip and the NaCl (100) surface. Consequently, a direct comparison can not be made. However, good agreement has been found between this work and another earlier computational study7 on the value of frictional forces per atom. In our model it is interesting to note that the average offaxial angles are little influenced by impulses for a given temperature, for example, 300K (see Table (1b)). Since the frictional force increases with the impulse velocity as shown in Fig. 3b, it follows that this velocity dependence is mainly due to energy dissipation at the edges of the tube (not due to offaxial motion), and the larger impulse velocity, the more energy dissipation at the edges of the tube. More importantly, we can extract the dynamic friction force due to offaxial wavy motion at 300K by simply extrapolating the data to zero impulse velocity. It yields a friction force of about 12 fN due to the offaxial wavy motion. For comparison, the contribution from the edge effects (158 fN) with an impulse velocity of 480m/s is about one order of magnitude larger than wavy motion at 300K with the same impulse velocity. This is supported by a recent independent study by Tangney et al25. } To conclude, we have carried out a quantitative characterization of dominant energy dissipation mechanisms for a carbon peapod nanooscillator that can be realized in the lab. Our molecular dynamics simulation results reveal significant effects of the temperature and the impulse velocity on friction. In particular, it has been shown that the edge effects are the main cause of the dynamic friction force. This novel nanooscillator 7 design proposed here holds great promise for applications in NEMS owing to its extremely low operating friction and easy adaption for impulse generation^{26. } Acknowledgement H.B.S. is grateful for kind assistances in coding from Y. Lansac, J. Dodson, and P. Meulbroek. The work at NTU is funded by NTUCoESUG under Grant No. M58070001. The facilities of the Materials and Process Simulation Center (MSC) used in these studies are funded by DURIP (ARO and ONR), NSF (CTS and MRI), and a SUR Grant from IBM. In addition, the MSC is funded by grants from AROMURI, NIH, ChevronTexaco, General Motors, SeikoEpson, the Beckman Institute, Asahi Kasei, and Toray Corp. 8 Reference: 1 K. E. Drexler, Nanosystems: Molecular Machinery, Manufacturing and Computation (Wiley, New York, 1992). 2 M. Roukes, Scientific American 285, 48 (2001); M. Roukes, Phys. World 14, 25 (2001). 3 J. Cumings and A. Zettl, Science 289, 602 (2000). 4 Q. S. Zheng and Q. Jiang, Physical Review Letters 88, 045503 (2002). 5 W. G. Cady, Piezoelectricity (McGrawHill Book Co. Inc., New York, 1946). 6 S. B. Legoas, V. R. Coluci, S. F. Braga, et al., Physical Review Letters 90, 055504 (2003). 7 Y. Zhao, C. C. Ma, G. H. Chen, et al., Physical Review Letters 91, 175504 (2003); C.C. Ma, Y. Zhao, C.Y. Yam, G.H. Chen, and Q. Jiang, Nanotechnology 16, 1253 (2005). 8 B. W. Smith and D. E. Luzzi, Chemical Physics Letters 321, 169 (2000). 9 M. Hodak and L. A. Girifalco, Physical Review B 68, 085405 (2003). 10 S. Okada, S. Saito, and A. Oshiyama, Physical Review Letters 86, 3835 (2001). 11 H. Ulbricht, G. Moos, and T. Hertel, Physical Review Letters 90, 095501 (2003). 12 Y. G. Yoon, M. S. C. Mazzoni, and S. G. Louie, Applied Physics Letters 83, 5217 (2003). 13 P. Jaroenapibal, S. B. Chikkannanavar, D. E. Luzzi, et al., Journal of Applied Physics 98 (2005). 14 D. Qian, W. K. Liu, and R. S. Ruoff, Journal of Physical Chemistry B 105, 10753 (2001). 15 P. Liu, Y. W. Zhang, and C. Lu, Journal of Applied Physics 97 (2005). 16 Y. J. Guo, N. Karasawa, and W. A. Goddard, Nature 351, 464 (1991). 17 M. Hodak and L. A. Girifalco, Chemical Physics Letters 363, 93 (2002). 18 W. Kohn, Y. Meir, and D. E. Makarov, Physical Review Letters 80, 4153 (1998). 19 B. N. J. Persson, Surface Science Reports 33, 85 (1999). 20 N. J. Brewer, B. D. Beake, and G. J. Leggett, Langmuir 17, 1970 (2001). 21 E. Gnecco, R. Bennewitz, T. Gyalog, et al., Physical Review Letters 84, 1172 (2000). 22 R. Wiesendanger, Scanning Probe Microscopy and Spectroscopy: Methods and Applications (Cambridge University Press, London, 1994). 23 D. Dowson, History of Tribology (American Soc. Mech. Engrs, New York, 1999). 24 J. Krim, Scientific American 275, 74 (1996); J. Krim, Am. J. Phys. 70, 890 (2002). 25 P. Tangney, S. G. Louie, and M. L. Cohen, Physical Review Letters 93, 065503 (2004). 26 Our proposed peapod oscillator can be modified by replacing the neutral C_{60 with NC59+. This heterofullerene ion is formed in the gap phase during fast atom bombardment of a clusteropened Nmethoxyethoxy methyl ketolactam27. It is possible to use same approach to put NC59+ inside SWNT as assembling C60 with SWNT. Then applying Efield can generate proper amount of momentum for NC59+ as initial impulse. For practical consideration, NC59+ can be stabilized in gas phase with nitrogen atmosphere. } 9 10 Table Table 1. Temperature and velocity effects on offaxial angles. In (1.a) the impulse is fixed as 480 m/s for all the temperatures. In (1.b) the temperature is kept to be 300 K while varying impulses. The offaxial angles are averaged for 80 periods during the oscillation. See text for detailed discussions.


+ نوشته شده در
ساعت توسط مجتبی گل بهار 

صفحه نخست پروفایل مدیر وبلاگ پست الکترونیک آرشیو وبلاگ عناوین مطالب وبلاگ 
درباره وبلاگ 
شیمی تجزیه

پیوندهای روزانه 
دانشگاه فنی انقلاب اسلامی آرشیو پیوندهای روزانه 
پیوندها 
**عکس مجانی پرینت کنید** 
RSS
