计算化学公社

标题: ORCA计算结合能，SCF总是不收敛怎么办？ [打印本页]

作者
Author: CayJ 时间: 2021-11-30 21:55
标题: ORCA计算结合能，SCF总是不收敛怎么办？

作者
Author: zjxitcc 时间: 2021-11-30 21:57
你需要ORCA 5.0.1，SCF大有改善。

作者
Author: wzkchem5 时间: 2021-11-30 22:34
这种超大幅度的SCF能量振荡，一般是基组线性相关问题。在%scf里把sthresh设大一些，比如设到1e-7或1e-6，再试

作者
Author: indec 时间: 2021-12-1 01:51
orca scf不收敛可以先试一下 slowconv 或者 veryslowconv

作者
Author: wzkchem5 时间: 2021-12-1 04:15

indec 发表于 2021-11-30 18:51
orca scf不收敛可以先试一下 slowconv 或者 veryslowconv

这两个关键词对于解决基组线性相关问题没用

作者
Author: CayJ 时间: 2021-12-1 11:08

wzkchem5 发表于 2021-11-30 22:34
这种超大幅度的SCF能量振荡，一般是基组线性相关问题。在%scf里把sthresh设大一些，比如设到1e-7或1e-6，再 ...

老师，您好！按照您的提示，我修改了输入文件：

%pal nprocs 4 end

! wB97M-V ma-def2-TZVP autoaux RIJCOSX strongSCF grid5 gridx5 noautostart miniprint nopop
%maxcore 2000

%scf
sthresh 1e-7
end

* xyz 0 1

结果还是报错（如下图），感觉是不是我的输入文件有问题？

作者
Author: hdhxx123 时间: 2021-12-1 12:45
本帖最后由 hdhxx123 于 2021-12-1 12:47 编辑

上ORCA5.0.1，直接tightscf就行。先从def2TZVP收敛好了再整madef2TZVP来读；另外，我记得sob说了MV不好收敛，用wB97XV也可以的。

作者
Author: CayJ 时间: 2021-12-1 14:36

hdhxx123 发表于 2021-12-1 12:45
上ORCA5.0.1，直接tightscf就行。先从def2TZVP收敛好了再整madef2TZVP来读；另外，我记得sob说了MV不好收敛 ...

谢谢您，我试一试

作者
Author: wzkchem5 时间: 2021-12-1 16:21

CayJ 发表于 2021-12-1 04:08
老师，您好！按照您的提示，我修改了输入文件：

%pal nprocs 4 end

那就是autoaux的问题，sthresh只能解决轨道基组的线性相关问题，不能解决辅助基组的线性相关问题。可以换ORCA 5.0.1试试，我记得autoaux有改善

作者
Author: shalene 时间: 2021-12-1 16:31
我不知道怎么解决，但能说这是正常情况。课题组合成的多个分子，原子数>70，就没有一个能用弥散基组成功收敛的，所以现在都用def2-tzvpp

作者
Author: CayJ 时间: 2021-12-1 17:03

shalene 发表于 2021-12-1 16:31
我不知道怎么解决，但能说这是正常情况。课题组合成的多个分子，原子数>70，就没有一个能用弥散基组成功收 ...

不用弥散，能量计算的结果不能用吧。

作者
Author: wzkchem5 时间: 2021-12-1 17:38

shalene 发表于 2021-12-1 09:31
我不知道怎么解决，但能说这是正常情况。课题组合成的多个分子，原子数>70，就没有一个能用弥散基组成功收 ...

主要是autoaux的问题，弥散函数本身不是问题，如果不用RIJCOSX（因此没有辅助基组），弥散函数收敛概率会高很多，尽管可能仍然比没有弥散函数要低。

作者
Author: shalene 时间: 2021-12-1 18:02

wzkchem5 发表于 2021-12-1 17:38
主要是autoaux的问题，弥散函数本身不是问题，如果不用RIJCOSX（因此没有辅助基组），弥散函数收敛概率会 ...

不用rijcosx，太慢了

作者
Author: shalene 时间: 2021-12-1 18:13

CayJ 发表于 2021-12-1 17:03
不用弥散，能量计算的结果不能用吧。

以前社长说过def2基组对弥散要求不高，加了色散矫正后算弱相互作用不用加弥散。你用def2-qzvpp，保证能收敛，我也不信有审稿人说你基组不行

作者
Author: hdhxx123 时间: 2021-12-1 18:37

shalene 发表于 2021-12-1 18:13
以前社长说过def2基组对弥散要求不高，加了色散矫正后算弱相互作用不用加弥散。你用def2-qzvpp，保证能收 ...

不过其实我觉得def2qzvp比madef2tzvp还是要大得多的

作者
Author: hdhxx123 时间: 2021-12-1 18:38

CayJ 发表于 2021-12-1 14:36
谢谢您，我试一试

你这个如果不优化的话tightscf就够了，如果优化记得用verytightscf

作者
Author: CayJ 时间: 2021-12-2 08:28

