计算化学公社

标题: 求助-VASP计算Ni 111表面吸附H的能量-可有偿咨询 [打印本页]

作者
Author: 爱上飞机的坦克 时间: 2023-8-30 09:03
标题: 求助-VASP计算Ni 111表面吸附H的能量-可有偿咨询
各位老师好，我用vasp计算Ni 111表面吸附H，能量经过了ZPE和TS的校正后。Ni顶位的△GH*是0.08eV （Ni-H键长1.46），桥位、hcp和fcc位的△GH*是-0.42eV左右（Ni-H键长1.7）。我看文献中Ni 111面对H吸附的△GH*在-0.30~-0.39 eV之间（Ni-H键长1.93）。附上INCAR，POSCAR和CONTCAR，请各位老师帮忙看看问题在哪里？一直卡在这个问题上，如果解决比价麻烦，可有偿咨询，谢谢！

作者
Author: 含光君 时间: 2023-8-30 13:01
如果确定建模、所用赝势、收敛标准、k点精度、截断能等都与原文无差别，应再检查你输入文件中的LREAL、是否需加溶剂模型，以及确认文献作者是否带着偶极校正优化（而不是直接优化后再加偶极校正算单点）

作者
Author: Stars 时间: 2023-12-27 22:05
你好，我想请问下，你说的计算吸附自由能为-0.39的，是哪篇？作者是不是用的PBE而非RPBE计算的？

作者
Author: 卡开发发 时间: 2023-12-28 01:51
看看原文用没用色散修正，有可能问题在色散修正上。

作者
Author: Stars 时间: 2023-12-28 10:23

卡开发发发表于 2023-12-28 01:51
看看原文用没用色散修正，有可能问题在色散修正上。

我记得在Norskov做表面吸附计算的年代，DFT-D的方法甚至还没提出来。他的结果是-0.30eV左右。问题是，即使是计算化学吸附，现在有了D3，也没理由不加吧？

作者
Author: gog 时间: 2023-12-28 10:43
吸附的是质子氢还是氢原子？

作者
Author: 卡开发发 时间: 2023-12-28 15:48

Stars 发表于 2023-12-28 10:23
我记得在Norskov做表面吸附计算的年代，DFT-D的方法甚至还没提出来。他的结果是-0.30eV左右。问题是，即 ...

D3对金属的sigma键有可能高估，这种情况PBE不加DFT-Dn可能都高估吸附能，如果RPBE不加色散修正反而不会，相关文献你可以找找。

作者
Author: 卡开发发 时间: 2023-12-28 15:59

Stars 发表于 2023-12-28 10:23
我记得在Norskov做表面吸附计算的年代，DFT-D的方法甚至还没提出来。他的结果是-0.30eV左右。问题是，即 ...

数据我可以给一个（该表格被裁剪过，原本包含更多），但是调研文献应该大家自己来。
(, 下载次数 Times of downloads: 14)

作者
Author: Stars 时间: 2023-12-28 18:05

卡开发发发表于 2023-12-28 15:59
数据我可以给一个（该表格被裁剪过，原本包含更多），但是调研文献应该大家自己来。

感谢分享，我找到了一篇使用DFT-D3计算吸附的数据，确实会高很多，下面分享上来，lz也可以看看
https://doi.org/10.1016/j.ijhydene.2023.10.032

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