计算化学公社

标题: 气相下计算T1能级 [打印本页]

作者
Author: yezhonghua 时间: 2016-3-3 16:21
标题: 气相下计算T1能级
  采用TD-DFT计算T1的能级，方法是这样的不知道对不对：
  1、对初始结构进行freq+opt优化，电荷、多重度0,1
  2、优化完的基态进行TD计算（图1）所示，log文件找到triplet 1。

以下是sober老师的方法（TD-DFT）：
1 优化T1结构：# B3LYP/6-31G* opt，电荷和自旋多重度为0 3。
2 取上一步优化得到的结构，用# B3LYP/6-31G* TD(triplet)，电荷和自旋多重度为0 1。第一激发态的激发能就是磷光发射能了。

这里的第一步直接用opt就可以么?不需要用freq么？另外我的计算（第一步）是否正确？
谢谢！

作者
Author: sobereva 时间: 2016-3-3 21:10
图中显然不能写两次nstates，冲突了

是否需要做freq来验证是否真正得到了基态极小点取决于你自己，你自己从优化完的结构上看有信心觉得这就是极小点就不用做freq验证是否有虚频，拿不准就做freq

作者
Author: yezhonghua 时间: 2016-3-4 09:32

sobereva 发表于 2016-3-3 21:10
图中显然不能写两次nstates，冲突了

是否需要做freq来验证是否真正得到了基态极小点取决于你自己，你自 ...

Soberl老师，按照您的建议,是否可以认为：
1.基态极小值点的确定通过freq优化，结果无虚拼的话，可以确定是基态极小值点。
2.我算T1的第一步时，是否是先进行电荷，多重度0，3的freq+opt优化；第二步，在优化完的基础上进行TD=（Triplet）电荷，多重度0，1的计算，log文件的excited state 1就算T1了？第一步，电荷，多重度0，3
第二步，电荷，多重度0，1
谢谢老师！

作者
Author: 小范范1989 时间: 2016-3-4 09:53
你是要找T1的能级？类似于文献中，比如说T1的能量为3.2eV？这样子吗？
我感觉，
1：可以优化基态，找到HF。
2：然后通过DFT （0,3）或者TDDFT（0，1）的方式优化T1，找到HF能量。当然，这两个方法算出来的HF可能会有一定的差距。
3：两者相减，就是这个绝热的T1的能量。
4：实验上很多时候，我看好像是得到了S1的能级，然后知道了S1-T1的gap，然后做差找到的T1能级。
以上是我的理解，不一定对哈，还望高手指点。

ps：我咨询一个问题：这个最小点，是不是没有虚频？极小点是不是也没有虚频？鞍点才有虚频？谢谢指点。

作者
Author: sobereva 时间: 2016-3-4 10:36

yezhonghua 发表于 2016-3-4 09:32
Soberl老师，按照您的建议,是否可以认为：
1.基态极小值点的确定通过freq优化，结果无虚拼的话，可以确 ...

1 是
2 是的

作者
Author: sobereva 时间: 2016-3-4 10:38

小范范1989 发表于 2016-3-4 09:53
你是要找T1的能级？类似于文献中，比如说T1的能量为3.2eV？这样子吗？
我感觉，
1：可以优化基态，找到HF ...

最小点或极小点，以及其附近一定区域内都没有虚频
鞍点以及其附近一定区域内都有虚频

作者
Author: 小范范1989 时间: 2016-3-4 10:38

sobereva 发表于 2016-3-4 10:38
最小点或极小点，以及其附近一定区域内都没有虚频
鞍点以及其附近一定区域内都有虚频，

恩，明白了，谢谢老师

作者
Author: yezhonghua 时间: 2016-3-4 15:44

sobereva 发表于 2016-3-3 21:10
图中显然不能写两次nstates，冲突了

是否需要做freq来验证是否真正得到了基态极小点取决于你自己，你自 ...

有2个states，但是计算结果里面有triplet能级。没有问题，计算没有报错！

作者
Author: sobereva 时间: 2016-3-4 21:42

yezhonghua 发表于 2016-3-4 15:44
有2个states，但是计算结果里面有triplet能级。没有问题，计算没有报错！

没有报错不代表合理、计算没有问题！！！！！！！！！
nstates只能写一次，你又写10又写20，你怎么知道程序到底会算10个态还是20个态？得程序算完了你才知道。如果你想算20个态，程序只给你算了10个态怎么办？又得重新算一次，白浪费了计算时间。如果发现程序算了20个态，那你写nstates=10又等于白写，增加了输入文件里无用的关键词的量（输入文件里无用的关键词越少说明用户的程序使用水平越高）

作者
Author: yezhonghua 时间: 2016-3-8 22:54

sobereva 发表于 2016-3-4 21:42
没有报错不代表合理、计算没有问题！！！！！！！！！
nstates只能写一次，你又写10又写20，你怎么知 ...

