让你变成回忆 发表于 2017-11-3 15:10 好的,谢谢 |
zhanghuize 发表于 2017-11-2 11:22 结果基本没差距,有些时候用0 3的结果反而会更好。 最主要的原因是用0 3优化的速度远快于opt td. |
yezhonghua 发表于 2016-4-18 09:13 你好,我想问一下,意思是你用的这两种方法算T1,只是方法不同,但结果是一样的,是吗? |
sobereva 发表于 2016-4-17 22:25 谢谢老师!我明白了,我一直以为UDFT也是计算S0-T1的吸收能,只是方法不一样,看了是我的理解错了! |
yezhonghua 发表于 2016-4-17 21:30 无论是TDDFT还是DFT都能算S0和T1的能量差,不同几何结构下能量差是不一样的。 你在基态平衡结构下去算能量差,得到的能量差是S0->T1垂直吸收能。你如果先用UDFT优化出T1的平衡结构,再去算能量差,得到的就是T1->S0荧光垂直发射能。 思路必须清楚。 |
本帖最后由 yezhonghua 于 2016-4-17 21:32 编辑 sobereva 发表于 2016-4-16 23:46 谢谢老师,我明白你的这两种方法。 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)基本一样。 |
yezhonghua 发表于 2016-4-16 22:44 #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多花很多时间,根本没意义 |
本帖最后由 yezhonghua 于 2016-4-16 22:46 编辑 sobereva 发表于 2016-4-15 03:16 老师您好! 按照您的方法,我对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 |
bomsaude 发表于 2016-4-15 01:11 这样只是对T1进行构象优化,并不能得到T1的能量。我试了一下! |
小范范1989 发表于 2016-4-15 08:06 谢谢! |
sobereva 发表于 2016-4-15 03:16 这些细节问题弄明白好一点,谢谢老师! |
| 对,楼上两位老师说的对,你看看sob老师的那个帖子先,不着急算,先理解一下,这样你后面会减少很多错误。我的感觉,不一定正确哈。 |
|
直接用2L的方法就完了,这里也写得很明确 Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法 http://sobereva.com/314 什么TD-SCF并不是标准用词,不要管gview里显示什么 |
|
优化T1用DFT的话,网上说过的,似乎是这样的: #p opt ub3lyp/6-31g(d) T1 opt 0 3 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-22 05:41 , Processed in 0.297717 second(s), 25 queries , Gzip On.