计算化学公社

标题: CP2K的SCF为什么在出现了更低的能量值后，收敛在了更高的能量值 [打印本页]

作者
Author: flyingchow 时间: 2023-2-16 01:03
标题: CP2K的SCF为什么在出现了更低的能量值后，收敛在了更高的能量值
我有一个问题关于CP2k的问题一直不明白，就是用CP2k跑一个scf的运算（有的时候的几何优化的一步）的时候，经常会出现在循环的过程中跑出一个低的能量值，但是最终收敛的能量值却高于这个值。
这个难道不和变分原理相违背吗？我用Quantum Espresso的时候从来也没有出现这个情况。

有老师能讲讲这是啥原因吗？谢谢。

作者
Author: jiangning198511 时间: 2023-2-16 08:58
在VASP里面也有这样的现象，其原因可能是使用了density-mixing, 就是迭代m部后，m+1部的密度是通过前m的线性组合得到的，这样就可能导致能量下降的过低

作者
Author: 北大-陶豫 时间: 2023-2-16 09:21

jiangning198511 发表于 2023-2-16 08:58
在VASP里面也有这样的现象，其原因可能是使用了density-mixing, 就是迭代m部后，m+1部的密度是通过前m的线 ...

所以为啥会出现比正解能量还低的密度，这和HK变分定理相违背吧？

作者
Author: sobereva 时间: 2023-2-16 10:10

北大-陶豫发表于 2023-2-16 09:21
所以为啥会出现比正解能量还低的密度，这和HK变分定理相违背吧？

HK第二定理的前提是电子密度满足N可表。不满足这条的虚假的密度可以违背HK第二定理。

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