研究单壁碳纳米管振动的分子力学方法.docVIP

研究单壁碳纳米管振动的分子力学方法 文章信息： 文章历史: 2010年1月29日收到投稿 2010年3月9日收到修订后的形式 2010年3月15日被收入 2010年4月10日允许网上发表 关键词： 振动分析 碳纳米管 分子力学 固有频率 振型 摘要 我们采用分子力学方法研究锯齿形的振动特性和扶手椅型单壁碳纳米管（CNT）。自振频率和他们相关的固有振动模式已经被获得。分别对四种类型的进行了锯齿形纳米管(5,0),(6,0)、(8,0)、(10 0)和三种类型的扶手椅纳米管(3),(4,4),(6 6)进行了模拟计算。通用力场势被应用于分子力学方法。前五个固有频率宽高比从5到20范围波动。结果表明,宽高比增加时固有频率降低。本文中结果与先前用结构力学的方法研究碳纳米管的结果有类似的趋势。 介绍 碳纳米管（CNT ）[1]的发现激发了对其卓越的机械性能[2,3]和在纳米技术领域的应用潜力[4,5]的广泛的研究。众所周知，碳纳米管具有非常高的弹性模量和低质量密度在纳米电子学纳米器件和NEMS中可以作为太赫兹纳米谐振器[6-9]。具体来说，特别是振动模式，例如，径向呼吸模式[10-14]，横向模式[15]和纵向模式[16]，对原子结构和碳纳米管的弹性性质提供了有价值的探索。另一方面，碳纳米管的振动通常改变了原先平直的单壁碳纳米管的几何结构（单壁碳纳米管）为圆形横截面，也改变了多壁碳纳米管的同心结构。这种变形可能导致几何对称性敏感碳纳米管的电气特性的突然变化[6,7,9,19–24]和连续介质力学理论[8,25-29]被用来预测声子色散[18-22]关系，获取振动的变形模式[27,28]和量化碳纳米管在几何尺寸[11,12]和分子结构[26]方面和振动频率的关系。这些研究，然而，主要集中在多壁碳纳米管.这些对单壁碳纳米管的振动分析主要集中在基本模式[6]，拉曼活性模式，包括径向呼吸模式和切向模式[10-13]。尽我们所知，对单壁碳纳米管的圆筒形结构的一般振动行为的全面研究和深刻理解都没有详尽的探究。尤其众所周知的是单壁碳纳米管具有各种原子的结构，即，扶手椅，手性和锯齿形，其手性角从0度上升到30度[30]。直至目前为止这种原子结构对单壁碳纳米管的振动行为的影响还没有得到足够重视，除了最近基于分子力学（ MM ）模型[31]对锯齿型碳纳米管的研究和一些对径向呼吸模式[13,14]的早期研究。 在这种情况下，本论文旨在通过计算单壁碳纳米管扶手椅和锯齿形结构的一般振动光谱，以及不同半径和长度与直径的纵横比，以填补此特定的空缺。相关振动模式也将被应用于单壁碳纳米管和用于了解它们的频率对单壁碳纳米管的几何尺寸的相关性。这里特别要注意单壁碳纳米管的原子结构对振动的影响。对连续模型忽略的已被使用在单壁碳纳米管的振动分析中的原子细节这个问题的研究有极大的兴趣。为了实现这些目标，一个占单壁碳纳米管的原子结构的MM模型已经被开发和利用去进行单壁碳纳米管的振动分析。本文以下面的方式组成。第2节的重点是对MM模型的开发和验证。数值结果与讨论将在第3节介绍。最后，本文将基于第3节的结果和分析的在第4节得出主要结论。 2.分子力学方法的简要概述 从分子力学的观点出发，总能量的一般表达式为能量由于价键合或相互作用和非键合相互作用的总和。 价相互作用包括键的伸缩和角度扭曲。角变形包括键角弯曲，二面角扭转和反转条件。在无粘合相互作用包括范德华力和静电条款。在这项研究中，我们使用了通用力场（UFF）模型[32]，其力场参数通过使用仅基于元素，其杂交及其连接的一般规则已经被估算。在本研究中使用的力场函数形式和参数都是按照Rappe等... 2.1 碳纳米管的优化 锯齿型碳纳米管（n,0）的最佳平衡构型是当n=5，6，8和10时。对于扶手椅碳纳米管（n,n）来说n=3,4和6，被认为是最佳的在这个振动分析中。几何优化均采用纳米管。C-C键的最佳键长和C-C键平均键角的优化的键长分别列于表1和表2。一些微小的数量变化可以在这两个表中看出。这种变化，通常忽略了结构力学和连续介质力学的方法，在被提出的方法中是考虑到的。 我们先从Hessian矩阵，其持有的电位E的二阶偏导数相对于直角坐标系（CAR）的原子位移开始 这是一个3N ×3N矩阵（N是原子的数目），其中用于在笛卡尔坐标的替换 指的衍生物采取的原子的平衡位置的事实。这些力常数随后被转换为质量加权笛卡尔坐标（MWC）。 其中等等。是对角化，产生一组3N特征向量和3N的特征值。 接下来的步骤是将质量中心到原点，并确定惯性矩和惯性积，以发现对角化的惯性矩张量的矩阵的目标。使用这个矩阵，我们可以找到对应的旋转和平移向量。一旦这些载体是已知的，我们知道，正常模式的其余部分是振动。质