求助葫芦脲与客体分子吸附问题的泛函基组推荐

DFT.. · 发表于 Post on 2025-7-8 22:27:51

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

Stardust0831 · 发表于 Post on 2025-7-8 22:46:07

本帖最后由 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结合量子化学程序方便地计算分子的各种热力学数据》

sobereva · 发表于 Post on 2025-7-9 05:02:10

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

DFT.. · 发表于 Post on 2025-7-9 14:31:38

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

感谢老师回复！

DFT.. · 发表于 Post on 2025-7-9 14:32:00

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

感谢老师回复！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 求助葫芦脲与客体分子吸附问题的泛函基组推荐

评分 Rate

浏览过的版块