计算化学公社

标题: 求助葫芦脲与客体分子吸附问题的泛函基组推荐 [打印本页]

作者
Author: DFT.. 时间: 2025-7-8 22:27
标题: 求助葫芦脲与客体分子吸附问题的泛函基组推荐
各位老师好，我最近在做葫芦脲与客体分子吸附问题，体系总共有160个分子左右，主要计算主客体间的弱相互作用，主客体分子仅涉及C、N、O、H、S这些元素，通过查阅文献打算采用wB97XD/6-311G**做结构优化及电子能量（为后续获得热力学数据）的计算，各位老师觉得这个泛函基组的组合够用吗？或者有其他更好的推荐吗？

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

作者
Author: Stardust0831 时间: 2025-7-8 22:46
本帖最后由 Stardust0831 于 2025-7-8 22:50 编辑

结构优化3zeta有点太大了，用wB97XD/6-31G**通常够用。《浅谈为什么优化和振动分析不需要用大基组》
如果目的是计算弱相互作用，可以考虑在wB97XD/6-31G**优化得到的基础上，用M06-2X/6-311G**算单点，记得考虑BSSE：《谈谈BSSE校正与Gaussian对它的处理》

我个人更推荐用ORCA的gCP来处理BSSE问题，gCP有解析梯度和hessian，基本上不会增加很多耗时。
如果换ORCA，我推荐RI-wB97M-V/def2-TZVP//RI-B3LYP-D3(BJ)/def2-SVP ,并且记得带上gCP

以ORCA 6.1.0为例，关键词可以是：
预优化：

!XTB2 opt

复制代码

结构优化：

! B3LYP D3 def2-SVP gCP(DFT/SVP) def2/J RIJCOSX tightopt freq verytightSCF noautostart miniprint defgrid3

复制代码

单点：

! wB97M-V def2-TZVP gCP(DFT/TZ) def2/J RIJCOSX verytightSCF noautostart miniprint defgrid3

复制代码

吉布斯自由能之类的热力学数据推荐用shermo:《使用Shermo结合量子化学程序方便地计算分子的各种热力学数据》

作者
Author: sobereva 时间: 2025-7-9 05:02
如果要算弱相互作用能，6-311G**结合counterpoise也不够用，应当用6-311+G(2d,p)结合counterpoise。可参考下文表格里不同基组计算弱相互作用能误差的对比
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

作者
Author: DFT.. 时间: 2025-7-9 14:31

sobereva 发表于 2025-7-9 05:02
如果要算弱相互作用能，6-311G**结合counterpoise也不够用，应当用6-311+G(2d,p)结合counterpoise。可参考 ...

感谢老师回复！

作者
Author: DFT.. 时间: 2025-7-9 14:32

Stardust0831 发表于 2025-7-8 22:46
结构优化3zeta有点太大了，用wB97XD/6-31G**通常够用。《浅谈为什么优化和振动分析不需要用大基组》
如果 ...

感谢老师回复！

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