| 查看: 1445 | 回复: 4 | |||
[求助]
请教FeO超胞收敛问题
|
|
系统是210原子的3*3*3FeO超胞，里面有一个defect结构，下面是算的单点能，输入如下： PREC = Accurate ENCUT = 600 NBANDS = 1024 NELM = 200 NELMIN = 10 EDIFF = 1E-05 ALGO = Fast ADDGRID =.TRUE. ISPIN = 2 MAGMOM = AFM 设置，由于太长，不列出来，总MAG=0 PSTRESS = 50 AMIX = 0.01 BMIX = 0.0001 LDAU = .TRUE. LDAUTYPE = 2 LDAUL = 2 -1 LDAUU = 5.0 0.0 LDAUJ = 0.0 0.0 LMAXMIX = 4 MAXMIX = 100 IBRION =-1 LORBIT = 11 ISMEAR = 0 (-1-Fermi, 1-Methfessel/Paxton) SIGMA = 0.05 系统能量无法收敛（NELM=200），我发现原因可能是电荷密度无法收敛，rms和rms(c)都不收敛，rms到1E-2就不减少了，rms(c)到1E-1也不减少。 我尝试了多种AMIX，BMIX组合，包括默认的，线性的（AMIX=0.1～0.2，BMIX=0.0001），也调整了MAXMIX。收敛依然不好，请大家指教，谢谢。 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP第一性原理计算 |
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有11人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有11人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
yuraining
木虫 (小有名气)
- 应助: 17 (小学生)
- 金币: 2187.2
- 帖子: 56
- 在线: 259.3小时
- 虫号: 1503708
- 注册: 2011-11-22
- 专业: 凝聚态物性 II ：电子结构
2楼2013-03-06 08:35:49
3楼2013-03-06 11:11:47
niuxiaoqi721
铜虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 138.7
- 散金: 2
- 红花: 1
- 帖子: 222
- 在线: 105小时
- 虫号: 1226906
- 注册: 2011-03-09
- 性别: GG
- 专业: 理论和计算化学
4楼2013-05-01 16:44:11
【答案】应助回帖
|
不知你的问题解决了没？ 关于收敛，我觉得调混合参数挺麻烦的，而且无法从实质上解决问题，给你几个方案，你可以挨个试一下：1，试试其他类型的PAW，比如PW91，PBE，PBEsol，甚至是LDA，收敛的情况和PAW类型有关系的；2、检查一下你的磁结构，是反铁磁，铁磁？反铁磁是什么类型的呢？层内铁磁，层内反铁磁？初始磁结构给对了，才能计算；3、另外，你加压了，有没有试过不加压的情况收敛否？如果不加压可以收敛，你不妨一点一点加压，比如先算不加压的，然后10,20，。。。4、可以先用一个小的截断能计算，如果收敛，再读入旧的波函数，用新的截断能计算；5、设定ISYM=-1，不限制对称性计算，看看能否好转；6、初始磁矩可以大点给出，比如Fe可以给出磁矩为9，防止体系出现低自旋！ 你可以先试试1、2、3、5 |
5楼2013-05-02 03:24:54
