高斯得到的相互作用能与ADF的差别很大

wangweihua

老师，我在高斯程序中用矫正后的总体系的能量分别减去两个片段的能量的到AB的相互作用能为31KJ/MoL,后来用ADF程序做能量分解时得到的bonding energy为46.6KJ/mol,已经考虑变形能了，但是差距还很大，请问怎么解决？

sobereva

Gaussian和ADF用的基函数明显不一样，前者用的gauss基函数，后者用的slater基函数，光是这点就让两个程序结果很难对比。而且用的数值上的设定还不一样（DFT积分格点精度、收敛限等），而且你还得确保两个程序用的泛函相同（但是ADF用杂化泛函的时候还有特殊性）。还有可能用的结构也不一样。
只要你在Gaussian里用了足够好的泛函，足够大的基组，甚至有余力情况下用counterpoise把BSSE也考虑了，Gaussian下的结果肯定是没错的。
算相互作用能的时候通常没必要考虑变形能。