24小时热门版块排行榜    

查看: 4413  |  回复: 6

ws19930819

新虫 (初入文坛)

[交流] lammps计算丢失原子错误 已有6人参与

计算时出现下面的错误:
LAMMPS (10 Aug 2015)
# ReaxFF potential for TATB system
# this run is equivalent to reax/in.reax.tatb,

units                real
boundary        p p p

atom_style        full
read_data            soot.data
  orthogonal box = (-0.029171 1.30916 -0.310918) to (99.9708 101.309 99.6891)
  1 by 2 by 5 MPI processor grid
  reading atoms ...
  3600 atoms
  scanning bonds ...
  3 = max bonds/atom
  reading bonds ...
  4500 bonds
  3 = max # of 1-2 neighbors
  6 = max # of 1-3 neighbors
  18 = max # of 1-4 neighbors
  18 = max # of special neighbors

pair_style      reax/c  NULL safezone 20 mincap 2000
pair_coeff      * * ffield.reax C H
Reading potential file ffield.reax with DATE: 2010-02-19

compute reax all pair reax/c

#variable eb           equal c_reax[1]
#variable ea           equal c_reax[2]
#variable elp          equal c_reax[3]
#variable emol          equal c_reax[4]
#variable ev          equal c_reax[5]
#variable epen          equal c_reax[6]
#variable ecoa          equal c_reax[7]
#variable ehb          equal c_reax[8]
#variable et          equal c_reax[9]
#variable eco          equal c_reax[10]
#variable ew          equal c_reax[11]
#variable ep          equal c_reax[12]
#variable efi          equal c_reax[13]
#variable eqeq          equal c_reax[14]

neighbor        0.5 bin
neigh_modify        delay 0 every 10 check no


fix                1 all nvt temp 300 800 2
fix             2 all qeq/reax 1 0.0 10.0 2.0e-5 reax/c
fix             3 all reax/c/species 1 1 100 species.tab
fix                   4 all reax/c/bonds 1000 bonds.reaxc

thermo                5
thermo_modify   lost warn
thermo_style         custom step temp epair etotal press
WARNING: New thermo_style command, previous thermo_modify settings will be lost (../output.cpp:665)
#               v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa
#                v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq

timestep        1

dump                1 all custom 10000 equilibrium.lammpstrj id type q x y z ix iy iz mass element
dump            2 all xyz 10000 equilibrium.xyz
dump_modify     1 element C H
#dump_modify     2 element C H
#dump                2 all image 5 image.*.jpg type type #                axes yes 0.8 0.02 view 60 -30
#dump_modify        2 pad 3

#dump                3 all movie 5 movie.mpg type type #                axes yes 0.8 0.02 view 60 -30
#dump_modify        3 pad 3
restart         1000 coro-oh.restart
run                2000
WARNING: Using a manybody potential with bonds/angles/dihedrals and special_bond exclusions (../pair.cpp:211)
Neighbor list info ...
  2 neighbor list requests
  update every 10 steps, delay 0 steps, check no
  master list distance cutoff = 10.5
  ghost atom cutoff = 10.5
Memory usage per processor = 467.982 Mbytes
Step Temp E_pair TotEng Press
       0            0   -2352.2438   -2352.2438   -274.44276
       5    185.08898   -2337.1083    -351.4857   -1631.3032
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 7 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 7 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 8 step (../fix_qeq_reax.cpp:713)
      10    4.6625616   -2631.9359   -2581.9162   -372.36604
      15    57.138078   -2678.1304   -2065.1568   -198.71646
      20    239.99433   -2848.4589   -273.81531   -259.82012
      25    455.78352   -2973.6212     1915.995    5.7903178
      30    160.43763   -3260.8257   -1539.6612   -1436.6448
      35    343.24441   -3298.7567     383.5471   -135.92057
      40    304.72779   -3244.3692     24.73085   -74.598508
      45    137.00317   -3620.7786   -2151.0175   -118.43819
      50    341.29056   -3769.6895   -108.34642   -2.1020943
      55    132.37224   -4031.7061   -2611.6253   -504.95925
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 58 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 59 step (../fix_qeq_reax.cpp:713)
      60    631.20217   -4053.3982    2718.0979    13.988525
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 61 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 61 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 62 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 62 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 63 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 63 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 64 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 64 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 65 step (../fix_qeq_reax.cpp:713)
      65    22.031636   -4886.7011   -4650.3471   -206.46782
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 68 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 68 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 70 step (../fix_qeq_reax.cpp:713)
      70    539.69265   -5033.8316    755.95624    59.795671
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 71 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 73 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 74 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 74 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 75 step (../fix_qeq_reax.cpp:713)
      75    399.87388   -5542.6245   -1252.8034   -38.885763
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 76 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 77 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 77 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 78 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 78 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 79 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 79 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 80 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 80 step (../fix_qeq_reax.cpp:713)
      80    265.67103   -5652.0085   -2801.9069   -81.313986
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 81 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 81 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 82 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 82 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 83 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 83 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 84 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 84 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 85 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 85 step (../fix_qeq_reax.cpp:713)
      85    290.71986   -5597.9779    -2479.154   -682.79084
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 86 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 86 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 87 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 87 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 88 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 88 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 89 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 89 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 90 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 90 step (../fix_qeq_reax.cpp:713)
      90     357.1032   -5614.4995   -1783.5194   -49.230798
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 91 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 91 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 92 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 92 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 93 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 93 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 94 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 94 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 95 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 95 step (../fix_qeq_reax.cpp:713)
      95            0    4495.0559    4495.0559   -71134.069
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 96 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 96 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 97 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 97 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 98 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 98 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 99 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 99 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 100 step (../fix_qeq_reax.cpp:713)
WARNING: Fix qeq/reax CG convergence failed after 200 iterations at 100 step (../fix_qeq_reax.cpp:713)
ERROR: Lost atoms: original 3600 current 3524 (../thermo.cpp:395)


谢谢哪位大神遇到过这种错误,求指教  谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张兴锋

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,请问问题解决了吗

» 本帖已获得的红花(最新10朵)

2楼2016-10-17 21:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfan

铁杆木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
可能是timestep取值问题,目前的1有些大。
XeTaL
3楼2016-11-19 02:20:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tina张

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
引用回帖:
2楼: Originally posted by 张兴锋 at 2016-10-17 21:35:35
你好,请问问题解决了吗

你好,你是不是也遇到同样的问题了?
不忘初心
4楼2016-11-26 11:11:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

eden_1936

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我也遇到这个问题,请问你是怎么解决的?
5楼2016-12-08 17:00:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

廖明俊

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我也遇到这个问题,请问你是怎么解决的?
6楼2016-12-09 16:48:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1378625

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by 张兴锋 at 2016-10-17 21:35:35
你好,请问问题解决了吗

您好,请问您这个问题是怎么解决的呀
7楼2022-05-27 11:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ws19930819 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上意见出来了 +11 黄鸟于飞Chao 2026-08-29 20/1000 2026-08-31 17:36 by dxcharlary
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见