计算化学公社

标题: 计算XeF6负离子结果莫名报错的问题 [打印本页]

作者
Author:
楚明城    时间: 7 hour ago
标题: 计算XeF6负离子结果莫名报错的问题
本帖最后由 楚明城 于 2025-8-12 12:52 编辑

各位老师好,我在计算XeF6负离子结构时Gaussian莫名报错,也没有看出错误原因在哪,请问一下这是什么情况出问题了呢
作者
Author:
tonganlhy    时间: 7 hour ago
输入文件也要贴上来。不过一眼就看出来自旋多重度错了
作者
Author:
楚明城    时间: 7 hour ago
tonganlhy 发表于 2025-8-12 12:41
输入文件也要贴上来。不过一眼就看出来自旋多重度错了

您好,自旋多重度确实输错了,但修改之后还是报错了
作者
Author:
tonganlhy    时间: 7 hour ago
所以说输入文件也要发上来。以后问报错输出输入文件都要放。
我估计你是赝势基组输错了
作者
Author:
楚明城    时间: 7 hour ago
tonganlhy 发表于 2025-8-12 12:55
所以说输入文件也要发上来。以后问报错输出输入文件都要放。
我估计你是赝势基组输错了

好的,已经放上来了,赝势是从BES上拷贝的,应该没啥问题?
作者
Author:
tonganlhy    时间: 7 hour ago
楚明城 发表于 2025-8-12 12:58
好的,已经放上来了,赝势是从BES上拷贝的,应该没啥问题?

小跑了一下。可以跑起来。不过第一轮SCF不收敛
参考此文解决
http://bbs.keinsci.com/thread-4829-1-1.html
作者
Author:
Uus/pMeC6H4-/キ    时间: 7 hour ago
这个赝势基组/赝势看起来是def2-SVP的数据?如果确实是,def2系列除镧系锕系外元素的基组在Gaussian已经内置而并不需要去BSE下载,可以在论坛搜电子版手册查一查(关键词写作没有横线的def2SVP)。

B3LYP泛函对这种超价体系可能无法正确优化VSEPR(外加电子对是否“立体活性”)预测的立体构型,这点在社长博文http://sobereva.com/557末尾和http://bbs.keinsci.com/thread-49406-1-1.html讨论的很相似的[IF6]-例子中有说明。氙的氟化物的高精度计算结果参考Inorg. Chem. 2010, 49, 1, 261–270这篇文章,注意不要尝试效仿此文试图在Gaussian拿CCSD(T)高精度后HF波函数方法做优化。(从上一帖来看是化竞生高考完了趁暑假入门计算化学?如果是这样,[XeF6]-这种体系稍显复杂了,不太适合尝试。)
作者
Author:
楚明城    时间: 5 hour ago
tonganlhy 发表于 2025-8-12 13:11
小跑了一下。可以跑起来。不过第一轮SCF不收敛
参考此文解决
http://bbs.keinsci.com/thread-4829-1-1. ...

好的,谢谢老师
作者
Author:
楚明城    时间: 5 hour ago
Uus/pMeC6H4-/キ 发表于 2025-8-12 13:12
这个赝势基组/赝势看起来是def2-SVP的数据?如果确实是,def2系列除镧系锕系外元素的基组在Gaussian已经内 ...

是的,其实就是想去算一算电子对在这种体系会产生什么影响,谢谢老师的推荐的文章和建议
作者
Author:
zjxitcc    时间: 5 hour ago
补充一点,如果是出于某种明确的特殊原因 不得不计算[XeF6]-阴离子,此时应当给所有原子使用弥散函数,否则阴离子所在轨道可能没有充分的基函数描述。例如给F原子使用ma-def2-SVP, 6-31+G(d), def2-SVPD中的任意一个基组,给Xe原子使用ma-def2-TZVP或def2-TZVPD。由于体系是二重态,还应当检验UB3LYP波函数稳定性,确保使用稳定波函数进行结构优化、频率计算和波函数分析。
作者
Author:
楚明城    时间: 5 hour ago
zjxitcc 发表于 2025-8-12 15:11
补充一点,如果是出于某种明确的特殊原因 不得不计算[XeF6]-阴离子,此时应当给所有原子使用弥散函数,否则 ...

感谢建议,不过这些还是超出了我目前的菜鸡水平,以后了解多了再用您的建议尝试尝试
作者
Author:
zjxitcc    时间: 4 hour ago
本帖最后由 zjxitcc 于 2025-8-12 15:31 编辑
楚明城 发表于 2025-8-12 15:15
感谢建议,不过这些还是超出了我目前的菜鸡水平,以后了解多了再用您的建议尝试尝试

这与菜鸡水平没有关系,我不是在说“当前计算级别已经合理,如何提升计算精度”,而是在说一些基本注意事项,并不是说“我写了一个gjf文件跑起来了,结果就是可以用的、合理的”,仔细阅读下文学习如何选用合适的基组
谈谈量子化学中基组的选择http://sobereva.com/336
不要自己吭哧算了一个月,以为好像合理,到时候到公共场合展示的时候被人说有低级错误和硬伤。
作者
Author:
Uus/pMeC6H4-/キ    时间: 4 hour ago
楚明城 发表于 2025-8-12 15:15
感谢建议,不过这些还是超出了我目前的菜鸡水平,以后了解多了再用您的建议尝试尝试

十分推荐先读一读http://sobereva.com/355以及底下一长串相关链接来入门哦,其中提到的科音的培训有很详细的指导和充足的计算案例。计算化学确实有很多语言与传统的化学不太一样,也确实在如何选取计算级别获得合理结果上很有讲究。“计算级别”既包括楼上提到的基组的问题也涵盖我说的泛函的问题,而关于泛函可以看http://sobereva.com/272这篇。

要说选择小型体系做练习的话,我觉得单是前三周期元素的简单分子就有很多可以算的了,更高周期的不一定更难但容易有种种影响准确性的细节。很多专业计算的benchmark文献都会提供充足的例子和信息(特别是补充材料里有完整的结构坐标文件的),适合作为建模与计算的参考(但计算级别是另一回事,以论坛里推荐的为准,拿不准还可以在论坛提问),比如GMTKN55数据集(原文Phys. Chem. Chem. Phys., 2017,19, 32184-32215,结构在课题组网站上)。像你前一帖计算的水分子辅助质子转移的烯醇化反应,实际上就是GMTKN55中牵涉反应能垒(过渡态搜索)的WCPT18子集的关心对象。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3