计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › CASTEP对Ni晶胞进行几何优化时，能量振荡不收敛

CASTEP对Ni晶胞进行几何优化时，能量振荡不收敛

查看数: 12678 | 评论数: 4 | 收藏 Add to favorites 5

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

kcd

发布时间: 2021-3-24 17:49

正文摘要:

请问各位老师，我在使用castep优化软件库里的Ni晶胞的时候，总是无法收敛，能量总是反复振荡的形式的，勾选了自旋极化，initial spin 设置为2，是我哪设置不对吗？

回复 Reply

卡开发发 发表于 Post on 2021-3-27 11:42:43

喵星大佬发表于 2021-3-27 02:09
不知道CASTEP里面的超软赝势质量怎么样，如果按照QE的赝势看，SSSP上pslibrary.1.0.0 US算Ni至少要到50-6 ...

QE的PSLib有些赝势标注的是50～60Ry，但相同的赝势有些体系实测40~45Ry就能收敛。PSlib给出的参考动能截断似乎是没啥价值的，建议好好做一下收敛性测试。反而CASTEP的默认动能截断偏低，比较可能的原因是这个参考值用的是相对能量的收敛而不是绝对能量的收敛，但是对于构型优化来说绝对能量收敛意义要更好一些（更严格应该看力收敛的情况）。另外CASTEP如果用的新版本就不建议用老的USPP了，尽量用OTFG，性能和精度能差出来一节。

喵星大佬 发表于 Post on 2021-3-27 02:09:10

本帖最后由喵星大佬于 2021-3-27 02:43 编辑

kcd 发表于 2021-3-25 09:17
谢谢，我把截断能提高到300 eV，SCF收敛限也提高了，但是能量还是反复振荡

不知道CASTEP里面的超软赝势质量怎么样，如果按照QE的赝势看，SSSP上pslibrary.1.0.0 US算Ni至少要到50-60Ry截断能之后才比较稳定，折算到eV大概是680-820，当然pslibrary.1.0.0的赝势不少元素是比较差劲的，如果算CASTEP里面的OTFG超软质量比较好的话至少也是500-600eV的截断能。

尝试了一下600eV截断能，k点888，开自旋极化默认初猜、开金属(smearing 0.1eV)其他默认，全优化，SCF收敛限-6，3步就收敛了
试了一下300eV截断能，没有问题能收敛，4步即可，然后优化出来的Ni-Ni距离偏长，进一步降低截断能晶格常数会更奇怪，感觉500eV不能再低了，但是并没出现收敛不了的情况
按照你给的参数，把charge mixing amplitude从0.5改到0.1会出现不收敛，用默认的0.5就好了

kcd 发表于 Post on 2021-3-25 09:17:34

喵星大佬发表于 2021-3-24 20:14
截断能提高，而且做优化把SCF收敛限设严一点

谢谢，我把截断能提高到300 eV，SCF收敛限也提高了，但是能量还是反复振荡

喵星大佬 发表于 Post on 2021-3-24 20:14:00

截断能提高，而且做优化把SCF收敛限设严一点

评分 Rate

参与人数 Participants 1	eV +2	收起理由 Reason
卡开发发	+ 2	我很赞同，精度过低会引入force和stress噪.

查看全部评分 View all ratings

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register