| 查看: 2108 | 回复: 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 |
» 猜你喜欢
职称评审没过，求安慰
已经有41人回复
-
回收溶剂求助
已经有7人回复
-
硝基苯如何除去
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
-
三无产品还有机会吗
已经有6人回复
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