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

第二步中，需要计算S1态的振动频率，计算方法参考：如何进行单重激发态的振动频率计算

S1态振动频率计算得到的S1_OptFreq.t21，那么在后面第三步的过程中，替换T1_OptFreq.t21即可。

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