计算化学公社

标题: 求助：GLOBAL＞PRINT_LEVEL为HIGH时，输出信息中关于原子力组分的若干问题 [打印本页]

作者
Author: 蟹黄宝 时间: 2024-5-4 23:54
标题: 求助：GLOBAL＞PRINT_LEVEL为HIGH时，输出信息中关于原子力组分的若干问题
本帖最后由蟹黄宝于 2024-5-4 23:56 编辑

在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参见附件

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

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