| 查看: 2749 | 回复: 9 | ||||
[交流]
【求助】对于难收敛的体系 已有6人参与
|
| 各位，对于vasp优化难收敛的体系，可不可以这样做，先用低的收敛标准使其收敛，然后读入波函数和电荷再采用高的收敛标准进行计算，可行否？多谢指教！ |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP优化不收敛 |
» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关高斯算晶体的
已经有7人回复
- vasp计算的半导体结构，bands带隙比dos里面的大，何故？
已经有33人回复
- 请各位大侠看下高斯出错
已经有3人回复
- opt=TS 循环步骤问题
已经有7人回复
- 【求助】K点取得越大越不收敛，什么原因！
已经有20人回复
- 【求助】结构优化不收敛
已经有9人回复
- 【求助】vasp体系过大“优化”难以收剑
已经有7人回复
- 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？
已经有24人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
2楼2010-09-19 12:47:42
fangyuan5687
捐助贵宾 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1092.5
- 散金: 2263
- 红花: 2
- 帖子: 481
- 在线: 93.3小时
- 虫号: 655357
- 注册: 2008-11-16
- 专业: 能源化工
3楼2010-09-19 17:12:57
chimegreen
木虫 (正式写手)
Dr.
- 应助: 8 (幼儿园)
- 金币: 3659.1
- 散金: 44
- 红花: 4
- 帖子: 521
- 在线: 165.3小时
- 虫号: 177844
- 注册: 2006-02-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-09-25 09:19:43
5楼2010-09-25 10:59:21
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
6楼2010-09-25 11:06:42
chimegreen
木虫 (正式写手)
Dr.
- 应助: 8 (幼儿园)
- 金币: 3659.1
- 散金: 44
- 红花: 4
- 帖子: 521
- 在线: 165.3小时
- 虫号: 177844
- 注册: 2006-02-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+3):谢谢。欢迎总结评论式的回复！ 2010-09-25 16:42:24
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+3):谢谢。欢迎总结评论式的回复！ 2010-09-25 16:42:24
|
如果已经开始弛豫了好几步了，就copy过来，如果第一步都没过去，那就不必了。 收敛参数的调整量子化学论坛以前有个老帖 你自己看看有没有你能用的 对于一些磁性体系、镧系和锕系元素及相关化合物的静态计算（电子迭代），经常会遇到“难收敛”的问题。 下面给出几个相关Flag及设置方法：（有些话从说明书上摘录的，不大好翻译，就copy了，请见谅） 1、LMAXMIX Default: LMAXMIX = 2 An additional flag controls up to which l quantum number the onsite PAW charge densities are passed through the charge density mixer. Higher l-quantum numbers are usually not handled by the mixer. In order to obtain fast convergence to the groundstate, you can try the following setting: LMAXMIX = 4 for d elements LMAXMIX = 6 for f elements 这个FLAG对于含d电子和f电子的体系是非常重要的，很大一部分体系的收敛问题可以通过设置合适的LMAXMIX值来解决。 2、ALGO, IALGO, LDIAG If the self-consistency loop does not converge within 40 steps, it will probably not converge at all. In this case you should reconsider the tags IALGO, LDIAG, and the mixing-parameters. 这是说明书上的建议。 一般情况下，或使用IALGO=48时遇到收敛问题的话，可以考虑设IALGO为38（4.5以前的版本可设为8），或设置ALGO=Normal or Fast (in VASP.4.5 and later versions)。 3、NELMDL NELMDL gives the number of non-selfconsistent steps at the beginning; if one initializes the wave functions randomly the initial wave functions are far from anything reasonable. The resulting charge density is also 'nonsense'. Therefore it makes sense to keep the initial Hamiltonian, which corresponds to the superposition of atomic charge densities, fixed during the first few steps. Choosing a 'delay' for starting the charge density update becomes essential in all cases where the SC-convergence is very bad (e.g. surfaces or molecules/clusters chains). Without setting a delay VASP will probably not converge or at least the convergence speed is slowed down. NELMDL might be positive or negative. A positive number means that a delay is applied after each ionic movement -- in general not a convenient option. A negative value results in a delay only for the start-configuration. 4、mixing-parameters 对于一些难收敛的体系，可以使用“linear mixing”,具体详见VASP说明书中的“Mixing-tags”。 For an initial linear mixing (BMIX ~ 0) an optimal setting for A(AMIX) can be found easily by setting Aopt=Acurrent*Γmean. For the Kerker scheme either A or q0(i.e. AMIX or BMIX) can be optimized, but we recommend to change only BMIX and keep AMIX fixed (you must decrease BMIX if the mean eigenvalue is larger than one, and increase BMIX if the mean eigenvalue is smaller than one). 尽管VASP说明书中给出了调节AMIX和BMIX的一些较为明确的建议，但是实际去调节的时候，还是挺难的，但原则上说，是可以通过调节这两个Flag来使得收敛问题得以解决的，只是得有耐心。 5、kmesh, SIGMA 收敛问题还跟kmesh及SIGMA(当使用ISMEAR不等于-5 和-4时)的设置有关。要达到同样的精度，较小的SIGMA则需要较大的kmesh；而且，当SIGMA较小时，若kpoints不够多，也会出现难收敛的情况。 P.S. 对于不同系统的计算，问题的原因不一定一样，因而可能解决之一问题的方法也不一定会一样。其实以上大部分都来自VASP的说明书，只是我在遇到收敛问题时试过这些方法. |
7楼2010-09-25 14:36:48
8楼2010-09-25 15:10:42
wuli8
荣誉版主 (知名作家)
…………
- 1ST强帖: 2
- 应助: 35 (小学生)
- 贵宾: 12.924
- 金币: 20192.4
- 散金: 15888
- 红花: 88
- 沙发: 4
- 帖子: 7840
- 在线: 1114.6小时
- 虫号: 465889
- 注册: 2007-11-23
- 专业: 物理学I
- 管辖: 计算模拟
9楼2013-08-22 16:07:58
10楼2013-08-28 10:06:28
