求助使用CP2K计算CeO2（111）切面晶胞单点能量与VASP计算所得能量差距过大的原因

Q157 · 发表于 Post on 2025-1-29 08:58:48

我是第一性原理计算的初学者，组内是用VASP的，我用VASP计算CeO2（111）切面晶胞的能量，使用DFT+U所得结果为-1148.5279 eV，INCAR和cif晶胞文件如下：
最近我在学习CP2K程序包的使用，对相同的cif文件，使用DFT+U，也考虑了色散校正，使用基组Ce DZV-MOLOPT-SR-GTH，赝势为GTH-PBE-q12，DFT+U中U-J设置为4.08，
具体输入文件在附件中，计算所得结果却可以达到-3404.3868071 hatree，换算后为-91918.4436 eV，请问怎么会有这么大的差距，会是哪里设置的错误？
请各位大佬解惑，祝大家新年快乐，万事如意！

CeO2.inp (12.01 KB, 下载次数 Times of downloads: 16)

yxdd98 · 发表于 Post on 2025-1-29 11:10:13

对比绝对能量无意义，不同程序的默认设置都不同

sobereva · 发表于 Post on 2025-1-29 16:11:46

正确拼写Hartree

基组、赝势、数值算法，特别是赝化的电子数，都有天壤之别，电子能量毫无可比性。只有计算出来的有化学意义的数据才可能能比，诸如功函数、原子电荷、势垒、反应能、吸附能，等等

Q157 · 发表于 Post on 2025-1-29 22:29:01

yxdd98 发表于 2025-1-29 11:10
对比绝对能量无意义，不同程序的默认设置都不同

好的，感谢解惑，我后面尝试计算吸附能看看

Q157 · 发表于 Post on 2025-1-29 22:29:53

sobereva 发表于 2025-1-29 16:11
正确拼写Hartree

基组、赝势、数值算法，特别是赝化的电子数，都有天壤之别，电子能量毫无可比性。只有 ...

感谢sob老师回复，以后会注意Hartree拼写，我准备计算吸附能比较看看

Jerry_Fan · 发表于 Post on 2025-4-13 19:50:40

Q157 发表于 2025-1-29 22:29
感谢sob老师回复，以后会注意Hartree拼写，我准备计算吸附能比较看看

你好，请问后续吸附能计算结果怎么样？我也是觉得cp2k计算的能量太大了，和lammps根本没有办法对比（上虽然确实对比绝对值不太应该）。对比起来，vasp计算的结果似乎还是能够在一个量级上比较的

Q157 · 发表于 Post on 2025-5-19 12:49:40

Jerry_Fan 发表于 2025-4-13 19:50
你好，请问后续吸附能计算结果怎么样？我也是觉得cp2k计算的能量太大了，和lammps根本没有办法对比（上虽 ...

你好，时间过去有点久了，我也挺久没继续算了，印象里计算吸附能再将Hartree单位换算为eV后，差别同VASP相比确实比较大，有需要的话可以参考看看有没有使用cp2k计算吸附能的文献

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

[CP2K] 求助使用CP2K计算CeO2（111）切面晶胞单点能量与VASP计算所得能量差距过大的原因

浏览过的版块