请教在优化激发态时波函数的稳定性问题

DUNING

老师好，我现在在做有关荧光发射的计算，在优化激发态S1的构型时，发现到最后几步第一激发态垂直发射能为负值，看了您之前的帖子，采用了stable=opt关键词寻找稳定基态波函数，刚开始输出文件中有“The wavefunction has an RHF -> UHF instability.”这句话，之后高斯可能又寻找了稳定的基态波函数，输出了“The wavefunction is stable under the perturbations considered.”,此时输出的<S**2>=0，所以后来我又用上一步产生的chk文件作为优化激发态时的波函数读入，发现在之后的激发态优化中，第一激发态后面的 <S**2>=2.000，也就是S(S+1)=2，而且方法RPEB1PBE全部变成了UPBE1PBE,难道说第一激发态的自旋多重度成了3？这不就成了磷光发射了吗？

sobereva

基态当成U来计算，这时候TDDFT没法指定要算哪个自旋多重度的激发态，恰好算出来的最低激发态就是T1。你可以从所有激发态里看<S**2>=0是第几个态，然后用root指定要求优化那个态