| 查看: 1365 | 回复: 5 | ||||
| 当前主题已经存档。 | ||||
[交流]
【求助】磁性计算时的收敛问题
|
||||
|
用VASP计算过渡金属表面吸附磁性分子的问题，体系比较大，总共30个原子，过渡金属和分子都要考虑磁性。在静态计算时，即只计算一个离子步，NSW=0，想得到体系能量，发现此体系很难收敛，电子步设置到200步，也不行。另外即使先计算非磁的波函数，也不收敛。 请问在INCAR文件里怎样设置参数，才能更快的收敛呢？ 谢谢！欢迎大家讨论、指教。 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有11人回复
-
限项规定
已经有9人回复
-
免疫学博士有名额，速联系
已经有4人回复
-
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有5人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
cbFeng200881
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 262.8
- 散金: 3414
- 红花: 2
- 沙发: 12
- 帖子: 2772
- 在线: 699.4小时
- 虫号: 524662
- 注册: 2008-03-14
- 专业: 凝聚态物性 II ：电子结构
2楼2009-10-23 23:34:18
3楼2009-10-23 23:44:03
这是以前论坛里就有的，帖出来你看看。
★ ★ ★ ★ ★
wlscxf(金币+3,VIP+0):谢谢，期待更多的讨论。 10-26 23:03
zxzj05(金币+2,VIP+0):给个红包，谢谢回帖交流 11-4 21:14
wlscxf(金币+3,VIP+0):谢谢，期待更多的讨论。 10-26 23:03
zxzj05(金币+2,VIP+0):给个红包，谢谢回帖交流 11-4 21:14
|
[VASP]电子迭代的收敛问题 提交者: avan 日期: 2006/11/5 13:16 阅读: 216 作者:avan 来源:量子化学网 摘要:对于一些磁性体系、镧系和锕系元素及相关化合物的静态计算（电子迭代），经常会遇到“难收敛”的问题。 下面给出几个相关Flag及设置方法：（有些话从说明书上摘录的，不大好翻译，就copy了，请见谅） Tag: VASP 软件经验 电子迭代 疑难解答 [作者注] 目前这文章看起来有点乱，过段时间我会再重新整理。 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的说明书，只是我在遇到收敛问题时试过这些方法. 写得比较仓促，以后再修改及补充。 希望有这方面经验的朋友也把自己的见解提出来一起讨论。 |
4楼2009-10-26 20:57:00
5楼2009-10-26 23:05:05
kongjzh
木虫 (正式写手)
- 应助: 139 (高中生)
- 贵宾: 0.11
- 金币: 1935.7
- 红花: 4
- 帖子: 657
- 在线: 127.3小时
- 虫号: 618224
- 注册: 2008-10-06
- 性别: GG
- 专业: 无机合成和制备化学
6楼2009-11-04 16:14:27