提问用molclus调用ORCA计算的结果正确与否和是否需要考虑BSSE

王涛

使用Molclus结合xtb做的动力学模拟对瑞德西韦(Remdesivir)做构象搜索
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=16255&fromuid=84268
(出处: 计算化学公社)
在根据老师的帖子做到用Molclus调用Gaussian和ORCA得到水环境下的高精度自由能这一步时，得到的orcaSP00001等文件都仅为77KB，具体文件已经以附件形式上传。
我所用的ORCA模板内容为：
! revDSD-PBEP86-D4/2021 def2-TZVPP def2/J def2-TZVPP/C RIJCOSX tightSCF noautostart miniprint nopop
%maxcore 9000
%pal nprocs 32 end
%cpcm
smd true
SMDsolvent "water"
end
* xyz -2 1
[GEOMETRY]
*
有几点想请教一下：
1.ORCA计算是否正确正常？
2.若正常，是否可以把pal nprocs提高然后把maxcore改小吗？因为总共用的也就77KB多一些。
3.可以Gaussian频率计算完成后，检查没有谱虚频后，再启用ORCA吗？主要是担心有虚频算出的单点能不能用，浪费计算资源。
4.使用ORCA计算弱相互作用体系单点能是否还要考虑BSSE（基组重叠误差）问题？如果需要应该如何处理？

sobereva

1 正常
2 为什么总有人不看文件内容而光盯着文件大小看？看文件大小有什么意义？
3 不能
4 完全取决于计算目的、用什么基组，以及什么方式考虑BSSE（counterpoise还是gCP）。不计算片段间弱相互作用能就牵扯不到counterpoise

王涛

sobereva 发表于 2025-11-13 01:41
1 正常
2 为什么总有人不看文件内容而光盯着文件大小看？看文件大小有什么意义？
3 不能

谢谢老师的解答。
在ORCA中做counterpoise校正并计算分子间结合能的例子
http://bbs.keinsci.com/forum.php ... s/logo-O35E636P.png
(出处: 计算化学公社)
我看您这个帖子，我想算的是A-B体系的结合能，请问
1.我构型搜索得到比如5个结构，都需要做counterpoise校正吗？还是可以选择能量最低的计算就可以呢？
2.molclus和ORCA中做counterpoise校正有没有联动呢？还是需要molclus做完结构优化和频率分析后，再单独做ORCA的counterpoise校正呢？
3.如果选择用gCP方法来校正，请问关键词应该如何写呢？

sobereva

王涛 发表于 2025-11-14 09:47
谢谢老师的解答。
在ORCA中做counterpoise校正并计算分子间结合能的例子
http://bbs.keinsci.com/forum ...

1 前面说了，对于算相互作用能活结合能问题，counterpoise的必要性首先取决于基组，不说这个最基本前提没法回答
倘若就是非要做counterpoise校正，显然对哪个构型算结合能就需要对哪个体系做

2 算结合能是自己的事，跟molclus干的事毫无直接关系。molclus做的仅仅是构型和构象搜索

3 RTFM