求助Gaussian09计算HOMO-LUMO 和 三重态,单重态 △E ST
各位虫友，本人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 ]
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
单体计算的话，你自己的机器还跑得动。
聚集体的话，理论来说我觉得也可以跑哈哈，但这七十个单体聚合，自己的电脑得跑爆炸。
我想尝试自己跑一下，但是对于高斯还不会用，所以想请教一下大家帮忙指点一下
各位大神帮帮忙
首先高斯一般情况下（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的电子能。二者相减就是 绝热能量差。
对于二聚体 三聚体 四聚体 计算方法基本一样。只是聚在一起可能会有不同构型。只要顺利找到优化出最稳定的构型。剩下的步骤就都一样，
非常有帮助，谢谢！另外，我想问一个问题，就是我的二聚体、三聚体、四聚体都是非共价相互作用聚集到一起的，我想知道的是怎么在高斯里导入这种非共价相互作用结合的分子，并且可以拿来计算优化？或者说是怎么在高斯里模拟这种结构进行计算优化？
高斯里也不区分共价还是非共价（相互作用比较弱就把各个单体相互摆的远一点）。你把初始结构设定好 算出来该是什么样子就就是什么样子。算出来之后 想知道它们之间是什么作用，可以进行进一步的波函数分析。但是算之前不需要考虑它们之间是什么相互作用。