计算化学公社

标题: 求助对单个过渡金属原子计算单点能（为了计算分子结合能） [打印本页]

作者
Author: JTDJZS 时间: 2025-10-18 22:40
标题: 求助对单个过渡金属原子计算单点能（为了计算分子结合能）
各位老师好，我想问一下，我是相对过渡金属内嵌C60体系进行结合能的计算，在对单个过渡金属原子进行单点计算时，出现502不收敛，尝试了SCF(novaracc,noincfock) int=acc2e=12，SCF=vshift=x，SCF=conver=N后仍然失败，于是采用XQC关键词，成功得到结果，但是在进行XQC计算时我加入了stable关键词，输出文件中却并未查找到关于波函数是否稳定的结果，想请问这是什么原因，以及这样的单点能结果是可靠的吗？附上我的输入输出文件如下

作者
Author: KazusaT 时间: 2025-10-19 04:00
你的输出文件写了“ The wavefunction has an RHF -> UHF instability.”
Ti原子基态肯定是三重态，你设置单重态难收敛很正常。

作者
Author: sobereva 时间: 2025-10-19 06:13
解决SCF不收敛问题的方法
http://sobereva.com/61

这都是博文里明确说过的

对于过渡金属配合物，自旋多重度设定得不符合实际基态（非基态的情况往往更难收敛）

并且博文也说了，纯泛函往往比杂化泛函更难收敛，HF成份高的泛函往往更容易收敛。色散校正的TPSSh、PBE0更推荐用于当前计算

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