wfx文件里的EDF是否影响能量密度？

beefly

做了个测试，Tl+阳离子，ECP有78个电子，计算产生wfx文件。然后用MultiWFN在坐标为(0,0,0.6)的点（原子单位）计算电子密度及各种性质

Density of all electrons:  0.9959744042E+01
Laplacian of electron density:  0.2677268998E+03
Lagrangian kinetic energy G(r):  0.1686141219E-02
Energy density E(r) or H(r):  0.1167802834E-02


删除EDF字段后，结果为

Density of all electrons:  0.2013599868E-04
Laplacian of electron density:  0.1141577621E-01
Lagrangian kinetic energy G(r):  0.1686141219E-02
Energy density E(r) or H(r):  0.1167802834E-02


可以看到，电子密度、Laplacian有显著变化，而G、H不变。这是程序bug，还是说EDF对G、H不重要呢？

sobereva

计算能量密度需要轨道波函数，所以EDF字段不影响能量密度
Multiwfn里EDF只影响电子密度、梯度、RDG、sign(lambda2)rho、电子密度拉普拉斯、电子密度Hessian等纯粹基于密度的。虽然静电势也纯粹基于密度，但目前EDF还不影响静电势结果。

beefly

明白了。谢谢！