搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

含缺陷碳纳米管及碳纳米豆荚静动力特性模拟研究

王磊 张冉冉 方炜

含缺陷碳纳米管及碳纳米豆荚静动力特性模拟研究

王磊, 张冉冉, 方炜
PDF
HTML
导出引用
  • 采用分子动力学方法, 对含双空位及多空位缺陷碳纳米管进行静动力特性模拟研究. 首先讨论了双原子空位缺陷以及多原子空位缺陷对碳纳米管的准静态力学性质的影响, 然后讨论了缺陷以及轴向预应力对碳纳米豆荚内C60分子振荡动力学的影响. 研究表明, 相对于无缺陷碳纳米管, 含不同类型双原子空位缺陷碳纳米管的极限应力、极限应变和弹性模量都大幅下降; 当碳纳米管缺陷原子较多, 缺陷连接在一起形成类似裂纹之后, 使得碳纳米管轴向抗压性能大幅降低, 裂纹沿周向发展相比于裂纹沿轴向发展, 其抗压能力下降得更多, 这类似于含裂纹的壳体模型结构抗压性能的下降; 缺陷碳纳米豆荚中C60分子的振荡频率受到缺失的碳原子数的影响, 单原子空位缺陷使得C60分子的振荡频率增大, 但随着空位数的增多, C60分子的振荡频率会逐渐减小; 当缺陷碳纳米豆荚存在轴向预应力时, C60分子的振荡不仅受到缺陷影响, 同时还受到轴向预应力的影响, 这使得C60分子振荡变得更为复杂.
      通信作者: 王磊, wangL@hhu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11472098)和中央高校基本科研业务费专项资金(批准号: 2018B57814)资助的课题.
    [1]

    Iijima S 1991 Nature 354 56

    [2]

    Krishnan A, Dujardin E, Ebbesen T W 1998 Phys. Rev. B 58 14013

    [3]

    Berber S, Kwon Y K, Tomanek D 2000 Phys. Rev. Lett. 84 4613

    [4]

    Rinzler A G, Hafner J H, Nikolaev P, Lou L, Kim S G, Tomanek D, Nordlander P, Colbert D T, Smalley R E 1995 Science 269 1550

    [5]

    Popov V N 2004 Mater. Sci. Eng. R 43 61

    [6]

    Harris P J F, Hernandez E, Yakobson B I 2004 Am. J. Phys. 72 413

    [7]

    Stone A J, Wales D J 1986 Chem. Phys. Lett. 128 501

    [8]

    Nardelli M B, Yakobson B I, Bernholc J 1998 Phys. Rev. B 57 4277

    [9]

    Hirai Y, Nishimaki S, Mori H, Kimoto Y, Akita S, Nakayama Y, Tanaka Y 2003 Jpn. J. Appl. Phys. 42 4120

    [10]

    辛浩, 韩强, 姚小虎 2008 物理学报 57 4391

    Xin H, Han Q, Yao X H 2008 Acta Phys. Sin. 57 4391

    [11]

    袁剑辉, 程玉民, 张振华 2009 物理学报 58 2578

    Yuan J H, Cheng Y M, Zhang Z H 2009 Acta Phys. Sin. 58 2578

    [12]

    Zhang Y Y, Xiang Y, Wang C M 2009 J. Appl. Phys. 106 113503

    [13]

    Kulathunga D D T K, Ang K K, Reddy J N 2010 J. Phys.: Condens. Matter 22 345301

    [14]

    王锋, 曾祥华, 徐秀莲 2002 物理学报 51 1778

    Wang F, Zeng X H, Xu X L 2002 Acta Phys. Sin. 51 1778

    [15]

    Liu P, Zhang Y W, Lu C 2005 J. Appl. Phys. 97 094313

    [16]

    Song H Y, Zha X W 2009 Phys. Lett. A 373 1058

    [17]

    崔柳, 冯妍卉, 檀鹏, 张欣欣 2015 科学通报 60 1414

    Cui L, Feng Y H, Tan P, Zhang X X 2015 Chin. Sci. Bull. 60 1414

    [18]

    方炜, 王磊 2018 材料导报 32 164

    Fang W, Wang L 2018 Mater. Rev. 32 164

    [19]

    Plimpton S 1995 J. Comput. Phys. 117 1

    [20]

    Brenner D W, Shenderova O A, Harrison J A, Stuart S J, Ni B, Sinnott S 2002 J. Phys.: Condens. Matter 14 783

    [21]

    Stuart S J, Tutein A B, Harrison J A 2000 J. Chem. Phys. 112 6472

    [22]

    Lennard-Jones J E, Dent B M 1926 Proc. R. Soc. A 112 230

    [23]

    Girifalco L A, Hodak M, Lee R S 2000 Phys. Rev. B 62 13104

    [24]

    Rafizadeh H A 1974 Physica 74 135

    [25]

    Cox B J, Thamwattana N, Hill J M 2007 Proc. R. Soc. A 463 461

    [26]

    Cox B J, Thamwattana N, Hill J M 2007 Proc. R. Soc. A 463 477

  • 图 1  (a) (10 10)无缺陷扶手椅型碳纳米管; (b) 六元环中碳原子编号

    Fig. 1.  (a) (10 10) Non-defective armchair carbon nanotube; (b) numbering of carbon atoms in six-membered ring.

    图 2  表1中编号1−6所对应的含双原子空位缺陷碳纳米管的局部原子构型图

    Fig. 2.  Local atomic configuration of carbon nanotubes containing diatomic vacancy defects corresponding to numbers 1−6 in Table 1

    图 3  含多原子空位缺陷碳纳米管模型 (a) 沿轴向分布的多原子空位缺陷; (b) 沿周向分布的多原子空位缺陷

    Fig. 3.  Carbon nanotubes with polyatomic vacancy defects: (a) Polyatomic vacancy defects distributed along the axial direction; (b) polyatomic vacancy defects distributed along the circumferential direction.

    图 4  无缺陷碳纳米豆荚模型

    Fig. 4.  Non-defective carbon nano-peapod.

    图 5  含缺陷碳纳米豆荚模型

    Fig. 5.  Defective carbon nano-peapod.

    图 6  含多原子空位缺陷碳纳米管在轴向荷载作用下屈曲失稳构型 (a) 沿轴向分布的含多原子空位缺陷管; (b) 沿周向分布的含多原子空位缺陷管

    Fig. 6.  Buckling instability configuration of carbon nanotubes with polyatomic vacancy defects under axial loading: (a) Carbon nanotubes with polyatomic vacancy defects distributed along the axial direction; (b) carbon nanotubes with polyatomic vacancy defects distributed along the circumferential direction.

    图 7  含多原子空位缺陷碳纳米管的弹性模量随原子缺失数的变化

    Fig. 7.  Variation of elastic modulus of carbon nanotubes with the number of missing atoms when polyatomic vacancy defects occurs.

    图 8  C60在一个运动周期内的受力随运动距离的变化

    Fig. 8.  The variation of force on C60 molecule with the distance of its motion in a cycle.

    图 9  C60分子的振荡频率随碳管中原子缺失数的变化

    Fig. 9.  Variation of the oscillation frequency of C60 molecule with the number of missing atoms in carbon tubes.

    图 10  缺陷碳纳米豆荚内C60分子的振荡频率随轴向预应力的变化

    Fig. 10.  Variation of the oscillation frequency of C60 molecules in the defective carbon nano-peapods with axial pre-stress.

    表 1  含双原子空位缺陷碳纳米管

    Table 1.  Carbon nanotubes containing diatomic vacancy defects.

    编号缺陷类型间隔碳原子数
    1缺失1, 2位置原子0
    2缺失1, 3位置原子1
    3缺失1, 4位置原子2
    4缺失1, 5位置原子1
    5缺失1, 6 位置原子0
    6缺失3, 6位置原子2
    下载: 导出CSV

    表 2  轴压下双原子空位缺陷对碳纳米管力学性能影响

    Table 2.  Effect of diatomic vacancy defects on mechanical properties of carbon nanotubes under axial compression.

    编号极限应力/GPa极限应变/%弹性模量/GPa
    127.963.76743.6
    228.264.00706.5
    328.284.00707.0
    429.064.24685.4
    531.214.20743.1
    630.344.30705.6
    740.864.85842.5
    下载: 导出CSV
  • [1]

    Iijima S 1991 Nature 354 56

    [2]

    Krishnan A, Dujardin E, Ebbesen T W 1998 Phys. Rev. B 58 14013

    [3]

    Berber S, Kwon Y K, Tomanek D 2000 Phys. Rev. Lett. 84 4613

    [4]

    Rinzler A G, Hafner J H, Nikolaev P, Lou L, Kim S G, Tomanek D, Nordlander P, Colbert D T, Smalley R E 1995 Science 269 1550

    [5]

    Popov V N 2004 Mater. Sci. Eng. R 43 61

    [6]

    Harris P J F, Hernandez E, Yakobson B I 2004 Am. J. Phys. 72 413

    [7]

    Stone A J, Wales D J 1986 Chem. Phys. Lett. 128 501

    [8]

    Nardelli M B, Yakobson B I, Bernholc J 1998 Phys. Rev. B 57 4277

    [9]

    Hirai Y, Nishimaki S, Mori H, Kimoto Y, Akita S, Nakayama Y, Tanaka Y 2003 Jpn. J. Appl. Phys. 42 4120

    [10]

    辛浩, 韩强, 姚小虎 2008 物理学报 57 4391

    Xin H, Han Q, Yao X H 2008 Acta Phys. Sin. 57 4391

    [11]

    袁剑辉, 程玉民, 张振华 2009 物理学报 58 2578

    Yuan J H, Cheng Y M, Zhang Z H 2009 Acta Phys. Sin. 58 2578

    [12]

    Zhang Y Y, Xiang Y, Wang C M 2009 J. Appl. Phys. 106 113503

    [13]

    Kulathunga D D T K, Ang K K, Reddy J N 2010 J. Phys.: Condens. Matter 22 345301

    [14]

    王锋, 曾祥华, 徐秀莲 2002 物理学报 51 1778

    Wang F, Zeng X H, Xu X L 2002 Acta Phys. Sin. 51 1778

    [15]

    Liu P, Zhang Y W, Lu C 2005 J. Appl. Phys. 97 094313

    [16]

    Song H Y, Zha X W 2009 Phys. Lett. A 373 1058

    [17]

    崔柳, 冯妍卉, 檀鹏, 张欣欣 2015 科学通报 60 1414

    Cui L, Feng Y H, Tan P, Zhang X X 2015 Chin. Sci. Bull. 60 1414

    [18]

    方炜, 王磊 2018 材料导报 32 164

    Fang W, Wang L 2018 Mater. Rev. 32 164

    [19]

    Plimpton S 1995 J. Comput. Phys. 117 1

    [20]

    Brenner D W, Shenderova O A, Harrison J A, Stuart S J, Ni B, Sinnott S 2002 J. Phys.: Condens. Matter 14 783

    [21]

    Stuart S J, Tutein A B, Harrison J A 2000 J. Chem. Phys. 112 6472

    [22]

    Lennard-Jones J E, Dent B M 1926 Proc. R. Soc. A 112 230

    [23]

    Girifalco L A, Hodak M, Lee R S 2000 Phys. Rev. B 62 13104

    [24]

    Rafizadeh H A 1974 Physica 74 135

    [25]

    Cox B J, Thamwattana N, Hill J M 2007 Proc. R. Soc. A 463 461

    [26]

    Cox B J, Thamwattana N, Hill J M 2007 Proc. R. Soc. A 463 477

  • [1] 张凯旺, 钟建新. 缺陷对单壁碳纳米管熔化与预熔化的影响. 物理学报, 2008, 57(6): 3679-3683. doi: 10.7498/aps.57.3679
    [2] 李威, 冯妍卉, 陈阳, 张欣欣. 碳纳米管中点缺陷对热导率影响的正交试验模拟分析. 物理学报, 2012, 61(13): 136102. doi: 10.7498/aps.61.136102
    [3] 辛 浩, 韩 强, 姚小虎. 单、双原子空位缺陷对扶手椅型单层碳纳米管屈曲性能的不同影响. 物理学报, 2008, 57(7): 4391-4396. doi: 10.7498/aps.57.4391
    [4] 李瑞, 孙丹海. 缺陷对碳纳米管摩擦与运动行为的影响. 物理学报, 2014, 63(5): 056101. doi: 10.7498/aps.63.056101
    [5] 曾永辉, 江五贵, Qin Qing-Hua. 螺旋上升对自激发锯齿型双壁碳纳米管振荡行为的影响. 物理学报, 2016, 65(14): 148802. doi: 10.7498/aps.65.148802
    [6] 金年庆, 滕玉永, 顾 斌, 曾祥华. 稀有气体原子注入缺陷性纳米碳管的分子动力学模拟. 物理学报, 2007, 56(3): 1494-1498. doi: 10.7498/aps.56.1494
    [7] 孟利军, 张凯旺, 钟建新. 硅纳米颗粒在碳纳米管表面生长的分子动力学模拟. 物理学报, 2007, 56(2): 1009-1013. doi: 10.7498/aps.56.1009
    [8] 李振武. 纳米CdS/碳纳米管复合材料的光电特性. 物理学报, 2012, 61(1): 016103. doi: 10.7498/aps.61.016103
    [9] 赵东林, 曾宪伟, 沈曾民. 碳纳米管/聚苯胺纳米复合管的制备及其微波介电特性研究. 物理学报, 2005, 54(8): 3878-3883. doi: 10.7498/aps.54.3878
    [10] 李 瑞, 胡元中, 王 慧, 张宇军. 单壁碳纳米管在石墨基底上运动的分子动力学模拟. 物理学报, 2006, 55(10): 5455-5459. doi: 10.7498/aps.55.5455
    [11] 徐葵, 王青松, 谭兵, 陈明璇, 缪灵, 江建军. 形变碳纳米管选择通过性的分子动力学研究. 物理学报, 2012, 61(9): 096101. doi: 10.7498/aps.61.096101
    [12] 杨成兵, 解辉, 刘朝. 锂离子进入碳纳米管端口速度的分子动力学模拟. 物理学报, 2014, 63(20): 200508. doi: 10.7498/aps.63.200508
    [13] 谢根全, 夏 平. 基于微极性弹性力学的碳纳米管中波的传播特性. 物理学报, 2007, 56(12): 7070-7077. doi: 10.7498/aps.56.7070
    [14] 廖庆亮, 齐俊杰, 黄运华, 邓战强, 高战军, 曹佳伟, 张 跃, 夏连胜. 丝网印刷制备碳纳米管阴极的强流脉冲发射特性研究. 物理学报, 2008, 57(4): 2328-2333. doi: 10.7498/aps.57.2328
    [15] 王畅, 王利光, 张鸿宇, Terence K. S. W.. 嵌入锂原子的zigzag型单壁碳纳米管的电子传导特性. 物理学报, 2010, 59(1): 536-540. doi: 10.7498/aps.59.536
    [16] 张忠强, 丁建宁, 刘珍, Y. Xue, 程广贵, 凌智勇. 碳纳米管-聚乙烯复合材料界面力学特性分析. 物理学报, 2012, 61(12): 126202. doi: 10.7498/aps.61.126202
    [17] 袁学松, 张宇, 孙利民, 黎晓云, 邓少芝, 许宁生, 鄢扬. 碳纳米管冷阴极脉冲发射特性及仿真模型研究. 物理学报, 2012, 61(21): 216101. doi: 10.7498/aps.61.216101
    [18] 麻华丽, 霍海波, 曾凡光, 向飞, 王淦平. 铜基碳纳米管薄膜的制备及其强流脉冲发射特性研究. 物理学报, 2013, 62(15): 158801. doi: 10.7498/aps.62.158801
    [19] 韩典荣, 王璐, 罗成林, 朱兴凤, 戴亚飞. (n, n)-(2n, 0)碳纳米管异质结的扭转力学特性. 物理学报, 2015, 64(10): 106102. doi: 10.7498/aps.64.106102
    [20] 蔡丹, 刘列, 巨金川, 王海涛, 赵雪龙, 王潇. 碳纳米管和天鹅绒阴极强流发射特性的对比研究. 物理学报, 2016, 65(4): 045202. doi: 10.7498/aps.65.045202
  • 引用本文:
    Citation:
