关于CT态

肖富贵

大家好，有个小问题，请讨论一下。

我用QM/MM方法计算蛋白中色氨酸的激发能，感兴趣的是两个态，S1态和CT态，因为色氨酸通过CT态迅速回到基态。

我首先用DFT优化基态结构，然后用cam-b3lyp做单点，找到CT态，最后用cam-b3lyp优化该CT态。（关于CT态的寻找，纯属感觉。比如某激发态中轨道a-->b的激发对应CT，并且轨道a-->b的系数最大，则把该激发态当成）。
  优化结束后，进行结果分析，发现最后优化到的态并非CT态。比如cam-b3lyp优化前，CT是第五激发态；优化结束后，第五激发态对应的不再是CT态，相反，另一个态，比如第三激发态，变成了CT态。
为了处理这种复杂的情况，我自己想到一种方法（可能很幼稚），也就是依然用cam-b3lyp优化CT态，得到结果后，用cam-b3lyp(nstates=20,singlet)计算多个激发态，然后用激发能乘以一个系数，这个系数是该激发态中CT激发的轨道所占的比例。举个例子，比如某激发态包含两个轨道激发，a-->c, b-->d。并且a-->c占40%，对应CT态；b-->d占60%，不是CT激发。那么用该激发能乘以40%，可以得到一个值。对这20个激发态每个激发态都做这样的运算，然后将得到的值最后相加，把最后得到的值定义为CT态的能量。


不好意思，非常冗余，但是也是想尽量把问题说清楚。

请有兴趣的朋友讨论一下，谢谢！

sobereva

这么算CT态的能量没什么物理意义。
优化第5激发态过程中，可以一直跟踪占主导成份的轨道跃迁的系数，看看怎么变化的。
第5激发态优化后，也可以在其结构基础上再优化第3激发态，看看结果最后它是否是CT态。