计算化学公社

标题: 关于相互作用能和能量分解 [打印本页]

作者
Author: captain 时间: 2017-4-18 16:27
标题: 关于相互作用能和能量分解
请问大神
两个单体形成复合物的体系，
先用wB97XD/6-311+G** 分别计算复合物和单体的单点能（考虑BSSE），再相减得到的相互作用能，研究相互作用能的大小
再用SAPT0/6-311+G** 进行能量分解，研究相互作用能中的各种能量成分
这样做合理么？两种不同的方法得到的相互作用能具有可比性吗？

比如说两个体系 A1-B1 和 A2-B2
wB97XD/6-311+G** 得到的相互作用能的相对大小与
SAPT0/6-311+G** 进行能量分解得到的相互作用能的相对大小
理论上来说应该是吻合的吧？
感谢解答！

作者
Author: sobereva 时间: 2017-4-18 20:00
两种方法算的总相互作用能的定量数值对不上，但用SAPT0讨论相互作用的内部成分比例还是可以的

不知道你用什么程序做的SAPT0，如果是molpro做的DFT-SAPT的话，算的相互作用能精度比wB97XD明显要高，都没必要报道wB97XD结果了。

作者
Author: 冰释之川 时间: 2017-4-18 20:09
本帖最后由冰释之川于 2017-4-18 20:16 编辑

sobereva 发表于 2017-4-18 20:00
两种方法算的总相互作用能的定量数值对不上，但用SAPT0讨论相互作用的内部成分比例还是可以的

不知道你 ...

那能不能通过SAPT0/6-311+G** 进行能量分解得到的成分和总能量与在wB97XD/6-311+G** 得到的相互作用能进行Scale一下，这样比例也可以在wB97XD/6-311+G** 得到的相互作用能下算出来了

作者
Author: sobereva 时间: 2017-4-18 20:15

冰释之川发表于 2017-4-18 20:09
那能不能通过SAPT0/6-311+G** 进行能量分解得到的成分和总能量与在B97XD/6-311+G** 得到的相互作用能进行 ...

可以
不过既然SAPT0算得动，还不如把DFT计算时候级别再弄好点。比如双杂化-D3结合ma-TZVP

作者
Author: captain 时间: 2017-4-18 22:17

sobereva 发表于 2017-4-18 20:00
两种方法算的总相互作用能的定量数值对不上，但用SAPT0讨论相互作用的内部成分比例还是可以的

不知道你 ...

大师感谢解惑
我是用PSI4算的SAPT0
两种方法比较相互作用能的绝对值肯定没有意义
但是定性来看
比如说wB97XD算出来A1-B1的作用能比A2-B2的作用能更负
SAPT0理论上应该也能得到相同的结果吧？

作者
Author: sobereva 时间: 2017-4-19 01:00

captain 发表于 2017-4-18 22:17
大师感谢解惑
我是用PSI4算的SAPT0
两种方法比较相互作用能的绝对值肯定没有意义

是的

作者
Author: captain 时间: 2017-4-19 09:13

sobereva 发表于 2017-4-19 01:00
是的

谢谢！

作者
Author: lao7 时间: 2017-4-19 14:49

sobereva 发表于 2017-4-18 20:00
两种方法算的总相互作用能的定量数值对不上，但用SAPT0讨论相互作用的内部成分比例还是可以的

不知道你 ...

他这个相互作用能分析，能否通过multiwfn分析实现？分析主要成键原子对配体贡献的能量百分比？谢谢！

作者
Author: sobereva 时间: 2017-4-19 21:53

lao7 发表于 2017-4-19 14:49
他这个相互作用能分析，能否通过multiwfn分析实现？分析主要成键原子对配体贡献的能量百分比？谢谢！

Multiwfn+Gaussian可以做简单的能量分解，但没有SAPT这么细
“成键原子对配体贡献的能量百分比”我不清楚什么意思

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