使用cp2k进行结构优化能量输出问题

yjb · 发表于 Post on 2022-5-19 13:22:55

使用cp2k对一个结构进行优化，Convergence check：四个收敛指标都YES。
图中 ENERGY| Total FORCE_EVAL ( QS ) energy [a.u.]: -132.146373010317205

这个数据是不是优化完结构的整体能量。

为什么在GEOMETRY OPTIMIZATION COMPLETED之后还会再算一步，又得到一个

ENERGY| Total FORCE_EVAL ( QS ) energy [a.u.]: -132.146373072488927

最终的能量应该取哪一个，谢谢各位老师。

elpa · 发表于 Post on 2022-5-19 14:51:20

输出坐标文件的最后两帧结构好像一模一样

丁越 · 发表于 Post on 2022-5-19 16:32:52

读取最后一个能量就行，这个和QE类似，做完结构优化还会最后做一次单点任务，能量从最后一步读取。

yjb · 发表于 Post on 2022-5-19 18:02:09

丁越发表于 2022-5-19 16:32
读取最后一个能量就行，这个和QE类似，做完结构优化还会最后做一次单点任务，能量从最后一步读取。

好的感谢

ChemG · 发表于 Post on 2022-5-19 22:07:41

在密度泛函理论（中文版）的第79页脚注中，对于晶胞体积变化的计算可能会存在Pulay stress，因此需要多做一步单点消除掉，我一直是这么理解的，不知是否有误。

yjb · 发表于 Post on 2022-5-19 22:30:11

ChemG 发表于 2022-5-19 22:07
在密度泛函理论（中文版）的第79页脚注中，对于晶胞体积变化的计算可能会存在Pulay stress，因此需要多做一 ...

好的谢谢

ma455173220 · 发表于 Post on 2022-5-26 11:02:53

我个人觉得这类似于VASP的单点计算，应该是为了在优化的结构基础上进一步算能量，所以最终能量要取用“Geometry optimization completed”下面给出的能量。

yjb · 发表于 Post on 2022-5-26 13:04:53

ma455173220 发表于 2022-5-26 11:02
我个人觉得这类似于VASP的单点计算，应该是为了在优化的结构基础上进一步算能量，所以最终能量要取用“Geom ...

好的，谢谢

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

[CP2K] 使用cp2k进行结构优化能量输出问题

浏览过的版块