计算化学公社

标题: 求教，CP2K跑AIMD出现中间几步不收敛的情况 [打印本页]

作者
Author: bentty_lee 时间: 2024-7-22 14:33
标题: 求教，CP2K跑AIMD出现中间几步不收敛的情况
各位老师：
刚开始上手使用CP2K，遇到了AIMD的SCF不收敛的情况，请教一下是什么问题，附上输入输出文件。
我的体系是GeSbTe，想跑一个非晶结构，所以需要跑一个2300K的heat的任务，根据参考文献，选择了TZVP-PBE方法，2fs的step，由于GST是窄带隙半导体，所以不能用OT。
体系是首先优化过，然后从300K heat到2300K，由于坐标文件太大了传不上来，截两张图。
刚开始我用的Multiwfn自动默认的SCF128步，发现中间有几步SCF不收敛，而且跑满128导致算的很慢，看了一下out文件，是震荡的趋势，感觉128的SCF其实有点多余了，反正都不收敛了，为了快点看到后续的MD的趋势，所以我杀掉任务，用50步的SCF rerun了，本来以为升温到2300K之后，体系后面会SCF收敛，但是仍然还是会出现偶尔中间有几步SCF不收敛的情况。
ps:我也曾经试过GEN1-xtb方法，改用过DZVP，尝试过改成step 0.1 fs, 都不行。

然后我用同样的输入文件，强制性抽取了一帧2300K的结构和速度，跑了个CSVR下800 TIMECON的cool退火，也是遇到了中间会有几步SCF不收敛的情况，感觉跟heat过程出现的问题一样。

想请问一下各位老师，这种情况的计算结果可信吗？该怎么修改参数？

作者
Author: bentty_lee 时间: 2024-7-22 14:46
更改一下问题，我之前test的时候用的是2fs，SCF不收敛后，改成了附件中小一点的1 fs 的step。

作者
Author: 吞木木 时间: 2024-7-22 18:38
本帖最后由吞木木于 2024-7-22 18:53 编辑

bentty_lee 发表于 2024-7-22 14:46
更改一下问题，我之前test的时候用的是2fs，SCF不收敛后，改成了附件中小一点的1 fs 的step。

高温用0.2fs就好了，没见过高温还跑步长2fs的。

另外打开
&OUTER_SCF  T
      EPS_SCF    1.0000000000000001E-05
      MAX_SCF  20
   &END OUTER_SCF
另外：强烈建议看完卢天老师的帖子：CP2K中遇到SCF难收敛时的解决方法 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)，里边有对角化不容易收敛的解决方式
你的inp文件里使用对角化算法，缺了好几处关键词，就比如K点，SMEAR等，看完再说吧

作者
Author: bentty_lee 时间: 2024-7-23 12:13

吞木木发表于 2024-7-22 18:38
高温用0.2fs就好了，没见过高温还跑步长2fs的。

另外打开

感谢回复。
我仔细阅读了一下卢天老师的帖子，改了smearing，现在已经能够收敛了。
1. 关于K点，按照我以往的经验，30/L取整，大概是够的，我的体系的lattice parameter是31 19 19，足够大了，所以用gamma模式是够的，我之前也测试过KPOINTS 1 2 2，跟gamma模式差不了多少。
2. 关于时间步长，这个是参考的参考文献里面的，文献里面取的2fs，感觉原子够重的话，应该是可以的，我之前也测试过0.1 fs，也是不收敛的，所以收敛性应该更多的是跟smearing相关的，跟时间步长关系不大。
3. 关于smearing，GST体系的带隙大概在0.4 eV左右，我现在用的是300K的电子温度,可能2300K下用这个电子温度会有点小，但是我后面要做cool down到室温，所以感觉可能可以先暂时用着300K？

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