使用reaxff力场，出现 ffield.reax.chon:46: Invalid force field file format报错

wangqian123

想请教各位大神，我在使用lammps做聚酰亚胺热裂解分子动力学研究时出现如下报错，请问如何解决ERROR on proc 0: ffield2016.reax.chon:46: Invalid force field file format (src/REAXFF/reaxff_ffield.cpp:584)
Last command: pair_coeff * * ffield2016.reax.chon n c' o' c5 cp h  o hn

data文件部分为
1970 atoms
   8 atom types
     4.641319048    50.569719048 xlo xhi
     3.988655190    49.917055190 ylo yhi
    -5.263177289    40.665222711 zlo zhi
Masses
   1  14.006700 # n
   2  12.011150 # c'
   3  15.999400 # o'
   4  12.011150 # c5
   5  12.011150 # cp
   6   1.007970 # h
   7  15.999400 # o
   8   1.007970 # hn
Atoms
      1      1   1 -0.110000    27.988040887    35.887421566    34.131797464    # n
      2      1   2  0.380000    29.162045629    35.975962615    34.847921842    # c'
      3      1   2  0.380000    27.848985595    37.034262261    33.380564516    # c'

in 文件部分语句：
read_data        reax.data
replicate 3 3 3
pair_style        reax/c NULL lgvdw yes
pair_coeff * * ffield2016.reax.chon n c' o' c5 cp h o hn
timestep        0.1

reaxff力场部分文件：

力场文件中45 46 47行内容如下
cov r3;Elp;Heat inc.;n.u.;n.u.;n.u.;n.u.
ov/un;val1;n.u.;val3,vval4
C 1.3742 4.0000 12.0000 1.9684 0.1723 0.8712
完整的发到附件了

Kmetsch

data文件中的元素不是力场支持的类型，你的pair_coeff行定义的名字也不是力场支持的名字。更改你的data使其只有CHON元素类型。

Kmetsch

另外，你在这一行启用了内层电子斥力和lg长程色散

pair_style        reax/c NULL lgvdw yes

，但你的力场并不是支持的类型。如果你的体系没有特殊关注的键的类型的话，可以去找一篇类似体系的reaxff/lg的文献的力场拿来用

25110031

你好请问你的问题解决了嘛

25110031

Kmetsch 发表于 2023-4-5 23:53
另外，你在这一行启用了内层电子斥力和lg长程色散，但你的力场并不是支持的类型。如果你的体系没有特殊关注 ...

你好能不能帮我看下我的问题

cyx98

in文件改成： pair_coeff * * ffield2016.reax.chon n c o c c h o h
再试一试

wangqian123

25110031 发表于 2023-4-19 09:02
你好请问你的问题解决了嘛

已经解决了，力场文件没选对，导出data文件也不错，改完就好了

wangqian123

Kmetsch 发表于 2023-4-5 23:49
data文件中的元素不是力场支持的类型，你的pair_coeff行定义的名字也不是力场支持的名字。更改你的data使其 ...

谢谢

耿无敌

请问您现在用的哪个力场，那个力场为什么不对啊？我看有的文献里就用的那个力场。

tiank

wangqian123 发表于 2023-4-21 10:20
已经解决了，力场文件没选对，导出data文件也不错，改完就好了

你好，我也是做热解的，跟你遇到了一样的情况，可以交流一下吗？q805698973

456321789

您好，我也是做聚酰亚胺材料的，也出现了这个问题，能问一下您用的是什么力场吗

linlink

您好，我也在做酰胺类热降解，请问您模拟产物后有没有出现产物的N原子个数大于原始模型的情况。

peekaboo

cyx98 发表于 2023-4-19 09:37
in文件改成： pair_coeff * * ffield2016.reax.chon n c o c c h o h
再试一试

您好能麻烦您帮我看一下我的问题吗