求助使用ORCA计算三个分子间的势能时RI、格点、counterpoise的问题

harper

各位老师好，最近使用orca计算三个分子间的势能，写了以下输入文件，但有几个问题：
1.设置以下关键字：
! RI-PWPB95 D4 RIJK def2-QZVPP def2/JK def2-QZVPP/C miniprint nopop tightscf
使用密度拟合方法加速时，泛函前一定要加RI-前缀吗：RI-PWPB95，不加会有什么影响；
这里没有出现RI-MP2加速方法，为什么还用到了def2-QZVPP/C 辅助基组；
ORCA5不能出现GRID关键字，那么我们该如何设置网格呢？

2.通过虚拟原子的方式计算的分子间势能与gaussian中通过counterpoise校正在效果上有什么区别
3.通过下方输入文件中的方式计算的势能
①分子体系整体的能量 E - 三个分子自身的能量 [E_A(ABC)+E_B(ABC)+E_C(ABC)] =三体势能 E_ABC(ABC) + 两两之间的二体势能 [E_AB(ABC)+E_BC(ABC)+E_AC(ABC)]
②分子体系整体的能量 E - 三个分子自身的能量 [E_A(ABC)+E_B(ABC)+E_C(ABC)] - 两两之间的二体势能 [E_AB(ABC)+E_BC(ABC)+E_AC(ABC)] =三体势能 E_ABC(ABC)
PS：我对势能的理解是ABC三个分子间的势能=三体势能+两两之间的二体势能，如有错误，感谢各位老师指正
PS：E_A(ABC)表示ABC基函数下A的能量，E_AB(ABC)表示ABC基函数下AB之间的二体势能，E_ABC(ABC)表示ABC基函数下ABC三个分子间的三体势能

输入文件
%pal nprocs   4 end

! RI-PWPB95 D4 RIJK def2-QZVPP def2/JK def2-QZVPP/C miniprint nopop tightscf
%maxcore  1000
* xyz   0   1
O                 -3.78900161    0.62830247   -0.08276004
H                 -2.82900161    0.62830247   -0.08276004
H                 -4.10945619    1.53323830   -0.08276004
O                 -1.23672207   -1.49094034   -0.07548478
H                 -0.27672207   -1.49094034   -0.07548478
H                 -1.55717665   -0.58600451   -0.07548478
O                 -1.08018734   -0.61603942   -2.14971258
H                 -0.12018734   -0.61603942   -2.14971258
H                 -1.40064193    0.28889641   -2.14971258
*

$new_job
! RI-PWPB95 D4 RIJK def2-QZVPP def2/JK def2-QZVPP/C miniprint nopop tightscf Pmodel
%maxcore  1000
* xyz   0   1
O                 -3.78900161    0.62830247   -0.08276004
H                 -2.82900161    0.62830247   -0.08276004
H                 -4.10945619    1.53323830   -0.08276004
O:                -1.23672207   -1.49094034   -0.07548478
H:                -0.27672207   -1.49094034   -0.07548478
H:                -1.55717665   -0.58600451   -0.07548478
O:                -1.08018734   -0.61603942   -2.14971258
H:                -0.12018734   -0.61603942   -2.14971258
H:                -1.40064193    0.28889641   -2.14971258
*

$new_job
! RI-PWPB95 D4 RIJK def2-QZVPP def2/JK def2-QZVPP/C miniprint nopop tightscf Pmodel
%maxcore  1000
* xyz   0   1
O:                -3.78900161    0.62830247   -0.08276004
H:                -2.82900161    0.62830247   -0.08276004
H:                -4.10945619    1.53323830   -0.08276004
O                 -1.23672207   -1.49094034   -0.07548478
H                 -0.27672207   -1.49094034   -0.07548478
H                 -1.55717665   -0.58600451   -0.07548478
O:                -1.08018734   -0.61603942   -2.14971258
H:                -0.12018734   -0.61603942   -2.14971258
H:                -1.40064193    0.28889641   -2.14971258
*

$new_job
! RI-PWPB95 D4 RIJK def2-QZVPP def2/JK def2-QZVPP/C miniprint nopop tightscf Pmodel
%maxcore  1000
* xyz   0   1
O:                -3.78900161    0.62830247   -0.08276004
H:                -2.82900161    0.62830247   -0.08276004
H:                -4.10945619    1.53323830   -0.08276004
O:                -1.23672207   -1.49094034   -0.07548478
H:                -0.27672207   -1.49094034   -0.07548478
H:                -1.55717665   -0.58600451   -0.07548478
O                 -1.08018734   -0.61603942   -2.14971258
H                 -0.12018734   -0.61603942   -2.14971258
H                 -1.40064193    0.28889641   -2.14971258
*

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “使用ORCA计算三个分子间的势能” 改了，以后务必注意！

sobereva

1 需。这种事对比便知
双杂化里有MP2部分，显然需要用/C
格点问题看下文
谈谈ORCA 5.0的新特性和改变
http://sobereva.com/604（http://bbs.keinsci.com/thread-23917-1-1.html）

2 本来Gaussian里也可以用Bq等效地实现counterpoise关键词的效果，仔细看下文搞清楚原理
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

3 你所谓的“势能”那叫“相互作用能”