请问为何质子化的能量变化在有/无溶剂模型下差别如此大？

sai77

110个原子的体系，含过渡金属，优化PBE0-D3BJ/6-31G*-SDD，单点PWPB95-D3BJ/def2TZVPP-SMD
其中一步反应是底物中的亚胺被酸质子化，
R-C=N-R'  +   TsOH   ->    RC=NH(+)-R'   +   TsO-
按照上述级别算下来这步反应自由能变化大概-60 kcal/mol，明显不对。
但是如果结构优化的时候加上溶剂模型IEFPCM，能量变化就是-6 kcal/mol。

想请问一下为什么结构优化时加不加溶剂模型对结果影响如此之大？谢谢！

sobereva

对比优化完的结构
溶剂环境可能严重影响局部带显著电荷体系的几何结构，下文说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）