计算化学公社

标题: HER论文计算部分的问题 [打印本页]

作者
Author: ilovelr 时间: 2022-12-1 15:43
标题: HER论文计算部分的问题
本帖最后由 ilovelr 于 2022-12-2 10:57 编辑

最近受隔壁组委托算合金体系的HER过程，发现建模很难达到他们预期想要的结果，但是很多文章里的VASP计算结果却很漂亮的与实验对上，于是我也看了很多文献，发现里面有很大的问题：
比如这篇文章https://doi.org/10.1016/j.jechem.2020.08.040
(, 下载次数 Times of downloads: 13)

文章计算了H在金属111表面的吸附自由能（文中说用PBE），根据文章示意图，4x4晶胞，1/16ML的覆盖率。
首先，这里Pt的吸附自由能大约是-0.09eV，很明显这是使用RPBE计算的结果（https://doi.org/10.1016/j.apsusc.2021.149627），而PBE计算的H吸附在Pt上的自由能应该是-0.2eV左右。
其次，图中NiCo的自由能是-0.58eV左右，而这个值经过测试，应该是用PBE算的，但是这个数值很明显是E(H*)而不是G(H*)。所以，我猜测应该是为了凑出Mo,W doped CoNi自由能在Pt和NiCo之间，而故意不给出NiCo的上的G(H*)，因为如果算自由能，NiCo应该是在W doped之上的。

再比如这篇文章https://doi.org/10.1016/j.ijhydene.2020.04.242,作者研究了MxNi1-x (M = Cr, Mo and W; x ≈ 0.2) 合金上的HER反应，在计算中使用了M1Ni3合金作为模型。
(, 下载次数 Times of downloads: 12)

这篇文章用了RPBE计算，而这里就有两个问题，一是为了凑出Ni上的自由能最低，只用E(H*)，按照作者的模型，Ni111上的G(H*)应该在-0.16 eV左右。
其二，有计算的文献中显示(https://doi.org/10.1016/j.apsusc.2021.149627),W1Ni3会有2个fcc点位，一个是-0.20eV, 一个是-0.04eV, 这里作者应该是采用了-0.04eV的点位的数据以支持他的实验。

综上，这些只是我看到文献中的一些问题，我说的不一定对，欢迎批评指正。
我也是初学者，通过看这些文献，尤其是实验里面的计算，同样的模型不同的文章却相反的结果，给我造成很大困惑。虽然实验论文的计算并不重要，但这样肆无忌惮的选数据真的让人很无语

作者
Author: hpq 时间: 2022-12-1 17:27
给实验算的吗随缘就行别较真

作者
Author: ilovelr 时间: 2022-12-1 22:12

hpq 发表于 2022-12-1 17:27
给实验算的吗随缘就行别较真

随缘算出来可对不上，只能选数据

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