计算化学公社

标题: 关于def2-TZVPPD基组搭配不同的辅助基组的简单测评 [打印本页]

作者
Author: 冰释之川 时间: 2022-7-28 08:38
标题: 关于def2-TZVPPD基组搭配不同的辅助基组的简单测评
本帖最后由冰释之川于 2022-7-28 13:09 编辑

据说，在ORCA里def2-TZVPPD没有配套的辅助基组用于RI-J加速，本文分别测试借用def2/J def2-TZVPPD/C辅助基组和autoaux自动拟合辅助基组两种基组方案精算单点能（分子结构保持一致），从而考察某个化学反应的自由能垒。

单点能方案1：PWPB95 D3 def2-TZVPPD def2/J def2-TZVPPD/C RIJCOSX tightSCF noautostart miniprint nopop SCFConvForced NOTRAH
单点能方案2：PWPB95 D3 def2-TZVPPD autoaux RIJCOSX tightSCF noautostart miniprint nopop SCFConvForced NOTRAH

考察的反应为：
(1)Alder-ene反应(中性体系)：
(, 下载次数 Times of downloads: 37)
(2)硫酯水解反应(阴离子体系)：
(, 下载次数 Times of downloads: 32)

计算结果：
(1)Alder-ene反应(中性体系)：
(, 下载次数 Times of downloads: 30)
(, 下载次数 Times of downloads: 43)
(2)硫酯水解反应(阴离子体系)：
(, 下载次数 Times of downloads: 35)

结论：
无论是中性还是阴离子体系，用def2-TZVPPD基组搭配def2/J def2-TZVPPD/C辅助基组与搭配autoaux自动拟合基组给出的能垒差异<0.03kcal/mol，基本上可以忽略。
所以，def2-TZVPPD基组搭配def2/J def2-TZVPPD/C辅助基组也还算可行？

友情提示：
用autoaux自动拟合基组方案似乎要比def2/J def2-TZVPPD/C辅助基组更加吃内存(所以附件里采用前者计算的时候，我只能降低核数，提升单核内存了，囧RZ。。)

单点能输出文件（给出一对ORCA输出文件）：
(, 下载次数 Times of downloads: 24)
(, 下载次数 Times of downloads: 9)

作者
Author: mfdsrax2 时间: 2022-7-28 09:16
加速效果有多少？

作者
Author: RAL 时间: 2022-7-28 10:02
本帖最后由 RAL 于 2022-7-28 20:10 编辑

也可能是弥散函数本身对这个体系不太重要？
补充：已看到阴离子体系的测试结果，辛苦楼主！

作者
Author: ghifi37 时间: 2022-7-28 10:37
没弄明白，基组是用def2-TZVPD还是def2-TZVPPD？

作者
Author: 冰释之川 时间: 2022-7-28 10:48

ghifi37 发表于 2022-7-28 10:37
没弄明白，基组是用def2-TZVPD还是def2-TZVPPD？

哈哈，笔误笔误，基组是def2-TZVPPD

作者
Author: ionexchangeC 时间: 2022-7-28 14:52
PSI4对于def2-nZVPD系列，默认用的辅助基组就是def2/J

作者
Author: GoldenBaby 时间: 2023-2-18 18:27
老师，我在ORCA 5.0.3的手册里看到了这种组合
! EA-EOM-DLPNO-CCSD NORMALPNO ma-def2-SVP RIJCOSX aug-cc-pVDZ/C def2/J
请问对于普通的DFT计算也可以这么组合吗？

作者
Author: wzkchem5 时间: 2023-2-18 18:34

GoldenBaby 发表于 2023-2-18 11:27
老师，我在ORCA 5.0.3的手册里看到了这种组合
! EA-EOM-DLPNO-CCSD NORMALPNO ma-def2-SVP RIJCOSX aug-cc ...

如果能找到文献证明这么组合对于EA-EOM-DLPNO-CCSD来说确实是合理的（光是在手册上看到，还不算特别有力的证据），那么可以用来证明对于DFT也可以这么组合

作者
Author: GoldenBaby 时间: 2023-2-19 10:46

wzkchem5 发表于 2023-2-18 18:34
如果能找到文献证明这么组合对于EA-EOM-DLPNO-CCSD来说确实是合理的（光是在手册上看到，还不算特别有力 ...

好的，感谢回复

欢迎光临计算化学公社 (http://bbs.keinsci.com/)