计算含有铜元素分子NTO时遇到三重态T1能量为负的情况

lmyqwq

首先非常感谢各位老师的指导，我在计算该分子的激发态时遇到了问题，对于该分子的优化和激发。我分别使用了下面的命令，
OPT：
%chk=1.chk
%mem=100GB
%nprocshared=24
#p opt=tight PBE1PBE/6-31g(d) em=GD3BJ

222.cdxml

0 1
C                 -7.70955900   -2.44450000   -0.61218100
C                 -8.61012900   -3.52916500   -1.06780800
C                 -6.38010200   -2.47786500   -1.00549500
下面省略。

计算完后进行频率优化
%chk=2.chk
%mem=100GB
%nprocshared=24
#p opt freq PBE1PBE/6-31g(d) em=GD3BJ

222.cdxml

0 1
C                 -7.70960800   -2.44449500   -0.61249100
C                 -8.61011900   -3.52915200   -1.06824900
C                 -6.38009500   -2.47786600   -1.00562000
C                 -8.17684300   -1.48495800    0.29803400


优化完频率后进行激发计算，
%chk=3.chk
%mem=100GB
%nprocshared=24
#p td(nstates=5,50-50) IOp(9/40=4) CAM-B3LYP/6-31g(d)

222.cdxml

0 1
C                 -7.70960800   -2.44449500   -0.61249100
C                 -8.61011900   -3.52915200   -1.06824900
C                 -6.38009500   -2.47786600   -1.00562000


但是得到的T1态能量为负值
Excited State   1:      Triplet-A      -0.7232 eV    -1714.40 nm  f=-0.0000  <S**2>=2.000
     173 -> 484       -0.00058
     173 -> 508        0.00012
     173 -> 511        0.00013
     173 -> 513        0.00016


优化结束后将其放入gview中的模型也并不符合预期。
恳请各位老师帮忙解答一下学生那里操作有误，不胜感激。

xsc6

确定你的电荷和自选多重度是否正确

wzkchem5

两个铜是反铁磁耦合吧？如果是，必须做对称性破缺处理，或者做spin flip TDDFT；如果不是，那么基态应该是三重态。按闭壳层单重态来算肯定是不对的

sobereva

老生常谈的问题
参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）的ppt：