计算化学公社

标题: ONIOM EmbedCharge [打印本页]

作者
Author: myyang 时间: 2017-6-20 16:01
标题: ONIOM EmbedCharge
Sob 老师，您好
最近在做ONIOM计算，处理能量是遇到疑惑：
当我的关键词为如下时：
oniom(m062x/6-31g**:amber=(softonly)) opt=(quadmac) freq nosymm scf(conver=6)

SCF Done:  E(RM062X) =  -883.437142081    A.U. after 1 cycles
ONIOM: gridpoint  1 method:  low system:  model energy:    0.598423197710
ONIOM: gridpoint  2 method:  high  system:  model energy:  -883.437142054235
ONIOM: gridpoint  3 method:  low system:  real  energy: -0.498654586223

ONIOM: extrapolated energy = -884.534219838169
SCF Done 的能量和 model system QM 能量相等

当关键词改变，加EmbedCharge 时：（我计算了另外一结构）
#p  oniom(m062x/6-31g**:amber=(softonly))=embedcharge opt=(quadmac)  nosymm scf(conver=6)
SCF Done:  E(RM062X) =  -915.130627735    A.U. after 2 cycles
ONIOM: calculating energy.
ONIOM: gridpoint  1 method:  low system:  model energy: -0.252072935201
ONIOM: gridpoint  2 method:  high  system:  model energy:  -883.779699054408
ONIOM: gridpoint  3 method:  low system:  real  energy: -1.092149799155
ONIOM: extrapolated energy = -884.619775918361

SCF Done 的能量和 model system QM 能量相差很大

我想知道加入EmbedCharge 关键词能量差别的原因以及此时SCF Done能量的意义
在最后能量判断时，我应该选取那个能量作为ONIOM 能量
谢谢

作者
Author: sobereva 时间: 2017-6-21 15:06
ONIOM: extrapolated energy 是最终要用的ONIOM单点能
EmbedCharge 是电子嵌入，此时MM区的原子电荷会对QM区产生静电势场，所以QM区能量明显变化

作者
Author: myyang 时间: 2017-6-24 09:28
谢谢老师！

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