计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2233|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] Lammps做反应力场报错ERROR on proc 0: step 19: hbondchk failed.

[复制链接 Copy URL]

33

帖子

0

威望

57

eV
积分
90

Level 2 能力者

各位大神!我做反引力场的时候报错ERROR on proc 0: step 19: hbondchk failed: H=79 end(H)=2476 str(H+1)=2470
大家知道这需要怎么该吗?
lammps结果
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
#Initialization
units           real

atom_style      charge

timestep        0.2

read_data       C:\Users\Lenovo\Desktop\poly1\2\2.data
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

pair_style    reax/c  NULL
pair_coeff           * * 2019.ff.data n c C o C h H N
fix             2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
neighbor        0.4 bin

neigh_modify    every 10 delay 0 check no

fix     1 all nve
velocity        all create 300.0 9999999 dist gaussian
fix     3 all temp/berendsen 500.0 500.0 100.0
dump    1 all atom 30 dump.reax.CHON

run     20000

CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE

Your simulation uses code contributions which should be cited:

- pair reaxff command:

@Article{Aktulga12,
author = {H. M. Aktulga, J. C. Fogarty, S. A. Pandit, A. Y. Grama},
title = {Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques},
journal = {Parallel Computing},
year =    2012,
volume =  38,
pages =   {245--259}
}

- fix qeq/reaxff command:

@Article{Aktulga12,
author = {H. M. Aktulga, J. C. Fogarty, S. A. Pandit, A. Y. Grama},
title = {Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques},
journal = {Parallel Computing},
year =    2012,
volume =  38,
pages =   {245--259}
}

CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE-CITE

Neighbor list info ...
  update every 10 steps, delay 0 steps, check no
  max neighbors/atom: 2000, page size: 100000
  master list distance cutoff = 10.4
  ghost atom cutoff = 10.4
  binsize = 5.2, bins = 6 6 6
  2 neighbor lists, perpetual/occasional/extra = 2 0 0
  (1) pair reax/c, perpetual
      attributes: half, newton off, ghost
      pair build: half/bin/newtoff/ghost
      stencil: full/ghost/bin/3d
      bin: standard
  (2) fix qeq/reax, perpetual, copy from (1)
      attributes: half, newton off
      pair build: copy
      stencil: none
      bin: none
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 0 (src/REAXFF/fix_qeq_reaxff.cpp:769)
Per MPI rank memory allocation (min/avg/max) = 152.2 | 152.2 | 152.2 Mbytes
   Step          Temp          E_pair         E_mol          TotEng         Press     
         0   300           -126403.21      0             -124642.44      118990.35   
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 1 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 2 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 3 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 4 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 5 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 6 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 7 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 8 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 8 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 9 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 9 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 10 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 11 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 11 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 12 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 13 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 14 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 15 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 15 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 16 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 16 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 17 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 18 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 19 (src/REAXFF/fix_qeq_reaxff.cpp:769)
WARNING: Fix qeq/reaxff CG convergence failed after 200 iterations at step 20 (src/REAXFF/fix_qeq_reaxff.cpp:769)
ERROR on proc 0: step 19: hbondchk failed: H=79 end(H)=2476 str(H+1)=2470
(src/REAXFF/reaxff_forces.cpp:119)
Last command: run     20000

我的in函数
#Initialization
units           real
atom_style      charge
timestep        0.2

read_data       C:\Users\Lenovo\Desktop\poly1\2\2.data
pair_style      reax/c  NULL
pair_coeff      * * 2019.ff.data n c C o C h H N

fix             2 all qeq/reax 1 0.0 10.0 1e-6 reax/c
neighbor        0.4 bin
neigh_modify    every 10 delay 0 check no

fix             1 all nve
fix             3 all temp/berendsen 500.0 500.0 100.0
dump            1 all atom 30 dump.reax.CHON

run             20000

模型文件和力场在附件里
求求各位大神了!帮我看看吧!


2.data

92.9 KB, 下载次数 Times of downloads: 5

模型

2019.ff.data

8.5 KB, 下载次数 Times of downloads: 10

势函数

1

帖子

0

威望

101

eV
积分
102

Level 2 能力者

2#
发表于 Post on 2023-12-25 15:28:29 | 只看该作者 Only view this author
请问您解决了吗,这个问题

1

帖子

0

威望

15

eV
积分
16

Level 1 能力者

3#
发表于 Post on 7 day ago | 只看该作者 Only view this author
请问这个问题您解决了吗,我现在也遇到了完全相同的问题

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 00:47 , Processed in 0.161075 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list