前几天使用CH.airebo势函数绘制了单层石墨烯拉伸的应力应变曲线，在查看轨迹文件后，发现拉伸过程中出现了不合理的多次强化阶段。

后来经查阅有关资料发现，出现该现象的原因是因为CH.airebo势函数第7行的参数rcmin_CC需要修改，默认是1.70。

在AIREBO势表达式中的键序项bij包含多组开关函数，通过cutoff distance来判断，拉伸键长的改变在断裂后原有的数值不能真实反映结果，需要对此进行修正。查阅了几篇文献发现该参数一般都是设置的2.0，我自己通过改变该参数做了个测试，测试结果如下：

拉伸的应变率为10e8/s,使用ppp边界+fix npt 来控制模拟环境，拉伸采用fix deform命令。如果文中有错误或者不合适的地方欢迎大家指正。