The radii and lengths of the modeled armchair CNTs
1. Introduction
Since Iijima reported MWCNTs in 1991, CNTs have captured the intensive attention of researchers worldwide due to the combination of their expected structural perfection, small size, low density, high stiffness, high strength, and excellent electronic properties. CNTs have been widely adopted as microscopic probing tips (Dai et al., 1996; Hafner et al., 2001), nanocomposites reinforcements (Bower et al., 1998; Jin et al., 1998), nanotweezers (Kim & Lieber, 1999), and nanoactuators (Baughman et al., 1999; Fennimore et al., 2003) due to their slender and high aspect ratio structures. Meanwhile, nanotubes are also highly susceptible to buckling under compression, which is a structural instability. Once the buckling of CNTs occurs, the load-carrying capability would suddenly reduce and lead to possible catastrophic failure of the nanotubes, which significantly limit the loading strengths of the probing tips and compressive strengths of nanocomposite structures. Even the physical properties such as conductance of carbon nanotube can be influenced by the occurrence of buckling (Postma et al., 2001). Hence, it is crucial to understand the mechanism of nanotube buckling and even predict the onset of buckling in order to improve the nanotube applications.
A review of the relevant literature shows that significant studies have employed both experimental (Falvo et al., 1997; Iijima et al., 1996; Thostenson & Chou, 2004; Waters et al., 2004) and theoretical (Ru, 2000; Yakobson & Avouris, 2001) approaches to investigate the bucking behaviors of CNTs. However due to the difficulties encountered at nanoscale, the experimental investigation of the buckling behaviors of CNTs remains a challenging problem and individual factors that affect buckling could not be easily identified. In theoretical study, the CNTs are commonly treated as beams or thin-shell tubes with certain wall thickness and elastic constants and, thus, it is difficult to consider the chirality and size effects on buckling behavior of CNTs because the continuum assumption disregards the discrete nature of atomic structures (Ru, 2000; Yakobson & Avouris, 2001). Some researchers attempted to introduce the atomic-continuum method combining the atomic detail in the continuum description and examine the various properties of CNTs (Chang, 2004; Guo et al., 2008; Li & Chou, 2003a, 2003b). The atomic-continuum method could shorten the computational time in larger atomic system.
As the fast development and rapid advancement of computers, molecular approaches have become important tools and are widely applied to study the factors that would influence the buckling of CNTs (Buehler et al., 2004; Cao & Chen, 2006a, 2006b; Huh & Huh, 2008; Liew et al., 2004; Ozaki et al., 2000). Although some researchers already discussed various aspects of the CNT buckling behavior, systematic analysis on the effect of geometry (i.e. radius, length and length-to-radius ratio) and chirality (i.e. armchair, zigzag and chiral) on the buckling mechanism is still lacking to the best of our knowledge. Consequently, the present study employs MD simulations based on the Tersoff many-body potential function (Tersoff, 1986, 1988, 1989) to perform a systematic and comprehensive investigation into the buckling behaviors of single-walled CNTs under uniaxial compressive displacement loading. Besides, the applicability of the continuum buckling theory, which has been well developed for thin tubes, on predicting the buckling behavior of the CNT will also be examined.
2. Methodology
Atomic model of the single-walled CNT with radius
During the simulations, the nanotube is compressed uniaxially and incrementally along
CNTs can be considered as a graphene sheet rolled into a cylindrical shell and three distinct types of nanotubes could be classified as armchair, zigzag and chiral according to the way of rolling-up. In order to assess the influence of geometry and chirality on the buckling behavior of the single-walled CNTs, various radii, lengths, length-to-radius ratios and helical types of CNTs would be simulated and elucidate the dependence of buckling strain. The length-to-radius ratio (or called slenderness ratio) is defined as
where
( |
(5,5) | (8,8) | (10,10) |
Radius(Å) | 3.39 | 5.424 | 6.78 |
Length(Å) / |
7.38 / 2.18 12.30 / 3.63 36.89 / 10.88 61.49 / 18.14 86.08 / 25.39 147.57 / 43.53 |
12.3 / 2.27 19.68 / 3.63 59.03 / 10.88 98.38 / 18.14 137.73 / 25.39 236.11 / 43.53 |
14.76 / 2.18 24.595 / 3.63 73.79 / 10.88 122.98 / 18.14 172.165 / 25.39 295.14 / 43.53 |
( |
(9,0) | (14,0) | (17,0) |
Radius(Å) | 3.523 | 5.48 | 6.655 |
Length(Å) / |
8.52 / 2.42 12.78 / 3.63 38.34 / 10.88 63.9 / 18.14 89.46 / 25.39 153.36 / 43.53 |
12.78 / 2.33 21.3 / 3.89 59.64 / 10.88 97.98 / 17.88 140.58 / 25.65 238.56 / 43.53 |
17.04 / 2.56 25.56 / 3.84 72.42 / 10.88 119.28 / 17.92 170.4 / 25.60 289.68 / 43.53 |
( |
(6,4) | (9,6) | (12,8) |
Radius(Å) | 3.413 | 5.119 | 6.825 |
Length(Å) / |
18.569 / 5.44 37.138 / 10.88 55.707 / 16.32 92.845 / 27.20 148.55 / 43.52 |
18.569 / 3.62 55.707 / 10.88 92.845 / 18.14 129.982 / 25.39 222.83 / 43.53 |
18.569 / 2.72 74.276 / 10.88 129.982 / 18.92 167.712 / 24.57 297.10 / 43.53 |
3. Results
MD approaches are utilized to simulate the uniaxial compression test of CNTs with different geometries and chiralities. The size and chirality effects on buckling behaviors are studied and systematically compared through the critical compressive strains. It is noticed that two distinct types of buckling configurations are commonly observed as shown in Fig. 3 (a) and (b). One is shell wall buckling with kinks on the wall with the tube’s centerline remaining straight, which is a local instability and the other is column buckling, which buckle as a whole. Moreover, there is another intriguing type of initial buckling configuration, which starts with a bump on the wall as shown in Fig. 3(c), and the CNT will eventually turn into column or shell wall buckling. It is also observed that the buckling kinks or bends appear mostly near the middle of the nanotubes, which indicates the loading mechanism does not impose any extra constraint on the two ends. Unlike the velocity controlled loading (Jeng et al, 2004), the drastic deformation starts from the imminent places near the ends due to Poisson’s effect. Hence, it is difficult to separate the boundary constraint effect from the other investigating factors. It is also noted that the radial distribution function before and after buckling is quite similar as shown in Fig. 4, which implies that the bond structures do not change due to buckling.
The buckling strains for different chirality and slenderness ratio of CNTs are shown in Fig. 5, 6 and 7, respectively. The filled marks indicate column buckling and the empty one indicate shell wall buckling. In general, the nanotube’s buckling behavior transits from shell wall type for a short tube to column type for a long tube of the same radius irrespective to the chirality of the CNTs. For slender nanotubes (
4. Analysis and discussion
From current simulation results, it is very difficult to deduce a clear relationship between the critical compressive buckling strain and the CNT geometry and helical type, not to mention predicting the buckling type. In continuum mechanics, there already exist theories to calculate the buckling strain of thin shell tube depending on the buckling type in macroscopic scale. However, the applicability of continuum buckling theory in nanoscale is still an unsolved question. Hence, the continuum predictions on the buckling strains and the simulation results will be compared in order to examine the applicability.
By using Euler-Bernoulli beam theory, the critical buckling strain of a tube with both ends clamped displaying the column-like buckling behavior can be described as (Landau et al., 1986)
where
It is noticed that the critical compressive buckling strain decreases as the inverse square of the slenderness ratio for column-type buckling. On the other hand, the critical compressive buckling strain of the tube displaying the shell-like buckling behavior is (Libai & Simmonds, 1998; Timoshenko & Gere, 1961)
The Poisson’s ratio, ν, of the CNTs is chosen as 0.19 (Yakobson et al. 1996). It is noted that the critical compressive strain of shell-like buckling depends only on the tube thickness and radius, but is independent of the tube length.
Fig. 8, 9 and 10 illustrate the comparisons between the continuum predictions on the buckling strains and the simulation results for respective set of similar radii. One of the uncertainties of the continuum model is the effective nanotube thickness. Here in this research, the effective nanotube thickness is chosen as 0.66Å (Yakobson et al. 1996), which is roughly the atomic radius of carbon, instead of the commonly used CNT wall thickness, 3.4 Å. Since the nanotube thickness is not negligible as compared to the radius, Eq. (1) is adopted in the comparison. The dash lines illustrate the continuum prediction of shell-like buckling while the solid lines depict the prediction based on Euler-Bernoulli beam theory. It is observed that the buckling resistance does not show obvious dependence on the CNT chirality under similar radii. The continuum prediction of shell-like buckling can serve as an upper bound for predicting the shell wall buckling but significantly overestimate the buckling strains for nanotubes with higher length-to-radius ratio irrespective to the nanotube chirality. As the slenderness ratio becomes higher, the continuum prediction based on Euler-Bernoulli beam theory could capture the trend of column type buckling but also overestimate the buckling strains for nanotubes with smaller length-to-radius ratio. It is observed that those buckling strain which deviate significantly from the continuum prediction are at the intersection region of the two continuum theories. It is speculated that the difference in buckling strain between the continuum theories and simulation results at the transition region could be due to the competition between two buckling mechanisms.
5. Conclusion
In this work, molecular dynamics approach is employed to study the buckling behaviors of single-walled carbon nanotubes with different geometric sizes and chiralities under room temperature. Based on the MD simulation results, it is observed that the nanotube’s buckling behavior transits from shell wall type for a short tube to column type for a long tube of the same radius irrespective to the chirality of the CNTs. Moreover, the buckling strain is getting smaller as the CNT becomes slender for most nanotubes, which implies that the slender nanotubes have lower buckling resistance. Under similar length-to-radius ratio, it is noted that the buckling strain decreases as the radius of the CNT increases especially for CNTs with smaller slenderness ratios. From the comparison with the prediction made by continuum buckling theories, it is concluded that the corresponding buckling strain and buckling type predicted by the continuum theory could agree reasonably well with MD simulations of the CNTs under compression except at the transition region. From the findings of this paper, it is suggested that the continuum buckling theory with proper choice of parameters, i.e., wall thickness and Poisson’s ratio, could capture the trend of the buckling strain on the length-to-radius ratio disregarding to the helical types (i.e., armchair, zigzag and chiral) of the nanotubes and, hence, could serve as a primitive guideline in predicting the buckling strain of the CNTs.
Acknowledgments
This research work is supported by National Science Council of Taiwan under the grants NSC97-2221-E-194-015 and NSC98-2221-E-194-012-MY2. The support of AFOSR under Contract No. FA2386-09-1-4152 AOARD 094152 is also acknowledged. We are also grateful to the Taiwan National Center for High-performance Computing for computer time and facilities.
References
- 1.
Baughman R. H. Cui C. X. Zakhidov A. A. Iqbal Z. Barisci J. N. Spinks G. M. Wallace G. G. Mazzoldi A. De Rossi D. Rinzler A. G. Jaschinski O. Roth S. Kertesz M. 1999 Carbon Nanotube Actuators,284 5418 (May 1999),1340 1344 ,0036-8075 - 2.
Bower C. Rosen R. Jin L. Han J. Zhou O. 1999 Deformation of Carbon Nanotubes in Nanotube-Polymer Composites, ,74 22 (May 1999),3317 3319 ,0003-6951 - 3.
Buehler M. Kong Y. Gao H. 2004 Deformation Mechanisms of Very Long Single-Wall Carbon Nanotubes Subject to Compressive Loading,,126 3 (July 2004),245 249 ,0094-4289 - 4.
Cao G. Chen X. 2006a Mechanisms of Nanoindentation on Single-Walled Carbon Nanotubes: The Effect of Nanotube Length, ,21 4 (April 2006),1048 1070 ,0884-2914 - 5.
Cao G. Chen X. 2006b Buckling Behavior of Single-walled Carbon Nanotubes and a Targeted Molecular Mechanics Approach, ,73 16 (October 2006),155435 1098-0121 - 6.
Chang T. 2004 Buckling of Single-Walled Carbon Nanotubes via a Hybrid Atomic/Continuum Approach, ,36 6 (November 2004),744 748 , - 7.
Dai H. Hafner J. H. Rinzler A. G. Colbert D. T. Smalley R. E. 1996 Nanotubes as Nanoprobes in Scanning Probe Microscopy, ,384 6605 (November 1996),147 150 ,0028-0836 - 8.
Falvo M. R. Clary G. J. Taylor R. M. Chi V. Brooks F. P. Washburn S. Superfine R. 1997 Bending and Buckling of Carbon Nanotubes under Large Strain, ,389 6651 ( October 1997 ),582 584 ,0028-0836 - 9.
Fennimore A. M. Yuzvinsky T. D. Han W. Q. Fuhrer M. S. Cumings J. Zettl A. 2003 Rotational Actuators based on Carbon Nanotubes,424 6947 (July 2003),408 410 ,0028-0836 - 10.
Guo X. Leung A. Y. T. He X. Q. Jiang H. Huang Y. 2008 Bending Buckling of Single-Walled Carbon Nanotubes by Atomic-Scale Finite Element, ,39 1 202 208 ,1359-8368 - 11.
Hafner J. H. Cheung C. L. Woolley A. T. Lieber C. M. 2001 Structural and Functional Imaging with Carbon Nanotube AFM Probes, ,77 1 73 110 ,0079-6107 - 12.
Huh J. Huh H. 2008 Effect of Helicity on the Buckling Behavior of Single-Wall Carbon Nanotubes, ,22 31-32 , (December 2008),5872 5877 ,0217-9792 - 13.
Iijima S. Brabec C. Maiti A. Bernholc J. 1996 Structural Flexibility of Carbon Nanotubes, ,104 5 (February 1996),2089 2092 ,0021-9606 - 14.
Jeng Y. R. Tsai P. C. Fang T. H. 2004 Effects of Temperature and Vacancy Defects on Tensile Deformation of Single-Walled Carbon Nanotubes, ,65 11 (November 2004),1849 1856 ,0022-3697 - 15.
Jin L. Bower C. Zhou O. 1998 Alignment of Carbon Nanotubes in a Polymer Matrix by Mechanical Stretching, ,73 9 ( August 1998),1197 1199 ,0003-6951 - 16.
Kim P. Lieber C. M. 1999 Nanotube Nanotweezers, ,286 5447 (December 1999),2148 2150 ,0036-8075 - 17.
Landau L. D. Pitaevskii L. P. Lifshitz E. M. Kosevich A. M. 1986 , Butterworth-Heinemann,978-0-75062-633-0 Oxford, UK - 18.
Li C. Chou T. W. 2003a A Structural Mechanics Approach for the Analysis of Carbon Nanotubes, ,40 10 (May 2003),2487 2499 ,0020-7683 - 19.
Li C. Chou T. W. 2003b Single-Walled Carbon Nanotubes as Ultrahigh Frequency Nanomechanical Resonators, ,68 7 (August 2003),073405 1098-0121 - 20.
Libai A. Simmonds J. G. 2005 Shells, Cambridge University Press,978-0-52101-976-7 Cambridge, UK - 21.
Liew K. M. Wong C. H. He X. Q. Tan M. J. Meguid S. A. 2004 Nanomechanics of Single and Multiwalled Carbon Nanotubes, ,69 11 (March 2004),115429 1098-0121 - 22.
Ozaki T. Iwasa Y. Mitani T. 2000 Stiffness of Single-Walled Carbon Nanotubes under Large Strain, ,84 8 (February 2000),1712 1715 ,0031-9007 - 23.
Postma H. W. C. Teepen T. Yao Z. Grigoni M. Dekker C. 2001 Carbon Nanotube Single-Electron Transistors at Room Temperature, ,293 5527 (July 2001),76 79 ,0036-8075 - 24.
Ru C. Q. 2000 Effective Bending Stiffness of Carbon Nanotubes, ,62 15 (October 2000),9973 9976 ,0163-1829 - 25.
Tersoff J. 1986 New Empirical-Model for the Structural-Properties of Silicon, ,56 6 (February 1986),632 635 ,0031-9007 - 26.
Tersoff J. 1988 New Empirical-Approach for the Structure and Energy of Covalent Systems, ,37 12 (April 1988),6991 7000 ,0163-1829 - 27.
Tersoff J. 1989 Modeling Solid-State Chemistry-Interatomic Potentials for Multicomponent Systems, ,39 8 (March 1989 ),5566 5568 ,1098-0121 - 28.
Thostenson E. T. Chou T. W. 2004 Nanotube Buckling in Aligned Multi-Wall Carbon Nanotube Composites, ,42 14 3015 3018 ,0008-6223 - 29.
Timoshenko S. P. Gere J. M. 1961 , McGraw-Hill,978-0-48647-207-2 New York, USA - 30.
Waters J. F. Riester L. Jouzi M. Guduru P. R. Xu J. M. 2004 Buckling Instabilities in Multiwalled Carbon Nanotubes under Uniaxial Compression, ,85 10 (September 2004),1787 1789 ,0003-6951 - 31.
Yakobson B. I. Avouris P. 2001 Mechanical Properties of Carbon Nanotubes, ,80 287 327 ,0303-4216 - 32.
Yakobson B. I. Brabec C. J. Bernholc J. 1996 Nanomechanics of Carbon Tubes: Instabilities beyond Linear Response, ,76 14 (April 1996),2511 2514 ,0031-9007