计算化学公社

标题: sobEDA计算时使用自定义基组时识别出错 [打印本页]

作者
Author: 西柚 时间: 2025-10-2 13:53
标题: sobEDA计算时使用自定义基组时识别出错
本帖最后由西柚于 2025-10-2 16:18 编辑

应审稿人要求按照一篇文章的方法补充EDA的分析，体系是MOF材料+Xe原子。我咨询了作者使用的基组，他回应说对于Xe最低要使用aug-def2-tzvpp （也叫def2-tzvppd），由于算力有限，打算使用def2-tzvpd代替，但该基组好像没有在Gaussian16中内置。我根据sob老师的文章（http://sobereva.com/60），自己学着写template.gjf文件，尝试了文中的两种方法（一种是采用混合基组，关键字为genecp，分别对CHON使用def2TZVP基组，对Xe使用def2-tzvpd，参数是在https://www.basissetexchange.org/上下载的，并在后面也写了对应的赝势的内容；另一种是下载了def2-tzvpd.gbs基组文件，通过@.. .gbs的方式让所有原子都使用该基组），都无法正常进行。希望各位老师指点一下（我不是专门做计算的，因反稿时间比较紧，一些基础薄弱的地方现在来不急慢慢补了）。用到的文件和结果上传在附件中。

作者
Author: sobereva 时间: 2025-10-2 18:10
aug-def2-tzvpp和def2-tzvppd完全是两码事，仔细看下文
给Ahlrichs的def2系列基组加弥散的方法
http://sobereva.com/340（http://bbs.keinsci.com/thread-3487-1-1.html）

根本没必要加弥散，def2-TZVP结合counterpoise对于计算Xe和MOF的相互作用能足矣，下文列出的sobEDA原文里的基组横向测试结果明确体现了
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

sobEDA教程里明确给了带自定义基组的例子，诸如(CO)5CrCH2

作者
Author: 西柚 时间: 2025-10-3 11:29
非常感谢老师！！

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