大侠好手段

sobereva 发表于 2016-12-15 21:55
你直接让Gaussian输出单重态基态到三重态激发态的跃迁密度矩阵，即使用density=transition=x IOp(6/8=3)肯 ...

太感谢sob老师的回复，一步步写的非常详细！

你直接让Gaussian输出单重态基态到三重态激发态的跃迁密度矩阵，即使用density=transition=x IOp(6/8=3)肯定行不通，因为波函数自旋部分正交，所以跃迁密度都是零，因此也没法用Multiwfn绘制。

对策就是让Multiwfn来产生S0->T的跃迁密度矩阵，Multiwfn比较聪明，会忽略掉自旋部分，而只考虑空间部分，因此对S0->T也可以给出跃迁密度矩阵，并进而绘图。

具体就是计算个体系，比如
%chk=c:\nozomi.chk
#p cam-b3lyp/6-31g(d) td=triplet IOp(9/40=4)

然后在Multiwfn里输入
nozomi.fch
18
1
nozomi.out
1  //假设对S0->T1感兴趣
4  //产生跃迁密度矩阵
y  //对跃迁密度矩阵做对称化
就在当前目录得到了tdmat.txt。

重新启动Multiwfn，输入
nozomi.fch
18
2
tdmat.txt
然后就可以照常绘制跃迁密度矩阵图了

brothers 发表于 2016-12-15 19:29
不考虑旋轨耦合的话S0-->Tn的跃迁密度矩阵当然都是0，跃迁偶极矩也是0

S0-->Tn的跃迁密度矩阵都是0,那么跃迁密度矩阵平面图肯定就得不到了，想分析s0-T1的电子跃迁所涉及到的原子以及原子间相干范围，不知道有什么有效的方法？

不考虑旋轨耦合的话S0-->Tn的跃迁密度矩阵当然都是0，跃迁偶极矩也是0

[已删除]没弄懂问题