计算XeF6负离子结果莫名报错的问题

楚明城

各位老师好，我在计算XeF6负离子结构时Gaussian莫名报错，也没有看出错误原因在哪，请问一下这是什么情况出问题了呢

tonganlhy

输入文件也要贴上来。不过一眼就看出来自旋多重度错了

楚明城

tonganlhy 发表于 2025-8-12 12:41
输入文件也要贴上来。不过一眼就看出来自旋多重度错了

您好，自旋多重度确实输错了，但修改之后还是报错了

tonganlhy

所以说输入文件也要发上来。以后问报错输出输入文件都要放。
我估计你是赝势基组输错了

楚明城

tonganlhy 发表于 2025-8-12 12:55
所以说输入文件也要发上来。以后问报错输出输入文件都要放。
我估计你是赝势基组输错了

好的，已经放上来了，赝势是从BES上拷贝的，应该没啥问题？

tonganlhy

楚明城 发表于 2025-8-12 12:58
好的，已经放上来了，赝势是从BES上拷贝的，应该没啥问题？

小跑了一下。可以跑起来。不过第一轮SCF不收敛
参考此文解决
http://bbs.keinsci.com/thread-4829-1-1.html

Uus/pMeC6H4-/キ

这个赝势基组/赝势看起来是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]-这种体系稍显复杂了，不太适合尝试。）

楚明城

tonganlhy 发表于 2025-8-12 13:11
小跑了一下。可以跑起来。不过第一轮SCF不收敛
参考此文解决
http://bbs.keinsci.com/thread-4829-1-1. ...

好的，谢谢老师

楚明城

Uus/pMeC6H4-/キ 发表于 2025-8-12 13:12
这个赝势基组/赝势看起来是def2-SVP的数据？如果确实是，def2系列除镧系锕系外元素的基组在Gaussian已经内 ...

是的，其实就是想去算一算电子对在这种体系会产生什么影响，谢谢老师的推荐的文章和建议

zjxitcc

补充一点，如果是出于某种明确的特殊原因 不得不计算[XeF6]-阴离子，此时应当给所有原子使用弥散函数，否则阴离子所在轨道可能没有充分的基函数描述。例如给F原子使用ma-def2-SVP, 6-31+G(d), def2-SVPD中的任意一个基组，给Xe原子使用ma-def2-TZVP或def2-TZVPD。由于体系是二重态，还应当检验UB3LYP波函数稳定性，确保使用稳定波函数进行结构优化、频率计算和波函数分析。

楚明城

zjxitcc 发表于 2025-8-12 15:11
补充一点，如果是出于某种明确的特殊原因 不得不计算[XeF6]-阴离子，此时应当给所有原子使用弥散函数，否则 ...

感谢建议，不过这些还是超出了我目前的菜鸡水平，以后了解多了再用您的建议尝试尝试

zjxitcc

楚明城 发表于 2025-8-12 15:15
感谢建议，不过这些还是超出了我目前的菜鸡水平，以后了解多了再用您的建议尝试尝试

这与菜鸡水平没有关系，我不是在说“当前计算级别已经合理，如何提升计算精度”，而是在说一些基本注意事项，并不是说“我写了一个gjf文件跑起来了，结果就是可以用的、合理的”，仔细阅读下文学习如何选用合适的基组
谈谈量子化学中基组的选择http://sobereva.com/336
不要自己吭哧算了一个月，以为好像合理，到时候到公共场合展示的时候被人说有低级错误和硬伤。

Uus/pMeC6H4-/キ

楚明城 发表于 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子集的关心对象。