计算化学公社

标题: 关于对于T1进行opt问题 [打印本页]

作者
Author: yezhonghua 时间: 2016-4-15 00:16
标题: 关于对于T1进行opt问题
文献上的意思对于S0、T1用DFT优化结构；对于S1用TD-DFT进行结构优化。
T1我按照文献进行计算完全一致，是在S0构型优化的基础上进行垂直吸收计算，td=triplet，电荷自由度0，1。
但是按照T1用DFT进行构型优化的话，method是选TD-SCF还是ground state？如果选TD-SCF那不就是TD-DFT了；如果选用ground state的话，这里的电荷自由度是0，1还是0，3？
这一点不太明白？
谢谢！

作者
Author: bomsaude 时间: 2016-4-15 01:11
优化T1用DFT的话，网上说过的，似乎是这样的：

#p opt ub3lyp/6-31g(d)

T1 opt

0 3

作者
Author: sobereva 时间: 2016-4-15 03:16
直接用2L的方法就完了，这里也写得很明确
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314
什么TD-SCF并不是标准用词，不要管gview里显示什么

作者
Author: 小范范1989 时间: 2016-4-15 08:06
对，楼上两位老师说的对，你看看sob老师的那个帖子先，不着急算，先理解一下，这样你后面会减少很多错误。我的感觉，不一定正确哈。

作者
Author: yezhonghua 时间: 2016-4-15 22:58

sobereva 发表于 2016-4-15 03:16
直接用2L的方法就完了，这里也写得很明确
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
htt ...

这些细节问题弄明白好一点，谢谢老师！

作者
Author: yezhonghua 时间: 2016-4-15 22:58

小范范1989 发表于 2016-4-15 08:06
对，楼上两位老师说的对，你看看sob老师的那个帖子先，不着急算，先理解一下，这样你后面会减少很多错误。 ...

谢谢！

作者
Author: yezhonghua 时间: 2016-4-16 22:33

bomsaude 发表于 2016-4-15 01:11
优化T1用DFT的话，网上说过的，似乎是这样的：

#p opt ub3lyp/6-31g(d)

这样只是对T1进行构象优化，并不能得到T1的能量。我试了一下！

作者
Author: yezhonghua 时间: 2016-4-16 22:44
本帖最后由 yezhonghua 于 2016-4-16 22:46 编辑

sobereva 发表于 2016-4-15 03:16
直接用2L的方法就完了，这里也写得很明确
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
htt ...

老师您好！
按照您的方法，我对TD-DFT计算得到的T1进行opt优化，由于想要和文献的结果进行对照，所以在T1的基础上采用了不同的方法;
1. #p opt b3lyp/6-31g(d)  0,3
2. #p opt td=(nstates=3) b3lyp/6-31g(d)  0,3
第一种方法其实是对T1进行opt优化，通过他的E-TDHF和S0 opt的E-TDHF得到的是 UDFT的能量。
第二种方法应该是T1的垂直吸收能进行opt优化得到其垂直发射T1能量。
第一种显然不是文献要的T1的垂直发射能，但是第二种方法计算结果如下：
Excited State 1:  3.009-A    0.7207 eV 1720.30 nm  f=0.0237  <S**2>=2.014
显然不对。
谢谢老师！
我准备试试第3种方法：
3. #p opt td=(nstates=3) b3lyp/6-31g(d)  0,1 td=triplet

作者
Author: sobereva 时间: 2016-4-16 23:46

yezhonghua 发表于 2016-4-16 22:44
老师您好！
按照您的方法，我对TD-DFT计算得到的T1进行opt优化，由于想要和文献的结果进行对照，所以 ...

#p opt td=(nstates=3) b3lyp/6-31g(d) 0,3这什么都不是

过程在我的帖子里写得很明确了
就是
1 #P opt B3LYP/6-31G* 0 3
2 #P TD(triplet) B3LYP/6-31G* 0 1
这么简单

你想用TDDFT优化T1只会比直接0 3多花很多时间，根本没意义

作者
Author: yezhonghua 时间: 2016-4-17 21:30
本帖最后由 yezhonghua 于 2016-4-17 21:32 编辑

sobereva 发表于 2016-4-16 23:46
#p opt td=(nstates=3) b3lyp/6-31g(d) 0,3这什么都不是

过程在我的帖子里写得很明确了

谢谢老师，我明白你的这两种方法。
1.是UDFT方法计算T1
2.是TD-DFT计算T1
两种方法应该都是计算的垂直吸收的T1，也就是EVA(T1），因为他们都是在构象不变的基础上计算得到的。但是我想要得到EVE（T1）也就是垂直激发的T1能量。文献先是通过TD-DFT算测EVA(T1），然后再在激发态T1的基础上进行opt算EVE（T1）。我通过#p opt td=(nstates=3) b3lyp/6-31g(d) 0,1 td=triplet 算出来额结果和文献差0.007ev。
请问一下老师，不知我这么认为对不对。还请老师指点一下，因为我通过UDFT计算的T1和TD-DFT计算的T1(EVA)基本一样。

作者
Author: sobereva 时间: 2016-4-17 22:25

yezhonghua 发表于 2016-4-17 21:30
谢谢老师，我明白你的这两种方法。
1.是UDFT方法计算T1
2.是TD-DFT计算T1

无论是TDDFT还是DFT都能算S0和T1的能量差，不同几何结构下能量差是不一样的。
你在基态平衡结构下去算能量差，得到的能量差是S0->T1垂直吸收能。你如果先用UDFT优化出T1的平衡结构，再去算能量差，得到的就是T1->S0荧光垂直发射能。
思路必须清楚。

作者
Author: yezhonghua 时间: 2016-4-18 09:13

sobereva 发表于 2016-4-17 22:25
无论是TDDFT还是DFT都能算S0和T1的能量差，不同几何结构下能量差是不一样的。
你在基态平衡结构下去算 ...

谢谢老师！我明白了，我一直以为UDFT也是计算S0-T1的吸收能，只是方法不一样，看了是我的理解错了！

作者
Author: zhanghuize 时间: 2017-11-2 11:22

yezhonghua 发表于 2016-4-18 09:13
谢谢老师！我明白了，我一直以为UDFT也是计算S0-T1的吸收能，只是方法不一样，看了是我的理解错了！

你好，我想问一下，意思是你用的这两种方法算T1，只是方法不同，但结果是一样的，是吗？

作者
Author: 让你变成回忆 时间: 2017-11-3 15:10

zhanghuize 发表于 2017-11-2 11:22
你好，我想问一下，意思是你用的这两种方法算T1，只是方法不同，但结果是一样的，是吗？

结果基本没差距，有些时候用0 3的结果反而会更好。
最主要的原因是用0 3优化的速度远快于opt td.

作者
Author: zhanghuize 时间: 2017-11-5 11:26

让你变成回忆发表于 2017-11-3 15:10
结果基本没差距，有些时候用0 3的结果反而会更好。
最主要的原因是用0 3优化的速度远快于opt td.

好的，谢谢

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