多谢老师指正！我用您的方法计算一下，和我的比较一下，对照文献看看哪种比较合适！

作者
Author: yezhonghua 时间: 2016-3-8 23:01

小范范1989 发表于 2016-3-4 09:53
你是要找T1的能级？类似于文献中，比如说T1的能量为3.2eV？这样子吗？
我感觉，
1：可以优化基态，找到HF ...

你算出来的应该是绝热能T1，但是一般实验上是垂直激发能和材料的peak 峰比较符合。
不知道是不是对的，多谢指教！

作者
Author: 小范范1989 时间: 2016-3-9 07:09

yezhonghua 发表于 2016-3-8 23:01
你算出来的应该是绝热能T1，但是一般实验上是垂直激发能和材料的peak 峰比较符合。
不知道是不是对的， ...

偶，谢谢你的指点，我是做理论的，所以都是想当然。
我咨询你点问题：你刚刚说的这个triplet的垂直激发能，我想问一下，磷光峰一般和td的这个计算数值对应？如果我想计算T1，T，T3的能级，来比较这个S1-T1gap，你说，我是通过分别优化S1，T1，T2，T3（也就是绝热能差）来讨论，还是基于基态结构做td来讨论？
非常感谢你的指点，谢谢

作者
Author: yezhonghua 时间: 2016-3-9 16:51

小范范1989 发表于 2016-3-9 07:09
偶，谢谢你的指点，我是做理论的，所以都是想当然。
我咨询你点问题：你刚刚说的这个triplet的垂直激发 ...

你的这个问题我也咨询过老师，S1的激发态频率计算好像时间很长，所以我都是用的基态做的垂直激发能。

作者
Author: yezhonghua 时间: 2016-3-9 19:58

小范范1989 发表于 2016-3-4 09:53
你是要找T1的能级？类似于文献中，比如说T1的能量为3.2eV？这样子吗？
我感觉，
1：可以优化基态，找到HF ...

采用DFT和TD-DFT计算T1还是有很大差距的；DFT算出来为2.69ev；TD-DFT算出来为2.48ev。这对于TADF材料来说计算△Est影响很大。不知道兄台采用哪种方法？
据说DFT计算得到的DFT好像更加接近真实的T1。
谢谢！

作者
Author: 小范范1989 时间: 2016-3-9 21:54

yezhonghua 发表于 2016-3-9 19:58
采用DFT和TD-DFT计算T1还是有很大差距的；DFT算出来为2.69ev；TD-DFT算出来为2.48ev。这对于TADF材料来说 ...

你好，你说：据说DFT计算得到的DFT好像更加接近真实的T1、也就是03的方式更加的接近实验数值？
太好了，你也是TADF方向的，你们主要是做计算，理论还是实验？
我之前计算能差，采用的是adachi教授，还有张老师的方法，optimal HF method。
后来，我也优化过激发态，找能查，我就是那个符合的好，就用那个。
望以后多多交流，共同进步，我懂得比较少，还要多和你学习。

作者
Author: sobereva 时间: 2016-3-9 21:58

yezhonghua 发表于 2016-3-9 19:58
采用DFT和TD-DFT计算T1还是有很大差距的；DFT算出来为2.69ev；TD-DFT算出来为2.48ev。这对于TADF材料来说 ...

原理上DFT方法更准一些，因为会考虑到T1时轨道的弛豫，而不像TDDFT那样是基于S0轨道来算的

作者
Author: yezhonghua 时间: 2016-3-9 23:15

小范范1989 发表于 2016-3-9 21:54
你好，你说：据说DFT计算得到的DFT好像更加接近真实的T1、也就是03的方式更加的接近实验数值？
太好了， ...

兄台过奖了，我刚刚接触高斯猜一个多月，更多东西还要向兄台请教，加个QQ459604470，还望不吝指教！
我们这边坐材料和器件!

作者
Author: yezhonghua 时间: 2016-3-9 23:40

小范范1989 发表于 2016-3-9 21:54
你好，你说：据说DFT计算得到的DFT好像更加接近真实的T1、也就是03的方式更加的接近实验数值？
太好了， ...

TADF的100多个原子的体系，S1的opt优化能好使，但是freq的优化能算动吗？谢谢！

作者
Author: 小范范1989 时间: 2016-3-10 07:13
本帖最后由小范范1989 于 2016-3-10 07:15 编辑

yezhonghua 发表于 2016-3-9 23:40
TADF的100多个原子的体系，S1的opt优化能好使，但是freq的优化能算动吗？谢谢！

freq如果用DFT算不懂，sob老师说过，现在高斯还不支持二阶解析导数，不过可以用CIS算算试试吧，你算激发态频率作甚？qq已添加，不过添加不上，需要一个验证。我的：502291243

作者
Author: 小范范1989 时间: 2016-3-10 07:45
我咨询你一个问题，这个实验上说的这个dual fluorescence是什么意思？谢谢指点。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)