计算化学公社

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

作者
Author: Q157 时间: 2025-1-29 08:58
标题: 求助使用CP2K计算CeO2（111）切面晶胞单点能量与VASP计算所得能量差距过大的原因
我是第一性原理计算的初学者，组内是用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，请问怎么会有这么大的差距，会是哪里设置的错误？
请各位大佬解惑，祝大家新年快乐，万事如意！
(, 下载次数 Times of downloads: 16)

(, 下载次数 Times of downloads: 10) (, 下载次数 Times of downloads: 9)

作者
Author: yxdd98 时间: 2025-1-29 11:10
对比绝对能量无意义，不同程序的默认设置都不同

作者
Author: sobereva 时间: 2025-1-29 16:11
正确拼写Hartree

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

作者
Author: Q157 时间: 2025-1-29 22:29

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

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

作者
Author: Q157 时间: 2025-1-29 22:29

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

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

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

作者
Author: Jerry_Fan 时间: 2025-4-13 19:50

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

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

作者
Author: Q157 时间: 2025-5-19 12:49

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

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

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