计算化学公社

标题: oniom中计算绝热激发能时S1态的能量统计哪一个？ [打印本页]

作者
Author: 暖空 时间: 2019-11-15 21:29
标题: oniom中计算绝热激发能时S1态的能量统计哪一个？
本帖最后由暖空于 2019-11-28 20:07 编辑

大佬们好！最近在学习用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)区域的能量相减呢？
谢谢指导！

作者
Author: sobereva 时间: 2019-11-16 12:58
绝热激发能是S1减S0，而不是S0减S1

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

作者
Author: 暖空 时间: 2019-11-16 16:38

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的输出已上传，麻烦您了。

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

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

作者
Author: 暖空 时间: 2019-11-17 16:14

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

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

作者
Author: 暖空 时间: 2019-11-28 19:56

暖空发表于 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文件我已压缩上传，麻烦老师有时间帮忙看一下，谢谢老师。

[attach]22646[/attach][attach]22647[/attach]

作者
Author: 暖空 时间: 2019-11-29 14:49
本帖最后由暖空于 2019-11-29 14:51 编辑

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

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

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

作者
Author: sobereva 时间: 2019-11-30 03:11

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

我已删除

作者
Author: 暖空 时间: 2019-11-30 12:07

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

嗯嗯，谢谢老师。

作者
Author: jiangyan 时间: 2024-7-1 22:03
请问最后的绝热激发能时S1态的能量统计哪一个呀？

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