计算化学公社

标题: 在笼中心添加额外基函数的方式求助 [打印本页]

作者
Author: 量化小白 时间: 2020-10-20 10:54
标题: 在笼中心添加额外基函数的方式求助
为了描述位于氟化碳笼中心的过量电子的扩散特性，我在笼的中心添加了一个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,如果不是应该怎样修改

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

作者
Author: 量化小白 时间: 2020-10-20 14:49

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

好的，非常感谢

作者
Author: 量化小白 时间: 2020-10-20 14:53

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

以后，会注意到的

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