计算化学公社

标题: vasp关于计算真空盒子中一个单分子的能量时，单分子总是解体的解决办法 [打印本页]

作者
Author: LZN 时间: 2024-5-7 20:43
标题: vasp关于计算真空盒子中一个单分子的能量时，单分子总是解体的解决办法
各位老师好，我用vasp计算NH在一个15*15*15的盒子中的能量时，结算结果总是出现NH解体，想请教一下怎么解决，我在算NH2、NH3时没有出现这个问题。下图是我的结果图和INCAR,Kpoints参数。

作者
Author: nullzy 时间: 2024-5-7 22:46
为啥算NH的能量，没意义吧，你要算NH吸附能吗，不是这样算的。

作者
Author: LZN 时间: 2024-5-8 15:04

nullzy 发表于 2024-5-7 22:46
为啥算NH的能量，没意义吧，你要算NH吸附能吗，不是这样算的。

你好老师，我是想算NH3每一步脱H的吸附能，我想着把整体、基底、吸附分子的能量算出来，这样不对吗，请问正确应该是怎么算呀

作者
Author: 白皛白 时间: 2024-5-9 10:56
本帖最后由白皛白于 2024-5-9 10:59 编辑

LZN 发表于 2024-5-8 15:04
你好老师，我是想算NH3每一步脱H的吸附能，我想着把整体、基底、吸附分子的能量算出来，这样不对吗，请问 ...

NHn物种的吸附自由能感觉可以参照OH的吸附自由能来算，类比这样的形式ΔGOH=GOH*-G*-(GH2O-1/2GH2)，替换成ΔGNH3=GNH2*-G*-(GNH3-1/2GH2)这样。
不过对于NOR或者NO3RR的研究来说，似乎不太会像ORR和OER这样计算每个中间体的吸附能，我看到的文章里都只算每一步的自由能变ΔG，类似NH3*→NH2*+(质子+电子)这一步，它的自由能变就参照下式ΔG=GNH2*-GNH3*+1/2GH2来处理了，在这一部分只需要计算NH3与H2的吉布斯自由能，而不需要计算NH2，NH，N这些的吉布斯自由能，实际上应该也没办法这样算。
也可以参考下https://onlinelibrary.wiley.com/doi/10.1002/adfm.202008533这篇文章的supporting information

作者
Author: LZN 时间: 2024-5-9 15:51

白皛白发表于 2024-5-9 10:56
NHn物种的吸附自由能感觉可以参照OH的吸附自由能来算，类比这样的形式ΔGOH=GOH*-G*-(GH2O-1/2GH2)，替换 ...

我明白了，谢谢老师

作者
Author: 卡开发发 时间: 2024-5-9 16:48

LZN 发表于 2024-5-8 15:04
你好老师，我是想算NH3每一步脱H的吸附能，我想着把整体、基底、吸附分子的能量算出来，这样不对吗，请问 ...

4#对如何计算每一步自由能变如何计算已经进行比较完整的说明，此处便不再赘述。需要补充的是，即便真的需要NH或NH2在气相下的能量，这里计算失败的原因有如下可能：
1、MAGMOM不太合理，对NH来说，可以N猜2，H猜0，可以从这个分子可能的电子排布进行猜测。
2、对称性应当被关闭，或者要把分子摆成不太对称的朝向，否则有可能程序会按照晶体空间群识别得到更高的对称性而得到不合理结果。
3、精度太低，ENCUT设置到350对N元素有可能会不够，至少应该测试一下。
4、初始结构是否合理，每一步的自洽场是否都已经达到收敛，这些问题暂时无从考证。

作者
Author: LZN 时间: 2024-5-10 09:23

卡开发发发表于 2024-5-9 16:48
4#对如何计算每一步自由能变如何计算已经进行比较完整的说明，此处便不再赘述。需要补充的是，即便真的需 ...

谢谢老师

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