计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8154|回复 Reply: 27
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] ORCA计算结合能,SCF总是不收敛怎么办?

[复制链接 Copy URL]

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

游客,本帖隐藏的内容需要积分高于 150 才可浏览,您当前积分为 0

上传的文件.zip

4.24 MB, 下载次数 Times of downloads: 14

opt和e文件

4290

帖子

4

威望

9558

eV
积分
13928

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2021-11-30 21:57:55 | 只看该作者 Only view this author
你需要ORCA 5.0.1,SCF大有改善。
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

3#
发表于 Post on 2021-11-30 22:34:34 | 只看该作者 Only view this author
这种超大幅度的SCF能量振荡,一般是基组线性相关问题。在%scf里把sthresh设大一些,比如设到1e-7或1e-6,再试
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

155

帖子

0

威望

3433

eV
积分
3588

Level 5 (御坂)

4#
发表于 Post on 2021-12-1 01:51:11 | 只看该作者 Only view this author
orca scf不收敛可以先试一下 slowconv 或者 veryslowconv

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

5#
发表于 Post on 2021-12-1 04:15:22 | 只看该作者 Only view this author
indec 发表于 2021-11-30 18:51
orca scf不收敛可以先试一下 slowconv 或者 veryslowconv

这两个关键词对于解决基组线性相关问题没用
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2021-12-1 11:08:04 | 只看该作者 Only view this author
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


结果还是报错(如下图),感觉是不是我的输入文件有问题?

qiuzhu.png (150.7 KB, 下载次数 Times of downloads: 36)

结果报错提示

结果报错提示

444

帖子

0

威望

2580

eV
积分
3024

Level 5 (御坂)

娃娃儿鱼

7#
发表于 Post on 2021-12-1 12:45:07 | 只看该作者 Only view this author
本帖最后由 hdhxx123 于 2021-12-1 12:47 编辑

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

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-12-1 14:36:21 | 只看该作者 Only view this author
hdhxx123 发表于 2021-12-1 12:45
上ORCA5.0.1,直接tightscf就行。先从def2TZVP收敛好了再整madef2TZVP来读;另外,我记得sob说了MV不好收敛 ...

谢谢您,我试一试

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

9#
发表于 Post on 2021-12-1 16:21:13 | 只看该作者 Only view this author
CayJ 发表于 2021-12-1 04:08
老师,您好!按照您的提示,我修改了输入文件:

%pal nprocs   4 end

那就是autoaux的问题,sthresh只能解决轨道基组的线性相关问题,不能解决辅助基组的线性相关问题。可以换ORCA 5.0.1试试,我记得autoaux有改善
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

227

帖子

0

威望

2800

eV
积分
3027

Level 5 (御坂)

10#
发表于 Post on 2021-12-1 16:31:23 | 只看该作者 Only view this author
我不知道怎么解决,但能说这是正常情况。课题组合成的多个分子,原子数>70,就没有一个能用弥散基组成功收敛的,所以现在都用def2-tzvpp

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2021-12-1 17:03:24 | 只看该作者 Only view this author
shalene 发表于 2021-12-1 16:31
我不知道怎么解决,但能说这是正常情况。课题组合成的多个分子,原子数>70,就没有一个能用弥散基组成功收 ...

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

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

12#
发表于 Post on 2021-12-1 17:38:29 | 只看该作者 Only view this author
shalene 发表于 2021-12-1 09:31
我不知道怎么解决,但能说这是正常情况。课题组合成的多个分子,原子数>70,就没有一个能用弥散基组成功收 ...

主要是autoaux的问题,弥散函数本身不是问题,如果不用RIJCOSX(因此没有辅助基组),弥散函数收敛概率会高很多,尽管可能仍然比没有弥散函数要低。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

227

帖子

0

威望

2800

eV
积分
3027

Level 5 (御坂)

13#
发表于 Post on 2021-12-1 18:02:41 | 只看该作者 Only view this author
wzkchem5 发表于 2021-12-1 17:38
主要是autoaux的问题,弥散函数本身不是问题,如果不用RIJCOSX(因此没有辅助基组),弥散函数收敛概率会 ...

不用rijcosx,太慢了

227

帖子

0

威望

2800

eV
积分
3027

Level 5 (御坂)

14#
发表于 Post on 2021-12-1 18:13:46 | 只看该作者 Only view this author
CayJ 发表于 2021-12-1 17:03
不用弥散,能量计算的结果不能用吧。

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

444

帖子

0

威望

2580

eV
积分
3024

Level 5 (御坂)

娃娃儿鱼

15#
发表于 Post on 2021-12-1 18:37:37 | 只看该作者 Only view this author
shalene 发表于 2021-12-1 18:13
以前社长说过def2基组对弥散要求不高,加了色散矫正后算弱相互作用不用加弥散。你用def2-qzvpp,保证能收 ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-27 13:06 , Processed in 0.170555 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list