如何往CFOUR的GENBAS基组库中添加基组

科研雪

使用CFOUR的时候想用def2-TZVP基组，但是GENBAS基组库中没有这个基组，我在CFOUR网站中找到了def2-TZVP基组对应的基组信息，并复制过去了，但是计算的时候报错了  @GTFLGS-W, Option Default :  
  @GTFLGS-W, Option def2-TZVP for keyword BASIS       not known.  Default used.
  @GTFLGS-I, Calculation type set to UHF.
@MKDECK-F, Premature end-of-file on MOL.
  @CHECKOUT-I, Total execution time :       0.0100 seconds.


  ERROR ERROR ERROR ERROR ERROR ERROR
  Job has terminated with error flag
  ERROR ERROR ERROR ERROR ERROR ERROR


--executable xjoda finished with status          256


请问有人往基组库添加过基组吗，应该注意些什么

wjc404

我试过可以先从basis set exchange网站上下载对应基组的cfour格式文本存为GENBAS文件，然后把里面每个元素开头位置的基组名称改为4-31G，最后在ZMAT中声明BASIS=4-31G，改过的GENBAS和ZMAT放一起可以运行xcfour算。可能有更好的方法。

zjxitcc

不用改CFOUR自带基组文件。fch2cfour小程序可解君愁（一键产生ZMAT, OLDMOS，GENBAS文件）。我用的是CFOUR-v2.1，先到CFOUR输出文件里复制一下Bohr为单位的直角坐标（因为fch2cfour目前无法预测CFOUR对分子的自动转向），写个gjf文件（不妨叫test.gjf）

1. %chk=test.chk
   
2. %mem=150GB
   
3. %nprocshared=28
   
4. #p UHF/def2TZVP nosymm int=nobasistransform stable=opt units=au
   
5. 
   
6. title
   
7. 
   
8. -1 2
   
9. B       -1.89884372     3.30923345     0.08175460
   
10. B       -4.44569510    -1.46862696     0.25661782
    
11. B        3.49387019     2.98229456     0.30279830
    
12. B        0.90348380    -4.50876009     0.27771596
    
13. B       -1.86366789     0.00000027    -0.55404169
    
14. B        0.87294323     1.59135024    -0.65182735
    
15. B        0.90348556     4.50875979     0.27771596
    
16. B        3.49386905    -2.98229585     0.30279830
    
17. B       -4.44569455     1.46862854     0.25661782
    
18. B       -1.89884497    -3.30923266     0.08175460
    
19. B        4.01215175    -0.00000078     0.01992303
    
20. B        0.87294265    -1.59135050    -0.65182735
    

复制代码

高斯算完后获得test.fch文件，运行

1. fch2cfour test.fch

复制代码

即生成ZMAT, OLDMOS，GENBAS文件。内坐标、基组数据都写好了，连初始轨道都提供了。自动生成的ZMAT文件里默认关键词是做UHF计算的，你想做CCSD(T)计算可以将文件末尾改成例如

1. *CFOUR(CALC=CCSD(T),REF=UHF,SYM=OFF,CHARGE=-1,MULTI=2,BASIS=PVTZ
   
2. CC_PROGRAM=ECC,ABCDTYPE=AOBASIS,MEMORY=400000000)

复制代码

此处PVTZ仅是借用字符串，真正基组信息在GENBAS里。提交给CFOUR，例如

1. xcfour >OUTPUT 2>&1 &

复制代码

可发现UHF 2圈收敛，立即进入CC计算步骤。省去了书写基组的麻烦，SCF时间也省了，波函数稳定性也在高斯中检查过了。fch2cfour小程序在开源程序MOKIT中，可以在线安装或下载预编译版解压即用。