在笼中心添加额外基函数的方式求助

量化小白

为了描述位于氟化碳笼中心的过量电子的扩散特性，我在笼的中心添加了一个H-Bq, 并针对笼中心的虚原子额外添加了4s4p4d,并将整个基组表示为6-311+G(d,p)+4s4p4d。其中s、p、d的指数来自于基组6-311+G(d,p)中碳的最低指数：0.0438,0.0438,0.626.具体输入文件如下：
%mem=1GB
%chk=E:\G092\12cage\12cageNBO.chk
#p b3lyp/6-311+g(d,p) pop=nbo IOp(3/24=1) extrabasis

12cageNBO

-1 2
C                 -0.60731700   -0.61763400   -0.16492400
C                  0.86284900   -0.31192500   -0.22564700
C                 -0.13380300    0.80768800   -0.12023800
C                 -0.91726200   -1.02972400    1.27844400
C                 -0.16325400    1.24114600    1.34977000
C                  1.42400700   -0.54259800    1.18173500
C                 -1.15878500    0.28792200    2.06408600
C                  1.18222100    0.77479800    1.96736400
C                  0.42846500   -1.49616800    1.89600000
C                 -0.59757700    0.05745600    3.47148800
C                  0.39863500   -1.06251500    3.36600900
C                  0.87269900    0.36263400    3.41072700
F                 -1.32265000   -1.15556400   -1.19856500
F                 -1.95366400   -1.96315200    1.39529900
F                 -2.48674600    0.73013700    2.04460600
F                  0.70282200   -2.85539500    1.70767000
F                 -1.38334800    0.28385900    4.56727300
F                  0.63013500   -1.97914900    4.35377800
F                  1.58844100    0.90008500    4.44432900
F                  2.21898600    1.70843700    1.85062600
F                 -0.43780200    2.60075800    1.53798800
F                  1.64857000   -0.53855600   -1.32146800
F                 -0.36561000    1.72428500   -1.10794900
F                  2.75190600   -0.98489700    1.20124600
H-Bq            0.07708265   -0.07815329    1.62530426

H 0
SP   1   1.00
      0.0438000              1.0000000              1.0000000
SP   1   1.00
      0.00438000              1.0000000              1.0000000
SP   1   1.00
      0.000438000              1.0000000              1.0000000
SP   1   1.00
      0.0000438000              1.0000000              1.0000000
D   1   1.00
      0.626                  1.000000
D   1   1.00
      0.1252                  1.000000
D   1   1.00
     0.02504           1.000000
D   1   1.00
    0.005008         1.000000
****

我想问一下我这样为H-Bq添加弥散函数，是否满足4s4p4d,如果不是应该怎样修改

sobereva

发帖时候注意恰当命名标题，当前的“为了描述位于氟化碳笼中心的过量电子的扩散特性，在6-311+G(d,p)基集中补充了额外...”明显不像个标题。我给你改了，以后注意。仔细阅读http://bbs.keinsci.com/thread-9348-1-1.html了解论坛里发帖要求。
当前写法不对。关键词里用了6-311+g(d,p)就已经给H定义了基函数了，后来文件末尾额外写的基函数是在此之上累加。自行用gfinput看一下实际用的计算便知。如果想让H只有那些基函数，用gen的方式定义。

量化小白

sobereva 发表于 2020-10-20 13:32
发帖时候注意恰当命名标题，当前的“为了描述位于氟化碳笼中心的过量电子的扩散特性，在6-311+G(d,p)基集中 ...

好的，非常感谢

量化小白

sobereva 发表于 2020-10-20 13:32
发帖时候注意恰当命名标题，当前的“为了描述位于氟化碳笼中心的过量电子的扩散特性，在6-311+G(d,p)基集中 ...

以后，会注意到的