计量
  • 文章访问数:  1132
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-23
  • 修回日期:  2019-06-10
  • 上网日期:  2019-08-01
  • 刊出日期:  2019-08-20

含缺陷碳纳米管及碳纳米豆荚静动力特性模拟研究

  • 河海大学力学与材料学院, 南京 211100
  • 通信作者: 王磊, wangL@hhu.edu.cn
    基金项目: 国家自然科学基金(批准号: 11472098)和中央高校基本科研业务费专项资金(批准号: 2018B57814)资助的课题.

摘要: 采用分子动力学方法, 对含双空位及多空位缺陷碳纳米管进行静动力特性模拟研究. 首先讨论了双原子空位缺陷以及多原子空位缺陷对碳纳米管的准静态力学性质的影响, 然后讨论了缺陷以及轴向预应力对碳纳米豆荚内C60分子振荡动力学的影响. 研究表明, 相对于无缺陷碳纳米管, 含不同类型双原子空位缺陷碳纳米管的极限应力、极限应变和弹性模量都大幅下降; 当碳纳米管缺陷原子较多, 缺陷连接在一起形成类似裂纹之后, 使得碳纳米管轴向抗压性能大幅降低, 裂纹沿周向发展相比于裂纹沿轴向发展, 其抗压能力下降得更多, 这类似于含裂纹的壳体模型结构抗压性能的下降; 缺陷碳纳米豆荚中C60分子的振荡频率受到缺失的碳原子数的影响, 单原子空位缺陷使得C60分子的振荡频率增大, 但随着空位数的增多, C60分子的振荡频率会逐渐减小; 当缺陷碳纳米豆荚存在轴向预应力时, C60分子的振荡不仅受到缺陷影响, 同时还受到轴向预应力的影响, 这使得C60分子振荡变得更为复杂.

