| 查看: 4181 | 回复: 14 | ||||
Illusionist银虫 (正式写手)
|
[求助]
激发态的homo和lumo
|
|||
|
小弟要计算激发态的homo和lumo 看了论坛的帖子 http://muchong.com/bbs/viewthread.php?tid=1459478 但是还是有些不明白 1、基态构型优化：# b3lyp/6-31g(d) opt freq 2、垂直激发计算：# b3lyp/6-31g(d) td=(singlets,nstates=6,root=1)。坐标部分采用1优化后的坐标 3、激发态结构优化：# rcis(root=1,singlets) opt freq。坐标部分采用2的坐标 4、再一次TD的话就可以找出荧光谱和磷光谱了。 不明的地方： 我的做法是： 1、基态构型优化：# b3lyp/6-31g(d) opt freq 2.然后用这个优化# opt cis=singlets/6-311g geom=connectivity 3.# cis=direct/6-311+g(d,p) td(singlets,nstates=6,root=1) 算的不对 4.再用# cis=(nstates=6,root=1,direct)/6-311+g(d,p) geom=connectivity 也不行 5.我用# b3lyp/6-31+g(d) td(singlets,nstates=6,root=1)结果差不多 为什么一定要用 “2、垂直激发计算：# b3lyp/6-31g(d) td=(singlets,nstates=6,root=1)。坐标部分采用1优化后的坐标” 这个计算。 不明白为什么cis算得就是不行。是我的设置还是顺序什么的有问题么？ |
回复此楼» 猜你喜欢
能源环境新材料研究团队 导师介绍 研究生招生
已经有0人回复
-
能源环境新材料研究团队 导师介绍 研究生招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有61人回复
-
能源环境新材料研究团队 导师介绍 研究生招生
已经有0人回复
-
348 求调剂
已经有5人回复
-
招收 材料、新能源、冶金、化工、环境 等工科类考生
已经有2人回复
-
招收 新能源、冶金、材料、化工、环境 等工科类考生
已经有0人回复
-
393求调剂推荐
已经有1人回复
-
能源环境新材料研究团队 招收 新能源、冶金、材料、化工、环境 等工科类考生
已经有0人回复
-
能源环境新材料研究团队 招收 新能源、冶金、材料、化工、环境 等工科类研究生
已经有0人回复
-
武汉纺织大学教育部重点实验室调剂
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算周期性体系激发态
已经有3人回复
- phosphorescent OLED material design-------high triplet level
已经有10人回复
- “文章故事”开版第一篇：基于TDDFT设计有机分子，以期实现太阳能的全谱吸收
已经有26人回复
- 基态的HOMO-LUMO能隙与吸收光谱规律不一致！有这种情况吗？
已经有12人回复
- 如何判断 intramolecular charge transfer (ICT)
已经有5人回复
- 【求助】homo和lumo的能级差到底可以说明什么？
已经有12人回复
- 【求助】用TDDFT优化出来激发态之后的布居分析
已经有20人回复
- 【求助】HOMO -LUMO 基态-激发态
已经有4人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13022.2
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
2楼2011-11-08 23:56:22
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
10楼2013-09-28 08:46:40
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13022.2
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
【答案】应助回帖
Illusionist(金币+10): 谢谢 2011-11-10 08:21:49
|
第一步应该用dft。这一步确定基态平衡结构 第二步在基态平衡结构计算激发能，根据Franck－Condon原理，近似对应电子吸收谱的能量 第三步用tddft优化第一激发态，最后得到的激发能近似对应这个电子态的电子发射谱能量。这里的的freq是数值方法（因为gaussian的tddft目前还没实现解析二阶导数），相当耗时，如果能肯定是能量极小点（例如，根据文献，或者没有使用对称性而且也没对结构加上限制），没必要做freq计算。 HF和DFT的标度远比post-HF低得多，因此如果体系不是太大（例如，800基函数以下），可以从头至尾都用最大的基组做DFT／TDDFT计算。这样可以省去第四步。 |
7楼2011-11-09 22:24:53
喻儿在线
金虫 (小有名气)
很不愿意努力但必须努力的大喵
- 应助: 6 (幼儿园)
- 金币: 1559.2
- 散金: 59
- 帖子: 187
- 在线: 100.1小时
- 虫号: 1341868
- 注册: 2011-07-10
- 性别: MM
- 专业: 化学计量学与化学信息学
送红花一朵 |
楼主，我还是想问一下您，对于激发态我确实有太多的问题，您在 1. # opt b3lyp/6-31+(d,p) geom=connectivity 是对基态优化，其相应的轨道计算也就是： # b3lyp/6-31+(d,p) pop=full geom=connectivity 而在 3. # opt td=(nstates=6,root=1)/6-31+g(d,p) geom=check guess=read 是对激发态进行优化，那么其对应的轨道计算也就是： # td=(nstates=6,root=1)/6-31+g(d,p) pop=full geom=check guess=read 我这样计算是对的吗？还是说激发态的前线轨道还有别的方法？ |
11楼2013-10-07 22:26:12
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
12楼2013-10-08 10:24:23
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
14楼2013-10-08 15:47:36
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
3楼2011-11-09 08:30:47
daiyulan85
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2573.7
- 散金: 85
- 红花: 1
- 帖子: 401
- 在线: 349.7小时
- 虫号: 1078645
- 注册: 2010-08-20
- 性别: MM
- 专业: 理论和计算化学
4楼2011-11-09 09:07:03
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
|
第一步 优化是指优化 基态结构吧，用这个算法应该够了吧 # opt freq hf/6-311g geom=connectivity 第二步td 第二部td的意义是什么。是为了计算激发态，通过td得到一个较粗糙的激发态模型 那么是否该这么写 # td=(nstates=6,root=1)/6-31+g(d,p) geom=connectivity 第三步 进行更精确的模型 # opt freq td=(nstates=6,root=1) b3lyp/6-311g geom=connectivity 第四步 # td=(nstates=6,root=1) b3lyp/6-311+g(d,p) geom=connectivity |
5楼2011-11-09 11:31:58
daiyulan85
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 2573.7
- 散金: 85
- 红花: 1
- 帖子: 401
- 在线: 349.7小时
- 虫号: 1078645
- 注册: 2010-08-20
- 性别: MM
- 专业: 理论和计算化学
6楼2011-11-09 17:24:11
Illusionist
银虫 (正式写手)
- 应助: 20 (小学生)
- 金币: 232.8
- 散金: 1001
- 红花: 15
- 帖子: 431
- 在线: 329.1小时
- 虫号: 1123329
- 注册: 2010-10-15
- 专业: 半导体晶体与薄膜材料
|
第一步 确定基态平衡结构 # opt freq b3lyp/6-311g geom=connectivity 第二部 在基态平衡结构计算激发能，根据Franck－Condon原理，近似对应电子吸收谱的能量。 # td=(nstates=6,root=1) b3lyp/6-311++g(d,p) geom=connectivity 第三步 用tddft优化第一激发态，最后得到的激发能近似对应这个电子态的电子发射谱能量。这里的的freq是数值方法（因为gaussian的tddft目前还没实现解析二阶导数），相当耗时，如果能肯定是能量极小点（例如，根据文献，或者没有使用对称性而且也没对结构加上限制），没必要做freq计算。 # opt freq td=(nstates=6,root=1) b3lyp/6-311++g(d,p) geom=connectivity 问一下 第二步的结果对第三步的影响有多大？ 不使用扩散函数的结果会不会受影响。即可否使用 # td=(nstates=6,root=1) b3lyp/6-311g geom=connectivity 这样比较节约时间 非常谢谢 |
8楼2011-11-10 08:21:40
喻儿在线
金虫 (小有名气)
很不愿意努力但必须努力的大喵
- 应助: 6 (幼儿园)
- 金币: 1559.2
- 散金: 59
- 帖子: 187
- 在线: 100.1小时
- 虫号: 1341868
- 注册: 2011-07-10
- 性别: MM
- 专业: 化学计量学与化学信息学
9楼2013-09-26 11:25:04
