CP2Kaimd过程中scf不收敛，计算达到指定迭代次数后直接卡住了

lalala_小当家

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

sobereva

尝试用更新的CP2K版本
尝试不用ELPA

gog

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

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

吞木木

可能我的处理方式不对，但是我用这种方式解决了我自己的百分之九十九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算法得到一个收敛的波函数，然后用对角化读取已经收敛的波函数，会起到一定的作用

sobereva

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

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

sobereva

当前问题和SCF是否收敛没直接关系，不管SCF是否收敛MD都会照常往下做。这是MD部分代码的问题。就算SCF收敛了，照样可能卡住

gog

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

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

sobereva

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而已

gog

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

赞

1138711019

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

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