计算化学公社

标题: 求助使用ORCA计算三个分子间的势能时RI、格点、counterpoise的问题 [打印本页]

作者
Author: harper 时间: 2024-4-6 11:25
标题: 求助使用ORCA计算三个分子间的势能时RI、格点、counterpoise的问题
各位老师好，最近使用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
*

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

作者
Author: sobereva 时间: 2024-4-7 05:39
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 你所谓的“势能”那叫“相互作用能”

欢迎光临计算化学公社 (http://bbs.keinsci.com/)