不同基组优化下水分子氢键形成自由能变化数据差异大

dididi

背景：本小白最近在计算有机反应过程中发现不同基组校正过后代算得到的各结构之间的自由能差值发生较大变化，但别人说应该是我数据算错了，基组的影响应该很小，为了简化问题，我设置了一个溶剂为甲苯，两个水分子形成氢键的体系来观察数据。拿不同基组进行了结构优化，并收集热力学数据。文件格式：
%chk=*******.chk
#p  B3LYP/gen   opt  optcyc=30  freq  nosymm  em=GD3BJ  scrf=(smd,solvent=toluene)  

********

0    1
8       -1.300742517      1.243072524     -3.268483408
1       -1.300742517      2.002409524     -2.672440408
1       -1.300742517      0.483735524     -2.672440408

H -O
6-31g**(6-311g**/6-311+g**/6-311++g(2df,2p)/def2tzvp)
****

数据如下（带Δ的数据列单位为kcal/mol）

由图可知，由不同基组计算得到的自由能差异很大，极差将近7 kcal/mol（最后一列），但加上零点校正之后的电子能之差却稳定在-3.几上下（倒数第三列），未加零点校正的电子能之差极差更是达到了10 kcal/mol左右（左数第四列）。这一巨大差异和我算自己的体系时观察到的结果一致，因为我会拿小基组做优化，大基组做单点后拿大基组的E和小基组的G校正量做加和，不同基组下E差值巨大就会对我G的数据影响很大。
是我其中用到的一些基组过于粗糙吗？是乱加弥散的结果吗？或许是我数据获取的地方、计算的方法有问题？
还请各位不吝赐教~

dididi

使用软件为guassian09

sobereva

dididi 发表于 2022-7-3 14:45
使用软件为guassian09

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

以后不要再拼错哪怕一次Gaussian

比较那么一大堆级别根本没意义。诸如def2-TZVP根本毫无测试的必要，有常识的人都知道3-zeta不带弥散函数而且不考虑BSSE问题算的弱相互作用结合能根本没法接受。6-31G**、6-311G**也都是瞎测。把下文仔细看了
谈谈量子化学研究中什么样的benchmark才有意义
http://sobereva.com/554（http://bbs.keinsci.com/thread-17826-1-1.html）
乱谈DFT-D
http://sobereva.com/83
谈谈弥散函数和“月份”基组
http://sobereva.com/119
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

通常算氢键相互作用能至少6-311+G(d,p)起步

dididi

抱歉沦为反例了，我一定吸取教训好好看第5条发帖规范，辛苦老师、谢谢老师了。我误以为氢键与弥散关系小加不加弥散无关系，希望看到帖子的人不要学我一样盲目猜测。