应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】NEB计算的一个stdout输出怎么解读?
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 2110  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[交流] 【求助成功】NEB计算的一个stdout输出怎么解读? 已有3人参与

放假前一天把vststool Code和vasp4.6一起编译了.在两个端点中间插入3个images.下面是中间02文件夹中的stdout,OSZICAR及INCAR.希望得到帮助.
问题:
(1)stdout中第一步后能量收敛就很快(很快到1E-7),说明什么?后面的quench提示是什么意思?
(2)OSZICAR中电子步中能量是-400多,结果出来F=-1400多.这个是输出错误吧?stdout中的离子步能量也是这样的.是不vststool与vasp编译的问题?OUTCAR中每离子步的能量与电子步收敛后的能量是相符的.
stdout:

vasp.4.6.28 25Jul05 complex
02/POSCAR found :  2 types and   55 ions
LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro
00/POSCAR found :  2 types and   55 ions
04/POSCAR found :  2 types and   55 ions
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: wrap around errors must be expected
FFT: planning ...           9
reading WAVECAR
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.554736476659E+04    0.55474E+04   -0.18260E+05  3200   0.138E+03
DAV:   2     0.154389023593E+04   -0.40035E+04   -0.37247E+04  3860   0.348E+02
..............................................................................................................
RMM:  31    -0.467372515298E+03   -0.41006E-03   -0.26808E-04  2111   0.242E-02    0.431E-02
RMM:  32    -0.467372569584E+03   -0.54286E-04   -0.10334E-04  1827   0.215E-02
   1 F= -.14034855E+04 E0= -.14032665E+04  d E =-.140349E+04
quench:  g(F)=  0.000E+00 g(S)=  0.000E+00 dE (1.order)= 0.000E+00
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.467372620179E+03   -0.46737E+03   -0.14539E-04  2198   0.171E-02    0.109E-02
RMM:   2    -0.467372627169E+03   0.69903E-05   -0.48155E-05 1624   0.143E-02    0.168E-02
RMM:   3    -0.467372628565E+03   -0.13963E-05   -0.16158E-05  1509   0.936E-03    0.156E-02
RMM:   4    -0.467372628187E+03    0.37851E-06   -0.53661E-06  1176   0.382E-03    0.389E-03
RMM:   5    -0.467372628897E+03   -0.71051E-06   -0.92475E-07   750   0.114E-03
2 F= -.14034856E+04 E0= -.14032666E+04  d E =-.140349E+04
quench:  g(F)=  0.000E+00 g(S)=  0.000E+00 dE (1.order)= 0.000E+00
bond charge predicted

OSZICAR:

       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.554736476659E+04    0.55474E+04   -0.18260E+05  3200   0.138E+03
DAV:   2     0.154389023593E+04   -0.40035E+04   -0.37247E+04  3860   0.348E+02
DAV:   3    -0.316744761608E+03   -0.18606E+04   -0.18065E+04  4256   0.194E+02
.............................................................................................
RMM:  30    -0.467372105235E+03   -0.17295E-02   -0.16898E-03  2655   0.635E-02    0.330E-02
RMM:  31    -0.467372515298E+03   -0.41006E-03   -0.26808E-04  2111   0.242E-02    0.431E-02
RMM:  32    -0.467372569584E+03   -0.54286E-04   -0.10334E-04  1827   0.215E-02
   1 F= -.14034855E+04 E0= -.14032665E+04  d E =-.140349E+04
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.467372620179E+03   -0.46737E+03   -0.14539E-04  2198   0.171E-02    0.109E-02
RMM:   2    -0.467372627169E+03   -0.69903E-05   -0.48155E-05  1624   0.143E-02    0.168E-02
RMM:   3    -0.467372628565E+03   -0.13963E-05   -0.16158E-05  1509   0.936E-03    0.156E-02
RMM:   4    -0.467372628187E+03    0.37851E-06   -0.53661E-06  1176   0.382E-03    0.389E-03
RMM:   5    -0.467372628897E+03   -0.71051E-06   -0.92475E-07   750   0.114E-03
   2 F= -.14034856E+04 E0= -.14032666E+04  d E =-.140349E+04
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.467372629463E+03   -0.46737E+03   -0.17888E-06  2112   0.391E-03    0.181E-03
RMM:   2    -0.467372629846E+03   -0.38340E-06   -0.27115E-09   505   0.434E-06    0.261E-03
RMM:   3    -0.467372629831E+03    0.14523E-07    0.00000E+00   504   0.000E+00    0.131E-03
RMM:   4    -0.467372629813E+03    0.18692E-07   -0.28324E-09   505   0.402E-06    0.181E-03
RMM:   5    -0.467372629420E+03    0.39265E-06   -0.48393E-07   635   0.531E-04
   3 F= -.14034856E+04 E0= -.14032665E+04  d E =-.140349E+04

INCAR:

ENCUT = 480
ISTART= 0 ;ICHARG= 2
NPAR = 2
NELMIN = 5
EDIFF = 1E-4
EDIFFG = -0.05
LREAL = Auto
NSW = 60
ALGO = Fast
LPLANE = .TRUE.
NSIM = 4
GGA = 91
POTIM = 0.0
ICHAIN = 0
LCLIMB = .TRUE.
LGLOBAL = .TRUE.
MAXMOVE = 0.2
INVCURV = 0.01
IMAGES = 3
SPRING = -5
IBRION = 3
ISMEAR=0
VOSKOWN=1
#PREC=Accurate

