cp2k做一个表面结构的结构优化100多s一个电子步，为什么这么慢

wjybit · 发表于 Post on 2023-12-7 23:15:21

最近在用cp2k做一个表面结构的结构优化，采用对角化算法，SCF收敛很小1E-3，因为1E-6收敛不了，请大佬帮忙看一下是不是我inp文件写的有问题，提前谢谢大佬。

reid · 发表于 Post on 2023-12-7 23:37:28

本帖最后由 reid 于 2023-12-7 23:47 编辑

对你这个模式不太熟，一直用Multiwfn生成输入文件。是不是你的ADDED_MOS 500有点过大？我们都是几十，如30、40。应该用的是对角化吧，150个重原子，又是UKS，肯定是要慢。你怎么没有k点？还有自旋多重度？
使用Multiwfn生成的，使用了对角化，设置了空轨道90，设置了3,3,1的k点。自旋多重度没设，你自己算算。再试试速度，估计也快不了，主要是使用了k点。但是不使用k点也不会收敛。

coord.inp (14.16 KB, 下载次数 Times of downloads: 6)

啊不错的飞过海 · 发表于 Post on 2023-12-8 19:29:20

删掉cholesky inverse，换elpa做对角化试试，一般能比scalapack快一半甚至更多。

wjybit · 发表于 Post on 2023-12-8 20:29:10

reid 发表于 2023-12-7 23:37
对你这个模式不太熟，一直用Multiwfn生成输入文件。是不是你的ADDED_MOS 500有点过大？我们都是几十，如30 ...

对的大佬，是用的对角化，没用k点是因为电脑加上k点内存就会爆，我加上k点和自旋多重度试试。

wjybit · 发表于 Post on 2023-12-8 20:34:01

啊不错的飞过海发表于 2023-12-8 19:29
删掉cholesky inverse，换elpa做对角化试试，一般能比scalapack快一半甚至更多。

谢谢大佬，alph是指的在MIXING中设置吗？类似这种
               &MIXING
      METHOD BROYDEN_MIXING
   ALPHA 0.1
   BETA 1.5
   NBROYDEN 8
   &END MIXING

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

[CP2K] cp2k做一个表面结构的结构优化100多s一个电子步，为什么这么慢