|
|
本帖最后由 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？
|
|
发表于 Post on 2024-12-5 11:16:38
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
发表于 Post on 2024-12-5 17:15:57
楼主 Author
。