偶氮苯一价负离子的激发态计算问题请教

VASP-DFT

大家好，想请教大家一个问题，我算了偶氮苯一价负离子的激发态，用的输入文件是
%chk=trans.chk                                                  
%mem=6GB                                                        
%nproc=20                                                      
# b3lyp/gen TD=(nstates=30)

trans

-1  2

得到的结果是本应该是<S**2>=0.770代表单重态，但是不知道为什么有的态出现了<S**2>=1.033以及<S**2>=1.890这种代表三重态的结果呢？正常来说偶氮苯一价负离子的激发态不可能是三重态的。
希望大家可以帮忙解答一下，非常感谢！

sobereva

基于非限制性参考态波函数没法指定算出来的激发态的自旋多重度，算出来是什么就是什么，虽然基态指定的是二重态，但实际算出来的某些激发态有的有严重的自旋污染，甚至从<S**2>看上去似乎是不可能出现的自旋多重度。这种激发态就别管它就完了