计算化学公社

标题: 求助在Gaussian中计算激发态加不加opt的问题 [打印本页]

作者
Author: bobo 时间: 2020-5-14 18:34
标题: 求助在Gaussian中计算激发态加不加opt的问题
各位老师好：

计算激发态：#p cam-b3lyp/gen scrf=(solvent=toluene) TD(singlets, nstates=13) geom=connectivity，最后优化的结果:
Excited State 1:    Singlet-A    1.9235 eV  644.59 nm  f=0.9168  <S**2>=0.000
   339 -> 341       0.49792
   340 -> 342    -0.49771
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -4689.97020560
Copying the excited state density for this state as the 1-particle RhoCI density.
是不是说明现在得到的分子结构就是S1态的分子结构。

那如果加入关键词opt: #p cam-b3lyp/gen scrf=(solvent=toluene) TD(singlets, nstates=13, root=1) opt geom=connectivity后，优化S1结构，这样得到的结果和上面有什么区别？

加不加opt两种计算,得到的最终分子结构和各种跃迁的结果有差别吗（如果其他关键词完全相同，只有opt的差别）？

感谢帮助。

作者
Author: Frank 时间: 2020-5-14 18:47
http://sobereva.com/314

作者
Author: bobo 时间: 2020-5-14 18:51

Frank 发表于 2020-5-14 18:47
http://sobereva.com/314

你好，谢谢回答。我就是按老师这篇来计算的。
我就是不大明白opt加不加的区别，如果不加opt得出的已经是优化的S1态了（out文件中显示：This state for optimization and/or second-order correction.），那加opt优化S1态的意义是为什么，是不是加不加没有区别。

作者
Author: 冰释之川 时间: 2020-5-14 19:02

bobo 发表于 2020-5-14 18:51
你好，谢谢回答。我就是按老师这篇来计算的。
我就是不大明白opt加不加的区别，如果不加opt得出的已经是 ...

不加opt的话，算出来的是激发能，加opt之后，将会优化激发态结构，对应算出来的是从该优化好的激发态回到基态的发射能

作者
Author: bobo 时间: 2020-5-14 19:16

冰释之川发表于 2020-5-14 19:02
不加opt的话，算出来的是激发能，加opt之后，将会优化激发态结构，对应算出来的是从该优化好的激发态回到 ...

非常感谢，我懂了。还想再多问一下：加opt之后，最终结构是激发态结构，total energy是激发态结构能量。不加opt那个分子结构其实还是基态结构，total energy 对应的是基态结构能量。我这样理解对吗？

作者
Author: 冰释之川 时间: 2020-5-14 19:46

bobo 发表于 2020-5-14 19:16
非常感谢，我懂了。还想再多问一下：加opt之后，最终结构是激发态结构，total energy是激发态结构能量 ...

Total Energy, E(TD-HF/TD-DFT) = 这块都是激发态能量，只不过一个是基态结构下的S1激发态能量，另一个是结构弛豫后的S1激发态能量

作者
Author: bobo 时间: 2020-5-14 19:50

冰释之川发表于 2020-5-14 19:46
Total Energy, E(TD-HF/TD-DFT) = 这块都是激发态能量，只不过一个是基态结构下的S1激发态能量，另一个 ...

我完全懂了，非常感谢！

作者
Author: bobo 时间: 2020-6-9 08:34

冰释之川发表于 2020-5-14 19:46
Total Energy, E(TD-HF/TD-DFT) = 这块都是激发态能量，只不过一个是基态结构下的S1激发态能量，另一个 ...

老师，您好。
   我最近计算S0->T1的绝热吸收最后得到的是负值，请问是哪里出错了。具体操作如下：因为还需要计算S0->T2的绝热吸收，所以使用的是TDDFT做法。
1. 优化基态# opt freq b3lyp/6-311++g(d,p)，从out文件里面找到最后的 SCF Done:  E(RB3LYP) =  -954.540240772    A.U. after 1 cycles，
2. 在1优化的基态结构基础上优化激发态T1几何结构 # b3lyp/6-311++g(d,p) TD(triplets,nstates=10,root=1) opt freq，从out文件里面中找到最后一次的 This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -954.541231313
3.将第二步得到的能量与第一步得到的能量求差值 -954.541231313-（-954.540240772 ）=−9.90541×10−4<0 最后得到值是小于0的

作者
Author: 冰释之川 时间: 2020-6-9 09:15

bobo 发表于 2020-6-9 08:34
老师，您好。
我最近计算S0->T1的绝热吸收最后得到的是负值，请问是哪里出错了。具体操作如下： ...

首先仔细检查读取的数据有无纰漏，基态结构是否为极小点结构，如果没问题，也许该分子T1才是真正的基态。

P.s.几何结构优化在绝大多数情况下不需要加弥散函数，你所使用的6-311++g(d,p)基组十分糟糕，还不如用6-311g(d,p)基组

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