请问在溶剂模型下怎么用ORCA计算柔性分子的结合能

子不语i · 发表于 Post on 2023-11-10 11:48:57

本帖最后由子不语i 于 2023-11-13 09:07 编辑

各位老师好，在sob老师的《在ORCA中做counterpoise校正并计算分子间结合能的例子》（http://sobereva.com/542）一文的例子中提到“注意按照下文例子算出来的实际上是复合物结构中单体间相互作用能。而一般意义的结合能在计算时是要考虑单体形成复合物过程中形变造成的能量改变（变形能）。由于下面的例子里单体是刚性的，变形能可以忽略不计，因此下文仍称作结合能。”请问如果想要计算柔性分子在溶剂模型下的结合能：
1.对于公式E(AB)-E(A)-E(B)+E_BSSE，其中的E(A)和E(B)是不是应该用初始未结合未变形的分子结构计算得到的能量呢？
2.在溶剂模型下计算结合能，是不是就像博文中提到的要在气相下计算E_BSSE，然后加到常规方式算的E(AB)-E(A)-E(B)上？
使用的方法，opt和freq是B97-3c，单点能计算是r2SCAN-3c。谢谢各位老师。

wzkchem5 · 发表于 Post on 2023-11-13 16:48:41

本帖最后由 wzkchem5 于 2023-11-13 09:51 编辑

r2SCAN-3c（以及几乎所有其他3c方法）自带gCP校正，不需要做counterpoise，做counterpoise反而会把BSSE校正算两遍（一遍是gCP，一遍是counterpoise），反而是错的。

子不语i · 发表于 Post on 2023-11-13 21:00:20

wzkchem5 发表于 2023-11-13 16:48
r2SCAN-3c（以及几乎所有其他3c方法）自带gCP校正，不需要做counterpoise，做counterpoise反而会把BSSE校正 ...

好的，谢谢老师

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

[ORCA] 请问在溶剂模型下怎么用ORCA计算柔性分子的结合能