石墨烯力学性能的分子动力学模拟.pdfVIP

纳米压痕试验模型中石墨烯薄膜为圆形，直径为22nm，原子数为 12678；压头采用5589 个原子组成的 半径为2nm 的金刚石球，其碳-碳键长取为0.155nm 。拉伸模型中固定薄膜一端一层碳原子，对石墨烯薄 [7] 膜施加均匀拉伸应变。纳米压痕试验模型固定边界原子，中央区域原子保持自由，如图 2 所示 。碳原 子质量取为12.01(原子质量单位) 。插图序号、标题及注释居中放在图的下方。 1.2 多层石墨烯 多层石墨烯的堆积方式为 Bernal (AB)堆积，石墨烯层之间的距离取为石墨晶体的层间距，即 0.335nm，模型尺寸为 19.6800nm ×19.4579nm。层石墨烯是由碳原子通过共价键结合而成的蜂窝状二维 晶体，其中碳-碳键长约为0.142nm 。本文采用的拉伸破坏模型为矩形，其尺寸为 19.6800nm ×19.8839nm， 图1 单层石墨烯拉伸模型 固定石墨烯原子 图2 纳米压痕试验模拟模型 2 模拟方法及过程 本文的研究基于分子动力学(MD)方法，其基本原理是从微观尺度建立一个原子系统，通过假定原子 间的相互作用势，求得每个原子所受到的力，在选定的时间步长、边界和初始条件下，建立原子的牛顿 运动方程，用数值方法求解，得到系统中原子的运动轨迹和速度分布，再应用统计理论得到系统的力学、 热学和动力学特性。 原子间势函数的选取是决定分子动力学模拟精度的关键。Tersoff 势函数是共价键体系使用最广泛的 势函数，它可以很好的模拟碳-碳共价键的各种特性，包括键长、键角、键能、键的断裂常数和晶格常数 等动态行为，能够较真实地反映碳元素所构成固体材料的物理特性。因此，本文对石墨烯层内碳-碳原子 间选用Tersoff 势函数。而石墨烯层间、以及石墨烯与金刚球压头的碳原子之间则采用Lennard-Jones 势 进行模拟。 Tersoff 势的表达式为: 2 1 E E V (1)   i 2 ij i ij   V f (r ) f (r ) b f (r ) (2) ij c ij R ij ij A ij 式中，E 是体系总能量，V 为i ,j 原子之间的成键能量，f 和f 分别是对势的吸引项和排斥项，f 是 ij A R c 光滑截断函数，b 是键序函数。其中键序函数b 用于表征