计算化学公社

标题: ERROR: Non-existent ReaxFF type 大神们帮我看看吧 lammps力场错误 [打印本页]

作者
Author: 456321789 时间: 2023-7-3 10:34
标题: ERROR: Non-existent ReaxFF type 大神们帮我看看吧 lammps力场错误
报错信息
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

力场在附件
(, 下载次数 Times of downloads: 17)

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

作者
Author: cyx98 时间: 2023-7-4 09:51
lammps的reaxff只区分元素类型，pair_coeff 里Cp改成C，改成对应的元素符号

作者
Author: 456321789 时间: 2023-7-4 11:46

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

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

作者
Author: 456321789 时间: 2023-7-4 15:38

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)

作者
Author: 456321789 时间: 2023-7-4 15:40
这是我的三个文件 data 力场 in函数

作者
Author: cyx98 时间: 2023-7-4 23:52
力场复制的有问题，或者整数问题。具体不太清楚，建议谷歌搜索。
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

作者
Author: 456321789 时间: 2023-7-5 12:52

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

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

作者
Author: evenn 时间: 2023-8-12 10:57

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

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

作者
Author: mengxiangidea 时间: 2023-8-15 12:45

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

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

作者
Author: evenn 时间: 2023-8-16 11:37

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

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

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