计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8665|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助:利用TDDFT单独计算S1及同时计算S1T1(50-50)的结果不一样

[复制链接 Copy URL]

12

帖子

0

威望

99

eV
积分
111

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 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到底是谁?

6万

帖子

99

威望

6万

eV
积分
125181

管理员

公社社长

2#
发表于 Post on 2019-12-13 23:33:15 | 只看该作者 Only view this author
早期G09的50-50有bug,用较新版本再试。如果还这样,以单独算的为准
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

12

帖子

0

威望

99

eV
积分
111

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2019-12-14 09:54:25 | 只看该作者 Only view this author
sobereva 发表于 2019-12-13 23:33
早期G09的50-50有bug,用较新版本再试。如果还这样,以单独算的为准

好的,谢谢社长!!

9

帖子

0

威望

37

eV
积分
46

Level 2 能力者

4#
发表于 Post on 2020-5-29 12:37:44 | 只看该作者 Only view this author
我想请问一个问题
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是什么意思,而且此值的正负代表什么?

105

帖子

0

威望

3323

eV
积分
3428

Level 5 (御坂)

Yukikaze

5#
发表于 Post on 2020-5-29 13:23:51 | 只看该作者 Only view this author
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

9

帖子

0

威望

37

eV
积分
46

Level 2 能力者

6#
发表于 Post on 2020-5-29 13:30:18 | 只看该作者 Only view this author
已明了,多谢!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-26 02:33 , Processed in 0.155411 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list