求助：cp2k优化Ni（111）面scf不收敛问题

时光的邂逅

最近在使用cp2k计算Ni（111）面吸附Cl离子的吸附能时，在结构优化Ni（111）面时scf收敛一直震荡，一直不收敛，应该是有些参数没有设置对，以下是我的输入文件，望各位老师不吝赐教！！！

ddddnight

Eudaimonia 发表于 2022-3-15 15:38
搞不懂为什么给了建议根本不听

我用你的cif使用multiwfn重新生成inp文件，设置XY方向周期性，k点设置为3 ...

看了老师您的输入文件，我有两点疑问1. 为什么inp文件里面没有开自旋极化（UKS）呢，α电子和β默认完全相等吗？
2. 卢老师说ANALYTIC很垃圾，为什么不用MT呢
我试了一下，开了自旋极化和换了MT，怎么换ALPHA和NBROYDEN都收敛不了

lhx-xwzj

丁越 发表于 2023-12-26 08:52
这种体系没必要去折腾收敛性了，换QE去做不存在这种比较稀奇古怪的收敛性问题

好的，谢谢丁老师，我这就去研究一下QE

丁越

lhx-xwzj 发表于 2023-12-25 19:09
丁老师您好，最近我用cp2k 2023.1版本，也遇到了Ni111结构优化不收敛的问题，根据帖子里的方案都进行了尝 ...

这种体系没必要去折腾收敛性了，换QE去做不存在这种比较稀奇古怪的收敛性问题

lhx-xwzj

丁越 发表于 2022-2-24 14:10
你得继续改小才行，另外你的结构不合理，用Gamma点计算你的晶胞太小了，而且还是d区金属，需要提前测试一 ...

丁老师您好，最近我用cp2k 2023.1版本，也遇到了Ni111结构优化不收敛的问题，根据帖子里的方案都进行了尝试，一百多个镍的表面怎么都不收敛，实在没办法了用您这个能收敛的输入文件直接计算，发现300次后还是不收敛，您那个是8.2版本，怎么版本越新越不收敛了呢？另外还想问问CP2K优化过渡金属是不是就是容易不收敛啊，您有什么建议吗？下面是输入与输出文件

Janson_Shi

Eudaimonia 发表于 2022-3-15 15:38
搞不懂为什么给了建议根本不听

我用你的cif使用multiwfn重新生成inp文件，设置XY方向周期性，k点设置为3 ...

我最近遇到相似的问题，按照这个帖子调试了很多参数，但就是收敛不了，麻烦指点一下http://bbs.keinsci.com/thread-38331-1-1.html FeP211 slab scf不收敛的问题，谢谢

Janson_Shi

卡开发发 发表于 2022-2-25 16:43
是应该调整alpha值，必要的情况也可以去调整下beta（虽然我个人经验绝大部分情况不需要）。

我最近遇到相似的问题，按照这个帖子调试了很多参数，但就是收敛不了，麻烦指点一下http://bbs.keinsci.com/thread-38331-1-1.html FeP211 slab scf不收敛的问题，谢谢

Eudaimonia

yy97500 发表于 2022-4-16 16:56
老师，您好。我也是遇到Ni的表面层优化难以收敛的问题。学习了您写的的Ni-111-331-spin3.inp，不过发现里 ...

这个情况我也不确定，你可以尝试对比一下加不加UKS的区别

yy97500

Eudaimonia 发表于 2022-3-15 15:38
搞不懂为什么给了建议根本不听

我用你的cif使用multiwfn重新生成inp文件，设置XY方向周期性，k点设置为3 ...

老师，您好。我也是遇到Ni的表面层优化难以收敛的问题。学习了您写的的Ni-111-331-spin3.inp，不过发现里面写的是    MULTIPLICITY    3 #Spin multiplicity，但是UKS 没有描述，UKS 默认是关闭的，这种情况下MULTIPLICITY    3 #Spin multiplicity能够生效吗？

时光的邂逅

Eudaimonia 发表于 2022-3-15 15:38
搞不懂为什么给了建议根本不听

我用你的cif使用multiwfn重新生成inp文件，设置XY方向周期性，k点设置为3 ...

非常感谢老师的热心解答，这里可能要跟老师您道个歉，当时我尝试过老师您给的建议，但是由于我电脑中安装的cp2k有点问题，只要一加k点就报错，所以只能用111算 ，然后scf一直不收敛，因此后面就没有用331，551做测试。然后当时这个计算，导师又催的比较急，所以就换了qe算。真心感觉自己有点对不起老师您的热心解答，希望老师您给的建议能为以后遇到相同问题的人提供帮助。

Eudaimonia

搞不懂为什么给了建议根本不听

我用你的cif使用multiwfn重新生成inp文件，设置XY方向周期性，k点设置为331，打开smearing，cp2k版本8.2
之后手动修改inp文件，CUTOFF 600，REL_CUTOFF 55和你输入文件精度保持一致，修改ALPHA = 0.2 NBROYDEN =16 smearing设置为 50 50，Ni磁矩设置为3
根据卡开发发老师的建议MULTIPLICITY    设置为73(体系36个Ni)
尝试了k点为331,551，均收敛，331还测试了MULTIPLICITY    3情况下的SCF收敛情况，与73的情况下基本一致！
k点为111时128步SCF后无法收敛但存在收敛趋势，根据我之前的测试 [url=加速SCF收敛可以尝试的一些手段 http://bbs.keinsci.com/forum.php ... 27630&fromuid=26742 (出处: 计算化学公社)]加速SCF收敛可以尝试的一些手段[/url]  修改NBROYDEN =24 ALPHA = 0.05 经过35步SCF后也收敛

jiangning198511

时光的邂逅 发表于 2022-3-15 00:41
没有呢，我放弃了，用qe算了

楼主 qe收敛效果如何？

时光的邂逅

bblovelp 发表于 2022-3-13 11:48
您好 ，请问下你的Ni(111)几何优化收敛了吗？

没有呢，我放弃了，用qe算了

bblovelp

您好 ，请问下你的Ni(111)几何优化收敛了吗？

卡开发发

丁越 发表于 2022-2-25 16:40
卡开发发老师，增大到这个值后SCF整不收敛啊，delta_E change 在10^-2到0之间跳跃。我去Google group上搜 ...

是应该调整alpha值，必要的情况也可以去调整下beta（虽然我个人经验绝大部分情况不需要）。