计算化学公社

标题: 求助：此SCF不收敛情况是否适用于增大SCF迭代次数上限 [打印本页]

作者
Author: Mr.zhen 时间: 2023-9-5 10:42
标题: 求助：此SCF不收敛情况是否适用于增大SCF迭代次数上限
本帖最后由 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版。

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

作者
Author: Mr.zhen 时间: 2023-9-5 15:27

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

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

作者
Author: ahxb 时间: 2023-9-5 21:22
在输出文件中有提示：

Extrapolation method: initial_guess

复制代码

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

作者
Author: Mr.zhen 时间: 2023-9-6 16:07

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

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

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

作者
Author: Mr.zhen 时间: 2023-9-6 16:16

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

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

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

作者
Author: sobereva 时间: 2023-9-8 01:17

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

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

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