计算化学公社

标题: 计算小分子（水、6个C左右的酯类）相互作用能时计算级别选择求助 [打印本页]

作者
Author: zbhyyds 时间: 5 day ago
标题: 计算小分子（水、6个C左右的酯类）相互作用能时计算级别选择求助
本帖最后由 zbhyyds 于 2025-8-7 22:27 编辑

我在计算水二聚体的时候发现B3LYP(D3)/6-311+G(d,p)级别优化的很好，然后打算对复合物的相互作用能进行计算（有进行BSSE校正）。选择计算能量的级别时一开始选择了M062X(D3)/def2-TZVP，对比文献结果误差在8%左右。然后又重新看了sob老师的基组选择帖子，选了弱相互作用能计算的更好一点的基组jul-cc-pVTZ，误差到3%以内了。但是在论坛里搜索发现用M062X(D3)/jul-cc-pVTZ的级别不多。我猜测jul-cc-pVTZ的弥散起作用，于是尝试了ma-def2-TZVP发现结果和def2-TZVP没变化。
1.看基组选择的帖子猜测，是不是jul-cc-pVTZ和双杂化泛函一起用的多一些？
2.想问问使用M062X(D3)/jul-cc-pVTZ计算可以嘛？不行的话是不是要考虑双杂化泛函？

作者
Author: commanderwilson 时间: 5 day ago
一般不，Dunning系列基组不是设计用来给DFT算的，明显不如def2合适，至于DFT的问题，如果想要精度很高的能量可以考虑用def2-QZVP，如果需要足够多的弥散函数，实在不行用自定义基组方法用def2-TZVPPD，水分子二聚体这个体系不大，真要算相互作用能不如直接CCSD(T)

作者
Author: sobereva 时间: 4 day ago
问诸如此类问题时，你先得说明文献里用的是什么级别。如果你就想精确重复文献数据，当然应该用文献里的级别，哪怕文献里用的级别并不理想。参看下文
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678（http://bbs.keinsci.com/thread-38911-1-1.html）

用def2-TZVP算弱相互作用能时必须考虑BSSE问题（通常以counterpoise方式）。看下文的表格中不同基组算弱相互作用能误差对比
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）
ma-def2-TZVP只有s、p角动量的弥散，然而一定程度弥散的d基函数对于准确计算弱相互作用能也是不可忽视的。不明确考虑BSSE校正就用这个基组算弱相互作用能说不过去。

虽说Dunning相关一致性基组结合DFT一般来说是不理想的选择，但对于算弱相互作用能来说，如果你不特别在乎性价比、懒得自定义基组，用jul-cc-pVTZ从结果精度上来说是没有问题的。Dunning相关一致性基组绝对不要结合DFT用来算势垒、反应能那些牵扯到化学键改变的能量相关问题。

需要说明双杂化具体用什么，不能带有任何含糊性。诸如B2PLYP和revDSD-PBEP86-D3(BJ)描述弱相互作用的能力有天壤之别，前者还不如M06-2X。

M062X(D3)应写为M06-2X-D3(0)。

假设就要用B3LYP-D3(BJ)或M06-2X-D3(0)，算弱相互作用能的基组最次6-311+G(2d,p)结合counterpoise，要求更好用def2-TZVP结合counterpoise（或者def2-QZVP，但Gaussian的杂化泛函不支持利用RI加速，用这个性价比太低），参看685博文里的讨论。如果想结果更好，改用revDSD-PBEP86-D3(BJ)，若能算得动且希望尽可能精确就用CCSD(T)/jul-cc-pVTZ结合counterpoise。

作者
Author: zbhyyds 时间: 4 day ago

commanderwilson 发表于 2025-8-7 23:02
一般不，Dunning系列基组不是设计用来给DFT算的，明显不如def2合适，至于DFT的问题，如果想要精度很高的能 ...

好的，谢谢！我考虑的是想在计算量和精度之间找一个平衡。那我再看看基组方面的帖子

作者
Author: zbhyyds 时间: 4 day ago

sobereva 发表于 2025-8-8 02:47
问诸如此类问题时，你先得说明文献里用的是什么级别。如果你就想精确重复文献数据，当然应该用文献里的级别 ...

感谢Sobl老师的回复！
1.抱歉老师，我的表述不太准确。我不打算精确重复文献的数据，文献用了CCSD(T)。我是想把文献的数据作为一个标准，然后想用低级别的方法达到小的误差。
2.抱歉老师，我的表述有点含糊。我考虑BSSE问题了，在试着计算时都加了counterpoise=2
3.嗯嗯，好的。我就打算计算一下相互作用，Dunning基组之前没接触过，我再看看相关内容。
4.好的老师，对于双杂化泛函我之前也不了解，我去细看一下。
5.抱歉老师，这里我的表述不专业，看帖子说M06-2X用零阻尼矫正，我计算时使用的是em=gd3。
6.好的，谢谢老师！

作者
Author: sobereva 时间: 4 day ago

zbhyyds 发表于 2025-8-8 10:12
感谢Sobl老师的回复！
1.抱歉老师，我的表述不太准确。我不打算精确重复文献的数据，文献用了CCSD(T)。 ...

对于已知文献是CCSD(T)，你想用更低级别获得差不多的结果，对于弱相互作用计算来说，很好的做法是北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）里讲“基组”部分时介绍的这个策略：

(, 下载次数 Times of downloads: 1)

其中的小基组可以用比如aug-cc-pVDZ，大基组可以是jul-cc-pVTZ。

作者
Author: zbhyyds 时间: 3 day ago

sobereva 发表于 2025-8-8 15:16
对于已知文献是CCSD(T)，你想用更低级别获得差不多的结果，对于弱相互作用计算来说，很好的做法是北京科 ...

好的，感谢sob老师！

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