透彻认识氢键本质、简单可靠地估计氢键强度：一篇2019年JCC上的重要研究文章介绍

wzkchem5

牧生 发表于 2022-4-10 08:37
重新进行了计算，使用orca得到输入文件。为了表明都是同一个级别，且为了照顾阴离子体系B，所以，都选用 ...

虽然这个跟你观察到的IGMH图的问题估计关系不大，但是如果用双杂化泛函，建议考虑MP2项对密度的贡献（%mp2 density relaxed end；如果算不动，把relaxed改成unrelaxed，但精度会差一点。作图时注意不要用SCF密度，要用MP2密度），否则算出来的密度的精度可能还不如普通杂化泛函，因为漏掉了一些correlation的贡献

sobereva

牧生 发表于 2022-4-10 15:37
重新进行了计算，使用orca得到输入文件。为了表明都是同一个级别，且为了照顾阴离子体系B，所以，都选用 ...

因为等值面不是无限薄，而函数也不是处处相同，因此两个方向看会有所不同。如果把等值面数值设大一点，等值面会有所收缩，两个方向看的差异会更小。反正怎么看中间都是蓝色，就没比要刻意区分了。

wB97X-2是双杂化泛函，默认情况下产生的波函数是它的杂化泛函部分的，然而这还不如用B3LYP等普通泛函算得好。要么产生双杂化泛函的波函数（需要额外时间），要么用B3LYP的波函数，建议后者。因为我的JCC 2019年的文章里的参数是在B3LYP/ma-TZVP下拟合的，套用预测公式时波函数也应当在那个级别下计算。

牧生

sobereva 发表于 2022-4-11 00:25
因为等值面不是无限薄，而函数也不是处处相同，因此两个方向看会有所不同。如果把等值面数值设大一点，等 ...

OK，问题解决，感谢解答
进行了B3LYP/ma-TZVP和wB97X-2 D3 ma-def2-QZVPP得到的结果对比，带电体系，得到的数值几乎完全精确相同，差异不超过0.1%；对于不带电的体系，大约有1~2%的差异。