关于MOF吸附乙烷乙烯的结合能疑问

红茶泡枸杞

各位老师好，我用CP2K对MOF吸附乙烷乙烯分子进行结合能的计算（代算），计算方法如下，PBE，DZVP-MOLOPT-SR-GTH，DFT-D3(BJ)，对角化算法，开启smearing，CUTOFF  500，REL_CUTOFF  60，ALPHA 0.2，NBROYDEN 12，用BSSE计算结合能，得到的结果为-6.32、-7.11kcal/mol，我看文献上的在-20几（文献大多采用Gaussian计算），请问我现在计算的方法是否适合当前体系？如果还可以的话，这个结合能和其他文献相差有点大，可以用来发文章吗？
我尝试过用Gaussian截取配体加金属来计算，b3lyp/6-311+g(d) em=gd3bj，counterpoise=2，结果是-6.10、-6.04kcal/mol。金属为铝，是否需要单独设置方法？

红茶泡枸杞

好的谢谢社长答复，我用TZVP-MOLOPT-GTH试一下。毕竟和文献相差了一个数量级，结果大概率是错的。

sobereva

基组太小，用counterpoise后依然不理想。不要不分场合地乱用DZVP-MOLOPT-SR-GTH

北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里我专门强调过选用的事



当前根本没必要开smearing



你这个体系的相互作用能不太可能达到-20 kJ/mol那么大，文献里的情况和你的估计没有可比性。