计算化学公社

标题: 多组份间相互作用能计算公式 [打印本页]

作者
Author: cgchen 时间: 2016-6-29 21:15
标题: 多组份间相互作用能计算公式
用MD模拟了含4种组份（A,B,C,D）的体系，想计算EA-B的相互作用能，不知下面的计算公式是否正确：

EA-B = Esys - EA - EB - EC - ED - (EA-C+EA-D+EB-C+EB-D+EC-D+EA-BCD+EB-ACD+EC-ABD+ED-ABC)

同一水平同一精度下

作者
Author: sobereva 时间: 2016-6-29 21:53
你搞得太复杂了。

如果已经跑出了轨迹，想研究轨迹中A-B的相互作用能，就直接计算E(A-B)就完了，用不着那么相减。因为一般用的分子力场都是对势，而不像量化那样需要考虑多体耦合。

作者
Author: cgchen 时间: 2016-6-29 22:25

sobereva 发表于 2016-6-29 21:53
你搞得太复杂了。

如果已经跑出了轨迹，想研究轨迹中A-B的相互作用能，就直接计算E(A-B)就完了，用不着 ...

基础啊！！！对势对势！！

这两种公式计算的能量之间有一定的差别但误差很小想知道误差愿意是什么呢

作者
Author: sobereva 时间: 2016-6-30 04:51

cgchen 发表于 2016-6-29 22:25
基础啊！！！对势对势！！

这两种公式计算的能量之间有一定的差别但误差很小想知道误差愿意是什 ...

这属于一些琐碎计算细节的事，没有具体情况不好说。反正你就直接计算A-B相互作用能就够了

作者
Author: cgchen 时间: 2016-6-30 08:34

sobereva 发表于 2016-6-30 04:51
这属于一些琐碎计算细节的事，没有具体情况不好说。反正你就直接计算A-B相互作用能就够了

哦了谢谢sob老师

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