计算化学公社

标题: bug report [打印本页]

作者
Author: beefly 时间: 2017-3-16 08:16
标题: bug report
版本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文件做类似的分析，未发现此错误。

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

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