求助：GLOBAL＞PRINT_LEVEL为HIGH时，输出信息中关于原子力组分的若干问题

蟹黄宝

在cp2k中将　GLOBAL＞PRINT_LEVEL设置为HIGH时，输出信息会打印每个原子的多个力组分，以BLYP-D3(BJ)为例，力组分信息如下：
｀｀｀
FORCES [a.u.]


  Atom    Kind        Component        X           Y           Z


     1       1          overlap     -0.102528   -0.109573    0.017021
     1       1     overlap_admm      0.000000    0.000000    0.000000
     1       1          kinetic     -0.214914   -0.223289    0.034852
     1       1          gth_ppl      0.034635    0.035671   -0.005697
     1       1         gth_nlcc      0.000000    0.000000    0.000000
     1       1         gth_ppnl     -0.067610   -0.067784    0.010520
     1       1     core_overlap     -0.000000   -0.000000    0.000000
     1       1         rho_core     -0.374422   -0.388901    0.047493
     1       1         rho_elec      0.720321    0.742719   -0.100837
     1       1     rho_lri_elec      0.000000    0.000000    0.000000
     1       1         ch_pulay      0.000000    0.000000    0.000000
     1       1       dispersion     -0.000250   -0.000123   -0.000138
     1       1              gCP      0.000000    0.000000    0.000000
     1       1            other      0.000000    0.000000    0.000000
     1       1          fock_4c      0.000000    0.000000    0.000000
     1       1        ehrenfest      0.000000    0.000000    0.000000
     1       1           efield      0.000000    0.000000    0.000000
     1       1              eev      0.000000    0.000000    0.000000
     1       1      mp2_non_sep      0.000000    0.000000    0.000000
     1       1            total     -0.004768   -0.011281    0.003216
｀｀｀

我有以下四个问题，请各位老师和前辈解答：
１．请问这些力组分分别是什么含义？
２．请问这些力组分信息是如何计算得到的？
3.　请问这些力组分是否与输出文件中的能量项（Core Hamiltonian energy，Hartree energy，Exchange-correlation energy等）存在对应关系？
４．我观察到上面`total`部分得到的结果，与后续｀ATOMIC FORCES in [a.u.]｀开头的结果正好输出值相反，请问这是为什么？
完整的输入文件和输出log参见附件

sobereva

这种事几句话说不清楚，需要看quickstep原文Computer Physics Communications 167 (2005) 103–128、GAPW原文Theor Chem Acc (1999) 103:124±140，搞懂了算法细节自然明白怎么回事