计算化学公社

标题: g16 6-31G*基组定义问题 [打印本页]

作者
Author: mutron 时间: 2019-7-26 20:31
标题: g16 6-31G*基组定义问题
在算一个有机分子的时候突然发现直接书写6-31G*和用gen输入6-31G*的基函数数目竟然不一样，有点奇怪，比如：

# hf/6-31G* GFPRINT

Title Card Required

0 1
C                0.00000000 0.00000000 0.00000000
H                0.00000000 0.00000000 1.07000000
H                1.00880579 0.00000000 -0.35666635
H                -0.50440312 -0.87365131 -0.35666635
H                -0.50440269 0.87365135 -0.35666686

Atom C1    Shell    4 D 1    bf 10 - 15       0.000000000000       -0.000000000000       0.000000000000
   0.8000000000D+00  0.1000000000D+01

和

# hf/gen GFPRINT

Title Card Required

0 1
C                0.00000000 0.00000000 0.00000000
H                0.00000000 0.00000000 1.07000000
H                1.00880579 0.00000000 -0.35666635
H                -0.50440312 -0.87365131 -0.35666635
H                -0.50440269 0.87365135 -0.35666686

C H 0
6-31G*
****

Atom C1    Shell    4 D 1    bf 10 - 14       0.000000000000       -0.000000000000       0.000000000000
   0.8000000000D+00  0.1000000000D+01

（其它输出都一致）

前者是23个基函数，后者是22个（当然能量也不一样）。我用ORCA算了一次，和后者相符。试过def2SVP、cc-pvdz、6-311G*、3-21G、6-31G，都没有问题，出问题的是：6-31G*、6-31+G*、6-31G**，似乎是6-31G*的D极化函数出问题？

又或者是我的g16出问题了？谁来帮忙验证一下（版本ES64L-G16RevA.03）

作者
Author: liyuanhe211 时间: 2019-7-26 20:48
用gen的时候默认5d，不用gen直接用6-31G(d)关键词时默认6d。输入里声明5d 7f就一致了。

作者
Author: mutron 时间: 2019-7-26 20:54

liyuanhe211 发表于 2019-7-26 20:48
用gen的时候默认5d，不用gen直接用6-31G(d)关键词时默认6d。输入里声明5d 7f就一致了。

ORCA uses pure d and f functions (5D and 7F instead of Cartesian 6D and 10F) for all basis sets.
This needs to be taken into account when results are compared with other programs, especially for
Pople-style basis sets that were optimized with Cartesian (6D) functions

之前一直没留意，谢谢老师

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