使用gCP方法考虑BSSE时如何计算弱相互作用能？

121089576

在用CP和gCP考虑BSSE计算弱相互作用能时，有几个感觉很迷惑人的概念，希望与大家共同交流交流，如下。

关于计算公式的问题——
使用CP方法时，相互作用能计算公式为E_interaction = E_AB + E_BSSE_AB - E(A) - E(B) ，
（1）那使用gCP方法时，计算公式也是E_interaction = E_AB + E_gCP_AB - E(A) - E(B) 吗？
（2）如果（1）的说法正确，那问题又来了：用CP方法是只对复合物(AB)进行BSSE校正，而Grimme的文章似乎说用gCP则需要对复合物(AB)和单体(A)(B)都进行校正，那这样的话计算公式就应该是E_interaction = E_AB + E_gCP_AB - [E(A) + E_gCP_A] - [E(B) + E_gCP_B]了。

关于结合DFT-D3的问题——
（3）比如用CP结合D3，其实就是在关键词行加上em=GD3对吗？那用gCP结合D3，也就是加上D3的关键词对吗？
（4）如果（3）的说法正确，那我是否可以认为 E_interaction = E_AB + E_BSSE_AB - E(A) - E(B) 这个计算公式并没有考虑D3？
        如果要考虑D3的话，那基于CP的计算公式就应该是E_interaction = E_AB + E_BSSE_AB + E_D3_AB - E(A) - E(B) 了。
        那考虑D3，基于gCP的计算公式就应该是E_interaction = E_AB + E_gCP_AB + E_D3_AB - [E(A) + E_gCP_A + E_D3_A] - [E(B) + E_gCP_B + E_D3_B]了。

以上说法是否正确？

sobereva

1 否。直接E_AB - E(A) - E(B)，其中每个部分的能量计算都带着gCP

3 无论是CP（counterpoise）还是gCP，跟带不带D3完全没关系。平时D3该怎么写还怎么写。带不带D3属于不同的理论方法，别当成一个对相互作用能的校正项看待。