求助量子化学与第一性原理结果的比较

SHANY

老师您好，

我的催化体系中既有均相的部分也有在材料表面反应的贡献。

但是我听说由于量化程序与第一性原理程序的泛函与基组的不同，其计算结果不能相互比较。

请问如果用cp2k使用相同泛函，基组与k点密度得出的结果，是不是就可以互相比较了？

谢谢

SHANY

多谢老师回答

sobereva

你说的不确切
首先，即便是两个不同的量子化学程序，都用相同的理论方法和基组，结果也未必能直接比较，详见
谈谈不同量子化学程序计算结果的差异问题
http://sobereva.com/573（http://bbs.keinsci.com/thread-20122-1-1.html）
其次，量子化学程序和第一性原理程序的结果未必不能比较。例如CP2K可以用GAPW形式做全电子DFT计算，在计算分子或团簇时，只要用的泛函和基组和量子化学程序相同，并且把两个程序的积分格点精度设得尽可能好、CP2K里真空区设得足够大、平面波截断能足够高、SCF收敛阈值足够严，等等，注意各种影响可比性的细节设置，结果的差异完全在可接受范围内。

大多数量子化学程序只支持孤立体系的计算，所以既然你提到“k点密度”这事，那必然是算的周期性体系，光凭这点通常就注定不可能比较。不过，Gaussian等少量量子化学程序倒也支持周期性计算并且可以考虑k点采样，若CP2K和Gaussian做周期性计算时用严格相同的计算级别，k点定义相同，并且注意前面提到的影响可比性的因素，结果还是基本能比的。注意由于Gaussian的周期性计算功能算三维周期性体系相当慢，CP2K常研究的体系规模用Gaussian一般根本算不动，所以实际中还是难以或无法对比。

所以简而言之，你的问题答案是“既可以比又不可以比”。