计算化学公社

标题: 如何往CFOUR的GENBAS基组库中添加基组 [打印本页]
作者
Author: 科研雪    时间: 2021-4-23 09:52
标题: 如何往CFOUR的GENBAS基组库中添加基组
本帖最后由 科研雪 于 2021-4-23 09:58 编辑

使用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


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



作者
Author: wjc404    时间: 2022-10-24 09:33
我试过可以先从basis set exchange网站上下载对应基组的cfour格式文本存为GENBAS文件，然后把里面每个元素开头位置的基组名称改为4-31G，最后在ZMAT中声明BASIS=4-31G，改过的GENBAS和ZMAT放一起可以运行xcfour算。可能有更好的方法。
作者
Author: zjxitcc    时间: 2023-2-18 00:27
本帖最后由 zjxitcc 于 2023-9-14 13:34 编辑

不用改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

  6. title

  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中，可以在线安装或下载预编译版解压即用。





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3