通过优化波函数稳定性解决SCF不收敛问题

mfdsrax2

之前优化某自由基过渡态结构时，SCF很难收敛，后来发现通过先检查并优化波函数稳定性，用stable=opt关键词，然后再读取上一步的波函数，可以有效解决SCF不收敛的问题。
当然，仅仅针对我个人计算的结构，样本空间较小，也没有处理过限制性优化方法，是否具有普适性还待进一步研究。
看了61号博文http://sobereva.com/61，没有发现介绍这种方法，可以供大家参考。

更新：输入输出文件已传附件，Direct是直接用opt=(calcall,ts,noeigen)计算，第一步SCF就不收敛。没有Direct后缀的是先用stable=opt guess=mix关键词优化波函数，然后再读取优化的波函数进行计算，可以看出能够优化出目标结构。

Freeman

我也碰到过这种情况，不过不是自由基而是高自旋的锕系元素。

sobereva

当前情况类似于博文中说的

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

zjxitcc

没有“优化波函数稳定性”这种说法，最好给出详细关键词，以便大家提供更严谨的描述和分类。

mfdsrax2

zjxitcc 发表于 2024-6-15 20:08
没有“优化波函数稳定性”这种说法，最好给出详细关键词，以便大家提供更严谨的描述和分类。

这个功能你们叫什么，重新优化波函数？

zjxitcc

mfdsrax2 发表于 2024-6-16 14:54
这个功能你们叫什么，重新优化波函数？

（1）stable是检验波函数稳定性，不论结论是否稳定，都终止程序正常结束。stable=opt是先检验波函数稳定性，若稳定则程序正常结束；若不稳定，按照检测结果给出的能量下降方向 继续优化波函数，收敛后再检验波函数稳定性，如此循环直至波函数稳定为止。
（2）二重态的非限制性计算 不需要写guess=mix。您给出的gjf文件不能支持您在帖子中的观点，毕竟一个写了guess=mix，另一个没写，有可能是guess=mix无心插柳柳成荫（当然，我们不能因为这点而把guess=mix当成帮助SCF收敛的一种手段）。
（3）输入文件不要有24核2GB这种古怪搭配，24核搭配至少24GB为宜。