计算化学公社

标题: 用ORCA做AIMD时SCF不收敛问题 [打印本页]

作者
Author: deltazheng 时间: 2021-12-11 20:56
标题: 用ORCA做AIMD时SCF不收敛问题
本帖最后由 deltazheng 于 2021-12-12 17:19 编辑

Sob老师，我仿照博文“使用ORCA做从头算动力学(AIMD)的简单例子（http://sobereva.com/576）”尝试做有机体系的AIMD计算。我的体系有145个原子（CHONP，无过渡金属），仿照例子用Multiwfn生成inp文件，MD信息如下：

%md
timestep 1.0_fs
initvel 298.15_K no_overwrite
thermostat berendsen 298.15_K timecon 30.0_fs
dump position stride 1 format xyz filename "pos.xyz"
constraint add center 0..144
run 2000
end

计算级别为B97-3c。经过多次尝试，在模拟进行到70-80步时都出现SCF不收敛，然后报错退出。out文件末尾信息如下：

      | Time |  t_SCF | t_Grad |    Temp |    E_Kin |       E_Pot |       E_Tot |    Cons.Qty | Drift
   Step | [fs] | | |    [K] | [Hartree] |    [Hartree] |    [Hartree] |    [Hartree] |    [K]

……
   70    70.0 114.7    21.2    286.24 0.195796 -4143.965805 -4143.770009 -4143.899443    1.32
   71    71.0 581.0    20.9    282.12 0.192981 -4144.936028 -4144.743047 -4144.872891 -705.33
Warning: The SCF did not converge. Trying to continue...
   72    72.0 590.4    21.1    278.89 0.190767 -4145.909189 -4145.718422 -4145.848753  -1413.73
Error: The SCF did again not converge. Aborting.

把模拟步长减小为0.7 fs，或者使用tightscf均无法改变这一现象。该用GFN2-xTB方法可以正常完成2000步模拟，且结果符合预期。请问B97-3c级别下问题可能的原因是什么，应该如何调整。谢谢~

作者
Author: sobereva 时间: 2021-12-12 05:05
为什么加了删除线？
说明ORCA版本

作者
Author: wzkchem5 时间: 2021-12-12 15:33
检查最后一步的几何结构是否合理

作者
Author: deltazheng 时间: 2021-12-12 17:20

sobereva 发表于 2021-12-12 05:05
为什么加了删除线？
说明ORCA版本

ORCA 4.2.1
删除线是误操作，已经去掉~

作者
Author: deltazheng 时间: 2021-12-12 17:24

wzkchem5 发表于 2021-12-12 15:33
检查最后一步的几何结构是否合理

最后一步结构合理。AIMD轨迹间附件文件。

作者
Author: wzkchem5 时间: 2021-12-12 17:51

deltazheng 发表于 2021-12-12 10:24
最后一步结构合理。AIMD轨迹间附件文件。

嗯，确实合理。我又看了一下你的输出文件，最后几步的能量是偏低而不是偏高，说明可能遇到了基组线性相关问题，可以考虑在%scf里增大sthresh重跑

作者
Author: deltazheng 时间: 2021-12-12 21:54

wzkchem5 发表于 2021-12-12 17:51
嗯，确实合理。我又看了一下你的输出文件，最后几步的能量是偏低而不是偏高，说明可能遇到了基组线性相关 ...

我把sthresh改成1e-7，重跑还是报错，有一个RETURN CODE : 65280，不知道是什么意思

Warning: The SCF did not converge. Trying to continue...
   72    72.0 406.7    20.9 2610.18 1.785446 -4143.410528 -4141.625082 -4141.681478 1611.38
   73    73.0 610.0    21.5 2490.31 1.703451 -4144.407348 -4142.703897 -4142.708618 865.76
Warning: The SCF did not converge. Trying to continue...
   74    74.0 593.5    21.2 4522.94 3.093826 -4145.375120 -4142.281294 -4142.186396 1244.85
Error: The SCF did again not converge. Aborting.

Errors occurred in the MD loop.
CSANExecutor::Execute(): Errors occurred during execution.
ERROR    : MD Program returns an error
            cannot continue the md simulation
COMMAND    : /public1/home/sc80404/package/ORCA421/orca_4_2_1_linux_x86-64_shared_openmpi314/orca_md INTER-0-major  orca
RETURN CODE : 65280
[file orca_main/mainmd.cpp, line 177]:

作者
Author: wzkchem5 时间: 2021-12-13 00:21

deltazheng 发表于 2021-12-12 14:54
我把sthresh改成1e-7，重跑还是报错，有一个RETURN CODE : 65280，不知道是什么意思

Warning: The SCF ...

这个具体要看scf部分的输出，不能看主输出文件

作者
Author: deltazheng 时间: 2021-12-14 07:52

wzkchem5 发表于 2021-12-13 00:21
这个具体要看scf部分的输出，不能看主输出文件

换用ORCA 5.0.2 就不再出现这个问题了~

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