Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法

sobereva · 发表于 Post on 2020-6-30 00:28:57

clayer 发表于 2020-6-29 22:17
哈哈，感谢sob，下次他们再问我，我就让他们去问guassion客服。

注意Gaussian的拼写

Believe · 发表于 Post on 2020-8-1 15:07:41

sob老师您好，如果基态是三重态（T0），计算激发态，会存在单重态激发态吗？还是只能做三重态的激发态。如果可以，如何用三重态基态计算单重态激发态呢？

喵星大佬 · 发表于 Post on 2020-8-1 15:47:32

Believe 发表于 2020-8-1 15:07
sob老师您好，如果基态是三重态（T0），计算激发态，会存在单重态激发态吗？还是只能做三重态的激发态。如 ...

那你觉得氧气怎么办

sobereva · 发表于 Post on 2020-8-2 04:10:45

Believe 发表于 2020-8-1 15:07
sob老师您好，如果基态是三重态（T0），计算激发态，会存在单重态激发态吗？还是只能做三重态的激发态。如 ...

诸如单重态氧气就是理论化学经常被研究的对象。
你直接设0 1，都不用TDDFT，就等于算S1了（这相当于delta-SCF方法计算）

Believe · 发表于 Post on 2020-8-8 14:54:24

sob老师您好，我还有一个问题没弄明白。就是直接设置0,1可以计算单重态第一激发态，那我想要计算单重态第二激发态第三激发态应该怎么办呢？是以单重态第一激发态01为基态进行TDDFT计算吗？还是以三重态为基态03进行TDDFT计算？还是这种情况只能计算单重态第一激发态呢？

sobereva · 发表于 Post on 2020-8-9 04:09:41

Believe 发表于 2020-8-8 14:54
sob老师您好，我还有一个问题没弄明白。就是直接设置0,1可以计算单重态第一激发态，那我想要计算单重态第二 ...

帖子里已经写得极为清楚了，参考态是单重态的情况（诸如设0 1），TD里nstates>1显然就会把更高的激发态也算出来

Believe · 发表于 Post on 2020-8-10 20:35:40

好的，谢谢sob老师，我明白您的意思了。我设置的是01，然后TD(nstates=4),但是得到的<S**2>并不是单重态，而是更接近于二重态。
是不是因为我的基态是开壳层的原因，所以只能<S**2>是什么态就是什么态呢？可是我输入文件设置的是01呀。是我哪里设置的不对吗？
Excited State 1:  2.437-A    0.0268 eV 46270.33 nm  f=0.0000  <S**2>=1.234
186A ->187A       0.80621
186B ->187B    -0.80318
186A <-187A       0.38438
186B <-187B    -0.38384
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2848.39703229
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:  1.655-A    0.1147 eV 10813.63 nm  f=0.0000  <S**2>=0.435
186A ->187A       0.72083
186B ->187B       0.72384
186A <-187A    -0.14784
186B <-187B    -0.14744

Excited State 3:  2.228-A    0.5154 eV 2405.79 nm  f=0.0000  <S**2>=0.991
184A ->187A       0.67362
185A ->187A       0.30865
184B ->187B       0.65118
185B ->187B    -0.15000

Excited State 4:  2.260-A    0.5170 eV 2398.25 nm  f=0.0002  <S**2>=1.027
185A ->187A       0.65922
185B ->187B       0.73824

sobereva · 发表于 Post on 2020-8-11 07:01:29

Believe 发表于 2020-8-10 20:35
好的，谢谢sob老师，我明白您的意思了。我设置的是01，然后TD(nstates=4),但是得到的并不是单重态，而是更 ...

要么是你写了U，要么是你用了guess=read，读取的之前开壳层的chk

Believe · 发表于 Post on 2020-8-12 14:25:12

嗯嗯，谢谢sob老师,问题已解决。我确实是加了U的。

Believe · 发表于 Post on 2020-10-19 21:44:46

"TDDFT计算精度严重依赖于泛函，必须根据体系特征、激发类型选择合适的泛函才能得到不错结果。至于基组，若泛函选择得当，对于大体系用6-31G*就能得到定性正确结果"---sob老师您好，请问有关于6-31G*适合大体系计算激发态的相关文献吗？我有用6-31G*基组进行TDDFT激发态计算，但是审稿人质疑此基组不适合计算激发态。认为基组过小，我没有找到合适的文献支撑去反驳他

sobereva · 发表于 Post on 2020-10-21 12:24:47

Believe 发表于 2020-10-19 21:44
"TDDFT计算精度严重依赖于泛函，必须根据体系特征、激发类型选择合适的泛函才能得到不错结果。至于基组，若 ...

你拿其它大把用6-31G*的文章怼他就完了，google学术随便一搜就一大堆
或者你也用6-311G*算，只算一个体系，少算一些态，和6-31G*对比，证明对结果影响小

例如TDDFT计算茜素染料分子，明显6-31G*是最小可接受级别，而继续加大也影响不显著，显然对于大体系耗时的计算是很适合的

Believe · 发表于 Post on 2020-10-23 20:58:38

谢谢sob老师，我明白了~

lyt0226 · 发表于 Post on 2021-4-24 11:54:11

各位老师大家好，

我在通过势能面扫描确定一个关于二价苯解离成两个一价碎片的路径，中间涉及一些质子转移的过渡态，这些过渡态和中间体的寻找用的方法和基组是m062x/def2TZVPP ，单点能计算用的是ccsd(t)/aug-cc-pvtz，
（1）现在想计算二价苯的激发态的能量，因为这些在一个反应路径里面，能量对基组和算法要求比较高，所以想用和其他过渡态一样的计算级别来计算它的激发态能量，也就是ccsd(t)/aug-cc-pvtz，请问各位老师有必要吗？
（2）如果想实现ccsd(t)/aug-cc-pvtz计算激发能可能吗？如果不能实现，那对于同一路径里面，初态的激发能和后面过渡态等能量的计算不是一个基组和方法有关系吗？
（3）应该用什么方法，td或者EOM，我现尝试用#p eomccsd ccsd(t)/aug-cc-pvtz进行计算，但是一直报错"Internal consistency failure #1 in ROv08"

谢谢各位老师！

sobereva · 发表于 Post on 2021-4-25 00:10:03

lyt0226 发表于 2021-4-24 11:54
各位老师大家好，

我在通过势能面扫描确定一个关于二价苯解离成两个一价碎片的路径，中间涉及一些质子转 ...

这种问题应当单独发主贴

wzkchem5 · 发表于 Post on 2021-4-25 01:38:20

lyt0226 发表于 2021-4-24 11:54
各位老师大家好，

我在通过势能面扫描确定一个关于二价苯解离成两个一价碎片的路径，中间涉及一些质子转 ...

支持在CCSD(T)级别做激发态（EOM-CCSD(T)）的程序很少，高斯是不支持的，这就是你报错的原因。
如果你要研究的激发态没有显著的多参考态性质，那么可以用DeltaSCF方法在CCSD(T)级别算；如果激发态是多参考的，可以用EOM-CCSD/aug-cc-pVTZ算，省略了(T)不是大问题，但是激发态的能量应该用CCSD(T)基态能量加上EOM-CCSD的激发能，不能直接用EOM-CCSD的激发态能量。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法