计算化学公社

标题: wfx文件里的EDF是否影响能量密度？ [打印本页]

作者
Author: beefly 时间: 2017-6-2 11:52
标题: wfx文件里的EDF是否影响能量密度？
做了个测试，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不重要呢？

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

作者
Author: beefly 时间: 2017-6-2 20:43
明白了。谢谢！

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