gaussian09计算ΔEst的疑问

qdzzzlj

各位前辈好：
      本人新手，想要计算化合物（化合物组成为C33H35I2N2O4+）的ΔEst，用的泛函是PBE1PBE，基组和赝势分别是6-31G**和SDD，用的方法是Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法 - 思想家公社的门口：量子化学·分子模拟·二次元网站中的4.2和4.3，在学习ΔEst计算方法的时候有一些疑问：

有的帖子中的计算方法：
delta_ST绝热能=S1的绝热能-T1的绝热能
S1的绝热能=E(TD-HF/TD-KS)+HF能量    td opt或tda opt计算
T1的绝热能=E(TD-HF/TD-KS)+HF能量    td(triplets) opt或tda(triplets) opt计算，或 opt使用0 3
对于HF/DFT计算  HF能量=SCF Done的能量


另一种方法：
0   1               
计算S1的能量，从输出文件中提取最后一轮的第一激发态的E(TD-HF/TD-KS)能量，记为Es1
0   3
提取最后的HF的能量，记为T1的能量，Et1
然后这个作为绝热的delta ST=Es1-Et1


因为这两种方法计算的结果差值比较大，不太清楚哪种方法更加准确，或者说有没有更准确的方法能指点我一下，想请教一下各位前辈。

wzkchem5

仔细阅读https://mattermodeling.stackexch ... e-using-o/6468#6468以及https://mattermodeling.stackexch ... -triplets/9357#9357

sobereva

E(TD-HF/TD-KS)+HF能量根本什么都不是。E(TD-HF/TD-KS)本身就已经是root指定的激发态的能量了，等于参考态的能量（SCF Done的能量）加上相应态的激发能。