通过溶质自由能计算结合能时高精度单点能的计算方法选择

xiewei

请问老师，我对两个分子进行结合能计算，先分别对两个单分子按照图中的流程进行了计算，在第二步计算真空中的高精度单点能用的B2PLYPD3/def2TZVP方法，但是按流程的第二步计算这两分子形成的复合物的高精度单点时由于原子数超过了50个，用B2PLYPD3/def2TZVP方法计算特别耗时间和内存，所以想改用M06-2X/def2TZVP，但是不知道和前面计算单分子的高精度单点采用了不同的计算方法对结合能计算是否由影响，我是否需要都统一使用M06-2X/def2TZVP来计算第二步的高精度单点能？还是在计算单分子时用B2PLYPD3/def2TZVP，计算复合物时用M06-2X/def2TZVP好一点呢？

sobereva

计算单分子和复合物时用的计算级别显然必须精确相同，否则求差毫无意义
def2-TZVP对于算化学键作用能可以，但用于弱相互作用能的计算太low，要么带上弥散函数，要么考虑counterpoise校正，二者都考虑更好，看下文里不同基组算相互作用能的对比表格
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

双杂化图快用ORCA。如果就用Gaussian算分子间相互作用能，M06-2X对于发表文章目的一般可以接受