|
|
最近有设计一些材料要初步判断一下能不能做TADF，读了论坛上各位老师的帖子之后，自己开始算了几个文件，把计算过程写在下面，也有几个疑问，希望能有老师指导一下~谢谢！
计算过程：
# opt freq td=singlets wb97xd/6-311g(d,p)
0 1
计算S1的能量，从输出文件中提取最后一轮的第一激发态的E(TD-HF/TD-KS)能量，记为Es1
然后，把分子的初始构型稍微调整下使对称性破坏，再计算：
# opt freq uwb97xd/6-311g(d,p)
0 3
提取最后的HF的能量，记为T1的能量，Et1
然后这个作为绝热的delta ST=Es1-Et1
问题：
1.上面的计算过程有没有什么问题？计算的是绝热的delta ST吗？另外，准确度是不是很低？
2.看过帖子和文献后觉得应该用校正w的泛函来求delta ST，因为分子太多并且每个的原子数也不少，每计算一个分子就要调节一次w很耗时间，就想先用wb97xd先计算一下看看，然后用LC-wPBE校正后再计算。
用校正的计算过程和上面一样，就是添加了两个IOP（···）.
之前看到有老师的文章和帖子里写，先用B3LYP或者CAM-B3LYP或者M062X优化基态和激发态结构，再用LC-wPBE计算激发态的单点能，以此能量为S1的能量。看到一个帖子里在计算单点能的时候关键词写的是#p b3lyp/6-31g(d) td scrf(solvent=toluene) nosymm
2.1 总感觉是在激发态的结构上又做了一次激发态，不太能理解这种做法，不写td的话，就感觉默认是基态？
2.2 加不加溶剂对这个delta ST影响很大吗？我的材料还没有合成，是不是计算的时候要猜测一下大概用什么可以溶解，在什么溶剂里测荧光，加上溶剂效应呢？想知道有没有必要
2.3 这个nosymm的意思是不保持对称性嘛，可是不加的话关键词里又没有opt，分子构型不是不会变的吗，单点能计算的就是当前结构下的能量，那为什么还要加这个关键词？
3.我在之前的计算过程是这样的：
s1计算和上面一样，
T1计算，上面只是第一步，然后又# td=triplets wb97xd/6-311g(d,p) 电荷和自旋多重度是 0 1，然后取这个激发态的能量（E(TD-HF/TD-KS)）和s1的能量相减，我觉得上面计算的s1的能量是发出荧光时候的能量，也就是s1的最低能量，这里计算的T1能量是发磷光时候的能量，是T1的最低能量，但是被说不对，说是用开头写的方法才对，也没有T0态，感觉有些晕，我还是不清楚，想知道对或者不对的原因。
4.
这个图，我为了理清思路画了下。
我的理解：当我们计算紫外吸收波长时，红色那一段是我们说的激发波长，荧光发射波长是计算的是蓝色的那一段，磷光发射是绿色那段。这样对不对？
4.2 红色那段与实验UV光谱的哪一个波长对应？起峰处，还是吸收强度最大处？
同理， 蓝色那段是荧光光谱的最大发射波长？（发射光强度最大处的波长？）绿色那段是磷光发射光谱的最大发射波长？
那小段黑色的是不是发射荧光以后，分子弛豫到基态s0的过程中的放出能量？棕色和绿色的差值那段，也是这个道理吗？
4.3 之前sob老师的博文中有写，激发态能量（E(TD-HF/TD-KS)）等于最后那个SCF Done能量加上激发能，我也验证过了。那在这个图里，比如S0->S1，这个SCF Done能量是A点的能量吗（对应到T1，SCF Done是B点的能量）？
4.4 还有横坐标，我看有好几种写法，应该是什么？
还请老师指点，谢谢啦！
|
|
发表于 Post on 2017-3-13 22:39:17
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author