| 查看: 4727 | 回复: 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算得就是不行。是我的设置还是顺序什么的有问题么？ |
回复此楼» 猜你喜欢
Journal of Alloys and Compounds投稿过程
已经有2人回复
-
澳大利亚皇家墨尔本理工大学（RMIT）理学院招收2027博士生
已经有62人回复
-
物理化学论文润色/翻译怎么收费?
已经有264人回复
-
为什么现在电催化越来越卷？从“只看效果”到“全周期低成本”选型转变
已经有0人回复
-
电催化氧化反应器现场运维，不可忽视的安全风险与防控措施
已经有0人回复
-
电池测试中dQ/dV怎么理解
已经有1人回复
-
精细化工浓水做零排放，前端到底该配什么氧化工艺？
已经有0人回复
-
大湾区大学姚瑶课题组2027届博士生招生 （与中山大学联合培养）
已经有24人回复
-
氯离子浓度对电催化处理效果的双重影响实测与参数优化
已经有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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13423
- 散金: 27307
- 红花: 302
- 帖子: 4688
- 在线: 929.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
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +459 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13423
- 散金: 27307
- 红花: 302
- 帖子: 4688
- 在线: 929.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
