| 查看: 3076 | 回复: 8 | |||
zzyuss铁杆木虫 (知名作家)
抗旱救灾委员会名誉会长
|
[求助]
求助Gaussian09计算HOMO-LUMO 和 三重态,单重态 △E ST 已有1人参与
|
|
各位虫友，本人gaussian计算小白，最近尝试用gaussian09计算一个化合物单体的homo-lumo，单线态能级分布，三线态能级分布，以及系间窜越的能差Δes╟t值。在此基础上继续进行这个化合物单体聚集形成二聚体、三聚体以及四聚体的homo-lumo，单线态能级分布，三线态能级分布，以及系间窜越的能差Δes╟t值。化合物单体含有一个碘代季铵盐，请问大家该如何进行计算？大概的步骤是什么？有没有大神能帮个忙啊？有偿的！ 另外，我的化合物单体能够形成70个有效相干长度的聚集体，请问各位能不能计算70个聚集体的homo-lumo，单线态能级分布，三线态能级分布，以及系间窜越的能差Δes╟t值呢？ [ Last edited by paramecium86 on 2020-8-10 at 16:26 ] |
回复此楼» 猜你喜欢
为什么现在电催化越来越卷？从“只看效果”到“全周期低成本”选型转变
已经有0人回复
-
电催化氧化反应器现场运维，不可忽视的安全风险与防控措施
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有139人回复
-
电池测试中dQ/dV怎么理解
已经有1人回复
-
精细化工浓水做零排放，前端到底该配什么氧化工艺？
已经有0人回复
-
大湾区大学姚瑶课题组2027届博士生招生 （与中山大学联合培养）
已经有26人回复
-
北京工业大学2026年公开招聘公告
已经有0人回复
-
氯离子浓度对电催化处理效果的双重影响实测与参数优化
已经有0人回复
-
阴极协同电催化技术处理高盐难降解化工废水的机理与工程应用研究
已经有0人回复
-
氧化铝浆液
已经有2人回复
zzyuss2楼2020-08-10 18:49:00
3楼2020-08-10 18:49:51
zzyuss
铁杆木虫 (知名作家)
抗旱救灾委员会名誉会长
- 应助: 124 (高中生)
- 金币: 5451.8
- 散金: 1771
- 红花: 27
- 沙发: 67
- 帖子: 5818
- 在线: 576.5小时
- 虫号: 3191629
- 注册: 2014-05-08
- 专业: 无神论
4楼2020-08-10 21:27:05
zzyuss
铁杆木虫 (知名作家)
抗旱救灾委员会名誉会长
- 应助: 124 (高中生)
- 金币: 5451.8
- 散金: 1771
- 红花: 27
- 沙发: 67
- 帖子: 5818
- 在线: 576.5小时
- 虫号: 3191629
- 注册: 2014-05-08
- 专业: 无神论
5楼2020-08-10 22:31:38
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26763.8
- 红花: 349
- 帖子: 2031
- 在线: 1465.1小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
zhou2009: 金币+10 2020-08-11 15:03:27
zzyuss: 金币+1000, ★★★★★最佳答案, 非常有帮助，谢谢！另外，我想问一个问题，就是我的二聚体、三聚体、四聚体都是非共价相互作用聚集到一起的，我想知道的是怎么在高斯里导入这种非共价相互作用结合的分子，并且可以拿来计算优化？ 2020-08-11 15:13:20
感谢参与，应助指数 +1
zhou2009: 金币+10 2020-08-11 15:03:27
zzyuss: 金币+1000, ★★★★★最佳答案, 非常有帮助，谢谢！另外，我想问一个问题，就是我的二聚体、三聚体、四聚体都是非共价相互作用聚集到一起的，我想知道的是怎么在高斯里导入这种非共价相互作用结合的分子，并且可以拿来计算优化？ 2020-08-11 15:13:20
|
首先高斯一般情况下（20-40核心左右的服务器，linux版的高斯）用普通dft搭配最最保底的基组对于基态结构的优化，最多可以算300-400个原子数目的体系。所以看你体系的原子数目可以大概判断出能不能算得了。而且因为可能要涉及到了激发态的优化，这样的话200个原子以上就基本算不动了。所以70个分子聚在一起是肯定算不动的。 单体如果比较大可能连四聚体都不那么好算。毕竟聚在一起多了可能还有不同的构型。 以下是计算的简单过程，我默认楼主知道如何用高斯进行，几何优化计算， 并且知道dft， tddft都是什么。知道激发态是什么。如果这些都不知道的话，最好先去入门一下。（具体如何入门该看什么资料。见我这个教程前三页https://muchong.com/t-14242735-1） 就单体而言，计算homo-lumo轨道其实就是对这个分子的基态进行结构优化之后，就可以查看所有分子轨道的信息了，包括各个轨道的能量或者要作图也可以看分子轨道的等值面图。或者可能还涉及到homo-lumo轨道重合度分析。就都是这一步就可以完成。 至于系间窜越就涉及到了激发态计算。现在比较流行的方法就是tddft了。（简单的介绍可看https://muchong.com/t-14242735-1教程第1.5.1节和1.5.4 节推荐的几篇博文。） 而涉及到系间窜越的速率k isc 一般用Marcus公式来计算窜越速率什么的。这个公式里的 △E ST 一般要用到垂直能量差（如果是TADF之类的要计算反系间窜越 k risc 有时候要用到绝热能量差） 垂直能差计算比较容易 大体步骤就是 优化好分子的基态S0结构之后。用优化好的结构计算tddft(50-50) 关键词来计算基态结构下的各个单重和三重激发态的能量。把S1 T1能量作差就是 垂直 △E ST 。在这个计算中因为无论哪个激发态，都是基于基态的结构，它们结构相同，所以叫做垂直能量差。 而绝热能差，既然叫做绝热，就需要 对S1 和 T1 分别都进行几何优化（显然S1 和 T1的最稳定结构肯定不一样）。这样S1 和T1 各自都处于不同结构下，这样计算出来的能量差就叫绝热能量差。 对S1 的优化可以使用 tddft opt 关键词进行优化。计算完成后可以得到s1的电子能。（用tddft计算后得到的激发态的电子能如何寻找请看我这个教程https://muchong.com/t-14242735-1最后一页） 对T1的优化同样可以使用tddft(triplet) opt来进行优化。计算完成后也可以得到t1的电子能。二者相减就是 绝热能量差。 对于二聚体 三聚体 四聚体 计算方法基本一样。只是聚在一起可能会有不同构型。只要顺利找到优化出最稳定的构型。剩下的步骤就都一样。 |
6楼2020-08-11 00:22:45
zzyuss
铁杆木虫 (知名作家)
抗旱救灾委员会名誉会长
- 应助: 124 (高中生)
- 金币: 5451.8
- 散金: 1771
- 红花: 27
- 沙发: 67
- 帖子: 5818
- 在线: 576.5小时
- 虫号: 3191629
- 注册: 2014-05-08
- 专业: 无神论
7楼2020-08-11 15:14:39
paramecium86
版主 (著名写手)
- QC强帖: 5
- 应助: 588 (博士)
- 贵宾: 0.01
- 金币: 26763.8
- 红花: 349
- 帖子: 2031
- 在线: 1465.1小时
- 虫号: 2466443
- 注册: 2013-05-16
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
8楼2020-08-11 15:22:31
zzyuss
铁杆木虫 (知名作家)
抗旱救灾委员会名誉会长
- 应助: 124 (高中生)
- 金币: 5451.8
- 散金: 1771
- 红花: 27
- 沙发: 67
- 帖子: 5818
- 在线: 576.5小时
- 虫号: 3191629
- 注册: 2014-05-08
- 专业: 无神论
9楼2020-08-11 15:50:32
