| 查看: 2219 | 回复: 3 | ||
| 【悬赏金币】回答本帖问题，作者花放迟将赠送您 60 个金币 | ||
[求助]
计算S1,T1以及Tn的计算过程，请求大佬看一下计算过程给予指导 已有1人参与
|
||
|
各位老师您好！ 我求助计算S1,T1以及Tn的计算过程，求各位i老师看一下我的计算过程，并给与我一些指导 首先我先进行了结构优化opt，指令为# b3lyp/6-31g(d,p) opt freq? ?opt??0 1得到HOMO和LUMO 然后取以上chk再优化激发态td，指令为# cam-b3lyp/6-311+g(d,p) td=(nstates=100) guess=read geom=check em=gd3bj opt??0 1 然后取以上chk再计算S1态能量，指令为# b3lyp/6-311+g(d,p) TD(nstates=3,root=1) guess=read geom=check opt=(maxcycle=30) freq em=gd3bj??opt??0 1得到S1能量 然后取opt优化后的chk计算T1态能量，指令为# cam-b3lyp/6-311+g(d,p) opt=(maxcycle=20) freq geom=check opt??0 3得到T1能量 然后取以上T1最后的chk后计算Tn，指令为# cam-b3lyp/6-311+g(d,p) TD(nstates=20,50-50) freq geom=check opt 0 1得到Tn能量 这样计算过程的到的计算结果可靠不？ 求各位老师指导，谢谢 @jiewei 发自小木虫手机客户端 |
回复此楼
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26428.6
- 红花: 349
- 帖子: 2026
- 在线: 1446.2小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
感谢参与，应助指数 +1
|
首先你还是没说清楚你算这些到底是要算什么。如果是计算光谱数据。 那么 整体来看 你要做的就是 1. 基态优化使用B3LYP/6-31G(d) 2. S1激发态优化使用 TD cam-B3LYP/6-311+G(d,p) 3. 优化T1结构 使用 0 3 cam-B3LYP 做普通的DFT优化。 4 用T1结构 做 TD-CAM-B3LYP做Tn 各个态的激发能。 整体流程上 你的第二步冗余了，可以不要。 而且nstate=100也没有什么实际意义 白白浪费时间。 第三步计算S1的能量使用B3LYP也跟其他步骤不一致。且还要单独计算freq比较怪。因为激发能指的都是电子能。单独算freq比较奇怪。建议S1的优化 就是用TD-cam-B3LYP的方式直接算opt freq得到激发能就可以了。 如果是要算绝热的自由能差（比如使用Marcus公式衍生出来的各种变体算各种速率）。那么保证所有T1 S1 S0 的结构优化使用的opt freq是同一个计算级别。（热力学校正值使用的方法一致）T1 S1 S0 的高精度电子能方式使用同一个方法比如 S1 T1 都用 TD-cam-b3lyp/6-311+G(d,p) S0用cam-B3LYP/6-311+G(d,p) |
» 本帖已获得的红花（最新10朵）
2楼2024-12-12 12:49:15
送红花一朵 |
老师您好，感谢您的细心回复，我说明一下我要计算的是荧光磷光的发射能，还有S1以及T1到Tn的能级，判断系间窜越是从那个Tn到的S1。老师您说第二步多余了，那我计算了是算错嘛，如果没错误是不是也可以，就是费时了点。另外，第三步老师你最后说就是加上Cam，然后去掉freq嘛，那后面算T的时候需要freq？ 发自小木虫手机客户端 |
3楼2024-12-14 09:26:01
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26428.6
- 红花: 349
- 帖子: 2026
- 在线: 1446.2小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
4楼2024-12-14 10:23:35
