如何手动编辑基组截去高角动量轨道？

alanzhou

  各位老师好！想请教一下关于基组截断的问题。比如相关一致基组cc-pVDZ计算氢原子如果不考虑p轨道函数进行计算，我理解就是在Gaussian的基组里把P轨道那一行删掉，不知道这样理解对不对？但是MRCC/CFOUR格式的基组表示形式不同，不知道该怎么理解，实在看不出来p轨道函数是怎么表示的。请问各位大佬MRCC/CFOUR格式的基组该怎么修改呢？
Gaussian16里的H原子cc-pVDZ基组：
! CC-pVDZ basis  5-Apr-2019
-H
S   3 1.00
0.1301000000D+02  0.1968500000D-01
0.1962000000D+01  0.1379770000D+00
0.4446000000D+00  0.4781480000D+00
S   1 1.00
0.1220000000D+00  0.1000000000D+01
P   1 1.00
0.7270000000D+00  0.1000000000D+01

MRCC 2022里的H原子cc-pVDZ基组：
H:cc-pVDZ               
Dunning's correlation consistent basis set

    2
    0    1
    2    1
    4    1

         13.010000000000          1.962000000000          0.444600000000          0.122000000000

          0.019685000000          0.000000000000
          0.137977000000          0.000000000000
          0.478148000000          0.000000000000
          0.501240000000          1.000000000000

          0.727000000000

          1.000000000000

zjxitcc

那不如直接使用3-21G，写基组名更方便

alanzhou

老师你好！其实是最近在看文献时发现一个名为W3X-L的组合方法，想应用一下来计算一个简单的自由基(DOI: 10.1021/acs.jctc.5b00135)。里面有一个能量校正量deltaQ是用CCSDT(Q)/VDZ(d)-CCSDT/VDZ(d)计算的。采用的VDZ(d)基组在计算H原子时就是用cc-pVDZ基组截断p轨道基函数做的（也就是只计算了s轨道，C原子等重原子只计算了sp轨道）。但是CCSDT(Q)方法我只知道MRCC程序可以支持，MRCC的基组定义是CFOUR格式的，所以想问一下如何通过自己定义基组来算能量。

sobereva

BSE上能给出CFOUR格式的基组定义，BSE上面也有6-31G*的，效仿着Gaussian格式和CFOUR格式的H的6-31G*的定义关系，去把Gaussian格式的去掉p的cc-pVDZ改写成CFOUR格式的即可。最后再用Gaussian和CFOUR做个HF计算检验一下能量对不对

alanzhou

谢谢zjxitcc与sobereva两位老师的指导，目前修改后的VDZ(d)基组计算验证已通过Molpro以及MRCC程序计算验证。基组如下：
H:cc-pVDZ(d)               
Dunning's correlation consistent basis set with truncated p

    1
    0   
    2   
    4   

         13.010000000000          1.962000000000          0.444600000000          0.122000000000

          0.019685000000          0.000000000000
          0.137977000000          0.000000000000
          0.478148000000          0.000000000000
          0.501240000000          1.000000000000
删去p轨道极化函数即可。以甲醇为例验证结果：
molpro(molpro可通过关键词‘s,H,VDZ;c;'实现截去p轨道函数，然后进行HF计算，再调用MRCC的CCSDT以及CCSDT(Q)方法实现；详见文献DOI: 10.1021/acs.jctc.5b00135附件):
cc-pVDZ:
       CCSDT(Q)           CCSDT          HF-SCF
   -115.42118046   -115.42047178   -115.04924300
cc-pVDZ(d):
        CCSDT(Q)           CCSDT          HF-SCF
   -115.37029536   -115.36956620   -115.02819585
MRCC通过手动编辑基组计算与molpro结果完全一致。
CFOUR格式基组定义：
http://slater.chemie.uni-mainz.d ... ntryInTheGENBASFile