计算化学公社

标题: 使用reaxff力场，出现 ffield.reax.chon:46: Invalid force field file format报错 [打印本页]

作者
Author: wangqian123 时间: 2023-4-5 21:06
标题: 使用reaxff力场，出现 ffield.reax.chon:46: Invalid force field file format报错
想请教各位大神，我在使用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力场部分文件：
(, 下载次数 Times of downloads: 36)
力场文件中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
完整的发到附件了

作者
Author: Kmetsch 时间: 2023-4-5 23:49
data文件中的元素不是力场支持的类型，你的pair_coeff行定义的名字也不是力场支持的名字。更改你的data使其只有CHON元素类型。

作者
Author: Kmetsch 时间: 2023-4-5 23:53
另外，你在这一行启用了内层电子斥力和lg长程色散

pair_style reax/c NULL lgvdw yes

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

作者
Author: 25110031 时间: 2023-4-19 09:02
你好请问你的问题解决了嘛

作者
Author: 25110031 时间: 2023-4-19 09:03

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

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

作者
Author: cyx98 时间: 2023-4-19 09:37
in文件改成： pair_coeff * * ffield2016.reax.chon n c o c c h o h
再试一试

作者
Author: wangqian123 时间: 2023-4-21 10:20

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

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

作者
Author: wangqian123 时间: 2023-4-21 10:21

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

谢谢

作者
Author: 耿无敌 时间: 2023-4-29 10:57
请问您现在用的哪个力场，那个力场为什么不对啊？我看有的文献里就用的那个力场。

作者
Author: tiank 时间: 2023-6-13 14:03

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

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

作者
Author: 456321789 时间: 2023-7-4 16:26
您好，我也是做聚酰亚胺材料的，也出现了这个问题，能问一下您用的是什么力场吗

作者
Author: linlink 时间: 2023-11-16 09:44
您好，我也在做酰胺类热降解，请问您模拟产物后有没有出现产物的N原子个数大于原始模型的情况。

作者
Author: peekaboo 时间: 2025-4-15 15:38

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)