计算化学公社

标题: 求助：UiO片段不收敛，已经给Zr簇补氢中和电荷 [打印本页]

作者
Author: lemon_electron 时间: 2025-2-21 19:35
标题: 求助：UiO片段不收敛，已经给Zr簇补氢中和电荷
各位老师好，最近在做一个研究，想看在Zr6O8团簇存在下，对一系列联苯衍生物各种亲电亲核指数的影响，以这种方式模拟UiO片段，在最初始的模型上就碰到了问题；

按照如图所示的结构进行计算，B3LYP理论级别下，给Zr用了def2-SV(P)基组，其它部分用了6-31G*，gjf和最终的报错文件已上传，希望各位老师指正

（ps* 关于Zr簇电荷的问题，对于kappa3-O已经补了一半成-OH，以中和Zr团簇的电荷）

作者
Author: sobereva 时间: 2025-2-22 09:01
甭写完全多余的noraman，仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

“给Zr用了def2-SV(P)基组，其它部分用了6-31G*”很莫名其妙，没有理由不直接都用def2-SV(P)

本来就SCF不收敛，还不用#P监控收敛趋势

其它方面下文该说的都说了
解决SCF不收敛问题的方法
http://sobereva.com/61

令Gaussian 16中SCF未收敛到默认收敛限也能继续做后续计算的方法
http://sobereva.com/625（http://bbs.keinsci.com/thread-26427-1-1.html）

作者
Author: lemon_electron 时间: 2025-2-22 13:55

sobereva 发表于 2025-2-22 09:01
甭写完全多余的noraman，仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http:/ ...

收到，谢谢sob老师，立马改掉坏习惯

因为之前给纯分子用了6-31G*，想着加上金属之后进行比较，所以还是保持了非金属元素的基组；还是直接都拿def2-SV(P)都算一遍吧

作者
Author: zjxitcc 时间: 2025-2-22 15:54
提醒两点基组细节，若您之前就注意到了则没关系；若没注意到，以后应避免踩坑：
（1）def2系列基组对Zr是有定义赝势的，若Gaussian输入文件使用自定基组/混合基组格式，应当写为
#p B3LYP/genecp ...
[省略]
Zr 0
def2SV
****
[空行]
Zr 0
def2
[空行]
[空行]

关键词、赝势信息要齐全，否则容易变成完全错误的计算。

（2）在计算A...B复合物时，往往用户也会计算孤立的分子A、孤立的分子B，这样才可以计算结合能。假设基组和赝势使用情况如下
i. 单独计算A，#p这一行里写def2SVP，则Gaussian默认5D 7F
ii. 单独计算B，#p这一行里写6-31G(d)，则Gaussian默认6D 7F
iii. 计算A...B复合物，#p这一行里写gen或genecp，坐标底下使用def2SVP和6-31G(d)混合基组/赝势，则Gaussian默认5D 7F

问题会出在ii上，需要用户在#p这一行加上关键词5D 7F（这两个关键词就是这么短），保证与i和iii计算级别相同。

作者
Author: lemon_electron 时间: 2025-2-22 17:18

zjxitcc 发表于 2025-2-22 15:54
提醒两点基组细节，若您之前就注意到了则没关系；若没注意到，以后应避免踩坑：
（1）def2系列基组对Zr是 ...

确实之前没注意到，太感谢了，邹老师，谢谢您的提醒！

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