multiwfn计算model文件输出电荷变化求助

user_fu

各位老师好，我的模型整体为硫酸亚铁超胞加显式水层，模拟发生在表面暴露的三配位二价Fe原子与水分子的反应，对铁原子开启了自旋极化，同时整体为反铁磁设置，设置了磁矩，表层铁为负方向磁矩，我用CP2K跑表面反应的AIMD，同时在这个过程中我每200步输出了model文件，通过model文件想对涉及反应的二价Fe原子的电荷变化进行输出，但是通过multiwfn计算model文件发现，二价Fe原子在吸附了两个水分子后，CM5电荷变化了0.13e，而二价Fe原子初始CM5电荷为0.69e，老师说我的二价铁原子电荷偏小，我想请教一下各位老师，这会是什么原因,二价铁的计算参数是DZVP-MOLOPT-SR-GTH-q16，整体xy周期性。