[CP2K]求助：结构优化时收敛标准判断及目前优化方式是否合理

本帖最后由 JillW 于 2024-12-5 11:16 编辑

各位老师好，最近刚刚开始接触CP2K，用multiwfn生成了输入文件并想要优化晶胞，因此打算首先加氢对体系进行结构优化，收敛后再变胞优化。目前结构优化发现总是达到两个yes（均为step size）后计算就自动停止了（如图1），上网查阅后尝试了一些解决方案：（图片太大了就附在了帖子最后，辛苦大家翻阅一下~）

①得知如果SCF收敛标准设置的宽泛其实不利于结构优化的收敛，因此从1*10E-5改回了1*10E-6，发现确实更接近收敛标准（图2），但是仍然没有达到4个yes，能量变化如图3，轨迹如附件所示；

②改用BFGS方法续算，但是收敛过程会变得很慢，甚至算了16、7个小时一个yes都没有；

③目前自己调了WANTED_PROJ_GRADIENT和WANTED_REL_F_ERROR之后续算，均由默认的E-16改成了E-20，目前还在计算过程中。

④如果还是没法达到4个yes，参考sob老师在其他帖子的建议可以用BFGS算法，然后减小步长，或者使用CG算法。

不知道自己上述的解决方案是否合理，抑或是其实目前的优化已经可以满足要求（因为看到部分帖子达到3个yes也可以认为是满足要求），并且还有些疑问如下，希望能得到各位大佬的不吝赐教

①CP2K结构优化时的收敛标准应该如何判断？是必须达到4个yes还是满足其中2~3个，其他接近收敛标准即可？是否是接近收敛标准又如何界定？

②用.restart文件进行续算时是从头开始优化了吗？如果不是，为何之前优化到了2个yes，续算之后又都是no并且远离了收敛标准呢？

③在优化过程中，4个收敛标准总是在不时震荡，例如RMS gradient逐渐接近收敛标准后又突然远离，然后再次接近，最后计算停止时可能还不如过程中的某一步接近收敛标准，请问这是正常的吗？

④晶胞优化时是否应该固定角度和对称性，即KEEP_ANGLES和KEEP_SYMMETRY？

202412051107548990..png

(76.36 KB, 下载次数 Times of downloads: 90)

下载附件 Download

2024-12-5 11:07 上传 Uploaded

图3

202412051106489349..png

(288.48 KB, 下载次数 Times of downloads: 64)

下载附件 Download

2024-12-5 11:06 上传 Uploaded

图2

202412051106239073..png

(306.05 KB, 下载次数 Times of downloads: 87)

下载附件 Download

2024-12-5 11:06 上传 Uploaded

图1

cp2k-pos-1.zip

2024-12-5 11:09 上传 Uploaded

点击下载Click to download

648.55 KB, 下载次数 Times of downloads: 18

xyz文件

cp2k.inp

2024-12-5 11:15 上传 Uploaded

点击下载Click to download

33.83 KB, 下载次数 Times of downloads: 37

最开始的输入文件