三重态铀酰配合物收敛失败

AndrewX

附件中1.gjf是为多重度为1的铀酰配合物阴离子结构，优化正常结束得到了稳定的波函数，将其存为输入并将多重度改为3计算(即1triplet.gjf)，便无法收敛。
基本上碰到L401(Inaccurate quadrature in CalDSu)和L9999(震荡才加，不震荡就继续让其往下算)就加int=ultrafine，碰到L502就加SCF=QC。
确实发现三重态的结构算着算着就和单重态的越来越偏离，但是就抱着不管结构偏离多少先得到稳定波函数的想法按照上述解决办法一路算。
但是最后的结果却是结构偏离到完全无法接受，且最终也仍未得到稳定的波函数。
同时还在算的一大批类似结构也出现了同样的问题，即用单重态的稳定结构去算三重态却无法收敛。
在网上检索没有得到靠谱的信息，求老师解答，感谢！

PS：势能面的起点是一个六价铀酰，中间会经历四个过渡态，最后得到的是五价铀酰，但势能面的能量不是很符合实验，因此考虑到自旋交叉，三重态的势能面可能在某处远低于单重态的势能面，因此将单重态的所有结构都改为三重态去优化。

sobereva

“碰到L502就加SCF=QC”是极度糟糕且危险的，是下下下策
仔细看
解决SCF不收敛问题的方法
http://sobereva.com/61
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

你都不做stable怎么确认波函数的稳定性，难道以为几何结构收敛了波函数就稳定了？

轻易用scf=QC只会严重提升收敛到不稳定波函数的几率。如果几何优化第一步的波函数都不稳定（至少对于当前算的自旋态来说），之后的优化都毫无意义