关于reaxff力场的“This force field file requires using 'lgvdw yes'报错

zzucaicai

各位大佬们好，小弟最近正在做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'，请问各位老师，大佬该如何解决呢？ 附上力场文件

XXXIN

请问大佬查找是哪一项问题的那个代码是在哪里运行的呢？

zzucaicai

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  ...

谢谢大佬！问题解决了

Kmetsch

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的第二行缺少空格

Kmetsch

        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，建议检查一下空格换行等格式问题