计算方法与emissionofpho一致，唯一差别：

第二步中，需要计算S1态的振动频率，计算方法参考：如何进行单重激发态的振动频率计算，对于原子个数较多的分子的激发态频率计算，参考：有没有特别高效的基于TDDFT的激发能计算方法？

假定S1态振动频率计算得到S1_OptFreq.t21

第三步，用S1_OptFreq.t21替换T1_OptFreq.t21即可。

第四步中的“S1振动基态”→“S0振动基态”跃迁能应该是S1态优化之后的S1激发能。