求助电场下二价金属离子与CE18膜的相互作用能的结果问题

海棠依旧cdut

请问各位大佬。我模拟的是在电场作用下，CE18膜对二价金属离子的渗透性行为研究。体系是水环境，模型是SPC/E，立场是opls，截断半径1.2nm。我做了一个离子与膜的相互作用的计算。计算代码是compute 1 Co group/group membrane kspace yes。最终计算出来的结果峰值太过于尖锐。并且该尖锐的峰值总是伴随着离子的渗透（即，离子穿过一个，该时间的相互作用能就会出现一个峰值）。由于本人之前没有进行过类似的研究，也没有其他类似的长这样的数据供参考。因此想问问各位大佬，interaction energy的图长这样是正确的吗？如果不正确是什么原因造成的呢（类似代码，参数有误等）？