计算化学公社

标题: PSI4进行SAPT能量分解得到的总能量与高斯计算得到的相互作用能不一致问题 [打印本页]

作者
Author: 爱喝茶 时间: 2025-5-28 15:54
标题: PSI4进行SAPT能量分解得到的总能量与高斯计算得到的相互作用能不一致问题
本帖最后由爱喝茶于 2025-5-28 16:07 编辑

研究一个弱相互作用体系，在高斯中使用revDSD-PBEP86-D3(BJ)/aug-cc-pVTZ计算了弱相互作用能，使用counterpoise=2考虑了BSSE，最后得到：
complexation energy =    -6.23 kcal/mole (raw)
complexation energy =    -5.59 kcal/mole (corrected)

使用psi4的SAPT2+(3)δMP2进行能量分解，已知SAPT本身没有BSSE问题，相互作用能结果为：
Total SAPT2+(3)dMP2             -9.78728672 [mEh]    -6.14161514 [kcal/mol]    -25.69651775 [kJ/mol]

可以看到，SAPT的结果与高斯未进行BSSE校正前的结果相近，那么哪一个才是真正的相互作用能？如果校正后的是准确的，那SAPT分解的各能量项结果还能用吗？
（看到sob老师http://sobereva.com/526中例1计算得到的相互作用能相当精确）

作者
Author: sobereva 时间: 2025-5-29 08:38
revDSD-PBEP86-D3(BJ)/aug-cc-pVTZ+CP的精度比不上SAPT2+(3)δMP2（假定SAPT2+(3)δMP2用的基组合适）。要质疑SAPT2+(3)δMP2的精度，得用CCSD(T)比

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