|
本帖最后由 751013040 于 2022-11-18 17:20 编辑
各位老师,我在计算羟基化针铁矿吸附磷酸盐时(PBE泛函),计算出的吸附能为-60 eV ,只能看出吸附前的slab模型和吸附后的total模型能量差异巨大,但不知是什么原因引起的,还请老师们多多指教!
(化学式:2≡FeOH+H3PO4→≡Fe2O2PO2H+2H2O,吸附能公式为Eads = Etotal-Eslab-EH3PO4+2EH2O)
目前可以确定的是①OUTCAR中所有原子都读进去了,不是前后原子数不一致的问题②测试了偶极校正对能量的影响,应该也不是这方面的问题
有一点是吸附前后体系的总磁矩发生了改变(0→6.52,尤其是Fe原子的总磁矩:0→5.48),因为体系是反铁磁结构,理论上Fe的总磁矩应该是0,但吸附后铁有净磁矩了,不知道是不是这个原因引起的。如果是的话,有没有什么办法可以固定Fe原子的磁矩呢?
附件中有吸附前后的CONTCAR、INCAR以及OUTCAR
|
-
2.png
(32.15 KB, 下载次数 Times of downloads: 15)
-
1.png
(155.92 KB, 下载次数 Times of downloads: 15)
-
图片2.png
(238.53 KB, 下载次数 Times of downloads: 15)
-
图片1.png
(251.99 KB, 下载次数 Times of downloads: 15)
-
-
file.rar
4.26 MB, 下载次数 Times of downloads: 18
|