计算化学公社

标题: 求助，cp2k计算MgO形成能问题 [打印本页]

作者
Author: Lzzaya 时间: 2025-6-7 15:47
标题: 求助，cp2k计算MgO形成能问题
如题，本人cp2k新手，最近准备先使用cp2k计算形成能练手。但是在计算MgO的形成能时结果始终无法和MP数据库中对应，cp2k计算结果为-5.7 eV/atom，计算结构获取自material project的mp-549706，其报道数值为-2.9 eV/atom，同时附上我所使用的.inp文件。请问我还需要怎么调整.inp文件？
(, 下载次数 Times of downloads: 7)
(, 下载次数 Times of downloads: 5)
(, 下载次数 Times of downloads: 10)

作者
Author: sobereva 时间: 2025-6-8 01:06
每个计算的输出文件都上传，把计算公式列出来，让别人知道你是怎么得到的-5.7 eV/atom

当前有显而易见的错误，氧气分子的基态是三重态而非单重态。而且孤立体系用什么k点，完全莫名其妙。D3完全没必要加

作者
Author: Lzzaya 时间: 2025-6-8 08:59

sobereva 发表于 2025-6-8 01:06
每个计算的输出文件都上传，把计算公式列出来，让别人知道你是怎么得到的-5.7 eV/atom

当前有显而易见的 ...

谢谢老师，因为我本身是计算机专业为了做力场想用cp2k算点数据，所以参数设置上有不专业的地方。
以下是我分别使用.inp文件计算得到的输出结果文件：
(, 下载次数 Times of downloads: 0)
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 1)

我使用的计算公式是：
(, 下载次数 Times of downloads: 13)

所有能量从结果文件的"ENERGY| Total FORCE_EVAL ( QS ) energy [hartree]"这一行提取，通过单位转换之后计算得到，总计算过程为(单位转换后)：
(-4324.5066 - (-1717.0121*2 + -433.8305*2) )/4 = -5.70535 eV/atom

作者
Author: sobereva 时间: 2025-6-9 16:58

Lzzaya 发表于 2025-6-8 08:59
谢谢老师，因为我本身是计算机专业为了做力场想用cp2k算点数据，所以参数设置上有不专业的地方。
以下是 ...

先把我说的自旋多重度问题解决了再说

另外，要弄清楚算的是什么。形成能一般是以每mol的formula unit (f.u.)给出的，比如MgO的生成焓是-601.6 kJ/mol（https://webbook.nist.gov/cgi/cbo ... its=SI&Mask=7EF），约为-6.23 eV/mol。

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