求助在Gaussian中计算激发态加不加opt的问题

bobo

各位老师好：

计算激发态：#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的差别）？

感谢帮助。

Frank

http://sobereva.com/314

bobo

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态的意义是为什么，是不是加不加没有区别。

冰释之川

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

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

bobo

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

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

冰释之川

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

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

bobo

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

我完全懂了，非常感谢！

bobo

冰释之川 发表于 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的

冰释之川

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

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

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