使用gaussian考虑溶剂效应时磷光发射计算的问题求助

白也

在sob老师写的Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法http://sobereva.com/314一文中写到了基于线性响应溶剂模型计算磷光发射能这里对于计算步骤有一些疑惑：
第一步是（1）# PBE1PBE/6-311G* scrf=solvent=ethanol opt结合电荷和自旋设置0 3，这一步是通过SCF方式直接算T1态并优化，可以节省时间，
得到T1态的极小点结构后，那么接下来的步骤应该是2a还是2b呢？

（2a）# PBE1PBE/6-311G* TD(triplet) scrf=solvent=ethanol  直接读取S0->T1激发能吗？这样是不是溶剂效应考虑反了？
因为sob老师提到一句，获得荧光发射能不要取S1极小点的结构再单独做一次电子激发计算得到S0->S1的能量，因为此时激发态处于非平衡溶剂，和发射光谱的溶剂实际特征恰相反。

（2b）# PBE1PBE/6-311G* TD(triplet) scrf=solvent=ethanol opt  直接取输出文件当中最后一次的S0->T1激发能？

麻烦各位大佬指点一二

wzkchem5

白也 发表于 2022-10-12 12:59
非常感谢，我明白了，是不是用eqsolv来明确当前T1激发态是平衡的，这样程序就知道s0->T1这个过程是非平衡 ...

对
有可能，所以先（1）再（2b）有可能比一上来就（2b）快，但是肯定不如先（1）再（2a加上eqsolv）快

白也

wzkchem5 发表于 2022-10-12 19:13
（2b）虽然可以，但是就不太能起到（1）的“节省时间”的目的了，等于你的T1态优化还是在TDDFT级别下做的。 ...

非常感谢，我明白了，是不是用eqsolv来明确当前T1激发态是平衡的，这样程序就知道s0->T1这个过程是非平衡溶剂到平衡溶剂的是吗？感觉通过这个关键词的方式更好一些，这样就可以用中等基组优化结构，大基组来只是算个能量。
我还有一个小疑问，如果说我用SCF的方式已经得到T1的极小点结构，那么我把这个结构作为初始构型用TD优化，因为本身初始就是极小点结构了，会不会很快没几步就优化收敛然后结束呢？

wzkchem5

（2b）虽然可以，但是就不太能起到（1）的“节省时间”的目的了，等于你的T1态优化还是在TDDFT级别下做的。
要想又用SCF优化T1，又正确地考虑溶剂效应，可以用td(triplet,eqsolv)，参见http://gaussian.com/td/。另一个方法是用（2b），但是读第一步结构优化的S0->T1激发能，不过为了不浪费机时，得在第一步结构优化的S0->T1激发能输出以后把任务杀掉。