计算化学公社

标题: 关于MOF吸附乙烷乙烯的结合能疑问 [打印本页]

作者
Author: 红茶泡枸杞 时间: 2024-9-18 12:51
标题: 关于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。金属为铝，是否需要单独设置方法？

作者
Author: sobereva 时间: 2024-9-18 16:48
基组太小，用counterpoise后依然不理想。不要不分场合地乱用DZVP-MOLOPT-SR-GTH

北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里我专门强调过选用的事
(, 下载次数 Times of downloads: 1)

当前根本没必要开smearing

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

作者
Author: 红茶泡枸杞 时间: 2024-9-18 20:59
好的谢谢社长答复，我用TZVP-MOLOPT-GTH试一下。毕竟和文献相差了一个数量级，结果大概率是错的。

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