求助：如何确定COF应该用什么势函数？相应的参数应该怎么确定呢？

Wenlin

各位老师，我刚接触lammps，想用它模拟一个有CHNO四种元素的COF与石墨烯组成异质结的动力学过程，但是在这里我不知道怎么确定势函数，对于石墨烯我用的airebo，石墨烯与COF之间准备用lj，然后我本来打算用reaxff去模拟COF的，但是好像反应势涉及到键的断裂和生成什么的需要平衡电荷，我觉得这里COF不考虑反应就放弃了，然后看到有推荐OPLSAA什么的还有UFF，不过要弄明白每个原子的类型对应什么参数，我不太清楚这里到底应该怎么确定这个原子类型，应该怎么弄出输入文件呢？请大家指点一下

sobereva

非常推荐用GAFF或GAFF2描述COF，用sobtop（http://sobereva.com/soft/Sobtop）能很方便地产生gromacs的拓扑文件，参考主页里的说明。gromacs跑这个速度比lammps快得多，用户多得多