计算化学公社

标题: 求助:利用TDDFT单独计算S1及同时计算S1T1(50-50)的结果不一样 [打印本页]

作者
Author:
q1w2e3    时间: 2019-12-11 13:30
标题: 求助:利用TDDFT单独计算S1及同时计算S1T1(50-50)的结果不一样
本帖最后由 q1w2e3 于 2019-12-11 13:32 编辑

各位老师好!想请教一下,我用TDDFT单独计算S1及同时计算S1T1(50-50)的S1结果不一样,他们都是用优化好的S0进行计算的,N都是5。

1、单独计算S1: #p td=(singlets,nstates=5,root=1) b3lyp/6-31g(d)

Excited State   1:   Singlet-A      3.2590 eV  380.44 nm  f=0.0810
     146 ->154        -0.10135
     153 ->154         0.57258
     153 ->155        -0.32178
This state for optimization and/or second-order correction.
Total Energy, E(RPA) =  -2169.48022981   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:   Singlet-A      3.3366 eV  371.59 nm  f=0.0027
     149 ->155         0.10448
     150 ->155        -0.10030
     151 ->154        -0.26769
     152 ->154         0.15719
     152 ->155         0.28622
     153 ->154         0.22372
     153 ->155         0.44302

Excited State   3:   Singlet-A      3.4660 eV  357.72 nm  f=0.1315
     151 ->154         0.29819
     152 ->154         0.26904
     152 ->155         0.47927
     153 ->155        -0.20738


2、同时计算S1T1: #p td=(50-50,nstates=5,root=1) b3lyp/6-31g(d)

Excited State   6:   Singlet-A      3.1133 eV  398.25 nm  f=0.0739
     153 ->154         0.55966
     153 ->155        -0.31449

Excited State   7:   Singlet-A      3.2590 eV  380.44 nm  f=0.0026
     149 ->155         0.10326
     151 ->154        -0.26459
     152 ->154         0.15535
     152 ->155         0.28289
     153 ->154         0.22113
     153 ->155         0.43784

Excited State   8:   Singlet-A      3.3366 eV  371.59 nm  f=0.1218
     151 ->154         0.29255
     152 ->154         0.26398
     152 ->155         0.47031
     153 ->155        -0.20351



更奇怪的是:同时计算S1TI中的S2对应单独计算S1中的S1,请问出现这种情况的原因是什么?真正的S1到底是谁?

作者
Author:
sobereva    时间: 2019-12-13 23:33
早期G09的50-50有bug,用较新版本再试。如果还这样,以单独算的为准
作者
Author:
q1w2e3    时间: 2019-12-14 09:54
sobereva 发表于 2019-12-13 23:33
早期G09的50-50有bug,用较新版本再试。如果还这样,以单独算的为准

好的,谢谢社长!!
作者
Author:
chmdong    时间: 2020-5-29 12:37
我想请问一个问题
Excited State   1:   Singlet-A      3.2590 eV  380.44 nm  f=0.0810
     146 ->154        -0.10135
     153 ->154         0.57258
比如从146到154轨道中,对应singlet-A下面的-0.10135是什么意思,而且此值的正负代表什么?
作者
Author:
xuhj199508    时间: 2020-5-29 13:23
chmdong 发表于 2020-5-29 12:37
我想请问一个问题
Excited State   1:   Singlet-A      3.2590 eV  380.44 nm  f=0.0810
     146 ->154 ...

请看《电子激发任务中轨道跃迁贡献的计算》http://sobereva.com/230
作者
Author:
chmdong    时间: 2020-5-29 13:30
已明了,多谢!




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3