| 查看: 2393 | 回复: 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 |
» 猜你喜欢
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
面上的水很深
已经有18人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
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
10