求助：此SCF不收敛情况是否适用于增大SCF迭代次数上限

Mr.zhen · 发表于 Post on 2023-9-5 10:42:40

本帖最后由 Mr.zhen 于 2023-9-5 10:46 编辑

在materialsproject上下载的CIF（Symmetrized）版本使用cp2k进行结构优化仍然长时间不收敛，打开out文件发现有些SCF循环Convergence能够从第1次的2.0降到第128次的0.0001，但是新的SCF循环开始后第1次仍然从2.0开始。本人新接触cp2k计算，不知这种情况是否正常，我是应该进行一些修改还是再等等，原以为从materialsproject上下载的结构不用结构优化或者很快就能优化完成。我有意将最大SCF迭代次数从128增加到更高的数字，但是浏览帖子后发现这似乎是最不推荐的做法，所以希望各位老师能够指点一二。计算的输入输出文件如下，out文件有些大压缩为了rar文件，cp2k的版本为2022.2的popt版。

wzkchem5 · 发表于 Post on 2023-9-5 15:15:58

应该是新的scf计算没有读初猜波函数，或者读取方式不对吧

Mr.zhen · 发表于 Post on 2023-9-5 15:27:01

wzkchem5 发表于 2023-9-5 15:15
应该是新的scf计算没有读初猜波函数，或者读取方式不对吧

那么我应当对输入文件进行一些修改吗？

ahxb · 发表于 Post on 2023-9-5 21:22:19

在输出文件中有提示：

Extrapolation method: initial_guess

复制代码

因此每次SCF都重新猜测波函数然后自洽。其实Multiwfn产生的输出文件已经提供了一些可调整部分，比如&QS字段里的EXTRAPOLATION行被注释掉了，可以尝试取消注释，使用上一步自洽产生的密度。也有其它选项，比如USE_PREV_RHO_R和LINEAR_P等，具体含义去翻cp2k manual吧，通常对SCF收敛步数减少有一定作用。

Mr.zhen · 发表于 Post on 2023-9-6 16:07:05

ahxb 发表于 2023-9-5 21:22
在输出文件中有提示：

因此每次SCF都重新猜测波函数然后自洽。其实Multiwfn产生的输出文件已经提供了一 ...

谢谢您，我去翻一下手册试试

Mr.zhen · 发表于 Post on 2023-9-6 16:16:47

ahxb 发表于 2023-9-5 21:22
在输出文件中有提示：

因此每次SCF都重新猜测波函数然后自洽。其实Multiwfn产生的输出文件已经提供了一 ...

不过EXTRAPOLATION USE_PREV_P看起来是明显有利于收敛的吧，但是sob老师却将其在Multiwfn的输入文件生成中默认注释掉了，是说明此种方法也有不适用的情况或缺点吗

sobereva · 发表于 Post on 2023-9-8 01:17:13

Mr.zhen 发表于 2023-9-6 16:16
不过EXTRAPOLATION USE_PREV_P看起来是明显有利于收敛的吧，但是sob老师却将其在Multiwfn的输入文件生成 ...

考虑k点时，USE_PREV_P对于>99%的时候都是建议用的，但个别情况可能引起报错，所以没作为默认设置

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助：此SCF不收敛情况是否适用于增大SCF迭代次数上限