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

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

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

谢谢老师回复，我知道了

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

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

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

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

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

不能用NPA，对静电势重现性很差，在《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）的静电势对比里有充分的体现

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