ONIOM EmbedCharge

myyang

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 能量
谢谢

sobereva

ONIOM: extrapolated energy 是最终要用的ONIOM单点能
EmbedCharge 是电子嵌入，此时MM区的原子电荷会对QM区产生静电势场，所以QM区能量明显变化

myyang

谢谢老师！