English Abstract

    • 自1991年日本科学家饭岛发现碳纳米管以来[1], 其已受到广泛研究和关注. 碳纳米管具有高杨氏模量、高导电性、结构完善和化学稳定性等优异的性质[2,3], 近年来已被广泛地应用于复合材料、电子科技、生物工程、环境保护等诸多领域[4,5].

      在碳纳米管的实际制备过程中, 完美的碳纳米管总是难以得到, 利用各种方法制得的碳纳米管存在着各种缺陷[6], 如原子空位缺陷、五-七环缺陷, 以及其他非拓扑结构缺陷[79]等. 缺陷的存在难以避免, 同时缺陷的存在又严重影响了碳纳米管的力学、电学等性能, 故对含缺陷碳纳米管进行模拟研究显得尤为重要. 2008年, 辛浩等[10]对含单、双原子空位缺陷的扶手椅型单壁碳纳米管进行了轴向压缩模拟, 研究结果表明, 管壁缺陷大大降低了碳纳米管的承载能力. 2009年, 袁剑辉等[11]探讨了单原子空位缺陷比率以及两个单原子空位缺陷分布位置对碳纳米管弹性性能的影响. 同年, Zhang等[12]运用分子动力学方法, 研究了在轴向压缩作用下, 单原子空位缺陷以及双原子空位缺陷对不同手性的单壁碳纳米管稳定性的影响. 2010年, Kulathunga等[13]进一步研究了对称空位缺陷以及不对称空位缺陷对单壁碳纳米管屈曲性能的影响. 可见, 研究者们已经对含原子空位缺陷碳纳米管进行了广泛的研究, 但对同一个六元环缺失任意两个碳原子的双原子空位缺陷和多原子空位缺陷连接在一起形成类似裂纹的研究则较少.

      在碳纳米管内部的空腔中填充富勒烯C60, 将构成一种新型纳米材料-碳纳米豆荚. 2002年, 王锋等[14]利用分子动力学方法, 发现了C60分子进入单壁碳纳米管形成的碳纳米豆荚的吸入和俘获机理, 并且揭示了吸入和俘获势垒只存在于管口区, 而在管内区, C60分子沿轴向运动几乎不受力. 2005年, Liu等[15]运用分子动力学方法, 研究了C60分子在范德瓦耳斯力的作用下, 在碳管中做衰减振荡运动, 这种衰减行为受到碳管直径及螺旋角的影响. 2009年, Song和Zha[16]研究了单原子空位缺陷和双原子空位缺陷对碳纳米豆荚振荡行为的影响, 结果表明缺陷的存在使得碳纳米豆荚的振荡稳定性变差, 并且会造成碳纳米豆荚的能量损耗. 2015年, 崔柳等[17]研究了碳纳米豆荚内C60分子的振荡行为, 结果表明C60分子的振荡行为受环境温度和碳管管壁层数的影响, 而C60分子填充个数的增加将会导致振荡不再发生. 2018年, 方炜和王磊[18]研究了轴向预应力对碳纳米豆荚中C60振荡行为的影响, 研究结果表明C60分子的振荡频率随轴向拉伸预应力的增加而单调减小; C60分子的振荡频率在压缩预应力的作用下, 呈分段线性衰减模式. 尽管国内外学者已开始关注碳纳米豆荚的振荡行为, 但大多集中于无缺陷碳纳米豆荚的研究, 针对含缺陷碳纳米豆荚的研究仍很不充分. 由于制备的碳纳米管不可避免地含有缺陷, 其与C60构成的碳纳米豆荚也会存在缺陷. 因此, 对含缺陷碳纳米豆荚的研究非常有必要, 并且研究缺陷碳纳米豆荚在不同工况下的振荡行为, 将对高频振荡器的制造提供理论参考.

      本文采用分子动力学方法, 模拟了含缺陷扶手椅型碳纳米管在单轴压缩下的力学行为, 以及含缺陷碳纳米豆荚中C60分子的振荡行为. 重点讨论了双原子空位缺陷和多原子空位缺陷对扶手椅型碳纳米管力学性能的影响, 以及缺失原子数和轴向预应力对碳纳米豆荚中C60分子振荡行为的影响.

    • 建立若干正常模态的扶手椅型碳纳米管, 如图1(a)所示. 碳纳米管的手性为(10 10), 直径和长度分别为1.356 nm, 6 nm. 删去该碳纳米管管壁上的部分原子, 形成含缺陷碳纳米管.

      图  1  (a) (10 10)无缺陷扶手椅型碳纳米管; (b) 六元环中碳原子编号

      Figure 1.  (a) (10 10) Non-defective armchair carbon nanotube; (b) numbering of carbon atoms in six-membered ring.

      图1(b)为在图1(a)中黑色方框中的一个六元环. 本文定义当一个六元环中缺失两个碳原子时, 则称为双原子空位缺陷. 从该六元环的六个碳原子位置中任取两个碳原子位置形成空缺, 以确立双原子空位缺陷类型. 考虑到(10 10)扶手椅型碳纳米管具有高度的对称性, 在轴压状态下, 15种双原子空位缺陷中部分缺陷类型受力情况相同, 故缺陷类型简化为6种, 如表1所列. 表1中, 编号1—5的双原子空位缺陷类型以1号位为起点, 另一点为终点, 编号6的缺陷类型以3号位为起点, 6号位为终点. 表1中编号1—6所对应的含双原子空位缺陷碳纳米管的局部原子构型图如图2所示.

      编号缺陷类型间隔碳原子数
      1缺失1, 2位置原子0
      2缺失1, 3位置原子1
      3缺失1, 4位置原子2
      4缺失1, 5位置原子1
      5缺失1, 6 位置原子0
      6缺失3, 6位置原子2

      表 1  含双原子空位缺陷碳纳米管

      Table 1.  Carbon nanotubes containing diatomic vacancy defects.

      图  2  表1中编号1−6所对应的含双原子空位缺陷碳纳米管的局部原子构型图

      Figure 2.  Local atomic configuration of carbon nanotubes containing diatomic vacancy defects corresponding to numbers 1−6 in Table 1

      图1(a)所示的碳纳米管的中心附近区域, 分别沿轴向和沿周向方向去掉一定数目的碳原子, 这些空位缺陷连接在一起形成裂纹, 最终表现为沿轴向分布和沿周向分布的含多原子空位缺陷碳纳米管, 如图3所示.

      图  3  含多原子空位缺陷碳纳米管模型 (a) 沿轴向分布的多原子空位缺陷; (b) 沿周向分布的多原子空位缺陷

      Figure 3.  Carbon nanotubes with polyatomic vacancy defects: (a) Polyatomic vacancy defects distributed along the axial direction; (b) polyatomic vacancy defects distributed along the circumferential direction.

    • 本文所研究的无缺陷碳纳米豆荚是由图1(a)所示的扶手椅型碳纳米管和一个C60分子构成, 如图4所示. C60分子的直径为0.72 nm. 在该碳纳米豆荚的中心附近区域去掉一定数目的碳原子, 形成含缺陷碳纳米豆荚, 如图5所示.

      图  4  无缺陷碳纳米豆荚模型

      Figure 4.  Non-defective carbon nano-peapod.

      图  5  含缺陷碳纳米豆荚模型

      Figure 5.  Defective carbon nano-peapod.

    • 本文采用开源并行模拟软件LAMMPS (large-scale atomic/molecular massively parallel simulator)[19]进行分子动力学模拟计算; 选用AIREBO (adaptive intermolecular reactive empirical bond order)[20, 21]势能函数来表述碳纳米管和C60分子中的C—C键的相互作用. 在分子动力学模拟中碳-碳原子或碳-氢原子间的相互作用可以被该势能函数准确地描述, 其表达式如下:

      $E = \frac{1}{2}\sum\limits_i {\sum\limits_{i \ne j} {\left[ {E_{ij}^{{\rm{REBO}}} + E_{ij}^{{\rm{LJ}}} + \sum\limits_{k \ne i,j} {\sum\limits_{l \ne i,j,k} {E_{kijl}^{{\rm{TORSION}}}} } } \right]} } ,$

      其中$E_{ij}^{{\rm{REBO}}}$为REBO势能函数项, $E_{ij}^{{\rm{LJ}}}$为长程相互作用项, $E_{kijl}^{{\rm{TORSION}}}$为依赖于二面角的四体势扭转项.

      选用Lennard Jones (L-J)势能函数[22]来表述碳管与C60分子的碳原子间的长程非键作用, 即范德瓦耳斯力, 其表达式如下:

      $U({r_{ij}}) = 4\varepsilon \left[ {{{\left( {\frac{\sigma }{{{r_{ij}}}}} \right)}^{12}} - {{\left( {\frac{\sigma }{{{r_{ij}}}}} \right)}^6}} \right], $

      其中, ${r_{ij}}$是相邻的成键原子i与原子j之间的间距; $\varepsilon $$\sigma $为势能参数, 能量参数$\varepsilon $ = 0.00239 eV[23], 长度参数$\sigma = 0.3345\;{\rm{nm}}$[24].

      在微正则系综(NVE)下进行整个模拟, 时间步长为0.001 ps. 模拟过程中, 先对缺陷碳纳米管或缺陷碳纳米豆荚进行无约束弛豫, 以保证其结构的能量最小化, 达到初始平衡态. 然后, 在研究双原子空位缺陷和多原子空位缺陷对碳纳米管力学行为的影响时, 在弛豫过的缺陷碳纳米管的底部和顶部各定义一部分边界原子. 底部边界中的碳原子在轴向方向固定, 其他方向不加约束. 对顶部边界中的碳原子施加轴向预应力, 加载方式为位移加载. 在研究轴向预应力对缺陷碳纳米豆荚中C60分子的振荡行为的影响时, 先对缺陷碳纳米豆荚中的缺陷碳纳米管进行相同的加载, 加载过程中, C60分子固定不动. 当碳纳米管加载完成后, 把C60分子从碳纳米管的底部管口处自由释放.

    • 采用分子动力学方法, 模拟在轴压下双原子空位缺陷对扶手椅型碳纳米管极限应力、极限应变以及弹性模量的影响, 并与无缺陷扶手椅型碳纳米管进行对比分析. 这里我们通过对轴压作用下缺陷碳纳米管应力应变曲线的分析得到缺陷碳纳米管的弹性模量. 在缺陷碳纳米管应力应变曲线呈线性变化阶段, 其应力与应变的比值可以看作是缺陷碳纳米管的弹性模量.

      表1中编号1—6的碳纳米管模型进行轴向压缩模拟, 结果如表2中编号1—6所示. 表2中编号7为图1(a)所示的无缺陷扶手椅型碳纳米管在轴压作用下模拟的结果. 从表2可以看出, 在轴压下, 双原子空位缺陷的出现使得碳纳米管的极限应力下降了24%—32%不等, 使得极限应变下降了11%—22%不等, 使得弹性模量下降了12%—19%不等. 空位的出现使得碳纳米管原本完美的六元环结构遭到破坏, 在承受轴向荷载时缺陷处会产生应力集中现象, 从而导致了碳纳米管在抗轴向压缩方面的能力的下降.

      编号极限应力/GPa极限应变/%弹性模量/GPa
      127.963.76743.6
      228.264.00706.5
      328.284.00707.0
      429.064.24685.4
      531.214.20743.1
      630.344.30705.6
      740.864.85842.5

      表 2  轴压下双原子空位缺陷对碳纳米管力学性能影响

      Table 2.  Effect of diatomic vacancy defects on mechanical properties of carbon nanotubes under axial compression.

      表2可以看出编号1和5缺陷碳纳米管的弹性模量明显大于编号2, 3, 4, 6的弹性模量. 一方面, 编号1和5缺陷碳纳米管中缺失的两个碳原子是相邻的, 产生了4个悬挂键; 编号2, 4中缺失的两个碳原子间相隔了1个碳原子, 产生了6个悬挂键; 编号3和6中缺失的两个碳原子之间相隔了2个碳原子, 也产生了6个悬挂键. 悬挂键的产生将引起碳纳米管结构的松弛, 从而导致碳纳米管弹性模量的下降. 另一方面, 从图2可以看出编号1和5缺陷碳纳米管的缺陷面积为4个六元环的面积之和, 而编号2, 3, 4, 6缺陷碳纳米管的缺陷面积为5个六元环的面积之和. 而缺陷面积的增大势必会使得碳纳米管弹性模量减小.

    • 多原子空位缺陷类型主要分为沿轴向分布和沿周向分布两种类型, 图3(a)为沿轴向分布的含多原子空位缺陷管, 图3(b)为沿周向分布的含多原子空位缺陷管. 通过改变缺失碳原子的个数来控制碳纳米管沿轴向分布或沿周向分布的裂纹长度, 从而探究多原子空位缺陷对碳纳米管弹性模量的影响.

      图6给出了沿轴向分布和沿周向分布的含多原子空位缺陷碳纳米管在轴向荷载作用下屈曲失稳的构型. 图7为含多原子空位缺陷碳纳米管的弹性模量随原子缺失数的变化曲线. 从图7可以看出, 无论缺失原子是沿轴向分布还是沿周向分布, 碳纳米管的弹性模量都随着缺失原子数的增加而减小. 当碳纳米管在沿轴向方向或沿周向方向缺失3个碳原子时, 其弹性模量较无缺陷碳纳米管的弹性模量分别下降了14%和17%, 当缺失原子数达到7个时, 其弹性模量分别下降了28%和44%. 可见当缺失相同的原子个数时, 沿周向分布比沿轴向分布对碳纳米管的弹性模量影响要大得多. 随着缺失原子数的增多(从3—7), 沿轴向分布的碳管的弹性模量下降了14%, 而沿周向分布的却下降了26%. 碳纳米管在轴向荷载作用下屈曲失稳时, 会在碳管轴线中心附近产生凹陷变形, 凹陷处出现应力集中现象, 让原结构遭到破坏, 使得碳管难以承受轴向压力而失稳. 而凹陷变形总是沿着周向, 所以对比二者的弹性模量, 在缺失相同原子数目时, 缺陷沿周向分布的碳管的弹性模量总是低于缺陷沿轴向分布的碳管, 且随着缺失原子数目的增加, 缺陷沿周向分布碳管的弹性模量的下降速率也大于缺陷沿轴向分布的碳管.

      图  6  含多原子空位缺陷碳纳米管在轴向荷载作用下屈曲失稳构型 (a) 沿轴向分布的含多原子空位缺陷管; (b) 沿周向分布的含多原子空位缺陷管

      Figure 6.  Buckling instability configuration of carbon nanotubes with polyatomic vacancy defects under axial loading: (a) Carbon nanotubes with polyatomic vacancy defects distributed along the axial direction; (b) carbon nanotubes with polyatomic vacancy defects distributed along the circumferential direction.

      图  7  含多原子空位缺陷碳纳米管的弹性模量随原子缺失数的变化

      Figure 7.  Variation of elastic modulus of carbon nanotubes with the number of missing atoms when polyatomic vacancy defects occurs.

    • 首先采用分子动力学方法, 模拟如图4所示的单根完好的碳纳米豆荚内C60分子的振荡. 结果表明长6 nm的(10 10)扶手椅型碳纳米豆荚内C60分子的振荡频率为32.73 GHz, 该结果与Cox等[25,26]利用机械原理建立方程求得的36.13 GHz的频率结果较为接近.

      以碳纳米豆荚为基础的高频振荡装置, 可能会因为各种原因出现管壁缺陷的情况, 而管壁的缺陷势必会对碳纳米豆荚内C60分子的振荡产生一定的影响. 本文建立了若干个中部缺失不同原子数的碳纳米豆荚模型, 并对其进行分子动力学模拟, 从而探究缺失原子个数对C60分子振荡的影响.

      图8给出了图4所示的无缺陷碳纳米豆荚和图5所示的含十二原子空位缺陷碳纳米豆荚中C60分子在一个运动周期(C60从碳管的底部运动到顶部, 再由顶部运动到底部)的受力情况, 取竖直向上为正. 从图8可以看出, 当C60分子在无缺陷管的底部由静止状态释放后, 会受到一个竖直向上的力, 这是由于受到管口碳原子的范德瓦耳斯力, C60分子开始由静止做加速运动. 大约在距管口六分之一位置时, C60分子所受的合力约为零, 这是由于C60分子所受的范德瓦耳斯力与滑动摩擦力几乎平衡, 此时C60分子速度达到最大, 并几乎保持不变, 开始做匀速运动. 当C60分子运动到碳管的另一管口时, C60分子受到竖直向下的作用力, 该力同样是管口碳原子的范德瓦耳斯力, 但方向相反, 将阻止C60分子离开碳管, C60分子开始做减速运动, 直至速度减为零, 随后又做反向加速运动, 返回到碳管中, 最终在管内做周期性振荡. 而缺陷的引入使得碳纳米管中C60分子在缺陷附近的范德瓦耳斯力和摩檫力都发生了改变. 当C60分子进入缺陷位置时, 因为碳管管壁缺失, 其所受摩檫力和范德瓦耳斯力均减小, 原本的平衡状态遭到破坏. 但相较于摩擦力来说, 范德瓦耳斯力减小得多, 此时, 摩擦力大于范德瓦耳斯力, C60开始做减速运动. 当C60分子离开缺陷位置时, 因为碳管缺陷消失, 其所受摩擦力和范德瓦耳斯力均增大. 但范德瓦耳斯力增大得多, 此时, 范德瓦耳斯力大于摩檫力, C60做加速运动. 离开缺陷位置后, C60经过短暂的速度调整, 又以某一速度匀速运动. 在C60分子通过缺陷的整个过程中, 其振荡频率主要受摩擦力和范德瓦耳斯力二者的共同影响, 碳纳米管与C60分子之间的摩擦力的减小使得C60分子的振荡频率增加, 而碳纳米管与C60分子之间的范德瓦耳斯力的减小使得C60分子的振荡频率降低.

      图  8  C60在一个运动周期内的受力随运动距离的变化

      Figure 8.  The variation of force on C60 molecule with the distance of its motion in a cycle.

      图9给出了C60分子的振荡频率随碳管原子缺失数的变化曲线. 可见, 缺陷对碳纳米豆荚中C60分子的振荡有很大影响, 随着碳管缺失原子数的增多, C60分子的振荡频率先增大后减小. 当缺失原子数为1个时, 范德瓦耳斯力对C60分子振荡频率的影响小于摩擦力对C60分子的影响, C60分子的振荡频率明显高于无缺陷管的. 但随着缺失原子数的增加, 范德瓦耳斯力对C60分子的振荡频率的影响越来越大, 其振荡频率也基本呈线性减小. 当缺失原子数增加到6个时, 缺陷豆荚内的C60分子的振荡频率与无缺陷的振荡频率基本持平. 当缺失原子超过6个时, 范德瓦耳斯力对C60分子振荡频率的影响大于摩擦力对C60分子的影响, 缺陷豆荚内C60分子的振荡频率低于无缺陷管内C60分子的振荡频率.

      图  9  C60分子的振荡频率随碳管中原子缺失数的变化

      Figure 9.  Variation of the oscillation frequency of C60 molecule with the number of missing atoms in carbon tubes.

    • 若将缺陷碳纳米豆荚运用到高频率运行的纳米器件中, 缺陷碳纳米豆荚难免会受到温度梯度或周围其他纳米器件对其的作用, 这样缺陷碳纳米豆荚将会被视为施加了作用力. 因此, 研究轴向预应力对缺陷碳纳米豆荚内C60分子的振荡的影响是非常有必要的. 本文通过模拟缺陷碳纳米豆荚在轴向预应力作用下内部C60分子的振荡, 从而探究轴向预应力对缺陷碳纳米豆荚内C60分子振荡的影响.

      缺陷碳纳米豆荚内C60分子振荡频率随轴向预应力的变化如图10所示, 其中正的预应力代表碳纳米管受拉, 而负的预应力代表碳纳米管受压. 由图10可以看出, 缺陷碳纳米豆荚内的C60分子的振荡频率随着轴向拉伸预应力的增大呈指数型减小. 首先, 在碳纳米管的顶部施加轴向拉伸预应力使得碳纳米管的长度增加, 当C60分子的运动速度不变时, C60分子的振荡频率必定会降低; 其次, 轴向拉伸应力使得C60分子在经过缺陷时所受的范德瓦耳斯力和摩檫力的改变, 导致了C60分子振荡频率的改变; 再者, 泊松效应将导致碳纳米管径向变形, 改变了碳纳米管与C60分之间的范德瓦耳斯力, 最终使得C60分子的振荡频率改变. 而轴向压缩预应力对于缺陷碳纳米豆荚的影响不像轴向拉伸预应力那样单调变化. 当轴向压缩预应力较小时, C60分子的振荡频率先单调减小后单调增大; 随着轴向压缩预应力的增大, 越过某一临界值后, C60分子的振荡频率迅速呈线性单调减小.

      图  10  缺陷碳纳米豆荚内C60分子的振荡频率随轴向预应力的变化

      Figure 10.  Variation of the oscillation frequency of C60 molecules in the defective carbon nano-peapods with axial pre-stress.

      分析可知, 在较小的轴向压缩预应力下, 轴向压缩预应力使得碳管发生了变形从而导致碳管与C60分子间范德瓦尔斯力发生改变, 与此同时轴向预应力同样改变了碳管的长度, 缺陷的出现使得C60分子的振荡频率先减小后增大; 但当轴向压缩预应力增大到一定程度后, 影响C60分子振荡频率的主要因素变为泊松效应引起的管口碳原子与C60分子间范德瓦耳斯力的改变, 管长和缺陷的影响相比之下, 对其影响甚小, 故此时无论缺陷存在与否, C60分子的振荡频率都急速下降.

    • 本文利用分子动力学方法模拟了含几种类型缺陷的碳纳米管的准静态力学性质及缺陷碳纳米豆荚中C60分子的振荡动力学行为. 研究表明: 1) 相对于无缺陷碳纳米管, 含各类型双原子空位缺陷的碳纳米管的极限应力、极限应变以及弹性模量均大幅度降低; 2) 多原子空位缺陷有沿轴向分布和沿周向分布两种类型, 但无论缺失原子沿轴向分布还是沿周向分布, 碳纳米管的弹性模量均随缺失原子数的增加而减小, 缺失原子数相同的前提下, 周向缺陷导致的模量减小更明显; 3) 碳纳米豆荚中C60分子的振荡频率主要由摩擦力和范德瓦耳斯力共同决定. 在有缺陷的碳纳米豆荚中, C60分子的振荡频率随着缺失原子数的增多先增大后减小; 4) 缺陷和轴向预应力均会对碳纳米豆荚内C60分子的振荡产生影响, 随着轴向拉伸预应力的增加, 缺陷碳纳米豆荚中C60分子的振荡频率呈指数型减小; 而当轴向压缩预应力较小时, C60分子的振荡频率先减小后增大; 随着轴向预应力的增大, 越过某一临界点后, C60分子的振荡频率呈线性单调减小.

参考文献 (26)

目录

    /

    返回文章
    返回