G16做激发态优化，反复震荡经多次调整无效

阿锋001

目的：得到某有机分子的发射能
过程：（1）采用G16计算，首先关键词为#p PBE1PBE/6-31g(d,p) td opt=cartesian freq scrf=(smd,solvent=thf)，然后发生震荡，
（2）关键词调整为#p PBE1PBE/6-31G(d,p) TD opt=(cartesian,maxcyc=300,maxstep=5,notrust,calcfc) freq scrf=smd继续优化，依然震荡，
（3）#p PBE1PBE/6-31G(d,p) TD opt=(gdiis,recalc=3,maxcyc=300,maxstep=3,notrust,tight) freq scrf=(smd,solvent=thf)，关键词越写越多，依然不行
（4）也更换过计算方法为M062X，不行
说明：第一次初始结构拿的是基态优化后的稳定结构，后面调整用的都是上一次震荡下能量最低的结构。
        每次错误基本都是 l9999，并往上找到了Optimization stopped.
                             -- Number of steps exceeded,  NStep= 168
                             -- Flag reset to prevent archiving.
       而且有类似的结构用类似的方法优化和调整，已经优化成功了，有几个就一直这样不行，所以想请教各位怎么办？谢谢！

pwzhou

你可能遇到了两个态存在交叉的情况，这种情况下TDDFT没啥好的办法，你可以用写字板打开输出文件，看看最后计算的激发能是不是都不是按顺序排列了，比如第一激发态380nm，第二激发态420nm这样。社长写过一篇博文，http://sobereva.com/468。你可以好好看看，然后仔细看看输出文件，看看是什么原因导致的，好对症下药，盲目乱试是解决不了问题的。

hunterpyj

如果发生势能交叉的情况只会导致最终优化的结构不是一开始想要算的结构，也就是态追踪失败的情况出现吧，为什么会出现优化震荡的问题呢？