计算化学公社

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

[ORCA] 用ORCA计算CCSD(T)/CC-pvqz出错

[复制链接 Copy URL]

4

帖子

0

威望

64

eV
积分
68

Level 2 能力者

求助贴,大家好!上次参加sob老量子化学培训班,sob老师推荐用orca计算能量比较快,所有我回来后准备采用orca算能量选取的 CCSD(T)/CC-PVQZ和CCSD(T)/CC-PVTZ方法.对于CCSD(T)/CC-PVTZ, 我所有的计算都能正常结束不报错. 但是用 CCSD(T)/CC-PVQZ方法计算的时候会出现问题,具体如下:

这是我的输入文件
! cc-pvqz CCSD(T) TightSCF
%pal nprocs 16 end
%maxcore 2500

* xyz 0 2
H -2.31345100 0.44361700 0.92904500
H -2.31351500 0.44358900 -0.92893700
C -0.69392900 -0.76268300 0.00001700
H -0.59099000 -1.38949300 -0.88579600
H -0.59094600 -1.38949400 0.88582400
N 0.60459400 0.06523500 -0.00001400
O 1.63100600 -0.57160700 -0.00001600
O 0.51664800 1.26879000 -0.00003500
C -1.90681800 0.07229500 0.00004600
*
用这样的输入文件,采用CCSD(T)/CC-PVQZ计算小分子的时候也( 例如C2H4,C2H3, NO2)也都能正常结束,但是计算大分子(比如 C2H4NO2, C2H4ONO2, CH2OCH2NO 和  NO2 + C2H4反应的过渡态)就报错. 报错信息如下:
Number of amplitudes to be optimized ... 43430812
Iter E(tot) E(Corr) Delta-E Residual Time <S|S>**1/2
0 --------------------------------------------------------------------------
mpirun noticed that process rank 4 with PID 0 on node localhost exited on signal 9 (Killed).
--------------------------------------------------------------------------
ORCA finished by error termination in MDCI
Calling Command: mpirun -np 16 --allow-run-as-root /home/orca/orca_4_0_0_linux_x86-64/orca_mdci_mpi C2H4NO2-TS-CH3CHNO2.mdciinp.tmp
[file orca_tools/qcmsg.cpp, line 432]:
.... aborting the run



请问要怎么解决这个问题???


谢谢

1万

帖子

0

威望

9851

eV
积分
22085

Level 6 (一方通行)

27#
发表于 Post on 2021-6-17 15:11:25 | 只看该作者 Only view this author
XZDZ 发表于 2021-6-17 06:22
老师,请问ORCA用CCSD(T)-F12时,cc-pVTZ基组不能加弥散吗?不能用aug-cc-pVTZ吗

CCSD(T)-F12计算应该用专门的F12基组,比如cc-pVTZ-F12。cc-pVTZ-F12和aug-cc-pVTZ差不多大,相当于已经自带弥散了(虽然名字里没有aug),所以不用额外加弥散
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
本团队长期招收研究生,有意者可私信联系

159

帖子

0

威望

1231

eV
积分
1390

Level 4 (黑子)

26#
发表于 Post on 2021-6-17 13:22:28 | 只看该作者 Only view this author
sobereva 发表于 2017-10-23 15:25
maxcore设小点再试,或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T ...

老师,请问ORCA用CCSD(T)-F12时,cc-pVTZ基组不能加弥散吗?不能用aug-cc-pVTZ吗

38

帖子

0

威望

591

eV
积分
629

Level 4 (黑子)

25#
发表于 Post on 2021-1-5 20:50:11 | 只看该作者 Only view this author
sobereva 发表于 2021-1-5 19:21
取自EXPLICITLY CORRELATED COUPLED-CLUSTER THEORY(from Recent Progress in Coupled Cluster Methods)
...

好的,这么说来确实是有依据的,感谢sob老师的指点

38

帖子

0

威望

591

eV
积分
629

Level 4 (黑子)

24#
发表于 Post on 2021-1-5 20:49:40 | 只看该作者 Only view this author
biogon 发表于 2021-1-5 19:19
那个benchmark我觉得不靠谱,之前看过的几个benchmark里面ccpvdzf12的精度确实是和ccpvqz相仿的

好的,谢谢您的回复

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

