ERROR: Non-existent ReaxFF type 大神们帮我看看吧 lammps力场错误

456321789

报错信息
Microsoft Windows [版本 10.0.19045.3086]
(c) Microsoft Corporation。保留所有权利。

C:\Users\Lenovo\Desktop\2>lmp -in in.lmp
LAMMPS (23 Jun 2022 - Update 3)
OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (src/comm.cpp:98)
  using 1 OpenMP thread(s) per MPI task
Loaded 1 plugins from F:\LAMMPS_MSMPI\LAMMPS 64-bit 23Jun2022-MSMPI with Python\plugins
Reading data file ...
  orthogonal box = (3.6313307 1.0774698 -0.7377933) to (32.653131 30.09927 28.284007)
  1 by 1 by 1 MPI processor grid
  reading atoms ...
  1970 atoms
  read_data CPU = 0.008 seconds
ERROR: Non-existent ReaxFF type (src/REAXFF/pair_reaxff.cpp:310)
Last command: pair_coeff           * * CHON-2019.ff n c c5 o cp h hn n3

C:\Users\Lenovo\Desktop\2>

in文件

#Initialization
units           real

atom_style      charge

timestep        0.2

read_data       C:\Users\Lenovo\Desktop\2\2.data

pair_style    reax/c  NULL
pair_coeff           * * CHON-2019.ff n c c5 o cp h hn n3

neighbor        0.4 bin
neigh_modify    every 10 one 10000

fix     1 all nve
fix     2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
fix     3 all temp/berendsen 500.0 500.0 100.0
dump    1 all atom 30 dump.reax.CHON

run     10000



力场在附件



大神们 帮我看看吧。我做的是反应力场 聚合物化学键的断裂和重组

cyx98

lammps的reaxff只区分元素类型，pair_coeff 里Cp改成C，改成对应的元素符号

456321789

cyx98 发表于 2023-7-4 09:51
lammps的reaxff只区分元素类型，pair_coeff 里Cp改成C，改成对应的元素符号

明白了，万分感谢。感谢感谢感谢感谢感谢

456321789

cyx98 发表于 2023-7-4 09:51
lammps的reaxff只区分元素类型，pair_coeff 里Cp改成C，改成对应的元素符号

你好，大佬还想请教您一个问题。我力场文件第一行就报错了，ERROR on proc 0: Not a valid integer number: '50.0000' (src/REAXFF/reaxff_ffield.cpp:584)。您知道是什么原因吗
40 ! Number of general parameters
50.0000 !p(boc1)
9.5469 !p(boc2)
26.5405 !p(coa2)
0.6863 !p(trip4)
2.7295 !p(trip3)
70.0000 !kc2
1.0588 !p(ovun6)
4.1262 !p(trip2)
12.1176 !p(ovun7)
13.3056 !p(ovun8)
-68.9784 !p(trip1)
0.0000 !Lower Taper-radius (swa)
10.0000 !Upper Taper-radius (swb)
0.0000 !not used
33.8667 !p(val7)
6.0891 !p(lp1)
1.0563 !p(val9)
2.0384 !p(val10)
6.1431 !not used
6.9290 !p(pen2)
0.3989 !p(pen3)
3.9954 !p(pen4)
0.0000 !not used
5.7796 !p(tor2)
10.0000 !p(tor3)
1.9487 !p(tor4)
0.0000 !not used
2.1645 !p(cot2)
1.5591 !p(vdW1)
0.1000 !Cutoff for bond order*100 (cutoff)
2.1365 !p(coa4)
0.6991 !p(ovun4)
50.0000 !p(ovun3)

456321789

这是我的三个文件 data 力场 in函数

cyx98

力场复制的有问题，或者整数问题。具体不太清楚，建议谷歌搜索。
https://zhuanlan.zhihu.com/p/610513062
https://matsci.org/t/error-of-not-a-valid-integer/40523/2
https://www.researchgate.net/pos ... his_error_in_LAMMPS

456321789

cyx98 发表于 2023-7-4 23:52
力场复制的有问题，或者整数问题。具体不太清楚，建议谷歌搜索。
https://zhuanlan.zhihu.com/p/610513062 ...

确实是力场问题 已经解决了感谢感谢

evenn

456321789 发表于 2023-7-5 12:52
确实是力场问题 已经解决了感谢感谢

你好，请问是力场哪里有问题呢？我也用这个立场运行的时候出现问题了

mengxiangidea

evenn 发表于 2023-8-12 10:57
你好，请问是力场哪里有问题呢？我也用这个立场运行的时候出现问题了

reaxff在读取的时候从第二行读取，你需要在第一行加上原有的注释话语，详细可见其他力场文件，第一行都有一句英文的话语

evenn

mengxiangidea 发表于 2023-8-15 12:45
reaxff在读取的时候从第二行读取，你需要在第一行加上原有的注释话语，详细可见其他力场文件，第一行都有 ...

好的，谢谢回复~确实是格式问题，已经解决啦