计算化学公社

标题: 求助用高斯计算反应的单点能两种方法得出结果为什么相差太大 [打印本页]

作者
Author: lhao9 时间: 2023-5-12 18:25
标题: 求助用高斯计算反应的单点能两种方法得出结果为什么相差太大
请问老师，用gaussian软件计算反应（反应引入了电场）过程中TS的零点能和单点能，零点能选用b3lyp/6-31g*计算，单点能分别选用了b3lyp/6-311++g** empiricaldispersion=gd3和CCSD/def2tzvp方法计算，将零点能和单点能相加后，两种方法计算出来得到的反应能量最大相差了16 kcal/mol，这是为什么呢？反应体系的原子主要有N、O、Si、C、H。

作者
Author: sobereva 时间: 2023-5-12 18:52
基组相差太悬殊，理论方法也相差太悬殊，根本没可比性。

如果体系是中性的，直接用CCSD(T)/cc-pVTZ，如果有余力就上到CCSD(T)/cc-pVQZ，算出来是什么就是什么（保证波函数稳定的前提下）。CCSD能算得动就毫无理由用B3LYP，如今用B3LYP算能量也极难发得出去文章，看
坚持使用B3LYP算有机体系能量的人的下场（http://bbs.keinsci.com/thread-12773-1-1.html）

作者
Author: lhao9 时间: 2023-5-13 15:11

sobereva 发表于 2023-5-12 18:52
基组相差太悬殊，理论方法也相差太悬殊，根本没可比性。

如果体系是中性的，直接用CCSD(T)/cc-pVTZ，如 ...

好的谢谢老师

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