vasp 计算双原子分子电子能量不准确吗

wzkchem5 · 发表于 Post on 2024-3-12 18:14:07

dingniu2 发表于 2024-3-12 10:18
请问VASP算不准O2与官网中“https://www.vasp.at/wiki/index.php/Available_PAW_potentials”提到的“Imp ...

泛函的SIE和赝势的误差是两个完全不同的因素。我不是太了解这里赝势的误差有多大，但我估计泛函的误差是主导的，因为反推O2能量需要的H2O、H2都有短键，如果赝势误差占主导，那么这种反推方法并不解决问题

13314875806 · 发表于 Post on 2024-3-13 19:59:47

wzkchem5 发表于 2024-3-12 16:52
差这么多的话，不是计算不准确的问题，而应该是计算哪里搞错了，例如该减的什么数没有减，或者反应没有配 ...

谢谢老师，非常感谢，我重新计算了一下，然后算出来画图，OOH,OH，O自由能数值为负数，老师请问这个跟什么有关，是我的结构不对吗？

wzkchem5 · 发表于 Post on 2024-3-13 20:16:08

13314875806 发表于 2024-3-13 12:59
谢谢老师，非常感谢，我重新计算了一下，然后算出来画图，OOH,OH，O自由能数值为负数，老师请问这个跟什 ...

OOH、OH、O的什么自由能是负数？吸附自由能？还是绝对自由能？如果是绝对自由能，是吸附物种的绝对自由能还是没有吸附的物种的绝对自由能？
你为什么认为这个自由能不应该是负数？只是因为个别几篇文献不是负数，还是有深层原因？

13314875806 · 发表于 Post on 2024-3-13 23:06:59

wzkchem5 发表于 2024-3-13 20:16
OOH、OH、O的什么自由能是负数？吸附自由能？还是绝对自由能？如果是绝对自由能，是吸附物种的绝对自由能 ...

对，我先找不同的原因，谢谢

Stars · 发表于 Post on 2024-6-19 14:54:43

wzkchem5 发表于 2024-3-12 18:14
泛函的SIE和赝势的误差是两个完全不同的因素。我不是太了解这里赝势的误差有多大，但我估计泛函的误差是 ...

你好，针对这个情况，我进行了测试，结果如下：
按照一般的不带后缀的赝势计算ENCUT为400，硬一点的O_h为700，模拟一般研究者的ENCUT卡着ENMAX选取的情况，能量分别为-9.884eV和-10.005eV
两者都取极高的ENCUT（1000eV），检测这两个赝势，能量分别为-9.900eV和-10.053eV。
需要注意的是，这两个赝势对绝对能量计算显然不可比。为了验证选用赝势对绝对能量的影响，我们还看了H2和H2O的能量，针对普通的（400）和带h的（700），结果分别是-6.759/-6.764和-14.221/-14.258，所以说对于一般的分子，不带_h的赝势一般误差也就在0.05eV的量级甚至更小，相当于O2这种的1/3不到，但对于O2，可以超过0.1eV。

wzkchem5 · 发表于 Post on 2024-6-19 17:23:10

Stars 发表于 2024-6-19 07:54
你好，针对这个情况，我进行了测试，结果如下：
按照一般的不带后缀的赝势计算ENCUT为400，硬一点的O_h ...

这么看来可能是多重键而非短键的影响。形成多重键的分子，电子密度偏离自由原子的程度比较大，而赝势必然是适合描述自由原子的，也往往是主要以自由原子数据为基础拟合的，所以可能导致不准。

Stars · 发表于 Post on 2024-6-19 17:50:06

wzkchem5 发表于 2024-6-19 17:23
这么看来可能是多重键而非短键的影响。形成多重键的分子，电子密度偏离自由原子的程度比较大，而赝势必然 ...

我觉得暂时还不能下定论，电子密度偏离自由原子的体系还有强离子性的体系，或许可以研究一下这样的体系。

小叶 · 发表于 Post on 2025-4-23 19:22:29

那计算过氧化氢的电子能量准确嘛，我想计算生成过氧化氢的能量，我计算出来是2.66eV,但是很多文献是1.4eV,感觉差别挺大的

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] vasp 计算双原子分子电子能量不准确吗