计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1361|回复 Reply: 7
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] vasp反铁磁羟基化针铁矿吸附磷酸吸附能 -60 eV?

[复制链接 Copy URL]

8

帖子

0

威望

91

eV
积分
99

Level 2 能力者

本帖最后由 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)

2.png

1.png (155.92 KB, 下载次数 Times of downloads: 15)

1.png

图片2.png (238.53 KB, 下载次数 Times of downloads: 15)

图片2.png

图片1.png (251.99 KB, 下载次数 Times of downloads: 15)

图片1.png

file.rar

4.26 MB, 下载次数 Times of downloads: 18

8

帖子

0

威望

91

eV
积分
99

Level 2 能力者

2#
 楼主 Author| 发表于 Post on 2022-11-18 11:04:57 | 只看该作者 Only view this author
吸附能那个是+2EH2O,这个打错了

185

帖子

0

威望

2345

eV
积分
2530

Level 5 (御坂)

3#
发表于 Post on 2022-11-18 11:09:26 | 只看该作者 Only view this author
ISPIN没有设置呐,U值测试过吗,一般看到的多的Fe是3.8-4.0eV的样子。

8

帖子

0

威望

91

eV
积分
99

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2022-11-18 11:25:10 | 只看该作者 Only view this author
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)

208

帖子

3

威望

1998

eV
积分
2266

Level 5 (御坂)

5#
发表于 Post on 2022-11-18 15:47:57 | 只看该作者 Only view this author
确定原子数,种类一致的话再检查基组电荷是不是都一致,这么大的值一般不太可能是构型不同引起的

8

帖子

0

威望

91

eV
积分
99

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2022-11-18 16:50:01 | 只看该作者 Only view this author
Eudaimonia 发表于 2022-11-18 15:47
确定原子数,种类一致的话再检查基组电荷是不是都一致,这么大的值一般不太可能是构型不同引起的

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

255

帖子

1

威望

2666

eV
积分
2941

Level 5 (御坂)

7#
发表于 Post on 2022-11-19 11:35:31 | 只看该作者 Only view this author
你这个严格意义上讲不是吸附能,而是发生了磷酸和表面的脱水反应,在最终形成的结构中,不光产生了新的化学键,还产生了几个内氢键,这都会导致体系能量降低,体现在所谓的吸附能中。更确切来说,现在算出的能量是真空条件下磷酸分子与羟基氧化铁表面脱水反应的能量降低值。
我唯一知道的就是我一无所知,但我是化学小迷弟

8

帖子

0

威望

91

eV
积分
99

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2022-11-19 19:12:29 | 只看该作者 Only view this author
ChemG 发表于 2022-11-19 11:35
你这个严格意义上讲不是吸附能,而是发生了磷酸和表面的脱水反应,在最终形成的结构中,不光产生了新的化学 ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-17 01:46 , Processed in 0.173925 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list