bug report

beefly

版本3.3.9。如下wfx文件，做电子密度拓扑分析。某临界点的G，K，V，H值：


Lagrangian kinetic energy G(r):  0.3771912731
Hamiltonian kinetic energy K(r):  0.2080016661
Potential energy density V(r): -0.5770091772
Energy density E(r) or H(r): -0.2080016661


发现H 不等于G + V，而是K直接取负值。换其他分子的wfx文件做类似的分析，未发现此错误。

sobereva

谢谢报告。
导致问题的原因是计算V的时候是根据V=(1/4)der2rho-2G的关系给出的，但是EDF字段只对der2rho有贡献而对G没有贡献。当不读取EDF字段的时候就没这个问题。
等过几天我会在官网上上传新的3.4(dev)版，把V的计算方式改成了V=-K-G，就确保H=V+G=-K了