计算化学公社

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

[Gaussian/gview] 气相下计算T1能级

[复制链接 Copy URL]

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

  采用TD-DFT计算T1的能级,方法是这样的不知道对不对:
  1、对初始结构进行freq+opt优化,电荷、多重度0,1
  2、优化完的基态进行TD计算(图1)所示,log文件找到triplet 1。

以下是sober老师的方法(TD-DFT):
1 优化T1结构:# B3LYP/6-31G* opt,电荷和自旋多重度为0 3
2 取上一步优化得到的结构,用# B3LYP/6-31G* TD(triplet),电荷和自旋多重度为0 1。第一激发态的激发能就是磷光发射能了。

这里的第一步直接用opt就可以么?不需要用freq么? 另外我的计算(第一步)是否正确?
谢谢!

图1.PNG (3.96 KB, 下载次数 Times of downloads: 43)

图1.PNG

1294

帖子

0

威望

6236

eV
积分
7530

Level 6 (一方通行)

20#
发表于 Post on 2016-3-10 07:45:19 | 只看该作者 Only view this author
我咨询你一个问题,这个实验上说的这个dual fluorescence是什么意思?谢谢指点。
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

1294

帖子

0

威望

6236

eV
积分
7530

Level 6 (一方通行)

19#
发表于 Post on 2016-3-10 07:13:29 | 只看该作者 Only view this author
本帖最后由 小范范1989 于 2016-3-10 07:15 编辑
yezhonghua 发表于 2016-3-9 23:40
TADF的100多个原子的体系,S1的opt优化能好使,但是freq的优化能算动吗?谢谢!

freq如果用DFT算不懂,sob老师说过,现在高斯还不支持二阶解析导数,不过可以用CIS算算试试吧,你算激发态频率作甚?qq已添加,不过添加不上,需要一个验证。我的:502291243
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2016-3-9 23:40:22 | 只看该作者 Only view this author
小范范1989 发表于 2016-3-9 21:54
你好,你说:据说DFT计算得到的DFT好像更加接近真实的T1、也就是03的方式更加的接近实验数值?
太好了, ...

TADF的100多个原子的体系,S1的opt优化能好使,但是freq的优化能算动吗?谢谢!

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

17#
 楼主 Author| 发表于 Post on 2016-3-9 23:15:44 | 只看该作者 Only view this author
小范范1989 发表于 2016-3-9 21:54
你好,你说:据说DFT计算得到的DFT好像更加接近真实的T1、也就是03的方式更加的接近实验数值?
太好了, ...

兄台过奖了,我刚刚接触高斯猜一个多月,更多东西还要向兄台请教,加个QQ459604470,还望不吝指教!
我们这边坐材料和器件!

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

16#
发表于 Post on 2016-3-9 21:58:58 | 只看该作者 Only view this author
yezhonghua 发表于 2016-3-9 19:58
采用DFT和TD-DFT计算T1还是有很大差距的;DFT算出来为2.69ev;TD-DFT算出来为2.48ev。这对于TADF材料来说 ...

原理上DFT方法更准一些,因为会考虑到T1时轨道的弛豫,而不像TDDFT那样是基于S0轨道来算的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

1294

帖子

0

威望

6236

eV
积分
7530

Level 6 (一方通行)

15#
发表于 Post on 2016-3-9 21:54:07 | 只看该作者 Only view this author
yezhonghua 发表于 2016-3-9 19:58
采用DFT和TD-DFT计算T1还是有很大差距的;DFT算出来为2.69ev;TD-DFT算出来为2.48ev。这对于TADF材料来说 ...

你好,你说:据说DFT计算得到的DFT好像更加接近真实的T1、也就是03的方式更加的接近实验数值?
太好了,你也是TADF方向的,你们主要是做计算,理论还是实验?
我之前计算能差,采用的是adachi教授,还有张老师的方法,optimal HF method。
后来,我也优化过激发态,找能查,我就是那个符合的好,就用那个。
望以后多多交流,共同进步,我懂得比较少,还要多和你学习。
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2016-3-9 19:58:08 | 只看该作者 Only view this author
小范范1989 发表于 2016-3-4 09:53
你是要找T1的能级?类似于文献中,比如说T1的能量为3.2eV?这样子吗?
我感觉,
1:可以优化基态,找到HF ...

