24小时热门版块排行榜    

查看: 2391  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chenchu

铜虫 (正式写手)

[求助] 困扰:计算激发态,使用TD-DFT,输出结果中激发态的多重度怎么确定? 已有2人参与

新手求助,刚开始学习计算B3S2-激发态,使用TD-DFT,先计算了B3S2-的稳定态,是三重态,再计算其激发态,输入文件如下:
#p ub3lyp/6-311+g(d,p) TD(Nstates=10) scf=tight

B3S2-design-b3lyp

0 2
B                  0.00000000    0.00000000    0.00000000
B                  0.00000000    0.00000000    1.56067600
B                  0.00000000    0.00000000   -1.56067600
S                  0.00000000    0.00000000    3.22385600
S                  0.00000000    0.00000000   -3.22385600
输出文件中激发态找到了。是这样的:

Excited state symmetry could not be determined.
Excited State   1:  1.994-?Sym    0.0075 eV 166094.52 nm  f=0.0000  <S**2>=0.744
     20A -> 25A        0.31398
     24A -> 25A        2.58499
     24A -> 37A        0.13714
     20A <- 25A        0.29554
     24A <- 25A        2.38674
     24A <- 37A        0.13050
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -870.978868190   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited state symmetry could not be determined.
Excited State   2:  3.435-?Sym    1.8924 eV  655.18 nm  f=0.0000  <S**2>=2.701
     21A -> 27A        0.16214
     23A -> 25A        0.75825
     20B -> 26B       -0.15358
     22B -> 24B       -0.61745
     23B -> 25B       -0.13630
     23A <- 25A        0.10672
     22B <- 24B       -0.10353
后面的省略,情况都差不多。
1.请教其他人说多重度应该是Excited State   1:  1.994-?Sym 中1.994这个位置上的数,但是别人的体系输出是类似Singlet等形式,我这个体系输出不是整数也不是英文,无法确定多重度。
2.在Excited State   1里还出现了 24A -> 25A        2.58499,2.588499正常吗?请教其他人说是跃迁的几率,不应该超过1,但是我这里还有负数的情况,到底问题出在哪?
请各位大神帮忙啦!求正解!!或者在哪里可以看到相关资料!万分感谢!悬赏20金币!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xieruobing

木虫 (著名写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by beefly at 2014-12-25 17:11:30
gaussian的开壳层tddft不考虑自旋反转(这从轨道跃迁也能看出来),因此自旋多重度和参考态相同。但是由于自旋污染,会有一定的偏差。

专家,小的可否加你QQ,或者邮件咨询啊!!
12楼2015-01-10 11:55:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

beefly

专家顾问 (职业作家)

地沟油冶炼专家

gaussian的开壳层tddft不考虑自旋反转(这从轨道跃迁也能看出来),因此自旋多重度和参考态相同。但是由于自旋污染,会有一定的偏差。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2014-12-25 17:11:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

更正:刚发现Gaussian09的开壳层TDDFT已经可以考虑自旋反转了。自旋多重度就是1.994-?Sym前面的数值取整,即自旋多重度=2。但是第二个态的自旋污染太严重了。

跃迁因子(注意:不是跃迁几率)平方和应当接近0.5(闭壳层TDDFT)或1.0(开壳层TDDFT)。但是从第一激发态的能量来看,它可能和基态简并,而且可能涉及两个电子的激发,或者参考态是多参考体系。TDDFT无法正确描述这样的电子态,因此会得到奇怪的跃迁因子。

计算结果说明TDDFT不适用这个体系,应当改用CASPT2一类的方法
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2014-12-25 17:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
beefly说的很对了。
我感觉你说的第二个问题,加起来的和不是1
这个,你应该看箭头的方向,   23B -> 25B这样的是平方相加。     22B <- 24B   这样的是,平方相减。
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
4楼2014-12-26 15:29:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 by 9i7PCSZyMzot
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:15 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 4pRyVGCVerAA 2026-09-19 4/200 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[教师之家] 现代”学阀”该如何界定 +3 iaeyuan 2026-09-19 3/150 2026-09-19 20:52 by 没昵称呀
[基金申请] 国社科系统bug了,是不是要放榜了? +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
信息提示
请填处理意见