|
Sob老师,在您的大作http://sobereva.com/46中提到:
BSSE校正后的真实的相互作用能这样计算:E_interaction = E_corrected - (E_A` + E_B`),其中E_A`和E_B`分别是A和B在孤立状态下经过优化后的能量。注意,相互作用能是复合物(考虑了BSSE时)减去单体在孤立状态时的能量差,单体在孤立状态的结构与在复合物中的结构并不相同,尤其单体结构呈柔性、分子间相互作用比较强时差异会较为明显。在实际中为了方便、省时往往做这样的近似:E_A`=E_A,E_B`=E_B,即是说计算单体能量时使用处在复合物状态下的结构,这个能量在计算E_BSSE时就已经顺便计算了,所以省得再单独算,此时相互作用能公式也有另一种等价写法E_interaction = E_AB - E_A,bAB - E_B,bAB,可见E_A和E_B并没被用到,其实可以不用在counterpoise任务中计算。
这段文字中的E_A,bAB和E_B,bAB所指的是不是分子A和分子B结合之前的孤立状态下的能量E_A和E_B减去其各自的BSSE啊?还是在结合后的结构下,单独计算A片段和B片段下的能量呢?
谢谢sob老师 |
|