计算化学公社

标题: 使用cp2k进行结构优化能量输出问题 [打印本页]

作者
Author: yjb 时间: 2022-5-19 13:22
标题: 使用cp2k进行结构优化能量输出问题
使用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

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

作者
Author: elpa 时间: 2022-5-19 14:51
输出坐标文件的最后两帧结构好像一模一样

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

作者
Author: yjb 时间: 2022-5-19 18:02

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

好的感谢

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

作者
Author: yjb 时间: 2022-5-19 22:30

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

好的谢谢

作者
Author: ma455173220 时间: 2022-5-26 11:02
我个人觉得这类似于VASP的单点计算，应该是为了在优化的结构基础上进一步算能量，所以最终能量要取用“Geometry optimization completed”下面给出的能量。

作者
Author: yjb 时间: 2022-5-26 13:04

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

好的，谢谢

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