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

kcd

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

喵星大佬

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

kcd

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

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

喵星大佬

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就好了

卡开发发

喵星大佬 发表于 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，性能和精度能差出来一节。