23#
发表于 Post on 2021-1-5 19:21:09 | 只看该作者 Only view this author
机智猴 发表于 2021-1-5 19:05
感谢您的回复。
我刚刚看了一下,可能是我回复错了帖子了还是怎么回事,往上翻没翻到sob老师的那段话。
...

取自EXPLICITLY CORRELATED COUPLED-CLUSTER THEORY(from Recent Progress in Coupled Cluster Methods)



北京科音自然科学研究中心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

1245

帖子

1

威望

3860

eV
积分
5125

Level 6 (一方通行)

22#
发表于 Post on 2021-1-5 19:19:44 | 只看该作者 Only view this author
机智猴 发表于 2021-1-5 19:05
感谢您的回复。
我刚刚看了一下,可能是我回复错了帖子了还是怎么回事,往上翻没翻到sob老师的那段话。
...

那个benchmark我觉得不靠谱,之前看过的几个benchmark里面ccpvdzf12的精度确实是和ccpvqz相仿的

38

帖子

0

威望

591

eV
积分
629

Level 4 (黑子)

21#
发表于 Post on 2021-1-5 19:05:31 | 只看该作者 Only view this author
wzkchem5 发表于 2021-1-5 18:21
这个没问题啊,如果都是CCSD(T)-F12计算,cc-pVDZ-F12的耗时肯定比aug-cc-pVTZ少,误差也大。cc-pVDZ-F12 ...

感谢您的回复。
我刚刚看了一下,可能是我回复错了帖子了还是怎么回事,往上翻没翻到sob老师的那段话。
原贴:http://bbs.keinsci.com/thread-7062-1-1.html。在这个帖子里,sob老师提到“用到QZ级别时建议用F12。写CCSD(T)-F12 cc-pvdz-F12 cc-pvdz-F12-CABS,精度和CCSD(T)/cc-pVQZ基本一样,耗时能省很多。”
因为在平时计算中我认为ccsd(t)结合aug-cc-pvtz或者结合cc-pvqz的耗时是差不多的,相当于一个级别的计算精度。然后看到sob老师这里讲到CCSD(T)-F12 cc-pvdz-F12 cc-pvdz-F12-CABS和CCSD(T)/cc-pvqz精度差不多,结合上述那篇文献里结论“cc-pvdz-f12比aug-cc-pvtz的精度差不少”,才有这个疑惑的。所以我认为,同样的CCSD(T)-F12级别计算,cc-pvdz-F12 cc-pvdz-F12-CABS的精度是比cc-pvqz低不少的。
感谢您的回复和指导。

1万

帖子

0

威望

9851

eV
积分
22085

Level 6 (一方通行)

