计算化学公社

标题: ONIOM模型通过QM/MM计算聚集诱导发光体系计算咨询 [打印本页]

作者
Author: 晴天 时间: 2018-10-10 10:52
标题: ONIOM模型通过QM/MM计算聚集诱导发光体系计算咨询
各位老师好，最近在学习用ONIOM模型通过QM/MM计算聚集诱导发光体系，主要用QM算单分子的发光性质，MM来考虑周围分子影响。在计算中有些疑惑，特来咨询，谢谢各位老师的指点。
其中优化基态和激发态的关键词设置如下（应该是对的吧）：
优化S0：#p opt freq oniom(b3lyp/6-31g(d):uff=qeq) geom=connectivity pop=CHELPG
优化s1：#p opt freq oniom(b3lyp/6-31g(d) td:uff=qeq) geom=connectivity pop=CHELPG

在计算单三态能级差时S1和T1的能量统计是Total Energy, E(TD-HF/TD-KS) =  -1795.04041106
还是ONIOM: gridpoint  1 method:  low system:  model energy:    0.342191713515
ONIOM: gridpoint  2 method:  high  system:  model energy: -1795.034155441291
ONIOM: gridpoint  3 method:  low system:  real  energy: 16.002803286222
ONIOM: extrapolated energy = -1779.373543868584里ONIOM: extrapolated energy的能量？
刚接触，有些不懂还请老师们多多指点，谢谢！

作者
Author: sobereva 时间: 2018-10-10 18:59
关键词不对。你要想表现低层原子对高层产生的静电相互作用，ONIOM后头得加=EmbedCharge。你写pop=CHELPG没有任何意义

ONIOM: extrapolated energy是整体的能量，E(TD-HF/TD-KS)是high层那一块的能量。对于当前目的取后者为宜

作者
Author: 晴天 时间: 2018-10-11 16:07

sobereva 发表于 2018-10-10 18:59
关键词不对。你要想表现低层原子对高层产生的静电相互作用，ONIOM后头得加=EmbedCharge。你写pop=CHELPG没 ...

刚开始摸索，有些疑问，非常感谢老师的指点。

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