计算化学公社

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

[ORCA] 关于ORCA4.2.1中使用PWPB95计算含碘化合物MP2能量错误的问题

[复制链接 Copy URL]

44

帖子

0

威望

589

eV
积分
633

Level 4 (黑子)

最近算了一个含碘元素的化合物,其结构如下:


计算流程(使用orca4.2.1版本,输入文件使用Multiwfn生成):
(1)使用B3LYP-D3优化结构并计算频率检验虚频,关键词如下
!B3LYP D3 def2-TZVP(-f) def2/J RIJCOSX opt freq tightSCF noautostart miniprint nopop
(2)上一步完成,没有虚频,基于最优结构采用双杂化泛函PWPB95计算电子能量
! PWPB95 D3 def2-QZVPP def2/J def2-QZVPP/C RIJCOSX tightSCF grid4 gridx4 noautostart miniprint nopop

现在的问题是,在第二步中,计算的MP2能量非常大,明显错误,简单列一下两部分能量和总能量(摘自输出文件)
DFT部分    Total Energy       :        -1270.13969511 Eh
MP2部分   MP2 TOTAL ENERGY:     -1162676939355587.500000000 Eh
色散校正部分   Dispersion correction           -0.023309235
总能量     FINAL SINGLE POINT ENERGY   -1162676939355587.500000000000

能否帮忙看一下是什么原因造成的?是因为含碘采用了赝势导致还是说其他原因?输入和输出文件已经上传。

7m_all_sp.inp (1.81 KB, 下载次数 Times of downloads: 4)

7m_all_sp.out (116.17 KB, 下载次数 Times of downloads: 4)



我见青山多妩媚,料青山见我应如是。

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

2#
发表于 Post on 2022-10-7 10:34:47 | 只看该作者 Only view this author
建议先尝试ORCA 5.0.3
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

44

帖子

0

威望

589

eV
积分
633

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-10-7 11:17:39 | 只看该作者 Only view this author
sobereva 发表于 2022-10-7 10:34
建议先尝试ORCA 5.0.3

刚用自己的电脑使用orca5.0.3跑了一遍,还真是orca版本的问题,现在正常了,谢谢社长。
我见青山多妩媚,料青山见我应如是。

10

帖子

0

威望

59

eV
积分
69

Level 2 能力者

4#
发表于 Post on 2024-3-7 17:19:58 | 只看该作者 Only view this author
sobereva 发表于 2022-10-7 10:34
建议先尝试ORCA 5.0.3

老师好,我也是一个含碘的化合物,请问在高斯09用b3lyp/6-311g(d,p)优化好的初始结构,经过multiwfn转化成ORCA(ORCA4.2.1)的输入文件计算更高精度的单点能时老是出现SCF不收敛的问题试过kdiis,调整积分精度都无法解决,请问有什么好的建议吗?
! PWPB95 D3 def2-TZVPP def2/J def2-TZVPP/C RIJCOSX grid5 gridx5 tightSCF noautostart kdiis miniprint nopop
%maxcore  3000
%pal nprocs   48 end
* xyz   0   3

1万

帖子

0

威望

7396

eV
积分
18151

Level 6 (一方通行)

5#
发表于 Post on 2024-3-7 17:42:30 | 只看该作者 Only view this author
chenjhust 发表于 2024-3-7 10:19
老师好,我也是一个含碘的化合物,请问在高斯09用b3lyp/6-311g(d,p)优化好的初始结构,经过multiwfn转化 ...

上传输出文件。如果涉及保密问题,至少上传输出文件的最后大约一百行
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

59

eV
积分
69

Level 2 能力者

6#
发表于 Post on 2024-3-7 17:46:59 | 只看该作者 Only view this author
本帖最后由 chenjhust 于 2024-3-7 18:17 编辑

复制代码

1万

帖子

0

威望

7396

eV
积分
18151

Level 6 (一方通行)

7#
发表于 Post on 2024-3-7 18:08:31 | 只看该作者 Only view this author

对于碘等较重主族元素的赝势计算,ORCA 4.2.1的初猜有bug,参见https://mattermodeling.stackexch ... e-failing/6312#6312
解决方法是改用最新版本的ORCA,如因为客观原因无法改用最新版,可以尝试用PAtom初猜
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

10

帖子

0

威望

59

eV
积分
69

Level 2 能力者

8#
发表于 Post on 2024-3-7 18:12:00 | 只看该作者 Only view this author
wzkchem5 发表于 2024-3-7 18:08
对于碘等较重主族元素的赝势计算,ORCA 4.2.1的初猜有bug,参见https://mattermodeling.stackexchange.co ...

好的老师,我在学习一下

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 11:06 , Processed in 0.170595 second(s), 23 queries , Gzip On.

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