计算化学公社

标题: CP2K计算REPEAT电荷时，SCF不收敛该如何解决 [打印本页]

作者
Author: a9471163 时间: 2021-10-2 23:11
标题: CP2K计算REPEAT电荷时，SCF不收敛该如何解决
本帖最后由 a9471163 于 2021-10-2 23:51 编辑

大家好，我用PBE+DZVP-MOLOPT-SR-GTH计算伊利石的REPEAT电荷，出现SCF不收敛的问题，我将.out文件和.inp文件，以及伊利石的.cif上传了。

.out中有警告“WARNING in qs_scf.F:576 :: SCF run NOT converged”，想请教各位，该警告是否影响计算结果呢，如何改进有望达到SCF收敛?

计算步骤：参考了http://sobereva.com/587中第3.7节，即首先生成.gjf，随后载入Multiwfn，
-4  //输出原子电荷
7  //REPEAT
8  //设置k点
1,1,1
0  //产生输入文件

（为了计算速度，k点选取1,1,1，而晶胞x,y方向的尺寸与k点乘积小于10，请问此时k点是否太小了）

(, 下载次数 Times of downloads: 12)
(, 下载次数 Times of downloads: 51)
(, 下载次数 Times of downloads: 21)

作者
Author: sobereva 时间: 2021-10-2 23:18
较长的输入输出文件别直接贴，应当作为附件上传

先把SCF不收敛问题解决（更何况当前距离收敛很遥远），否则算的原子电荷没有意义

作者
Author: a9471163 时间: 2021-10-2 23:25
本帖最后由 a9471163 于 2021-10-2 23:47 编辑

sobereva 发表于 2021-10-2 23:18
较长的输入输出文件别直接贴，应当作为附件上传

先把SCF不收敛问题解决（更何况当前距离收敛很遥远）， ...

非常感谢您的提醒，下次发帖我会注意的，本帖已编辑以便缩短篇幅。
我记得您发过专门的文章，以便解决gaussian计算时SCF不收敛，论坛有人写过解决orca计算时SCF不收敛的问题

不过cp2k的关键词我目前不懂，那些关键词跟gaussian、orca也不一样，目前看着inp输入文件，不太清楚加什么语句可能奏效

想请问一下可以在Multiwfn创建inp时，通过Multiwfn对inp文件做些调整吗（我现在是会修改inp，增加SCF的迭代次数了，不过我记得SCF盲目增加迭代次数，不一定是好的方式；收敛限可以放宽，不过我个人的经验不足，不知道收敛限最低可以用到多少，担心我自己改完以后收敛限太宽松了）

作者
Author: a9471163 时间: 2021-10-9 19:19
问题解决了，需要开启smearing，这是SCF收敛的关键

作者
Author: Siling_S 时间: 2021-12-24 10:29

a9471163 发表于 2021-10-9 19:19
问题解决了，需要开启smearing，这是SCF收敛的关键

感激感激，我试了同样的做法也收敛了~~

作者
Author: han- 时间: 2022-12-5 08:55

a9471163 发表于 2021-10-9 19:19
问题解决了，需要开启smearing，这是SCF收敛的关键

老师请问cp2k如何开启smearing？

作者
Author: zongxupeng 时间: 2023-1-6 17:13

han- 发表于 2022-12-5 08:55
老师请问cp2k如何开启smearing？

   &SMEAR
      METHOD FERMI_DIRAC
      ELECTRONIC_TEMPERATURE 300
   &END SMEAR

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