有机分子三重态激发态计算求助

zoey

老师好，刚接触一个涉及激发态的课题，需要计算有机分子三重态时激发态的能量，第一次算激发态有很多不懂的地方。在计算上文献提到：“We studied the energetics of 3a(singlet) and 3a* (triplet) computationally using the B3LYP/6-311+G(2d,p) method” ，不太明白3a的singlet和3a*的triplet优化计算怎么处理的
3a的singlet优化是 # opt  B3LYP/6- 311+G(2d,p) freq  电荷0 自旋多重度1
3a*的triplet优化是# opt uB3LYP/6- 311+G(2d,p) td freq    电荷0 自旋多重度3
如上设置关键词合适吗？3a*的triplet优化B3LTP能加U吗？

其次，文献提到“Thecomputed energy gap (58.5 kcal mol−1) was lower than thetriplet energy of Ir*[dF(CF3)ppy]2(dtbbpy)+ (60.1 kcalmol−1), suggesting that energy transfer from the catalyst tothe substrate was feasible” energy gap应该是指上面3a的singlet和3a*的triplet的差，是不是3a的singlet优化时需要加td才能和3a*的triplet相比较？

喵星大佬

优化激发态结构不需要写td，只要你自旋写的不是1程序会默认用udft所以写不写ub3lyp无所谓
不用加弥散函数，想要更高精度可以统一用def2tzvp
再就是优化激发态结构的时候B3LYP不算稳妥，PBE0(Gaussian里写pbe1pbe，直接写pbe0是个巨坑)，M06-2X或者wB97XD比较好
输入文件写作
#p opt pbe1pbe 6-311G**

Title

0  3
坐标略

这样就可以了
用TDDFT方法优化过渡态结构也可以要用td(triplet)关键字
#p opt pbe1pbe  6-311G**  td(triplet)

Title

0  1
坐标略
就可以了，但是结果和上面的没区别，耗时倒是会增加几倍，属于比较2的做法。

liyuanhe211

这儿都有答案：http://sobereva.com/314
另外2L说的有几点错误，比如“但是（TD）结果和上面的（△SCF）没区别”是错误的，具体看上文就知道了。

sobereva

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“激发态计算”改了，以后务必注意

td结合0 3优化的根本不是T1，得到的是什么根本没法预料
怎么优化T1态这里都写明了，仔细阅读，别乱用关键词
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

把文献里的段落粘过来的时候注意缺失空格的问题，这造成别人阅读费劲。

zoey

喵星大佬 发表于 2020-6-27 21:12
优化激发态结构不需要写td，只要你自旋写的不是1程序会默认用udft所以写不写ub3lyp无所谓
不用加弥散函数 ...

多谢指导

zoey

liyuanhe211 发表于 2020-6-28 05:27
这儿都有答案：http://sobereva.com/314
另外2L说的有几点错误，比如“但是（TD）结果和上面的（△SCF）没 ...

谢谢 提供的资料

zoey

sobereva 发表于 2020-6-28 19:16
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反 ...

好的，标题和粘贴问题以后会注意

谢谢Sob老师

yiwind

我得到了 三重态的结构，现在想计算三重态的单点能，自选多重度是设置 “0 1”还是设置为“0 3”呀？

zjxitcc

yiwind 发表于 2023-2-11 14:27
我得到了 三重态的结构，现在想计算三重态的单点能，自选多重度是设置 “0 1”还是设置为“0 3”呀？

（1）是自旋，不是自选；
（2）说清楚要用什么理论方法算，以及基态是否为闭壳层单重态

柚子

zjxitcc 发表于 2023-2-11 15:44
（1）是自旋，不是自选；
（2）说清楚要用什么理论方法算，以及基态是否为闭壳层单重态

老师您好，请问如果基态为闭壳层单重态，DFT，若要计算T1态的单点能，自旋多重度是“0 1”还是“0 3”呢？

zjxitcc

柚子 发表于 2023-4-18 16:10
老师您好，请问如果基态为闭壳层单重态，DFT，若要计算T1态的单点能，自旋多重度是“0 1”还是“0 3”呢 ...

有两种计算方式，以高斯为例：
（1）UDFT计算，坐标上面一行的自旋多重度写3，不要搭配TD关键词；
（2）TDDFT计算，坐标上面一行的自旋多重度写1，搭配TD(triplet)关键词。

柚子

zjxitcc 发表于 2023-4-18 16:31
有两种计算方式，以高斯为例：
（1）UDFT计算，坐标上面一行的自旋多重度写3，不要搭配TD关键词；
（2 ...

感谢老师的回答。
这两种方法我都做了，但是两个能量相差较大。研究对象是一个TADF分子，计算得到的S1态能量为40kcal/mol，我是以基态单重态的能量为基准，用UDFT得到的能量为33kcal/mol，用TDDFT得到的能量为40kcal/mol。
请问老师这种情况下怎么选择呢？

wzkchem5

柚子 发表于 2023-4-18 10:09
感谢老师的回答。
这两种方法我都做了，但是两个能量相差较大。研究对象是一个TADF分子，计算得到的S1态 ...

建议用高精度方法算单点能，例如双杂化泛函、CC2、ADC(2)、DLPNO-STEOM-CCSD、NEVPT2等

柚子

wzkchem5 发表于 2023-4-18 17:15
建议用高精度方法算单点能，例如双杂化泛函、CC2、ADC(2)、DLPNO-STEOM-CCSD、NEVPT2等

谢谢老师的建议，我尝试一下。

symbolhone

wzkchem5 发表于 2023-4-18 17:15
建议用高精度方法算单点能，例如双杂化泛函、CC2、ADC(2)、DLPNO-STEOM-CCSD、NEVPT2等

wzkchem5老师，您好。我最近在用双杂化泛函计算TADF分子ΔEST，我是先用gaussian做的opt，然后再用ORCA做的双杂话泛函计算，其中一个分子，它的out里显示这样：
TD-DFT/TDA EXCITED STATES (SINGLETS)
STATE  1:  E=   0.099942 au      2.720 eV    21934.7 cm**-1 <S**2> =   0.000000
STATE  2:  E=   0.116849 au      3.180 eV    25645.5 cm**-1 <S**2> =   0.000000
STATE  3:  E=   0.126277 au      3.436 eV    27714.6 cm**-1 <S**2> =   0.000000

TD-DFT/TDA EXCITED STATES (TRIPLETS)
STATE  1:  E=   0.102962 au      2.802 eV    22597.6 cm**-1 <S**2> =   2.000000
STATE  2:  E=   0.096927 au      2.638 eV    21273.0 cm**-1 <S**2> =   2.000000
STATE  3:  E=   0.114811 au      3.124 eV    25198.2 cm**-1 <S**2> =   2.000000

其中STATE  1 (TRIPLETS)的能量是2.802eV， STATE  2 (TRIPLETS)则是2.638eV，这个是什么原因导致的呢，出现这种情况是不是意味着我这个计算任务就需要重新进行呢。