基于分子力场的能量分解和IGMH分析遇到的问题

Jiang152

老师们好，模拟的两组片段：第一组片段对应的IGMH和EDA-FF分别为A1，A2；第二组片段对应的IGMH和EDA-FF分别为B1，B2。
    关于各个项的正负值有如下两个问题：
1. 在EDA-FF中分为静电作用，起排斥作用的交换互斥项以及起吸引作用的色散项。当静电值为正值时，是否代表此时静电作用表现为排斥力？反之，负值代表此时静电作用表现为排斥力？
2. A1中蓝色氢键作用力突出，所以A2的静电力为负值而且很大；但是B1中显示的蓝色氢键作用力大于斥力，为什么B2中静电力为正值？

sobereva

绝对不要光给出没什么用的散点图而不给出非常直观的IGMH填色等值面图

Jiang152

sobereva 发表于 2026-1-5 00:04
绝对不要光给出没什么用的散点图而不给出非常直观的IGMH填色等值面图

实在抱歉，sob老师。我在原来问题上把IGNH等值面图加上了。图B1等值面图中显示蓝色吸引力居多，但是图B2中静电力为正值，我不太理解。（针对这个现象，我再一次计算了第二组片段的EDA-FF，结果还是这样。）

sobereva

那叫IGMH
1 EDA-FF里静电项为负说明静电作用对结合是有利的，为正说明不利。顺带一提，如果你不打算分解成原子、原子对贡献的话，我更建议用下文的sobEDAw方法，比EDA-FF精度明显更高使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

2 我怀疑对此结构你的EDA-FF计算方式不对，连总相互作用能都成了正值

仔细检查EDA-FF分析给出的各个原子以及原子对对相互作用能的贡献，总能弄清楚原因

Jiang152

sobereva 发表于 2026-1-6 00:42
那叫IGMH
1 EDA-FF里静电项为负说明静电作用对结合是有利的，为正说明不利。顺带一提，如果你不打算分解成 ...

谢谢sob老师的专业指导，我再尝试一下。（ps: 不好意思，IGMH粗心打错成IGNH啦）

Jiang152

Jiang152 发表于 2026-1-6 09:22
谢谢sob老师的专业指导，我再尝试一下。（ps: 不好意思，IGMH粗心打错成IGNH啦）

原因找到啦，由于没严格按照操作过程进行，实际原子序数跟txt文档中的序数不一致导致。