计算化学公社

标题: CP2Kaimd过程中scf不收敛，计算达到指定迭代次数后直接卡住了 [打印本页]

作者
Author: lalala_小当家 时间: 2023-4-1 09:56
标题: CP2Kaimd过程中scf不收敛，计算达到指定迭代次数后直接卡住了
各位老师，我在做aimd模拟过程中我的scf始终不收敛（c尝试改了许多参数）此外，计算达到指定迭代次数后直接卡住了，out文件不在有任何信息输出，但计算没有终止，想请大家帮我参谋下这两个问题。

作者
Author: sobereva 时间: 2023-4-1 15:04
尝试用更新的CP2K版本
尝试不用ELPA

作者
Author: gog 时间: 2023-4-2 11:43

sobereva 发表于 2023-4-1 15:04
尝试用更新的CP2K版本
尝试不用ELPA

不用elpa，计算速度会变慢吧。这个库还是很重要的。还是尽可能编译好软件。

作者
Author: 吞木木 时间: 2023-4-2 16:23
本帖最后由吞木木于 2023-4-2 16:26 编辑

可能我的处理方式不对，但是我用这种方式解决了我自己的百分之九十九scf不收敛的问题，谨慎使用。在你的out文件里，可以看到如下信息：
95 Broy./Diag. 0.10E+00 24.3    0.12546865 -32556.7677472373 -6.13E+00
96 Broy./Diag. 0.10E+00 24.3    0.03637298 -32552.3503186206  4.42E+00
97 Broy./Diag. 0.10E+00 24.2    0.06718901 -32555.2395700660 -2.89E+00
98 Broy./Diag. 0.10E+00 24.2    0.08131985 -32565.1681747740 -9.93E+00
99 Broy./Diag. 0.10E+00 24.2    0.11434304 -32564.4504850039  7.18E-01
100 Broy./Diag. 0.10E+00 24.5    0.04457700 -32554.3426249676  1.01E+01

也就是判断scf收敛的这一列参数一直处于一种震荡的状态，现在增加scf的最大循环圈数已经是无效的了。
在输入文件里，scf的收敛算法你选择用的是对角化，我个人觉得对于稍微大的体系，你不需要输出体系的能级相关信息的时候，选用OT算法（关键词MINIMIZER选择CG会更快的加速收敛）能很快地达到scf收敛的状态
或者你必须使用对角化算法的时候，可以先用OT算法得到一个收敛的波函数，然后用对角化读取已经收敛的波函数，会起到一定的作用

作者
Author: sobereva 时间: 2023-4-2 16:53

gog 发表于 2023-4-2 11:43
不用elpa，计算速度会变慢吧。这个库还是很重要的。还是尽可能编译好软件。

若不是跨节点大规模并行，ELPA影响甚微。个人双路服务器上跑，基本体现不出来和scalapack的性能差异

作者
Author: sobereva 时间: 2023-4-2 16:54
当前问题和SCF是否收敛没直接关系，不管SCF是否收敛MD都会照常往下做。这是MD部分代码的问题。就算SCF收敛了，照样可能卡住

作者
Author: gog 时间: 2023-4-3 12:28

sobereva 发表于 2023-4-2 16:53
若不是跨节点大规模并行，ELPA影响甚微。个人双路服务器上跑，基本体现不出来和scalapack的性能差异

ELPA是算电子结构时候用的。

作者
Author: sobereva 时间: 2023-4-4 05:38

gog 发表于 2023-4-3 12:28
ELPA是算电子结构时候用的。

不仅仅和电子结构部分计算有关
ELPA是大规模并行的对角化的库，和对角化有关的情况都可能涉及。例如BFGS几何优化过程用ELPA可能卡住已经是已知事件，如https://groups.google.com/g/cp2k/c/BAvvW_qGG2I/m/OoskF9ePBAAJ、https://groups.google.com/g/cp2k/c/qKcf49JFj0w。
虽然MD表面上并不直接涉及到对角化，但不排除一些细节过程上利用了ELPA做对角化。作为已知可能导致卡住的因素之一，至少值得尝试把ELPA改成scalapack一试，只不过输入文件&GLOBAL里加一行PREFERRED_DIAG_LIBRARY SL而已

作者
Author: gog 时间: 2023-4-4 08:26

sobereva 发表于 2023-4-4 05:38
不仅仅和电子结构部分计算有关
ELPA是大规模并行的对角化的库，和对角化有关的情况都可能涉及。例如BFGS ...

作者
Author: 1138711019 时间: 2023-4-4 15:21

sobereva 发表于 2023-4-4 05:38
不仅仅和电子结构部分计算有关
ELPA是大规模并行的对角化的库，和对角化有关的情况都可能涉及。例如BFGS ...

确实，我之前那个问题就是ELPA导致的，就算SCF收敛也是会卡住。现在换成SL了，感觉变慢了不少

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