计算化学公社

标题: 请问用高斯计算，描述氢键比较好的基组有哪些 [打印本页]

作者
Author: Yoghurt 时间: 2019-3-30 16:18
标题: 请问用高斯计算，描述氢键比较好的基组有哪些
如题，现用高斯计算基态和激发态下的有机小分子构型优化和红外谱、势能面，其中分子内含有氢键作用，请问目前比较公认的方法和基组有哪些？描述氢键作用比较好的基组推荐一下吧，谢谢

作者
Author: alonewolfyang 时间: 2019-3-30 19:27
http://sobereva.com/336

作者
Author: sobereva 时间: 2019-3-31 04:38
含氢键的体系opt freq目的用6-311G**或def-TZVP就足够了

至于“势能面”指代不明。如果你想说计算单点，看实际情况而定，仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

其它相关信息看
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387（http://bbs.keinsci.com/thread-6600-1-1.html）
乱谈DFT-D
http://sobereva.com/83

作者
Author: wuzhiyi 时间: 2019-4-2 00:55
本帖最后由 wuzhiyi 于 2019-4-2 04:39 编辑

sobereva老师好,我想对一个有一定共价性的O-H...O的氢键做单点结合能计算. 系统有90个原子左右,系统中只有碳氧氮氢. 但体系中有酪氨酸和組氨酸的pi-pi stacking.
系统已经用D3LYP-D3(BJ)/def2-TZVP优化完毕因为要算键级所以没有加弥散原计划是用MP2/jul-cc-pVQZ加counterpoise算结合能
但因为集群设置的问题,只有50TB的scratch space可以用.jul-cc-pVQZ好像略微有一点大.现在想保持MP2往下减基组.
我想问一下是往def2-TZVP上加弥散函数比较明智还
还是把jul-cc-pVQZ保持弥散函数砍到3zeta
还是保持4zeta砍月份比较明智?谢谢

作者
Author: sobereva 时间: 2019-4-2 08:44

wuzhiyi 发表于 2019-4-2 00:55
sobereva老师好,我想对一个有一定共价性的O-H...O的氢键做单点结合能计算. 系统有90个原子左右,系统中只有 ...

别考虑MP2
像这样较大体系，MP2或双杂化用Gaussian算极其困难，如果要结果很准确，强烈建议用ORCA做PWPB95-D3(BJ)计算，基组jul-cc-pVTZ，一定要利用RI节约耗时。

如果不用双杂化，用Gaussian靠M06-2X-D3/ma-TZVP + 0.5*CP也可以得到足够合理的结果，内行审稿人绝对不会挑刺。

作者
Author: wuzhiyi 时间: 2019-4-2 16:33
本帖最后由 wuzhiyi 于 2019-4-2 16:39 编辑

sobereva 发表于 2019-4-2 08:44
别考虑MP2
像这样较大体系，MP2或双杂化用Gaussian算极其困难，如果要结果很准确，强烈建议用ORCA做PWPB ...

谢谢sobereva老师，想问一下是用ma-TZVP而不是更新的ma-def2-TZVP吗？还有ma是may的意思吗？
我还是新手不太明白0.5×CP是什么意思，CP是指counterpoise么？0.5×CP是指0.5×BSSE吗？谢谢

作者
Author: sobereva 时间: 2019-4-2 17:40

wuzhiyi 发表于 2019-4-2 16:33
谢谢sobereva老师，想问一下是用ma-TZVP而不是更新的ma-def2-TZVP吗？还有ma是may的意思吗？
我还是新手 ...

仔细看
给ahlrichs的def2系列基组加弥散的方法
http://sobereva.com/340（http://bbs.keinsci.com/thread-3487-1-1.html）

0.5*CP就是用一半的counterpoise校正能，这里说了
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

作者
Author: wuzhiyi 时间: 2019-4-2 18:02

sobereva 发表于 2019-4-2 17:40
仔细看
给ahlrichs的def2系列基组加弥散的方法
http://sobereva.com/340（http://bbs.keinsci.com/thre ...

谢谢sob老师的答复我明白了

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