hdhxx123 发表于 2021-12-1 12:45
上ORCA5.0.1，直接tightscf就行。先从def2TZVP收敛好了再整madef2TZVP来读；另外，我记得sob说了MV不好收敛 ...

老师，您好！您说的“先从def2TZVP收敛好了再整madef2TZVP来读”
按照您的方法，我用ORCA4.2.0计算结合能，用的是def2TZVP，很快就收敛了。
然后是在收敛的gbw文件上，建立计算结合能的输入文件，加上弥散，继续计算结合能吗？

作者
Author: hdhxx123 时间: 2021-12-2 09:36

CayJ 发表于 2021-12-2 08:28
老师，您好！您说的“先从def2TZVP收敛好了再整madef2TZVP来读”
按照您的方法，我用ORCA4.2.0计算结合 ...

是的，如果你没法换版本的话，那就4.2.0下读gbw换基组为madef2tzvp就行了，两种方法：1.inp文件和gbw文件不换名字，加个autostart关键词就完了；2.gbw换个名字，然后在inp文件里写moread和%moinp=XXX.gbw（可能不是这样写的，好久没写忘了，你看下手册）这样读，这个方法gbw不用担心坏了得重新搞，比较推荐

作者
Author: sobereva 时间: 2021-12-3 03:53
直接用def2-QZVPP(-g)就完了，别跟弥散函数+autoaux较劲，来回折腾白浪费时间

作者
Author: wzkchem5 时间: 2021-12-3 04:56

sobereva 发表于 2021-12-2 20:53
直接用def2-QZVPP(-g)就完了，别跟弥散函数+autoaux较劲，来回折腾白浪费时间

话说不知道有没有人做过benchmark看看def2/J能不能做ma-类的基组的辅助基组。。。虽然没有benchmark的情况下大家都不敢这么用，但是搞不好其实是可以这么用的也不好说

作者
Author: sobereva 时间: 2021-12-3 05:20

wzkchem5 发表于 2021-12-3 04:56
话说不知道有没有人做过benchmark看看def2/J能不能做ma-类的基组的辅助基组。。。虽然没有benchmark的情 ...

我一直觉得应该专门搞个适合ma-的辅助基组，既划算又没有太严重的线性相关问题的。这对于算弱相互作用必要性极高
直接用def2/J我觉得悬...有个benchmark也好

作者
Author: CayJ 时间: 2021-12-3 08:09

sobereva 发表于 2021-12-3 03:53
直接用def2-QZVPP(-g)就完了，别跟弥散函数+autoaux较劲，来回折腾白浪费时间

谢谢老师的解答，但是我用如下的输入文件计算结合能，ORCA5.0.1报错说“def2-QZVPP(-g)”有问题

%pal nprocs 4 end

! wB97M-V def2-QZVPP(-g) RIJCOSX tightSCF noautostart miniprint nopop
%maxcore 2000
* xyz 0 1

我查了手册，没有找到(-g)的半弥散的表达，是不是在ORCA5.0.1里没有这样的表达？

作者
Author: sobereva 时间: 2021-12-3 08:39

CayJ 发表于 2021-12-3 08:09
谢谢老师的解答，但是我用如下的输入文件计算结合能，ORCA5.0.1报错说“def2-QZVPP(-g)”有问题

%pal ...

自定义
5.0开始已经没法直接用了

作者
Author: wzkchem5 时间: 2021-12-3 15:36

sobereva 发表于 2021-12-2 22:20
我一直觉得应该专门搞个适合ma-的辅助基组，既划算又没有太严重的线性相关问题的。这对于算弱相互作用必 ...

同意，其实可能在def2/J基础上加几个弥散函数就行了

作者
Author: CayJ 时间: 2021-12-4 20:00

sobereva 发表于 2021-12-3 08:39
自定义
5.0开始已经没法直接用了

sob老师，ORCA4.2.0也法直接用def2-QZVPP(-g)
请问这个自定义该怎么写输入文件呢？请老师稍微指点一下。

作者
Author: wzkchem5 时间: 2021-12-4 20:08

CayJ 发表于 2021-12-4 13:00
sob老师，ORCA4.2.0也法直接用def2-QZVPP(-g)
请问这个自定义该怎么写输入文件呢？请老师稍微指点一下。

就用def2-QZVP或def2-QZVPP也行，差别不大

作者
Author: CayJ 时间: 2021-12-4 20:21

wzkchem5 发表于 2021-12-4 20:08
就用def2-QZVP或def2-QZVPP也行，差别不大

好的，谢谢老师

作者
Author: sobereva 时间: 2021-12-4 23:24

CayJ 发表于 2021-12-4 20:00
sob老师，ORCA4.2.0也法直接用def2-QZVPP(-g)
请问这个自定义该怎么写输入文件呢？请老师稍微指点一下。

def2-QZVPP能算得动就直接用，算不动就根据手册说明自定义def2-QZVPP(-g)
去掉g极化能节约一些时间，对于普通泛函计算g极化不会有什么影响

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