更新补充一点我自己的理解。

先前一直困扰我的一个问题是，从表达式上看，如果E(r)=-K(r)，那岂不是说，体系的能量就简单地等于动能的相反数？
近期重新看Multiwfn的手册，发现E(r)=-K(r)的前提是：“Virial field”，也就是满足维里定理；

根据维里定理，当原子核受力为0时，若波函数是精确的，则势能和动能的绝对值之比（-V/T）为2。
于是，在维里场内，有E=V+T=-2T+T=-T；
但实际情况下的维里比值较之2有一定的偏差，当考虑高斯实际输出的（-V/T）的值，确实可以仅根据高斯输出文件中的动能值算出体系总能量。

附CHEMICAL MODELS AND THE LAPLACIAN一书中的相关内容：
https://chempedia.info/page/1751 ... 230002187148068252/

sobereva 发表于 2019-2-20 17:55
1 对K(r)和G(r)全空间积分的结果完全一样，和量化程序直接输出的电子动能完全相同。微小不同是来自数值积分 ...

明白了，非常感谢Sob老师！！！

1 对K(r)和G(r)全空间积分的结果完全一样，和量化程序直接输出的电子动能完全相同。微小不同是来自数值积分精度问题，增加积分格点数目就可以解决。
动能密度的定义有极大任意性，只要满足全空间积分等于实际以电子动能即可。然而电子密度拉普拉斯函数全空间积分为0，因此对动能密度可以任意加上一个 [常数]*电子密度拉普拉斯函数 的项。K(r)和G(r)也正是相差1/4*电子密度拉普拉斯函数。

2 能。就是Gaussian输出的电子能量（不包含核互斥能）。更具体来说，就是Gaussian输出的单点能减去输出的核互斥能。