| 查看: 2179 | 回复: 5 | |||
dfjk-123新虫 (小有名气)
|
[求助]
OSZICAR中数据显示电子自洽和离子弛豫收敛标准都已经达到，为何OUTCAR还在继续算？？
|
|
已设定EDIFF=0.1E-02，EDIFFG=0.1E-01，按理电子自洽前后两步的dE小于 0.1E-02就会收敛，而离子弛豫dE小于0.1E-01就会收敛，OSZICAR已经显示这两个标准达到了呀，为何OUTCAR还在继续算呢？？ N E dE d eps ncg rms rms(c) RMM: 1 -0.329376635237E+03 -0.36277E-02 -0.21313E-01 6954 0.791E-01 0.511E-01 RMM: 2 -0.329476459066E+03 -0.99824E-01 -0.29802E-02 7173 0.371E-01 0.432E+00 RMM: 3 -0.329363777178E+03 0.11268E+00 -0.13767E-02 6313 0.236E-01 0.142E+00 RMM: 4 -0.329370731742E+03 -0.69546E-02 -0.50608E-03 5370 0.182E-01 0.579E-01 RMM: 5 -0.329371599246E+03 -0.86750E-03 -0.12849E-03 4117 0.962E-02 82 F= -.32937160E+03 E0= -.32937160E+03 d E =0.143016E-02 BRION: g(F)= 0.281E-01 g(S)= 0.000E+00 retain N= 3 mean eig= 1.73 eig: 4.931 0.034 0.225 bond charge predicted N E dE d eps ncg rms rms(c) RMM: 1 -0.329403251181E+03 -0.32519E-01 -0.15638E+00 7158 0.206E+00 0.105E+00 RMM: 2 -0.329640404592E+03 -0.23715E+00 -0.13164E-01 7308 0.887E-01 0.668E+00 RMM: 3 -0.329381072940E+03 0.25933E+00 -0.40045E-02 7035 0.372E-01 0.260E+00 RMM: 4 -0.329379776035E+03 0.12969E-02 -0.12272E-02 6227 0.269E-01 0.125E+00 RMM: 5 -0.329381911263E+03 -0.21352E-02 -0.50947E-03 5617 0.164E-01 0.957E-01 RMM: 6 -0.329378903883E+03 0.30074E-02 -0.14167E-03 4264 0.108E-01 0.569E-01 RMM: 7 -0.329380210311E+03 -0.13064E-02 -0.59560E-04 3737 0.624E-02 0.466E-01 RMM: 8 -0.329380494345E+03 -0.28403E-03 -0.30049E-04 3591 0.500E-02 83 F= -.32938049E+03 E0= -.32938049E+03 d E =-.889510E-02 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
一志愿华中师范071000，325求调剂
已经有6人回复
-
070303一志愿西北大学学硕310找调剂
已经有6人回复
-
311求调剂
已经有4人回复
-
0703化学调剂
已经有8人回复
-
070303 总分349求调剂
已经有5人回复
-
309求调剂
已经有5人回复
-
070300化学学硕求调剂
已经有5人回复
-
277材料科学与工程080500求调剂
已经有3人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有5人回复
-
环境工程调剂
已经有4人回复
dhc198611
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 2488.4
- 散金: 20
- 红花: 2
- 帖子: 280
- 在线: 573.9小时
- 虫号: 780288
- 注册: 2009-05-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-01-24 22:15:02
dfjk-123
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 29.9
- 散金: 757
- 帖子: 119
- 在线: 24.5小时
- 虫号: 665485
- 注册: 2008-12-01
- 性别: MM
- 专业: 理论和计算化学
3楼2013-01-25 10:57:13
dhc198611
木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 2488.4
- 散金: 20
- 红花: 2
- 帖子: 280
- 在线: 573.9小时
- 虫号: 780288
- 注册: 2009-05-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2013-01-26 15:12:37
原来是梦8804
木虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 1388.3
- 红花: 1
- 帖子: 51
- 在线: 99.1小时
- 虫号: 1134938
- 注册: 2010-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
franch: 金币+5, 鼓励交流，， 2013-01-28 22:25:57
感谢参与，应助指数 +1
franch: 金币+5, 鼓励交流，， 2013-01-28 22:25:57
|
我来说下自己的理解，不对的地方请大家多多指正。 我在一开始研究到底使用能量收敛还是力的收敛标准的时候测试研究过，发现了这么一个有趣的事情：在OSZICAR大循环后的dE并不一定是当前大循环的收敛能量减去前一步大循环的收敛能量。所以我去研究了一下IBRION=2用的CG算法是怎么回事。 在VASP说明书里http://cms.mpi.univie.ac.at/vasp/vasp/IBRION_2.html，LZ可以去看一下，大概的意思是CG算法经过以下步骤：1 初始结构；2 猜测步骤；3 修正步骤；4 猜测步骤；5 修正步骤...在对猜测步骤进行修正之后，会进行力的计算，不管你是否采用力的收敛标准，（见说明书：iii) After the corrector step the forces and energy are recalculated and it is checked whether the forces contain a significant component parallel to the previous search direction. If this is the case, the line minimization is improved by further corrector steps using a variant of Brent's algorithm[28].）以便决定是否对猜测步骤再次进行修正，也就是说修正步骤可能有多步，这就解决了我自己的问题，就是当前大循环的dE不一定是当前大循环能量减去上一步大循环的能量，因为如果当前是第二次修正的话，dE应该是当前能量减去上上一步的能量（dE计算的是修正步骤的能量与猜测步骤的能量）。当然，这也同时解决了LZ的问题，就是你采用的即使是能量的收敛标准，在对猜测步骤进行修正时，仍要计算力来决定是否要进行第二次修正（前一次修正可能不正确）。 LZ可以用第84步大循环的能量分别减去83和82步的能量，不出意外的话，结果应该是84步的dE=84步的能量减去82步的能量，而不是83步的能量。也就是说，82步是猜测步骤，83步做了修正，然后计算发现83步的修正不正确（即使能量已经收敛，但是力不正确），所以要进行84步的修正，也可能还有85步的修正... |
5楼2013-01-27 03:05:47
原来是梦8804
木虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 1388.3
- 红花: 1
- 帖子: 51
- 在线: 99.1小时
- 虫号: 1134938
- 注册: 2010-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2013-01-27 03:09:10
20