基于分子力场进行能量分解的问题

段家康

大家好，我在基于分子力场做能量分解时，得到的结合能是-756.5KJ/mol，我自己计算得到的结合能是-58KJ/mol。但是能量分析的结果和预想的一样。这个结果能用吗，另外我想请问下结合能差距过大是因为我计算的时候加入了隐式溶剂模型吗，这是我计算的体系以及基组：b3lyp/6-31g(d,p) scrf=(solvent=water) em=gd3bj。我的本意是想通过能量分析，得到氢键的键能和氢键在结合能中的占比

wzkchem5

应该就是隐式溶剂模型的原因。可以想一下有电介质时的库伦能表达式U = kQq/(epsilon*R)，水的介电常数epsilon=80，所以阴阳离子离得足够远时，考虑水的结合能是不考虑水的结合能的1/80。当然实际你的计算结果没有这么夸张，一方面是因为阴阳离子是直接接触的，没有间隔水层，另一方面是有色散作用的贡献。

sobereva

应该跟Gaussian真空下算的对比。溶剂效应对相互作用的影响可以单作为一项。

另外，还得确保Multiwfn做能量分解时你设的原子电荷合理，原子类型分配的得当

段家康

wzkchem5 发表于 2023-2-13 23:09
应该就是隐式溶剂模型的原因。可以想一下有电介质时的库伦能表达式U = kQq/(epsilon*R)，水的介电常数epsil ...

谢谢老师，我明白了

段家康

sobereva 发表于 2023-2-14 09:25
应该跟Gaussian真空下算的对比。溶剂效应对相互作用的影响可以单作为一项。

另外，还得确保Multiwfn做能 ...

谢谢卢老师。