带电周期slab吸附分子算出来相互作用能特别高怎么办

Aridea

各位大佬，我最近用cp2k MT solver算了一套XY周期带电slab(+2e, -2e这种)吸附分子（A）的体系，计算后用相互作用能ΔE=E(slab_A+2e)-E(slab+2e)-E(A)衡量相互作用大小。中性体系还算正常，0.8-3eV的作用能都有（一个解离吸附5eV）；
奇怪的是带电体系作用能比较离谱：-4到-11eV的相互作用能都出来了，ρ(slab_A+2e)-ρ(slab+2e)-ρ(A)查看荷电体系电子转移曲线明显比中性来的猛烈。其中一个分子在与slab之间”电子转移后“相互排斥分离后被排斥到盒子中间，此时计算ΔE还有-3.6eV。这里太奇怪了不知怎么解释，感觉直接这么计算带电体系似乎不太合理（slab电势降低可能比较夸张），不过考虑算了这么久，而且算都算完了，有没有可能在不换软件前提下。换一种相互作用能计算方法或有他办法补救一下嘛，还请大佬们给点意见建议呀~ 不胜感激