[ Last edited by cenwanglai on 2010-7-1 at 20:45 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-06-17 09:29:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813276

铜虫 (小有名气)
cenwanglai(金币+6):谢谢.汗!都没检查出来少了IOPT. 2010-06-17 14:18:40
要加IOPT=1 不加的话就没优化了 qench是指IBRION=3的关系 没事的
那个输出问题不知道怎么回事
赞一下(1人) 回复此楼
2楼2010-06-17 12:31:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
引用回帖:
Originally posted by gump_813276 at 2010-06-17 12:31:31:
要加IOPT=1 不加的话就没优化了 qench是指IBRION=3的关系 没事的
那个输出问题不知道怎么回事

请问OUTCAR文件中:

(1)POSITION          TOTAL-FORCE
这个是每个原子的坐标及受力.对于非端点的image,所考察的移动原子在不同image中受力是不受EDIFFG=-xE-y(如-1E-2)约束的吧?根据NEB的原理,垂直与chain方向受力应该->0,但是沿着chain方向是受力的.请问EDIFFG对哪些力有约束呢?

(2)TANGET          CHAIN-FORCE
跟(1)有对应的六个数,表示什么意思呢?分别表示切向受力和垂直与CHAIN方向受力?

(3)前两着之间有如下信息:
energy of chain is (eV)  0.0000 for this image  0.000000
tangential force (eV/A)  -2.644444
left and right image  0.6379  0.6379 A
这行说明中的文字部分是什么意思?数字表示什么物理意义?
赞一下 回复此楼
3楼2010-06-17 14:35:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813276

铜虫 (小有名气)
★ ★
hedaors(金币+2):谢谢回复 2010-06-17 23:21:13
cenwanglai(金币+6):呵呵,谢谢回复! 2010-06-18 08:42:12
cenwanglai(金币+1):呵呵,有空请看看下面一条. 2010-06-18 09:52:40
(1)EDIFFG应该是垂直方向上的力的收敛条件
(2)这个部分我算的输出和你的有点不一样 可能是版本更新了吧
     前三个数应该是TANGENT方向 CHAIN-FORCE不知道具体指什么
(3) 我也没看到这部分输出
      energy of chain 不知道什么东西 而且是0 有点奇怪
      tangential force就是切向力
      left and right image是和旁边image的距离 这个距离应该差不多等距 不过climbing image两边是不等距离的


另外, 你可以在网页上下载vtstscripts里面有nebef.pl这个script 可以直接输出每个image上的力和能量
这个力就是收敛标准所用的力。 里面还有很多和neb相关的script
赞一下(1人) 回复此楼
4楼2010-06-17 23:18:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
算过度态,怎么选择固定原子呢?比如可不可以把slab全部固定,只留考察的目标移动原子?
赞一下 回复此楼
5楼2010-06-18 09:51:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813276

铜虫 (小有名气)
★
cenwanglai(金币+5):呵呵，我另外有个Au吸附O胡帖子，还没有人回复。 2010-06-18 10:51:37
zzy870720z(金币+1):感谢交流 2010-06-18 11:23:37
可以啊 不过rigid的slab算出的能类不准
方法就是selective dynamics 不动的就F F F
赞一下(1人) 回复此楼
6楼2010-06-18 10:48:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by cenwanglai at 2010-06-17 14:35:54
请问OUTCAR文件中:

(1)POSITION          TOTAL-FORCE
这个是每个原子的坐标及受力.对于非端点的image,所考察的移动原子在不同image中受力是不受EDIFFG=-xE-y(如-1E-2)约束的吧?根据NEB的原理,垂直与chain ...

请问区长，关于这个问题
>energy of chain is (eV)  0.0000 for this image  0.000000
tangential force (eV/A)  -2.644444
left and right image  0.6379  0.6379 A
这行说明中的文字部分是什么意思?数字表示什么物理意义?
>>您现在知道什么意思不？
赞一下 回复此楼
慢慢的思考。。。
7楼2014-04-04 09:03:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿大连理工大学材料求调剂 +5 Gymno 2026-03-30 5/250 2026-03-31 04:18 by fmesaito
[考研] 一志愿085600中科院宁波所276分求调剂 +4 材料学257求调剂 2026-03-28 5/250 2026-03-31 03:57 by fmesaito
[考研] 抱歉 +4 田洪有 2026-03-30 4/200 2026-03-30 21:26 by mumin1990
[考研] 085602化工求调剂（331分） +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 一志愿华东师范大学有机化学专业，初试351分，复试被刷求调剂! +7 真名有冰 2026-03-29 8/400 2026-03-30 21:08 by maaj02
[考研] 材料工程专硕求调剂 +8 hyl3153942 2026-03-29 8/400 2026-03-30 20:45 by dophin1985
[考研] 085600，专业课化工原理，320分求调剂 +5 大馋小子 2026-03-29 5/250 2026-03-30 19:02 by 源_2020
[考研] 材料学硕333求调剂 +14 北道巷 2026-03-24 14/700 2026-03-30 18:59 by 源_2020
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +5 贝呗钡钡 2026-03-30 5/250 2026-03-30 18:16 by 无际的草原
[考研] 317分 一志愿南理工材料工程 本科湖工大 求调剂 +12 芋泥小铃铛 2026-03-28 12/600 2026-03-30 17:06 by wangjy2002
[考研] 071010 323 分求调剂 +3 Baekzhy 2026-03-27 3/150 2026-03-30 14:24 by andresqi
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 311求调剂 +10 lin0039 2026-03-26 10/500 2026-03-30 10:26 by herarysara
[考研] 317求调剂 +10 蛋黄咸肉粽 2026-03-26 10/500 2026-03-30 09:45 by longlotian
[考研] 085602 307分 求调剂 +7 不知道叫什么！ 2026-03-26 7/350 2026-03-28 09:57 by 神马都不懂
[考研] 一志愿上海理工能源动力（085800）310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 300分，材料，求调剂，英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
信息提示
关闭
请填处理意见
关闭