计算化学公社

标题: EDA-FF能否用NPA替换MK电荷 [打印本页]

作者
Author: 七尺贱 时间: 2022-10-14 19:14
标题: EDA-FF能否用NPA替换MK电荷
本帖最后由七尺贱于 2022-10-14 21:17 编辑

老师您好，我的体系是碳环与金属富勒烯的复合，在用EDF-FF进行能量分解时，发现静电相互作用非常小，占比几乎为0。用ADF对类似的体系进行过能量分解，如下图所示，色散一项(灰色)占据绝对的主导，但静电(绿色)一项的贡献也有20%左右，所以我将MK电荷替换为了NPA电荷，结果有了一定的改进，不知这种方法是否可行？

作者
Author: sobereva 时间: 2022-10-14 20:01
不能用NPA，对静电势重现性很差，在《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）的静电势对比里有充分的体现

本身碳环的原子电荷就全都为零，就算受到外环境影响而稍微被极化，也不可能静电相互作用能有多大。

作者
Author: 七尺贱 时间: 2022-10-14 21:23

sobereva 发表于 2022-10-14 20:01
不能用NPA，对静电势重现性很差，在《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/ ...

谢谢老师，我又上传了一张图，我在计算B碳环形成的复合物时，算出来的elec是正的，但是我计算了二者的电荷转移，B碳环向富勒烯转移的电荷数是多于A碳环的，这感觉解释不了，所以我换成了NPA电荷，结果和预期是一致的
Interaction energy components between all fragments:（from MK charge）
                     Electrostatic Repulsion Dispersion    Total
Frag 1 -- Frag 2:       0.12    223.38    -388.37    -164.87

Interaction energy components between all fragments:（from NPA charge）
                     Electrostatic Repulsion Dispersion    Total
Frag 1 -- Frag 2:    -17.95    245.26    -406.60    -179.29

作者
Author: sobereva 时间: 2022-10-14 21:25

七尺贱发表于 2022-10-14 21:23
谢谢老师，我又上传了一张图，我在计算B碳环形成的复合物时，算出来的elec是正的，但是我计算了二者的电 ...

我不知道你具体拿什么结构算的原子电荷
始终记住要拿单体的结构计算，不要用复合物的。EDA-FF也并不体现出微量电子转移对结合能的轻微影响

作者
Author: 七尺贱 时间: 2022-10-14 21:26

sobereva 发表于 2022-10-14 20:01
不能用NPA，对静电势重现性很差，在《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/ ...

就是我想表达的意思是静电相互作用很微弱，但是最起码是起积极贡献的，但是起到了负贡献有点不太正常了

作者
Author: 七尺贱 时间: 2022-10-14 21:32

sobereva 发表于 2022-10-14 21:25
我不知道你具体拿什么结构算的原子电荷
始终记住要拿单体的结构计算，不要用复合物的。EDA-FF也并不体现 ...

谢谢老师回复，我知道了

作者
Author: sobereva 时间: 2022-10-14 21:50

七尺贱发表于 2022-10-14 21:26
就是我想表达的意思是静电相互作用很微弱，但是最起码是起积极贡献的，但是起到了负贡献有点不太正常了

理解清楚原子电荷计算原理，以及基于原子电荷的静电相互作用计算公式
也要明白，理想的结合状态是让整体势能最小，并不意味着此时片段间静电相互作用一定是最强的。只有静电作用对结合有显著贡献的情况，静电相互作用能才理应算出来是负值

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