计算化学公社

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

[VASP] RuO2(110)OER计算结果与文献差异较大

[复制链接 Copy URL]

114

帖子

0

威望

865

eV
积分
979

Level 4 (黑子)

参考《Breaking the Activity-Stability Trade-Off of RuO2 via Metallic Ru Bilateral Regulation for Acidic Oxygen Evolution Reaction》文献支撑信息部分,搭建RuO2(110)计算酸性情况下OER过程,超胞大小p(2*2),输入文献如下所示。当K点111时,△G*OH=0.52 eV;当K点421时,△G*OH=-0.11 eV。两者均与文献结构差异较大(△G*OH=0.8 eV),想请问是哪里出现了问题


文献计算细节为:

The DFTcalculations were performed using the projector-augmented wave (PAW) method, asimplemented in the Vienna Ab initio Simulation Package (VASP).48, 49 The Perdew–Burke–Ernzerhof (PBE) exchange-correlation functional withinthe generalized gradient approximation (GGA) was employed. For the optimization of the bulkRuO2 (space group: P42/mnm), a uniform 6×6×8 k-meshgrid was used for sampling the Brillouin zone. The RuO2slab model consisted of fivelayers with (110) surface exposed. A vacuum layer of ~16 Å wasintroduced along the z-direction to prevent inter-slab interactions. The unit cell of Ru adopts a hexagonal close-packed (hcp) structure with P63/mmcsymmetry, while RuO2 belongs to the P42/mnmsymmetry group. The Ru/RuO2 interfacemodel was constructed by combining Ru (100) and RuO2 (110),[50] as these facets are commonlyused due to their well-defined atomic arrangement, moderate surface energy, andstructural resemblance. The resulting model exhibited a lattice mismatch errorof less than 5%, ensuring a stable and physically meaningful interface. This interfacial model had fivelayers, with the slabs separated by ~16 Å of the vacuum layer. A 4×2×1 Monkhorst-Pack k-pointgrid was used for Brillouin zone sampling. The kinetic energy cutoff for thewavefunctions was set at 450 eV for the bulk and slab structures without andwith molecular adsorption. The geometry optimizations were consideredconvergent when the force on each atom was less than 0.02 eV/Å.


数据计算过程.png (37.89 KB, 下载次数 Times of downloads: 5)

数据计算过程.png

INCAR

346 Bytes, 下载次数 Times of downloads: 0

CONTCAR

15.34 KB, 下载次数 Times of downloads: 0

K111时OH吸附

CONTCAR

15.09 KB, 下载次数 Times of downloads: 0

K111时slab

CONTCAR

15.09 KB, 下载次数 Times of downloads: 0

K421时slab

CONTCAR

15.34 KB, 下载次数 Times of downloads: 0

K421时OH吸附

44

帖子

0

威望

653

eV
积分
697

Level 4 (黑子)

2#
发表于 Post on 前天 10:52 | 只看该作者 Only view this author
它文中只是提及了基础的模型数据和最基础的计算数据。你像是否打开溶剂化,磁矩设置是多少,是否用 DFT+U ,这里面都没说,这下都会对结果有比较大的影响。同时,INCAR 中 AMIX BMIX 这些计算参数也会有影响。
顺便,我计算的 110 面的 OH 吸附(开溶剂化,用 DFT+U)的结果是,0.325

114

帖子

0

威望

865

eV
积分
979

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on yesterday 10:44 | 只看该作者 Only view this author
shinnashiori 发表于 2026-3-10 10:52
它文中只是提及了基础的模型数据和最基础的计算数据。你像是否打开溶剂化,磁矩设置是多少,是否用 DFT+U  ...

1、参考《Eliminating over-oxidation of ruthenium oxides by niobium for highly stable electrocatalytic oxygen evolution in acidic media》和《Breaking the Activity-Stability Trade-Off of RuO2 via Metallic Ru Bilateral Regulation for Acidic Oxygen Evolution Reaction》,△G*OH约为0.8 eV。

2、参考《Dual-Descriptor-Guided Screening of Stable Metal-Doped RuO2 Catalysts for Acidic Oxygen Evolution》(U−J) value is 2.4。

44

帖子

0

威望

653

eV
积分
697

Level 4 (黑子)

4#
发表于 Post on yesterday 16:58 | 只看该作者 Only view this author
飞翔的猪 发表于 2026-3-11 10:44
1、参考《Eliminating over-oxidation of ruthenium oxides by niobium for highly stable electrocataly ...

额,你参考的参数是来自多篇文献的?然后想复现其中一篇的结果?你想复现的话肯定要用复现的那篇文章的计算参数去做计算啊。
如果计算参数都一致那可以考虑是否是收敛到能量最低情况(结构优化可能没有到最好的情况)

本版积分规则 Credits rule

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

GMT+8, 2026-3-12 20:53 , Processed in 0.171580 second(s), 23 queries , Gzip On.

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