| 查看: 2561 | 回复: 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 ] |
回复此楼» 猜你喜欢
PbS量子点如何确定浓度
已经有11人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有86人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）
已经有22人回复
-
PbS量子点紫外吸收
已经有0人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
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
- 应助: 587 (博士)
- 贵宾: 0.01
- 金币: 26383.1
- 红花: 349
- 帖子: 2025
- 在线: 1442小时
- 虫号: 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
- 应助: 587 (博士)
- 贵宾: 0.01
- 金币: 26383.1
- 红花: 349
- 帖子: 2025
- 在线: 1442小时
- 虫号: 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
