求助，cp2k计算MgO形成能问题

Lzzaya · 发表于 Post on 2025-6-7 15:47:14

如题，本人cp2k新手，最近准备先使用cp2k计算形成能练手。但是在计算MgO的形成能时结果始终无法和MP数据库中对应，cp2k计算结果为-5.7 eV/atom，计算结构获取自material project的mp-549706，其报道数值为-2.9 eV/atom，同时附上我所使用的.inp文件。请问我还需要怎么调整.inp文件？

Mg.inp (3.01 KB, 下载次数 Times of downloads: 7)

O.inp (2.69 KB, 下载次数 Times of downloads: 5)

static.inp (3.15 KB, 下载次数 Times of downloads: 10)

sobereva · 发表于 Post on 2025-6-8 01:06:05

每个计算的输出文件都上传，把计算公式列出来，让别人知道你是怎么得到的-5.7 eV/atom

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

Lzzaya · 发表于 Post on 2025-6-8 08:59:12

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

当前有显而易见的 ...

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

Mg.out (24.85 KB, 下载次数 Times of downloads: 0)

O.out (35.3 KB, 下载次数 Times of downloads: 1)

static.out (24.88 KB, 下载次数 Times of downloads: 1)

我使用的计算公式是：

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

sobereva · 发表于 Post on 2025-6-9 16:58:41

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。

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

[CP2K] 求助，cp2k计算MgO形成能问题

浏览过的版块