计算化学公社

标题: 求助三重态T1到Tn的吸收怎么计算-高斯 [打印本页]

作者
Author: xinren 时间: 2019-9-5 20:34
标题: 求助三重态T1到Tn的吸收怎么计算-高斯
本帖最后由 xinren 于 2019-9-5 20:36 编辑

最近在做三重态T1到Tn三重态之间的吸收，不知道怎么设输入文件，之前想的是T1作为三重态的基态Tn是三重态的激发态，输入文件是
第一步：优化T1:opt freq 0 3
第二步：td=(nstates=6) 0 3
输出文件：
Excited State 1:  3.020-A    1.1960 eV 1036.70 nm  f=0.0000  <S**2>=2.029
   38B -> 39B       0.98422
   38B -> 40B       0.11768
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -498.243571305
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:  3.019-A    1.6811 eV  737.54 nm  f=0.0021  <S**2>=2.029
   40A -> 41A       0.46357
   40A -> 43A    -0.15651
   37B -> 39B       0.86550
不知道文件输入对不对，请各位指导一下！另外还需要看T1到Tn的跃迁轨道，不知道具体操作是什么?

作者
Author: sobereva 时间: 2019-9-6 00:31
第二步关键词明显是错的。你当前的基态（作为参考态）是单重态，显然第二步必须写0 1

Gaussian只能给出激发态之间的能量差，但是激发态之间的跃迁偶极矩Gaussian没法直接给出来，如果要考虑这点，需要用Multiwfn
使用Multiwfn计算激发态间的跃迁偶极矩和各个激发态的偶极矩
http://sobereva.com/227

没有所谓的激发态之间的跃迁轨道一说，每个激发态都是基于相同的参考态的分子轨道来描述的

作者
Author: xinren 时间: 2019-9-6 15:34
好的谢谢老师，那是不是把第二步的0 3 改成0 1求得的能量差就是Tn与T1之间的能差吗？之前看到求磷光的方法是
1 优化T1结构：# B3LYP/6-31G* opt，电荷和自旋多重度为0 3。
2 取上一步优化得到的结构，用# B3LYP/6-31G* TD(triplet)，电荷和自旋多重度为0 1。第一激发态的激发能就是磷光发射能了。
请问如果改成0 1,那求出的能差不是就表示第一激发态是T1与S0的能差，第二激发态就是表示T2与S0的能差吗？另外是用singlet还是用triplet呢？
另外用Multiwfn求Tn与T1的跃迁矩的输入文件是B3LYP/6-31+G(d) TD(triplet，nstates=10) IOp(9/40=4) 0 1吗？

作者
Author: sobereva 时间: 2019-9-7 03:05

xinren 发表于 2019-9-6 15:34
好的谢谢老师，那是不是把第二步的0 3 改成0 1求得的能量差就是Tn与T1之间的能差吗？之前看到求磷光的方法 ...

每个结构下都有各个Tn与T1的能量差。显然你把当前结构下算出来的三重态激发能之间求差就能得到当前结构下的Tn与T1能量差

你要算三重态激发态之间的能量差，显然得写triplet，否则算出的激发态都是单重态的有什么用？

是

作者
Author: xinren 时间: 2019-9-22 17:14
好的，谢谢老师！之前Tn与T1的能量差算出来了，想要再问一下老师，相应的Tn与T1的分子轨道该怎么计算呢？

作者
Author: sobereva 时间: 2019-9-23 08:16

xinren 发表于 2019-9-22 17:14
好的，谢谢老师！之前Tn与T1的能量差算出来了，想要再问一下老师，相应的Tn与T1的分子轨道该怎么计算呢？

对于TDDFT，根本没法说什么T1、Tn的轨道。顶多是说在相应激发态优化的极小点结构下算的基态轨道
(, 下载次数 Times of downloads: 54)

作者
Author: xinren 时间: 2019-9-30 10:31
好的，谢谢！！！！！！

作者
Author: xinren 时间: 2019-10-2 11:24
老师，实验测得的三重态吸收峰有两个，一个大约是在400nm处，一个大约是在450nm处，但是高斯计算出来的振子强度比较强的两个峰值都比较集中在400nm区域，分别是390nm和410nm左右，请问会不会是计算方法引起的误差？

作者
Author: sobereva 时间: 2019-10-3 04:28

xinren 发表于 2019-10-2 11:24
老师，实验测得的三重态吸收峰有两个，一个大约是在400nm处，一个大约是在450nm处，但是高斯计算出来的振子 ...

可能

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