计算弱相互作用的问题

molx79 · 发表于 Post on 2015-5-9 11:27:09

各位老师：我是做弱相互作用的。我现在手动计算BBSE。
1、复合物的初始构型 opt freq（得到不含BSSE能量矫正）
2、上步得到的复合物的稳定构型 counter=2，得到BSSE能量。
那么在做NBO、AIM、构型（键长键角），EDA能量分析，是否都可以用第一步计算得到的构型（即不含counter=N的关键词）？

此外，计算弱相互作用含不含counter=N，复合物的构型会有差别，甚至差别很大（不是一种构型），那么在讨论NBO、AIM、构型（键长键角），EDA能量分析，用哪一种构型更合理？

恳请各位老师指导！

sobereva · 发表于 Post on 2015-5-9 16:36:29

NBO、AIM、构型都可以用第一步的。
EDA取决于具体怎么搞，如果是LMOEDA，程序直接就会输出考虑和不考虑BSSE情况下的能量分解结果。

几何优化的时候用的基组有足够的弥散函数时就不需要考虑BSSE对几何结构的影响，影响会甚微。不建议在几何优化时用counterpoise考虑BSSE，因为此时没有解析梯度，会导致几何优化耗时增加一个数量级甚至更多。所以几何优化时应当带弥散函数，仅在最后算能量的时候用counterpoise考虑BSSE校正能，这样最划算。
另外，也可以在几何优化时用ORCA支持的gCP方法考虑BSSE，此时基组没弥散函数也无妨，而且不会增加额外增加任何耗时，这使得def2-SVP结合gCP可以又便宜又好地优化含弱相互作用的大体系。

molx79 · 发表于 Post on 2015-5-10 14:36:18

谢谢sob老师。我以SCl2和HCHO为例尝试了一下做LMOEDA,用含ounterpoise的构型计算结果比不含ounterpoise关键词的Eint能量的绝对值要小。如图1所示。
另外您所说的gamess程序直接就会输出考虑和不考虑BSSE情况下的能量分解结果。如何看呢？我每次都是看OWN BASIS SET 的结果，其中的TOTAL INTERACTION ENERGY HF OR DFT 是什么意思？ALL BASIS SET是什么意思，两者的差是EBSSE吗？谢谢sob老师。

sobereva · 发表于 Post on 2015-5-10 16:10:13

ALL BASIS SET就是考虑了BSSE情况下的能量成分和总相互作用能。和OWN BASIS SET输出的总相互作用能的差值大抵就是BSSE能，你可以自行验证一下。

做MP2前会先做HF，所以你在MP2下得到的能量分解的结果会顺带着输出HF级别下的结果。

molx79 · 发表于 Post on 2015-5-10 16:32:09

sobereva 发表于 2015-5-10 16:10
ALL BASIS SET就是考虑了BSSE情况下的能量成分和总相互作用能。和OWN BASIS SET输出的总相互作用能的差值大 ...

万分感谢Sob老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 计算弱相互作用的问题