标题: 关于reaxff力场的“This force field file requires using 'lgvdw yes'报错 [打印本页] 作者Author: zzucaicai 时间: 2023-4-24 16:06 标题: 关于reaxff力场的“This force field file requires using 'lgvdw yes'报错 各位大佬们好，小弟最近正在做lammps关于铁铬合金的氧化，先简单做了一个Cr原子和氧气的反应，in文件如下：units real
dimension 3
boundary p p p
atom_style charge
read_data Cr.data
pair_style reax/c NULL
pair_coeff * * ffield.reax.Fe_O_Cr O Cr
neighbor 2.0 bin
neigh_modify every 10 delay 0 check no
fix 1 all nvt temp 1150.0 1150.0 10.0
fix 2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
timestep 1.0
thermo_style custom step temp etotal ke pe etotal dt time
thermo 100
dump 1 all custom 1000 CrO.xyz id type x y z fx fy fz
run 10000
但是跑了之后报错ERROR on proc 0: ffield.reax.Fe_O_Cr:70: This force field file requires using 'lgvdw yes'，我在网上查了说是要使用长程作用力，之后将pair_style 改为pair_style reax/c NULL lgvdw yes，报错ERROR on proc 0: Not a valid floating-point number: 'H'，请问各位老师，大佬该如何解决呢？ 附上力场文件
作者Author: Kmetsch 时间: 2023-4-24 17:42
for (i = 0; i < ntypes; ++i) {
// line one
values = reader.next_values(0);
++lineno;
if ((values.count() < 8) && !lgflag)
THROW_ERROR("This force field file requires using 'lgvdw yes'");
你的力场文件并不是lg形式的，产生这个报错的原因是你的力场格式出现了问题，存在某一行可被读取的参数小于8，建议检查一下空格换行等格式问题作者Author: Kmetsch 时间: 2023-4-24 17:45
X -0.1000 2.0000 1.0080 2.0000 0.0000 0.0100 -0.1000 6.0000
10.0000 2.5000 4.0000 0.0000 0.0000 5.00009999.9999 0.0000
-0.1000 0.0000 -2.3700 8.7410 13.3640 0.6690 0.9745 0.0000
-11.0000 2.7466 1.0338 2.0000 2.8793 0.0000 0.0000 0.0000
我帮你找出错误了，X的第二行缺少空格作者Author: zzucaicai 时间: 2023-4-24 18:24