计算荧光与磷光时涉及到三重态的结构优化以及吸收能和发射能求助

zhouyulin

学习了sob老师《Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法》，里面讲解了两种方法：
“方法一：
1 优化T1结构：# B3LYP/6-31G* opt，电荷和自旋多重度为0 3。
2 取上一步优化得到的结构，用# B3LYP/6-31G* TD(triplet)，电荷和自旋多重度为0 1。第一激发态的激发能就是磷光发射能了。

方法二：
1 优化T1结构：# B3LYP/6-31G* opt，电荷和自旋多重度为0 3。取最后一步的能量，作为T1结构下T1能量。
2 取上一步优化得到的结构，0 1下用# B3LYP/6-31G*计算单点能，得到T1结构下S0能量。
3 将第1步的能量减去第2步的能量即是磷光发射能”


我想请教一些相关的衍生问题是否是这样：
（1）s0基态几何优化：# B3LYP/6-31g* opt ，电荷和自旋多重度为0 1 ；
（2）取上一步s0优化的结构，# B3LYP/6-31g* td opt，电荷和自旋多重度为0 1，得到的结果就是s1优化的结构
问题1：依旧取s0优化的结构，分别设置关键词# B3LYP/6-31g* td(triplet) opt和# B3LYP/6-31g* td(triplet,root=n) opt，电荷和自旋多重度为0 1，这样计算出来的就算T1和Tn优化好的结构？
问题2：如果问题1的情况成立，那这样的话对T1优化的结构# B3LYP/6-31g* td(triplet) ，电荷和自旋多重度为0 1，这样得到的是不是T1的磷光发射能。同理，如果按照这样的方法对Tn优化的结构也进行# B3LYP/6-31g* td(triplet) ，电荷和自旋多重度为0 1的操作，是不是得到Tn和S0之间的磷光发射能，Tn到T1之间的激发能是不是直接相减就可以了？
问题3：用BDF计算soc的话，是不是直接在s0优化的结构上进行就可以了？
问题3：按照这个情况，假设需要用momap计算s1-T2的系间窜越速率，可以直接用上述过程已经得到的s1激发态优化好的结构和t2激发态优化好的结构输入条件，分别用BDF在s0优化的结构上计算得到s1-t2的soc以及用momap计算得到s1-t2的evc.cart.dat文件，最后再通过momap计算得到s1-t2的输入文件。想看看这样是否合理
问题4：之前在查找有关基态和激发态偶极矩的时候，提到了density这个关键词，那么在s0优化好的结构上进行# B3LYP/6-311+g* td density计算是不是计算的是s1的激发态偶极矩，同理如果是# B3LYP/6-31g* td(root=n) density是不是计算的是Sn的激发态偶极矩，# B3LYP/6-31g* td(triplet，root=n) density是不是计算的就算Tn的三重激发态偶极矩

sobereva

1 是

2 “B3LYP/6-31g* td(triplet) ，电荷和自旋多重度为0 1”最后一次输出的第一激发态的激发能是T1的垂直磷光发射能

“B3LYP/6-31g* td(triplet,root=n) ，电荷和自旋多重度为0 1”最后一次输出的第n激发态的激发能是Tn的垂直磷光发射能

优化较高激发态容易牵扯到势能面交叉，建议参考下文
谈谈势能面交叉对激发态优化的影响
http://sobereva.com/468（http://bbs.keinsci.com/thread-12318-1-1.html）

3 算SOC用什么结构取决于具体目的。论坛首页全文搜索框搜，论坛里早就有大量讨论

4 是