| 查看: 2119 | 回复: 6 | |||
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 ] |
回复此楼» 猜你喜欢
（调剂）一志愿报考哈尔滨工业大学0857资源与环境专业378分考生
已经有6人回复
-
262求调剂
已经有7人回复
-
材料专硕322分
已经有10人回复
-
材料调剂
已经有8人回复
-
085600，320分求调剂
已经有8人回复
-
一志愿 江南大学 085602 化工专硕 338分求调剂
已经有11人回复
-
求调剂
已经有6人回复
-
086000调剂
已经有3人回复
-
求调剂
已经有4人回复
-
0703调剂，一志愿天津大学319分
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 新手请教用NEB方法计算energy barrier.
已经有15人回复
- vasp用neb计算迁移势垒
已经有11人回复
- vasp NEB 如何选取原胞大小 ？
已经有10人回复
- 【求助】请教CNEB
已经有15人回复
- 【求助】NEB计算outcar中无各个镜像的能量输出
已经有6人回复
- 【求助】关于CI-NEB
已经有21人回复
- 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析
已经有30人回复
- 【求助】关于vasp NEB
已经有8人回复
- 【其他】vasp neb 计算讨论
已经有20人回复
- 【求助成功】怎么样提高NEB计算速度?
已经有19人回复
- 【求助/交流】请教NEB的pmeI的酶切体系
已经有7人回复
- 【讨论】NEB寻找过渡态专题讨论
已经有49人回复
- 【求助】vasp的NEB计算,输出需要哪些参数?
已经有7人回复
2楼2010-06-17 12:31:31
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
|
请问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
★ ★
hedaors(金币+2):谢谢回复 2010-06-17 23:21:13
cenwanglai(金币+6):呵呵,谢谢回复! 2010-06-18 08:42:12
cenwanglai(金币+1):呵呵,有空请看看下面一条. 2010-06-18 09:52:40
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 |
4楼2010-06-17 23:18:36
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
5楼2010-06-18 09:51:58
6楼2010-06-18 10:48:59
刘小伟110
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1635.7
- 散金: 51
- 帖子: 256
- 在线: 120.5小时
- 虫号: 2096954
- 注册: 2012-10-30
- 专业: 理论和计算化学
7楼2014-04-04 09:03:08
