计算化学公社

标题: vasp反铁磁羟基化针铁矿吸附磷酸吸附能 -60 eV? [打印本页]

作者
Author: 751013040 时间: 2022-11-18 11:03
标题: vasp反铁磁羟基化针铁矿吸附磷酸吸附能 -60 eV?
本帖最后由 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

作者
Author: 751013040 时间: 2022-11-18 11:04
吸附能那个是+2EH2O，这个打错了

作者
Author: 876449830 时间: 2022-11-18 11:09
ISPIN没有设置呐，U值测试过吗，一般看到的多的Fe是3.8-4.0eV的样子。

作者
Author: 751013040 时间: 2022-11-18 11:25

876449830 发表于 2022-11-18 11:09
ISPIN没有设置呐，U值测试过吗，一般看到的多的Fe是3.8-4.0eV的样子。

首先感谢回复~
对于ISPIN，先没开自旋极化预优化，然后再开ISPIN并设置磁矩为反铁磁。对于铁的U值我没有测试，是参照的一篇2012年用vasp研究U值对bulk针铁矿影响的（doi: 10.1016/j.matchemphys.2012.11.017）

作者
Author: Eudaimonia 时间: 2022-11-18 15:47
确定原子数，种类一致的话再检查基组电荷是不是都一致，这么大的值一般不太可能是构型不同引起的

作者
Author: 751013040 时间: 2022-11-18 16:50

Eudaimonia 发表于 2022-11-18 15:47
确定原子数，种类一致的话再检查基组电荷是不是都一致，这么大的值一般不太可能是构型不同引起的

感谢回复~
吸附前后原子数和种类是一致的，然后对于基组电荷我不太懂你的意思，我是检索了吸附前后的OUTCAR的NELECT（电子的总数）是对的上的，应该不是这方面的问题

作者
Author: ChemG 时间: 2022-11-19 11:35
你这个严格意义上讲不是吸附能，而是发生了磷酸和表面的脱水反应，在最终形成的结构中，不光产生了新的化学键，还产生了几个内氢键，这都会导致体系能量降低，体现在所谓的吸附能中。更确切来说，现在算出的能量是真空条件下磷酸分子与羟基氧化铁表面脱水反应的能量降低值。

作者
Author: 751013040 时间: 2022-11-19 19:12

ChemG 发表于 2022-11-19 11:35
你这个严格意义上讲不是吸附能，而是发生了磷酸和表面的脱水反应，在最终形成的结构中，不光产生了新的化学 ...

哦哦，好的好的，你这么一说我也感觉应该不能叫做吸附能了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)