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

白也

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


计算流程（使用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)

2022-10-7 10:05 上传 Uploaded

点击下载
Click to download

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

2022-10-7 10:05 上传 Uploaded

点击下载
Click to download

chenjhust

wzkchem5 发表于 2024-3-7 18:08
对于碘等较重主族元素的赝势计算，ORCA 4.2.1的初猜有bug，参见https://mattermodeling.stackexchange.co ...

好的老师，我在学习一下

wzkchem5

chenjhust 发表于 2024-3-7 10:46
好的，谢谢

对于碘等较重主族元素的赝势计算，ORCA 4.2.1的初猜有bug，参见https://mattermodeling.stackexch ... e-failing/6312#6312
解决方法是改用最新版本的ORCA，如因为客观原因无法改用最新版，可以尝试用PAtom初猜

chenjhust

1. 
   

复制代码

wzkchem5

chenjhust 发表于 2024-3-7 10:19
老师好，我也是一个含碘的化合物，请问在高斯09用b3lyp/6-311g(d,p)优化好的初始结构，经过multiwfn转化 ...

上传输出文件。如果涉及保密问题，至少上传输出文件的最后大约一百行

chenjhust

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

白也

sobereva 发表于 2022-10-7 10:34
建议先尝试ORCA 5.0.3

刚用自己的电脑使用orca5.0.3跑了一遍，还真是orca版本的问题，现在正常了，谢谢社长。

sobereva

建议先尝试ORCA 5.0.3