计算化学公社

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

[Gaussian/gview] 求助三重态T1到Tn的吸收怎么计算-高斯

[复制链接 Copy URL]

35

帖子

0

威望

143

eV
积分
178

Level 3 能力者

本帖最后由 xinren 于 2019-9-5 20:36 编辑

最近在做三重态T1到Tn三重态之间的吸收,不知道怎么设输入文件,之前想的是T1作为三重态的基态Tn是三重态的激发态,输入文件是
第一步:优化T1:opt freq 0 3
第二步:td=(nstates=6) 0 3
输出文件:
Excited State   1:  3.020-A      1.1960 eV 1036.70 nm  f=0.0000  <S**2>=2.029
     38B -> 39B        0.98422
     38B -> 40B        0.11768
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -498.243571305   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:  3.019-A      1.6811 eV  737.54 nm  f=0.0021  <S**2>=2.029
     40A -> 41A        0.46357
     40A -> 43A       -0.15651
     37B -> 39B        0.86550
不知道文件输入对不对,请各位指导一下!另外还需要看T1到Tn的跃迁轨道,不知道具体操作是什么?

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

2#
发表于 Post on 2019-9-6 00:31:32 | 只看该作者 Only view this author
第二步关键词明显是错的。你当前的基态(作为参考态)是单重态,显然第二步必须写0 1

Gaussian只能给出激发态之间的能量差,但是激发态之间的跃迁偶极矩Gaussian没法直接给出来,如果要考虑这点,需要用Multiwfn
使用Multiwfn计算激发态间的跃迁偶极矩和各个激发态的偶极矩
http://sobereva.com/227

没有所谓的激发态之间的跃迁轨道一说,每个激发态都是基于相同的参考态的分子轨道来描述的
北京科音自然科学研究中心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

35

帖子

0

威望

143

eV
积分
178

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2019-9-6 15:34:44 | 只看该作者 Only view this author
好的谢谢老师,那是不是把第二步的0 3 改成0 1求得的能量差就是Tn与T1之间的能差吗?之前看到求磷光的方法是
1 优化T1结构:# B3LYP/6-31G* opt,电荷和自旋多重度为0 3。
2 取上一步优化得到的结构,用# B3LYP/6-31G* TD(triplet),电荷和自旋多重度为0 1。第一激发态的激发能就是磷光发射能了。
请问如果改成0 1,那求出的能差不是就表示第一激发态是T1与S0的能差,第二激发态就是表示T2与S0的能差吗?另外是用singlet还是用triplet呢?
另外用Multiwfn求Tn与T1的跃迁矩的输入文件是B3LYP/6-31+G(d) TD(triplet,nstates=10) IOp(9/40=4) 0 1吗?

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

4#
发表于 Post on 2019-9-7 03:05:45 | 只看该作者 Only view this author
xinren 发表于 2019-9-6 15:34
好的谢谢老师,那是不是把第二步的0 3 改成0 1求得的能量差就是Tn与T1之间的能差吗?之前看到求磷光的方法 ...

每个结构下都有各个Tn与T1的能量差。显然你把当前结构下算出来的三重态激发能之间求差就能得到当前结构下的Tn与T1能量差

你要算三重态激发态之间的能量差,显然得写triplet,否则算出的激发态都是单重态的有什么用?

北京科音自然科学研究中心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

35

帖子

0

威望

143

eV
积分
178

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2019-9-22 17:14:43 | 只看该作者 Only view this author
好的,谢谢老师!之前Tn与T1的能量差算出来了,想要再问一下老师,相应的Tn与T1的分子轨道该怎么计算呢?

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

6#
发表于 Post on 2019-9-23 08:16:58 | 只看该作者 Only view this author
xinren 发表于 2019-9-22 17:14
好的,谢谢老师!之前Tn与T1的能量差算出来了,想要再问一下老师,相应的Tn与T1的分子轨道该怎么计算呢?

对于TDDFT,根本没法说什么T1、Tn的轨道。顶多是说在相应激发态优化的极小点结构下算的基态轨道

北京科音自然科学研究中心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

35

帖子

0

威望

143

eV
积分
178

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2019-9-30 10:31:36 | 只看该作者 Only view this author
好的,谢谢!!!!!!

35

帖子

0

威望

143

eV
积分
178

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2019-10-2 11:24:22 | 只看该作者 Only view this author
老师,实验测得的三重态吸收峰有两个,一个大约是在400nm处,一个大约是在450nm处,但是高斯计算出来的振子强度比较强的两个峰值都比较集中在400nm区域,分别是390nm和410nm左右,请问会不会是计算方法引起的误差?

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

9#
发表于 Post on 2019-10-3 04:28:43 | 只看该作者 Only view this author
xinren 发表于 2019-10-2 11:24
老师,实验测得的三重态吸收峰有两个,一个大约是在400nm处,一个大约是在450nm处,但是高斯计算出来的振子 ...

可能
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 06:32 , Processed in 0.154839 second(s), 23 queries , Gzip On.

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