关于弱相互作用中Mayer能量分解的数据分析

molx79

各位老师：

    大家好。在学习了sobereva老师的博文《Mayer能量分解及APEX4程序使用简介

》。尝试用Mayer方案分析弱相互作用。以SHCl…SHCl二体中的S…S键为例。

处理过程：

Step1: MP2/6-311G**  opt+freq 时产生chk文件，转成Fchk

Step2：用APEX4程序计算

部分结果如下：

"Exact" energy decomposition matrix

              1  S      2 H         3 Cl         4 S       5 H        6 Cl

   1 S -397.4180   -0.1457  -0.1569    0.0027      0.0015     -0.0011

   2 H   -0.1457   -0.4278    0.0120    0.0058      0.0014    -0.0033

   3Cl   -0.1569    0.0120  -459.4713   -0.0022   -0.0014    0.0030

   4 S    0.0027    0.0058   -0.0022   -397.4229  -0.1443   -0.1561

   5 H    0.0015    0.0014   -0.0014   -0.1443     -0.4288    0.0126

   6Cl   -0.0011   -0.0033    0.0030   -0.1561     0.0126    -459.4676

问题1：如何计算SHCl与SHCl形成二体的作用能Eint?（博文里提到“如果想获得片段间作用能将相应原子对儿作用能加和即可”，不太理解这句话）

问题2：博文中提到“每一对儿原子间的能量可以被分解为"动能"+"静电"+"交换"+"重叠"项”，如果想把第一步得到的两个片段间的作用能分解为：动能"+"静电"+"交换"+"重叠"，如何从这些数据ELECTROSTATIC CONTRIBUTIONS 、EXCHANGECONTRIBUTIONS、OVERLAP CONTRIBUTIONS、KINETIC ENERGY (CECA/T) 中获取。

    谢谢老师们赐教！

sobereva

1  ∑[i∈片段1]∑[j∈片段2] E_i,j。当然这根实际相互作用能（整体减片段）得到的会差很多。所以Mayer能量分解并不能准确分解分子间相互作用能。准确的分解需要EDA、SAPT等方式，但这些方法又没法分解成原子-原子之间的作用项。

2 也是按照如上方式，只不过E_i,j此时就不是原子间的总作用能了，而是相应能量成份