采用DFT和TD-DFT计算T1还是有很大差距的;DFT算出来为2.69ev;TD-DFT算出来为2.48ev。这对于TADF材料来说计算△Est影响很大。不知道兄台采用哪种方法?
据说DFT计算得到的DFT好像更加接近真实的T1。
谢谢!

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2016-3-9 16:51:44 | 只看该作者 Only view this author
小范范1989 发表于 2016-3-9 07:09
偶,谢谢你的指点,我是做理论的,所以都是想当然。
我咨询你点问题:你刚刚说的这个triplet的垂直激发 ...

你的这个问题我也咨询过老师,S1的激发态频率计算好像时间很长,所以我都是用的基态做的垂直激发能。

1294

帖子

0

威望

6236

eV
积分
7530

Level 6 (一方通行)

12#
发表于 Post on 2016-3-9 07:09:41 | 只看该作者 Only view this author
yezhonghua 发表于 2016-3-8 23:01
你算出来的应该是绝热能T1,但是一般实验上是垂直激发能和材料的peak 峰比较符合。
不知道是不是对的, ...

偶,谢谢你的指点,我是做理论的,所以都是想当然。
我咨询你点问题:你刚刚说的这个triplet的垂直激发能,我想问一下,磷光峰一般和td的这个计算数值对应?如果我想计算T1,T,T3的能级,来比较这个S1-T1gap,你说,我是通过分别优化S1,T1,T2,T3(也就是绝热能差)来讨论,还是基于基态结构做td来讨论?
非常感谢你的指点,谢谢
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2016-3-8 23:01:14 | 只看该作者 Only view this author
小范范1989 发表于 2016-3-4 09:53
你是要找T1的能级?类似于文献中,比如说T1的能量为3.2eV?这样子吗?
我感觉,
1:可以优化基态,找到HF ...

你算出来的应该是绝热能T1,但是一般实验上是垂直激发能和材料的peak 峰比较符合。
不知道是不是对的,多谢指教!

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1 没错

查看全部评分 View all ratings

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2016-3-8 22:54:48 | 只看该作者 Only view this author
sobereva 发表于 2016-3-4 21:42
没有报错不代表合理、计算没有问题!!!!!!!!!
nstates只能写一次,你又写10又写20,你怎么知 ...

多谢老师指正!我用您的方法计算一下,和我的比较一下,对照文献看看哪种比较合适!

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

9#
发表于 Post on 2016-3-4 21:42:56 | 只看该作者 Only view this author
yezhonghua 发表于 2016-3-4 15:44
有2个states,但是计算结果里面有triplet能级。没有问题,计算没有报错!


没有报错不代表合理、计算没有问题!!!!!!!!!
nstates只能写一次,你又写10又写20,你怎么知道程序到底会算10个态还是20个态?得程序算完了你才知道。如果你想算20个态,程序只给你算了10个态怎么办?又得重新算一次,白浪费了计算时间。如果发现程序算了20个态,那你写nstates=10又等于白写,增加了输入文件里无用的关键词的量(输入文件里无用的关键词越少说明用户的程序使用水平越高)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

430

帖子

0

威望

783

eV
积分
1213

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2016-3-4 15:44:58 | 只看该作者 Only view this author
sobereva 发表于 2016-3-3 21:10
图中显然不能写两次nstates,冲突了

是否需要做freq来验证是否真正得到了基态极小点取决于你自己,你自 ...

有2个states,但是计算结果里面有triplet能级。没有问题,计算没有报错!

1294

帖子

0

威望

6236

eV
积分
7530

Level 6 (一方通行)

7#
发表于 Post on 2016-3-4 10:38:57 | 只看该作者 Only view this author
sobereva 发表于 2016-3-4 10:38
最小点或极小点,以及其附近一定区域内都没有虚频
鞍点以及其附近一定区域内都有虚频,

恩,明白了,谢谢老师
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

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

GMT+8, 2024-11-24 22:00 , Processed in 0.369689 second(s), 25 queries , Gzip On.

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