计算化学公社

标题: 如何处理阴离子添加弥散函数后导致计算量增加的问题 [打印本页]

作者
Author: bianwenbo 时间: 2016-6-6 10:01
标题: 如何处理阴离子添加弥散函数后导致计算量增加的问题
Sob老师，您在《大体系弱相互作用计算的解决之道》中说过“对于分子间相互作用的计算，弥散函数起到的作用和CP实际上是相仿佛的”。我的研究对象是水溶剂中的阴离子团簇，在进行量化计算时必须添加弥散函数，这就导致在ORCA中使用AotuAux辅助基组从而计算量大增。因此为了降低计算量，我想要向您请教三个问题：

是否可以使用文中的DFT-D3 + gCP的方法替代原来必须添加的弥散函数？
如果不可以，那么使用最外层基函数是半弥散的def2-QZVPP(-g,-f)基组代替常规的ma-def-XXXX系列基组用于阴离子是否可行？
如果必须添加弥散函数，能否使用诸如jul-cc-pVTZ或may-cc-pVQZ的“月份”基组配合PWPB95-D3双杂化泛函处理阴离子体系？

谢谢您指导。

作者
Author: sobereva 时间: 2016-6-6 10:20
1 不可以。gCP只是解决缺乏弥散函数时计算分子间相互作用时的BSSE问题，而弥散函数的用处不仅限于此，诸如计算阴离子、里德堡态，电子分布都很弥散，没有弥散函数来描述肯定不行。
2 对于阴离子、里德堡态，QZ级别的半弥散和全弥散也依然有差距。
3 可以。阴离子体系有sp弥散，至多再增加个d弥散就足够了。

作者
Author: bianwenbo 时间: 2016-6-6 12:51
本帖最后由 bianwenbo 于 2016-6-6 13:48 编辑

sobereva 发表于 2016-6-6 10:20
1 不可以。gCP只是解决缺乏弥散函数时计算分子间相互作用时的BSSE问题，而弥散函数的用处不仅限于此，诸如 ...

感谢您的回答，我想再深入探讨一下。

这个我明白了，谢谢。
关于全弥散和弥散，您在《谈谈弥散函数和“月份”基组》中提到“另外，提高zeta数对结果的改进比起进一步增加弥散函数的最高角动量明显得多，在弥散函数上花费过多不值得。然而，需要弥散函数的任务中，给cc-pVnZ哪怕仅加入maug-级别的弥散对结果的改进都要比提高zeta数要显着得多”，我想这针对Dunning系列基组的推论应该也可以用在Ahlrichs系列基组上，即针对阴离子而言，全弥散ma-def2-SVP是否比半弥散def2-QZVPP(-g,-f)更准确（或更合适）？

3.1 在ORCA手册3.2.4 Basis sets的部分指出，软件支持的含有弥散函数的辅助基组仅有/C后缀的aug-cc-pVDZ/C和aug-cc-pVTZ/C。因此能否说使用 `RI-PWPB95 D3 aug-cc-pVTZ aug-cc-pVTZ/C RIJCOSX ...` 比 `RI-PWPB95 D3 ma-def2-QZVPP(-g,-f) AutoAux RIJCOSX ...` 来计算阴离子团簇是更好的选择？（因为后者需要AutoAux辅助基组）
3.2 在ORCA手册3.2.4，支持的基组中没有aug系列“月份”基组，所以我们想要在ORCA中使用的话就需要自行修改基组文件么？

谢谢。

作者
Author: sobereva 时间: 2016-6-6 13:06

bianwenbo 发表于 2016-6-6 12:51
感谢您的回答，我想再深入探讨一下。

2 是
3.1 是。
3.2 是。密度拟合基组你还用aug的，普通基组把高角动量弥散砍掉

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