|
需要你在gmx文件夹下的/share/gromacs/top/oplsaa.ff/文件夹下修改atomname2type.n2t文件内容,添加: C opls_145 0 12.011 3 C 0.142 C 0.142 C 0.142 C opls_145 0 12.011 2 C 0.142 C 0.142 C opls_145 0 12.011 1 C 0.142 这样就会把石墨烯的碳原子定义为OPLS-AA力场内的芳香sp2碳,以上是我根据培训讲义内容中的碳纳米管模拟部分(讲义对应GROMOS 54A7)拓展到OPLS-AA力场下的石墨烯体系,修改完后执行echo 16|gmx x2top -f model.gro -o md.top -ff select -nopbc -name md,然后删除前22行与最后8行后将md.top修改为md.itp即可 |
aymy 发表于 2026-1-25 16:11 n2t文档内容写的不对,网上有很多范例可以查到n2t文档应该怎么写 |
sobereva 发表于 2026-1-25 16:05 老师您好,因为我正在复现一篇文献里的内容,文献作者提到了使用OPLS-AA,所以发帖求助,也可以学习,明白自己的错误。 那我先继续用sobtop继续模拟着,感谢 |
|
无论是GAFF还是OPLS-AA,都没有专门给石墨烯搞的参数,显然也不存在非要用OPLS-AA不可的理由 但凡OPLS-AA能描述的体系,GAFF或GAFF2几乎都能合理描述;但凡GAFF或GAFF2无法描述好的情况,尝试用OPLS-AA也基本徒劳 用GAFF可以享受用sobtop创建拓扑文件的极度便利,而且又有极为丰富的与之兼容的力场,而且原子类型名字又远远比OPLS-AA十分抽象的数字清楚得多、友好得多 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-16 20:18 , Processed in 0.195723 second(s), 25 queries , Gzip On.