计算化学公社

标题: SAPT分解得出能量与高斯得出能量不一致问题 [打印本页]

作者
Author: 陌上花开小伟 时间: 2020-1-3 15:59
标题: SAPT分解得出能量与高斯得出能量不一致问题
老师您好想请教一个能量分解问题，我看Sob老师一篇描述氢键的文章里面SAPT分解出来的能量和高斯计算出的相互作用能有很好的一致性。但是我用SAPT分解出来的能量和高斯中计算出来的能量差距较大。我是计算了SAPT0下的6-311+g(d,p) 的能量同时高斯里面计算能量的基组用的是B3LYP/6-311+G(d,p) ，计算的结构是一个两性型的氨基酸和一个咪唑阳离子，其中在SAPT分解下的总能量为-119.61kJ/mol，高斯计算的能量为-76.77kJ/mol，整体相差42.84kJ/mol 。想问一下老师是因为我使用的基组不够精确所以才出现能量不相等，还是存在别的原因？还想请问一下这种能量不相等的情况下可否将用高斯算出的能量比较总相互作用能，而用SAPT能量比较经典，色散，诱导等能量的占比？

作者
Author: sobereva 时间: 2020-1-3 17:41
计算级别完全不一样，而且B3LYP算弱相互作用完全是垃圾，结果当然大相径庭。
仔细看下文了解弱相互作用计算常识，不要做无意义的计算
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413
DFT-D色散校正的使用
http://sobereva.com/210
乱谈DFT-D
http://sobereva.com/83
大体系弱相互作用计算的解决之道
http://sobereva.com/214

我下面那篇文章里用的SAPT2+(3)dMP2和CCSD(T)都是高精度级别，因此结果必然能很好吻合
透彻认识氢键本质、简单可靠地估计氢键强度：一篇2019年JCC上的重要研究文章介绍
http://sobereva.com/513（http://bbs.keinsci.com/thread-14600-1-1.html）

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