计算化学公社

标题: cp2k做一个表面结构的结构优化100多s一个电子步，为什么这么慢 [打印本页]

作者
Author: wjybit 时间: 2023-12-7 23:15
标题: cp2k做一个表面结构的结构优化100多s一个电子步，为什么这么慢
最近在用cp2k做一个表面结构的结构优化，采用对角化算法，SCF收敛很小1E-3，因为1E-6收敛不了，请大佬帮忙看一下是不是我inp文件写的有问题，提前谢谢大佬。

作者
Author: reid 时间: 2023-12-7 23:37
本帖最后由 reid 于 2023-12-7 23:47 编辑

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

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

作者
Author: wjybit 时间: 2023-12-8 20:29

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

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

作者
Author: wjybit 时间: 2023-12-8 20:34

啊不错的飞过海发表于 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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)