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

卡开发发 · 发表于 Post on 2022-2-24 22:16:10

丁越发表于 2022-2-24 21:23
之前看到sob老师说结合smear时默认不会保持alpha和beta电子数的差值，设置自旋多重度没有意义了

尝试 ...

可能还不够大，这个多重度指的是per cell的，MULTIPLICITY=2*Natoms+1试试。

丁越 · 发表于 Post on 2022-2-25 16:40:34

本帖最后由丁越于 2022-2-25 16:42 编辑

卡开发发发表于 2022-2-24 22:16
可能还不够大，这个多重度指的是per cell的，MULTIPLICITY=2*Natoms+1试试。

卡开发发老师，增大到这个值后SCF整不收敛啊，delta_E change 在10^-2到0之间跳跃。我去Google group上搜了一下，有一个邮件就是问NI111表面优化不收敛，给的建议也是减小ALPHA值。现在我也不知道怎么搞了

卡开发发 · 发表于 Post on 2022-2-25 16:43:37

本帖最后由卡开发发于 2022-2-25 16:54 编辑

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

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

bblovelp · 发表于 Post on 2022-3-13 11:48:43

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

时光的邂逅 · 发表于 Post on 2022-3-15 00:41:37

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

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

jiangning198511 · 发表于 Post on 2022-3-15 09:21:11

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

楼主 qe收敛效果如何？

Eudaimonia · 发表于 Post on 2022-3-15 15:38:54

本帖最后由 Eudaimonia 于 2022-3-15 15:39 编辑

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

我用你的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后也收敛

时光的邂逅 · 发表于 Post on 2022-3-16 23:51:01

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

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

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

yy97500 · 发表于 Post on 2022-4-16 16:56:11

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 · 发表于 Post on 2022-4-19 16:26:33

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

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

Janson_Shi · 发表于 Post on 2023-7-2 16:02:02

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

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

Janson_Shi · 发表于 Post on 2023-7-2 16:02:39

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

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

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

lhx-xwzj · 发表于 Post on 2023-12-25 19:09:50

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

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

丁越 · 发表于 Post on 2023-12-26 08:52:31

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

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

lhx-xwzj · 发表于 Post on 2023-12-26 12:46:46

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

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

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

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

评分 Rate