模拟硫酸亚铁超胞加显式水层，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周期性。

sobereva

你是指molden文件？

原子电荷和氧化态完全是两码事，要算氧化态用下文的方法
使用Multiwfn结合CP2K计算晶体中原子的氧化态
http://sobereva.com/711（http://bbs.keinsci.com/thread-45876-1-1.html）


要考察电子转移量就用原子电荷说事，甭管绝对大小、就看变化量，并且别像你老师那样人为臆测而无根据地否定计算结果

看下文了解原子电荷和氧化态的常识，二者在绝对数值上天差地别
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表！
http://sobereva.com/714（http://bbs.keinsci.com/thread-46067-1-1.html）
使用Multiwfn通过LOBA方法计算氧化态
http://sobereva.com/362（http://bbs.keinsci.com/thread-5118-1-1.html）

user_fu

sobereva 发表于 2026-6-9 05:04
你是指molden文件？

原子电荷和氧化态完全是两码事，要算氧化态用下文的方法

好的，谢谢sob老师