oniom中计算绝热激发能时S1态的能量统计哪一个？

暖空

大佬们好！最近在学习用oniom模拟固体发光，在统计QM的绝热激发能时有点儿困惑，还望指教。


S0能量找的是最后一个extrapolated energy上面QM区的能量：
ONIOM: gridpoint  1 method:  low   system:  model energy:     0.271088170508
ONIOM: gridpoint  2 method:  high  system:  model energy: -1685.253804385772
ONIOM: gridpoint  3 method:  low   system:  real  energy:    11.178913478074
ONIOM: extrapolated energy =   -1674.345979078206                           

优化S1时只对QM区域做了td，我现在找的有两个能量
（1）和S0一样找的最后一个extrapolated energy上面QM区域的能量：
ONIOM: gridpoint  1 method:  low   system:  model energy:     0.255063841593  
ONIOM: gridpoint  2 method:  high  system:  model energy: -1685.146374708537  
ONIOM: gridpoint  3 method:  low   system:  real  energy: 50038.496166870114  
ONIOM: extrapolated energy =   48353.094728319986                             

（2）最后一个excited state 1下面的能量：
Excited State   1:      Singlet-A      2.3802 eV  520.90 nm  f=0.0049  <S**2>=0.000
     144 ->145         0.69561
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -1687.68914004   

我用S1态统计的第二种能量减去S0的能量，绝热能差为-2.43a.u. 感觉不对。
请问要是做S1和S0之间的绝热能差是不是用的都是整体能量extrapolated energy上面的QM区域的能量相减？
还有，如果要做S1和T1之间的绝热能差用的是extrapolated energy上面的QM区域的能量还是都用的是E(TD-HF/TD-DFT)区域的能量相减呢？
谢谢指导！

sobereva

绝热激发能是S1减S0，而不是S0减S1

如果你用ONIOM优化了S1，把整个输出文件压缩后上传。或者搜搜输出文件，弄清楚-1685.146374708537是从哪来的。

暖空

sobereva 发表于 2019-11-16 12:58
绝热激发能是S1减S0，而不是S0减S1

如果你用ONIOM优化了S1，把整个输出文件压缩后上传。或者搜搜输出文 ...

谢谢老师指点， 帖子中已更正。
我看不太懂S1的-1685.146374708537这个能量具体是怎么来的，我的理解这一块是高斯手册中说的QM层的能量，但是又感觉S1整体的extrapolated energy和上面的ONIOM: gridpoint  3 method的的能量都很奇怪，压缩的S0和S1的输出已上传，麻烦您了。

sobereva

拿S1优化的最后一个结构算个单点再看看。怀疑其中一个数据是没有考虑电子嵌入时的结果。也可以对比一下用不用电子嵌入时的结果试图搞明白原因。
文件太大，浏览速度超慢，建议用#T

另外，QEq电荷对静电势重现性非常差，此时并不能很好描述周围环境产生的静电作用

暖空

sobereva 发表于 2019-11-17 15:38
拿S1优化的最后一个结构算个单点再看看。怀疑其中一个数据是没有考虑电子嵌入时的结果。也可以对比一下用不 ...

恩恩，谢谢老师指导！
我用RESP电荷作为拟合静电势电荷算了其他的分子，我再对比对比，谢谢老师。

暖空

暖空 发表于 2019-11-17 16:14
恩恩，谢谢老师指导！
我用RESP电荷作为拟合静电势电荷算了其他的分子，我再对比对比，谢谢老师。

sob老师好！
在优化S1态的时候，我分别做了加不加电子嵌入两种计算，发现无论加不加最后的extrapolated energy和上面的三种能量都差不多（差小数点后几位），而且extrapolated energy和real energy都是正的4、5千，应该是做S1激发的时候，虽然我加了embedcharge关键词，但是计算的时候没考虑。请问这种是什么原因呢？
#T opt freq oniom(bmk/6-31g(d) td(nstates=7):uff=qeq)=embedcharge geom=connectivity
ONIOM: gridpoint  1 method:  low   system:  model energy:     0.255063567794
ONIOM: gridpoint  2 method:  high  system:  model energy: -1685.146374745875
ONIOM: gridpoint  3 method:  low   system:  real  energy: 50038.496166608551
ONIOM: extrapolated energy =   48353.094728294884

#T opt freq oniom(bmk/6-31g(d) td(nstates=7):uff=qeq) geom=connectivity
ONIOM: gridpoint  1 method:  low   system:  model energy:     0.258937528198
ONIOM: gridpoint  2 method:  high  system:  model energy: -1685.139792394517
ONIOM: gridpoint  3 method:  low   system:  real  energy: 50038.496227715761
ONIOM: extrapolated energy =   48353.097497793045

我也用了RESP电荷作为拟合静电式电荷做了S1，能量这一块的输出和qeq的类似。

相关的log文件我已压缩上传，麻烦老师有时间帮忙看一下，谢谢老师。

暖空

sobereva 发表于 2019-11-16 12:58
绝热激发能是S1减S0，而不是S0减S1

如果你用ONIOM优化了S1，把整个输出文件压缩后上传。或者搜搜输出文 ...

sob老师，您不用帮我看啦！
我应该已经找到原因了，现在正在测试，可能是我构型中C和H之间的连接信息丢了，gview打开C和H之间没有单键。
弱弱问一下，怎么把自己上传的附件删掉呢？

sobereva

暖空 发表于 2019-11-29 14:49
sob老师，您不用帮我看啦！
我应该已经找到原因了，现在正在测试，可能是我构型中C和H之间的连接信息丢 ...

我已删除

暖空

sobereva 发表于 2019-11-30 03:11
我已删除

嗯嗯，谢谢老师。

jiangyan

请问最后的绝热激发能时S1态的能量统计哪一个呀？