计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8155|回复 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文件

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

28#
发表于 Post on 2021-12-4 23:24:05 | 只看该作者 Only view this author
CayJ 发表于 2021-12-4 20:00
sob老师,ORCA4.2.0也法直接用def2-QZVPP(-g)
请问这个自定义该怎么写输入文件呢?请老师稍微指点一下。

def2-QZVPP能算得动就直接用,算不动就根据手册说明自定义def2-QZVPP(-g)
去掉g极化能节约一些时间,对于普通泛函计算g极化不会有什么影响
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

27#
 楼主 Author| 发表于 Post on 2021-12-4 20:21:41 | 只看该作者 Only view this author
wzkchem5 发表于 2021-12-4 20:08
就用def2-QZVP或def2-QZVPP也行,差别不大

好的,谢谢老师

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

26#
发表于 Post on 2021-12-4 20:08:51 | 只看该作者 Only view this author
CayJ 发表于 2021-12-4 13:00
sob老师,ORCA4.2.0也法直接用def2-QZVPP(-g)
请问这个自定义该怎么写输入文件呢?请老师稍微指点一下。

就用def2-QZVP或def2-QZVPP也行,差别不大
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 (黑子)

25#
 楼主 Author| 发表于 Post on 2021-12-4 20:00:46 | 只看该作者 Only view this author
sobereva 发表于 2021-12-3 08:39
自定义
5.0开始已经没法直接用了

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

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

24#
发表于 Post on 2021-12-3 15:36:55 | 只看该作者 Only view this author
sobereva 发表于 2021-12-2 22:20
我一直觉得应该专门搞个适合ma-的辅助基组,既划算又没有太严重的线性相关问题的。这对于算弱相互作用必 ...

同意,其实可能在def2/J基础上加几个弥散函数就行了
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

23#
发表于 Post on 2021-12-3 08:39:25 | 只看该作者 Only view this author
CayJ 发表于 2021-12-3 08:09
谢谢老师的解答,但是我用如下的输入文件计算结合能,ORCA5.0.1报错说“def2-QZVPP(-g)”有问题

%pal  ...

自定义
5.0开始已经没法直接用了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

22#
 楼主 Author| 发表于 Post on 2021-12-3 08:09:41 | 只看该作者 Only view this author
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里没有这样的表达?

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

21#
发表于 Post on 2021-12-3 05:20:37 | 只看该作者 Only view this author
wzkchem5 发表于 2021-12-3 04:56
话说不知道有没有人做过benchmark看看def2/J能不能做ma-类的基组的辅助基组。。。虽然没有benchmark的情 ...

我一直觉得应该专门搞个适合ma-的辅助基组,既划算又没有太严重的线性相关问题的。这对于算弱相互作用必要性极高
直接用def2/J我觉得悬...有个benchmark也好
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

20#
发表于 Post on 2021-12-3 04:56:03 | 只看该作者 Only view this author
sobereva 发表于 2021-12-2 20:53
直接用def2-QZVPP(-g)就完了,别跟弥散函数+autoaux较劲,来回折腾白浪费时间

话说不知道有没有人做过benchmark看看def2/J能不能做ma-类的基组的辅助基组。。。虽然没有benchmark的情况下大家都不敢这么用,但是搞不好其实是可以这么用的也不好说
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
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125208

管理员

公社社长

19#
发表于 Post on 2021-12-3 03:53:58 | 只看该作者 Only view this author
直接用def2-QZVPP(-g)就完了,别跟弥散函数+autoaux较劲,来回折腾白浪费时间

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
StormSpirts + 5 哈哈,原来如此!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

444

帖子

0

威望

2580

eV
积分
3024

Level 5 (御坂)

娃娃儿鱼

18#
发表于 Post on 2021-12-2 09:36:42 | 只看该作者 Only view this author
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不用担心坏了得重新搞,比较推荐
真·探

75

帖子

0

威望

632

eV
积分
707

Level 4 (黑子)

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

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

444

帖子

0

威望

2580

eV
积分
3024

Level 5 (御坂)

娃娃儿鱼

16#
发表于 Post on 2021-12-1 18:38:47 | 只看该作者 Only view this author
CayJ 发表于 2021-12-1 14:36
谢谢您,我试一试

你这个如果不优化的话tightscf就够了,如果优化记得用verytightscf
真·探

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还是要大得多的
真·探

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

GMT+8, 2026-2-27 14:33 , Processed in 0.218870 second(s), 25 queries , Gzip On.

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