求助力场计算结果和量化计算结果相差较大

panMD

          在研究过程中，将CO2分子分别摆放到染料分子的各个官能团位置，通过量子化学计算结果与基于分子力场的能量分解分析（EDA-FF）结果结合讨论分析相互作用能。CO2在其他官能团位置处计算出的结果差异不大，例如，摆放到羰基位置处时，量化结果：-12.48 kj/mol  EDA-FF结果：-13.27 kj/mol。

          然而，二氧化碳在靠近氯原子及周围官能团位置处，EDA-FF计算出的交换互斥作用显著增强，导致计算出的总结合能与量化计算结果差异较大。例如，摆放到氨基位置时，二氧化碳中原子与氯原子的最短距离为3.529埃（图1），此时量化结果：--9.31kj/mol  EDA-FF结果：-0.18 kj/mol. 观察原子贡献，氯原子贡献的交换互斥作用较大（图2）。为了进一步说明EDA-FF计算有误，进行了IGMH分析，结果展现为色散吸引主导，这与EDA-FF计算出的交换互斥主导结论不一致。检查了原子类型等也没有发现错误，所以想请教一下问题出在了哪里，为什么会高估与氯原子间的交换互斥作用。

wudihua

量化计算有没有加色散呢？可以把输入文件放上来看一下

sobereva

这么小体系显然不该用EDA-FF，而应当用sobEDA，仔细看
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

panMD

sobereva 发表于 2024-10-8 03:37
这么小体系显然不该用EDA-FF，而应当用sobEDA，仔细看
使用sobEDA和sobEDAw方法做非常准确、快速、方便、 ...

因为研究中不仅限于对片段间相互作用能进行分解，反而更关注各个原子的贡献，所以优先选择的EDA-FF，没有考虑用您提出的sobEDA和sobEDAw能量分解。

panMD

wudihua 发表于 2024-10-7 17:19
量化计算有没有加色散呢？可以把输入文件放上来看一下

是有加的，输入文件如下：
%nprocshared=10
%mem=12GB
%chk=D:\Gaussian\gausssian data\ran-CO2\4TZVP-NH2(DUI).chk
#p B3LYP/def2tzvp opt=(calcall,ts,noeigen,gdiis) em=gd3bj

Generated by Multiwfn

0 1
。。。
。。。

@C:\ma-TZVP.txt

sobereva

可以尝试先用力场做几何优化再用EDA-FF

另注意基于力场准确描述CO2的静电作用必须考虑原子多极矩，光靠原子电荷算不准

panMD

感谢社长指点，我去学习。