计算化学公社

标题: 关于弱相互作用中Mayer能量分解的数据分析 [打印本页]

作者
Author: molx79 时间: 2016-7-7 18:27
标题: 关于弱相互作用中Mayer能量分解的数据分析
本帖最后由 molx79 于 2016-7-7 18:30 编辑

各位老师：

大家好。在学习了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) 中获取。

谢谢老师们赐教！

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

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

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