|
|
老师们好,我是按照http://bbs.keinsci.com/thread-38799-1-1.html这个例子来做的,在https://www.scm.com/doc/ReaxFF/Included_Forcefields.html找的关于C H O Fe Cl的reaxff力场文件,data文件通过ms导出pdb文件,然后通过VMD topo writelammpsdata <filename> charge 命令导出了data文件。但是每次运行lammps都会提示ERROR: Incorrect args for pair coefficients (src/REAXFF/pair_reaxff.cpp:280)Last command: pair_coeff * * /Users/tianzhongxun/lammps/lammpscalculation/ffield.reax C H O Fe Cl。因为我的体系中没有O,无论我把O甚至为NULL还是把O设置进去都会出现这个错误,是不是必须元素全部吻合才可以呢?或者我在M$中用forcite的cvff力场去计算了一下能量然后导出了car格式的文件,导致了data文件中有33个atom_types,是不是不需要用力场去assign,计算能量呢?我的文件在附件中
|
|