求助：SMD水模型下计算的相互作用能过小

陌上花开小伟

各位老师：     我在计算咪唑阳离子和醋酸根阴离子的相互作用能，我在SMD水模型下计算出来的相互作用能为-16.72KJ/mol。 我觉得过小，然后我又在这个结构基础上计算得到气相条件下相互作用能为-395.69KJ/mol。我觉得是我隐式水模型计算出了问题，想请各位老师帮忙分析一下是什么问题。
我采用的相互作用能计算公式：E(AB)-E(A)-E(B) ，直接从稳定的AB结构输出文件中剪出单独的A B 结构并直接计算单点能
结构优化使用的关键词：# opt=gdiis freq b3lyp/6-311+g(d,p) scrf=(smd,solvent=water) pop=nbo em=gd3 int=acc2e=12 scf=novaracc
能量计算使用关键词 b3lyp/6-311+g(d,p) scrf=(smd,solvent=water) pop=nbo em=gd3  或者 b3lyp/6-311+g(d,p) pop=nbo em=gd3
因为用了隐式水模型后不太容易收敛，我就加了gdiis。

sobereva

没有什么问题
极性溶剂会把阴阳离子间静电作用极大幅度屏蔽掉