计算化学公社

标题: 如何手动编辑基组截去高角动量轨道？ [打印本页]

作者
Author: alanzhou 时间: 2022-11-2 20:21
标题: 如何手动编辑基组截去高角动量轨道？
  各位老师好！想请教一下关于基组截断的问题。比如相关一致基组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

作者
Author: zjxitcc 时间: 2022-11-2 20:28
那不如直接使用3-21G，写基组名更方便

作者
Author: alanzhou 时间: 2022-11-2 22:19
老师你好！其实是最近在看文献时发现一个名为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格式的，所以想问一下如何通过自己定义基组来算能量。

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

作者
Author: alanzhou 时间: 2022-11-3 12:05
谢谢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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)