激发态计算中,如何得到T1,T2激发态的相应的轨道信息?
已经通过对结构进行优化(# B3LYP/6-31G* opt,电荷和自旋多重度为0 3), 取优化得到的结构,用# B3LYP/6-31G* TD(triplet,nstate=5),电荷和自旋多重度为0 1, 进行了计算。
然后有一些问题不是特别明白,也查阅了资料,结果不太确认,希望各位前辈能多多帮助下。
第一个问题:就是表格中得到的T1, 和T2中的dominant excitation 的计算结果是不是就是通过刚刚计算的第二步td结果中的log文件对应的 excited state 1 和excited state 2 的数值进行相应换算得到的?
Excited State 1: Triplet-A 1.5734 eV 788.01 nm f=0.0000 <S**2>=2.000
274 -> 279 0.19704
278 -> 279 0.63751
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) = -3290.38532782
Copying the excited state density for this state as the 1-particle RhoCI density.
Excited State 2: Triplet-A 1.8788 eV 659.91 nm f=0.0000 <S**2>=2.000
274 -> 279 0.15985
275 -> 279 -0.20843
276 -> 279 0.15620
276 -> 282 0.11142
277 -> 279 0.36316
277 -> 281 -0.18816
278 -> 279 -0.19542
278 -> 280 -0.34433
278 -> 281 -0.11651
Excited State 3: Triplet-A 1.9378 eV 639.81 nm f=0.0000 <S**2>=2.000
275 -> 279 -0.11025
276 -> 279 0.16483
276 -> 281 -0.16490
276 -> 282 -0.11299
277 -> 279 0.17993
277 -> 280 -0.25444
277 -> 281 -0.12976
277 -> 282 0.15435
278 -> 280 0.46294
278 -> 281 -0.11587
第二个问题:关于图中的T1, T2 激发态的character和assignment 是怎么分析得到的?是不是需要在T1态的计算下,添加些关键词然后需要分析些数据得到,比如pop=full ,还是说是其他的什么关键词? 同样的,是否还需要对T2态进行一些相应的计算才能得到结果?
返回小木虫查看更多
第一个问题第二个问题一起解决 是的 加pop=full 你要通过高斯sum做轨道成分分析 查看跃迁轨道的成分 统计大于10%的homo lumo 得到跃迁属性 关于mlct成分 可以做nbo得到
非常谢谢你的回答。还有几个细节方面的小问题,需要请教下。就是对于统计大于10%的H/L的轨道的,这个我明白。
就是后面的激发类型,电子由pi激发到反键轨道pi*的等类似的分析,以及后面判断属于所属的电荷转移激发类型,是不是通过上面的第二步计算结果得到的H/L.cube图形分析得到?还是我的理解是不对的?而MCLT的成分,向你说的,需要继续做个NBO的分析,我的最开始计算中没有添加这个关键词,是不是要继续重新添加了之后算?
这样跟你举个例子 在你的轨道成分图中HOMO的成分为(L1(25)+M(50)+L2(25)) LUMO的成分为(L1(5)+M5)+L2(90)) 则H->L为π(L1)+M+π(L2)->π*(L2)
至于NBO你需要做新的计算加关键词 当然在http://bbs.keinsci.com/thread-6355-1-3.html这里有一些其他很好的方法 做片段电荷 我们可以好好交流一下
,
第一个逗号之前就错了
额? 第一个逗号? 你指的是优化T1结构的输入有问题?