关于使用stable测试波函数稳定性

小强

背景: 一反应体系推测其机理由三重态双自由基(A)经ISC过程形成单重态自由基(B)，随后通过B进一步反应生成产物(C)。

计算流程: 在实际计算过程中，我是通过使用关键词 guess=mix 并结合U(UB3LYP/6-31G*)进行的，并且顺利找到了该过渡态结构；随后通过stable关键词在UB3LYP/6-31G*计算水平下对该结构进行波函数稳定性测试，结果显示稳定；在B3LYP/6-31G*下测试时，提示internal instablity。

问题：  （1）上述结果是不是说明当前的计算级别是合理的，结果能正常使用？对于后续高精度单点的计算是不是也应该通过guess=mix结合U来实现？
            （2）如果stable测试在B3LYP/6-31G*下提示稳定，那是不是说明当前计算级别的能量数据不能用，需要读取稳定波函数的数据重新在闭壳层下寻找过渡态结构？

(该问题在群里问过，但不是太理解老师给的回复，故重新发一遍，谢谢大家)

sobereva

1 如果你在用UB3LYP/6-31G*做stable的时候，初始波函数是对称破缺的，而且直接收敛到的波函数也是对称破缺的，此时stable提示稳定，就说明此时对称破缺波函数就是当前结构下基态波函数
2 能用。本身是闭壳层的体系，你按对称破缺方式来算，结果自然而然就收敛到了和闭壳层计算相同的结果。

小强

sobereva 发表于 2018-4-23 11:18
1 如果你在用UB3LYP/6-31G*做stable的时候，初始波函数是对称破缺的，而且直接收敛到的波函数也是对称破缺 ...

懂了，谢谢Sob老师。