|
|
你好,我也是自己看了一点资料,了解的也不是很深。说一下我的理解。
1:这个问题,就是计算的问题了。关于荧光和磷光一些概念还有计算,可以先看看 sob老师的这个帖子,非常非常的好 http://sobereva.com/314
这里如果是要计算垂直能差,那就是基于S0的结构,然后做td(50-50)、至于这个50-50就是singlet和triplet一起计算,出来如下的形式:
Excited State 1: Triplet-A 3.1058 eV 399.20 nm f=0.0000 <S**2>=2.000
121 ->131 0.25084
121 ->142 0.11098
123 ->131 -0.17376
124 ->131 0.47164
128 ->132 -0.17700
129 ->131 -0.27270
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) = -1837.69264810
Copying the excited state density for this state as the 1-particle RhoCI density.
Excited State 2: Triplet-A 3.2394 eV 382.74 nm f=0.0000 <S**2>=2.000
130 ->131 0.65932
130 ->133 -0.11968
130 ->134 -0.17054
Excited State 3: Singlet-A 3.2637 eV 379.89 nm f=0.0063 <S**2>=0.000
130 ->131 0.66857
130 ->133 -0.11552
130 ->134 -0.16898
所以这个S1-T1的vertical gap就是:0.1579eV。
2:关于垂直能查,就是上面的这种方式计算。然后绝热能查的话,就是分别优化S1和T1。sob老师帖子中应该有。关于能量如何读取,
1:你可以把优化之后的log文件,直接放到高斯view中。查看summary里面有,我记得。
2:或者搜寻exicted state
倒着搜索,就是最后一次。就是这个能量。这个能量=最后的HF能量+td能量。你可以手动算一下。
Excited State 1: Singlet-A 2.7310 eV 453.99 nm f=0.02 <S**2>=0.000
130 ->131 0.68673
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) = -1837.69717112
Copying the excited state density for this state as the 1-particle RhoCI density.
#p BMK/6-31g* opt td 这是我优化S1的表头。
关于T1,表头这样写#p BMK/6-31g* opt td(triplet,root=1) PS:我一般是采用01的方式哈。sob老师中帖子也说了,有 03 的方式,你看看。我采用01主要是为了和S1对应起来。
计算结果如下:
Excited State 1: Triplet-A 2.8122 eV 440.88 nm f=0.0000 <S**2>=2.000
122 ->131 0.22005
123 ->131 0.24389
124 ->131 0.10107
130 ->131 0.56246
130 ->133 -0.14508
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) = -1837.70029912
所以这个S1-T1的绝热能查为:(-0.69717112+0.70029912)*27.2116=0.085eV
3:最近我也看到了jacs那个文献,说的是严格的一样。最近也打算看看量化里面的公式推导。之前我问过一个帖子http://bbs.keinsci.com/forum.php?mod=viewthread&tid=239&extra=
这里一般这个T1是比较乱的,所以很多时候采用NTO的方式。
同时你看他的按个jacs也说了,不光与这个能差有关,还有SOC有关。
我感觉,小的能查是实现TADF的一个条件,小能查有利于实现TADF,并不是所有小gap就有TADF,这不是还有SOC吗。但是有的理论计算文章,就是通过能差衡量。比如有的计算文章,计算了实验的分子,准确计算了能查,然后基于实验分子,设计分子,算出来的能查比实验还小。(我之前几篇小文章就这个思路,档次不高)。
4:第四个问题,我第二个已经说了,应该能理解吧。
那就说一下这个调节w的lc-wpbe方法。孙老师和钟老师都在sob老师的群里,去北京参加培训的时候还碰到了两位老师。我当时咨询过孙老师这个方法。
你仔细看看文献,他没用这个方法来优化激发态,lc-wpbe的方法不支持激发态的优化。里面所有调节w的。都是单点计算。
我记得他激发态是通过cam-b3lyp优化的,要是没记错的话。
然后用lc-wpbe的方法来计算单点。
这种lc-wpbe的方法,以及adachi和张其胜老师的OHF方法,以及黄维老师的NTO方法,都是基于基态计算的。这个在论坛有讨论的帖子。
lc-blyp的方法,调节w,这个支持激发态的优化。顺带我当时组会的ppt。
个人见解,不一定准确
|
-
组会-忠.png
(387.71 KB, 下载次数 Times of downloads: 54)
-
组会-1.png
(445.28 KB, 下载次数 Times of downloads: 49)
|