20#
发表于 Post on 2021-1-5 18:21:57 | 只看该作者 Only view this author
机智猴 发表于 2021-1-5 17:22
老师好,我看文献里(标题:Assessment of basis sets for F12 explicitly correlated molecular electro ...

这个没问题啊,如果都是CCSD(T)-F12计算,cc-pVDZ-F12的耗时肯定比aug-cc-pVTZ少,误差也大。cc-pVDZ-F12基组的大小类似于aug-cc-pVDZ。
但是如果你比CCSD(T)-F12/cc-pVDZ-F12和不用F12的CCSD(T)/aug-cc-pVTZ,那么大概率是前者更好
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
本团队长期招收研究生,有意者可私信联系

38

帖子

0

威望

591

eV
积分
629

Level 4 (黑子)

19#
发表于 Post on 2021-1-5 17:22:06 | 只看该作者 Only view this author
sobereva 发表于 2017-10-23 15:25
maxcore设小点再试,或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T ...

老师好,我看文献里(标题:Assessment of basis sets for F12 explicitly correlated molecular electronic-structure methods)评测的cc-pvdz-f12基组的误差明显大于aug-cc-pvtz(我认为aug-cc-pvtz和cc-pvqz的精度差不多)。所以我认为应该是cc-pvtz-f12和cc-pvqz差不多。这里您是笔误了还是说我自己哪部分搞错了?

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

18#
发表于 Post on 2018-4-30 18:16:31 | 只看该作者 Only view this author
长颈鹿先森 发表于 2018-4-30 14:52
sob老师,我使用orca  ! CCSD(T) cc-pVQZ tightSCF pal4 计算一个大分子单点能出现错误,错误结果如图, ...

尝试linux版
北京科音自然科学研究中心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

103

帖子

0

威望

780

eV
积分
883

Level 4 (黑子)

17#
发表于 Post on 2018-4-30 14:52:15 | 只看该作者 Only view this author
sobereva 发表于 2017-10-23 15:25
maxcore设小点再试,或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T ...

sob老师,我使用orca  ! CCSD(T) cc-pVQZ tightSCF pal4 计算一个大分子单点能出现错误,错误结果如图,后来就试了试小分子CF4,结果还是图片所示;我调小了设置的内存和并行核数还是不行,结果还是那样。输入文件和输出文件放在这里,您帮我看看,我不知道怎么回事。我又试了一个原子F的,结果也是这样。

1525070780(1).png (11.08 KB, 下载次数 Times of downloads: 75)

1525070780(1).png

SPE input.txt

353 Bytes, 阅读权限: 30, 下载次数 Times of downloads: 4

SPE output.txt

47.77 KB, 阅读权限: 30, 下载次数 Times of downloads: 2

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

16#
发表于 Post on 2018-1-9 14:21:10 | 只看该作者 Only view this author
尚艳磊 发表于 2018-1-9 13:55
好的,老师!
还有个问题,Auxiliary basis set, /C 这个关键字 是什么意思呢?
比如 cc-pVDZ/C 和 cc- ...

cc-pVDZ/C和cc-pVTZ/C分别是适合cc-pVDZ和cc-pVTZ用的电子相关计算时的辅助基组。对辅助基组不了解的话看此文
大体系弱相互作用计算的解决之道
http://sobereva.com/214
北京科音自然科学研究中心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

47

帖子

0

威望

478

eV
积分
525

Level 4 (黑子)

15#
发表于 Post on 2018-1-9 13:55:05 | 只看该作者 Only view this author
本帖最后由 尚艳磊 于 2018-1-9 14:03 编辑
sobereva 发表于 2018-1-7 18:54
CCSD(T)别写RI-,速度不会有提升
只有耦合簇结合F12的时候才有利用RI的意义

好的,老师!
还有个问题,Auxiliary basis set, <basis set>/C 这个关键字 是什么意思呢?
比如 cc-pVDZ/C 和 cc-pVTZ/C 有什么区别呢?我看了manual 也没有解释(也可能我看了没找到。。),我也试了一下,算出来的能量差别也不是很大。Auxiliart basis sets for correlation calculations

6万

帖子

99

威望

6万

eV
积分
125122

管理员

公社社长

14#
发表于 Post on 2018-1-7 18:54:25 | 只看该作者 Only view this author
尚艳磊 发表于 2018-1-7 18:40
sob老师,RI加速只支持RHF体系吗?我选择了RHF,但还是出现了问题。我的输入文件如下:
|  1> ! RHF RI- ...


CCSD(T)别写RI-,速度不会有提升
只有耦合簇结合F12的时候才有利用RI的意义
北京科音自然科学研究中心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

47

帖子

0

威望

478

eV
积分
525

Level 4 (黑子)

13#
发表于 Post on 2018-1-7 18:40:04 | 只看该作者 Only view this author
sobereva 发表于 2018-1-6 07:10
44核才64GB内存太小了
1.4GB*44=61.6GB,本身ORCA实际对内存使用量往往会高于maxcore一些,再加上系统和 ...

sob老师,RI加速只支持RHF体系吗?我选择了RHF,但还是出现了问题。我的输入文件如下:
|  1> ! RHF RI-CCSD(T) cc-pVTZ cc-pVTZ/C
|  2> # Test for NO2
|  3>
|  4> *xyz 0 2
|  5>  N                  0.00000000    0.00000000    0.31500700
|  6>  O                  0.00000000    1.09057300   -0.13781600
|  7>  O                  0.00000000   -1.09057300   -0.13781600
|  8> **
而且它提示的也是这样:WARNING: your system is open-shell and RHF/RKS was chosen
  ===> : WILL SWITCH to UHF/UKS

但是结果却说:[file orca_mdci/mdci_uhf_sigma.cpp, line 3388]: Error (ORCA_MDCI): K(C)-RI not implemented in UHF branch
请问老师,这该咋么解决啊?谢谢啦

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

GMT+8, 2026-2-17 18:17 , Processed in 0.201961 second(s), 30 queries , Gzip On.

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