高斯CCSD(T)/CBS算不动，有计算精度相近的方法基组吗？

林守龙

审稿人说我用CBS-QB3计算单点能量精度不够，推荐我用CCSD(T)/CBS，结果我的体系算了三天也没算出来，最终超过服务器的时间限制。请问各位老师，有没有其他和CCSD(T)/CBS相近精度的方法基组，耗时又不是太高，用来和CBS-QB3的计算结果做一个对比验证。

sobereva

DLPNO-CCSD(T)，或者结合MP2做外推
北京科音基础（中级）量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）ppt：



前不久刚回复过几乎完全一样的问题，提问前先用首页google框搜论坛再问

林守龙

sobereva 发表于 2024-3-7 16:21
DLPNO-CCSD(T)，或者结合MP2做外推
北京科音基础（中级）量子化学培训班（http://www.keinsci.com/worksho ...

好的谢谢老师

林守龙

sobereva 发表于 2024-3-7 16:21
DLPNO-CCSD(T)，或者结合MP2做外推
北京科音基础（中级）量子化学培训班（http://www.keinsci.com/worksho ...

老师，不会ORCA，只会使用高斯，您的QUICKCBS可以结合MP2做外推吗

zjxitcc

林守龙 发表于 2024-3-7 16:52
老师，不会ORCA，只会使用高斯，您的QUICKCBS可以结合MP2做外推吗

可使用MOKIT的autosr自动做DLPNO-CCSD(T)计算，不用学ORCA输入文件语法就可以现在开始计算（当然学也不难）。以水分子为例，gjf文件如下

1. %mem=20GB
   
2. %nprocshared=4
   
3. #p DLPNO-CCSD(T)/cc-pVTZ
   
4. 
   
5. mokit{}
   
6. 
   
7. 0 1
   
8. O     0.000000    0.000000    0.062007
   
9. H     0.000000   -0.783976   -0.492052
   
10. H     0.000000    0.783976   -0.492052

复制代码

完全符合Gaussian用户习惯，程序会自动分配辅助基。提交任务命令（假定已安装好Gaussian和ORCA，没装看ORCA安装教程）

1. autosr h2o.gjf >h2o.out 2>&1

复制代码

输出内容有

1. E(ref)  =       -76.05651770 a.u.
   
2. E(corr) =        -0.27471406 a.u.
   
3. E(DLPNO-CCSD(T)) =       -76.33123176 a.u.
   
4. T1_diag =    0.00645
   
5. Leave subroutine do_cc at Thu Mar  7 16:52:05 2024

复制代码

若想进一步提高精度，可将基组提升至cc-pVQZ，也可以将方法提升至DLPNO-CCSD(T1)。

yflchx

林守龙 发表于 2024-3-7 16:52
老师，不会ORCA，只会使用高斯，您的QUICKCBS可以结合MP2做外推吗

ORCA自学不难，Gaussian是起跳板，不能只局限于一个量子化学程序。

wzkchem5

林守龙 发表于 2024-3-7 09:52
老师，不会ORCA，只会使用高斯，您的QUICKCBS可以结合MP2做外推吗

审稿人可不管你会不会ORCA。
如果审稿人要求你做DLPNO-CCSD(T)计算，你是没有任何理由不算的，不可能以不会为借口不算，必须学。毕竟ORCA是免费的，如果是收费软件，还可以以没买为借口不用。

sobereva

林守龙 发表于 2024-3-7 16:52
老师，不会ORCA，只会使用高斯，您的QUICKCBS可以结合MP2做外推吗

可以

quickCBS手册里都给明了例子




相对于CCSD(T)，MP2外推就太轻松了