计算化学公社

标题: Fe，C原子的6-31+g(d)基组形式的问题 [打印本页]

作者
Author: chenxin199261 时间: 2014-12-6 16:03
标题: Fe，C原子的6-31+g(d)基组形式的问题
本帖最后由 chenxin199261 于 2014-12-6 18:11 编辑

我在gaussian中，用iop(3/24=10) 把基组打印了出来
Fe（设定为6-31+g(d)）的基组如下：
S 6 1.00    0.000000000000
   0.6113262000D+05  0.1766110976D-02
   0.9179342000D+04  0.1353037982D-01
   0.2090857000D+04  0.6673127910D-01
   0.5892479000D+03  0.2314822969D+00
   0.1887543000D+03  0.4797057935D+00
   0.6444629000D+02  0.3501975953D+00
S 1 1.00    0.000000000000
   0.1256700000D-01  0.1000000000D+01
SP 6 1.00    0.000000000000
   0.1259980000D+04  0.2438014027D-02  0.4028018665D-02
   0.2998761000D+03  0.3224048035D-01  0.3144646739D-01
   0.9684917000D+02  0.1265724014D+00  0.1368316886D+00
   0.3631020000D+02 -0.3139902035D-01  0.3487235710D+00
   0.1472996000D+02 -0.6207593068D+00  0.4617930616D+00
   0.6066075000D+01 -0.4502914050D+00  0.2043057830D+00
SP 6 1.00    0.000000000000
   0.5043485000D+02 -0.3873255984D-02 -0.7017127880D-02
   0.1683929000D+02  0.7196597971D-01 -0.2877659951D-01
   0.7192086000D+01  0.2556590990D+00  0.6181382895D-01
   0.3053420000D+01 -0.2882836988D+00  0.3954945933D+00
   0.1273643000D+01 -0.7342821970D+00  0.4989058915D+00
   0.5040910000D+00 -0.2049352992D+00  0.1791250969D+00
SP 3 1.00    0.000000000000
   0.1950316000D+01  0.5694869031D-01 -0.4593796163D+00
   0.7367210000D+00  0.2882915015D+00  0.2852139102D+00
   0.1141770000D+00 -0.1138159006D+01  0.9076485323D+00
SP 1 1.00    0.000000000000
   0.4114800000D-01  0.1000000000D+01  0.1000000000D+01
P 1 1.00    0.000000000000
   0.1349150000D+00  0.1000000000D+01
P 1 1.00    0.000000000000
   0.4184300000D-01  0.1000000000D+01
D 3 1.00    0.000000000000
   0.2314994000D+02  0.8876935479D-01
   0.6122368000D+01  0.3896319210D+00
   0.1846601000D+01  0.7014816379D+00
D 1 1.00    0.000000000000
   0.5043610000D+00  0.1000000000D+01
D 1 1.00    0.000000000000
   0.1133000000D+00  0.1000000000D+01
F 1 1.00    0.000000000000
   0.8000000000D+00  0.1000000000D+01
============================================================
我想知道P型的高斯函数是什么，他在gamess里面的符号又是什么呢 ?

C（设定为6-31+g(d)）的基组如下：
S 6 1.00    0.000000000000
   0.3047524880D+04  0.1834737132D-02
   0.4573695180D+03  0.1403732281D-01
   0.1039486850D+03  0.6884262226D-01
   0.2921015530D+02  0.2321844432D+00
   0.9286662960D+01  0.4679413484D+00
   0.3163926960D+01  0.3623119853D+00
SP 3 1.00    0.000000000000
   0.7868272350D+01 -0.1193324198D+00  0.6899906659D-01
   0.1881288540D+01 -0.1608541517D+00  0.3164239610D+00
   0.5442492580D+00  0.1143456438D+01  0.7443082909D+00
SP 1 1.00    0.000000000000
   0.1687144782D+00  0.1000000000D+01  0.1000000000D+01
SP 1 1.00    0.000000000000
   0.4380000000D-01  0.1000000000D+01  0.1000000000D+01
SP 1 1.00    0.000000000000
   0.1687144782D+00  0.1000000000D+01  0.1000000000D+01
D 1 1.00    0.000000000000
   0.8000000000D+00  0.1000000000D+01
======================================
与EMSL数据库给出的基组进行比对，发现多了一个
SP 1 1.00    0.000000000000
   0.1687144782D+00  0.1000000000D+01  0.1000000000D+01
请问，这个是什么原因呢？为什么C有这个现象，而氧原子的基组跟EMSL给出的基组一样呢？

谢谢各位了!!!!

作者
Author: sobereva 时间: 2014-12-6 18:22
先说C，你的输入文件大抵有问题，比如
#p 6-31+G* gfinput

Kein

0 1
C    0. 0. 0.

输出出来的是
S 6 1.00    0.000000000000
   0.3047524880D+04  0.1834737132D-02
   0.4573695180D+03  0.1403732281D-01
   0.1039486850D+03  0.6884262226D-01
   0.2921015530D+02  0.2321844432D+00
   0.9286662960D+01  0.4679413484D+00
   0.3163926960D+01  0.3623119853D+00
SP 3 1.00    0.000000000000
   0.7868272350D+01 -0.1193324198D+00  0.6899906659D-01
   0.1881288540D+01 -0.1608541517D+00  0.3164239610D+00
   0.5442492580D+00  0.1143456438D+01  0.7443082909D+00
SP 1 1.00    0.000000000000
   0.1687144782D+00  0.1000000000D+01  0.1000000000D+01
SP 1 1.00    0.000000000000
   0.4380000000D-01  0.1000000000D+01  0.1000000000D+01
D 1 1.00    0.000000000000
   0.8000000000D+00  0.1000000000D+01
你的情况多出来一个SP怀疑是你用了其它关键词额外加的。

作者
Author: sobereva 时间: 2014-12-6 18:32
至于Fe，你用6-31G* 输出一下，一比较就能看出，一个S和两个P都是因为加了弥散函数所致
注意这种加弥散函数的方式是不合理的，因为明显看出，外加的两个P和6-31G的已有的SP指数很接近，因此是无意义的。
本来6-31+G*对于Fe就是没有定义的，非要这么写，高斯虽不报错，但不代表加的弥散函数是有意义的。

作者
Author: chenxin199261 时间: 2014-12-6 18:40

sobereva 发表于 2014-12-6 18:22
先说C，你的输入文件大抵有问题，比如
#p 6-31+G* gfinput

是的，我拿单个的碳原子进行测试了，基组是正常的。
我在仔细的check一下我的输入文件，谢谢你了
我用的是混合基组，很乱的。

谢谢你了

作者
Author: chenxin199261 时间: 2014-12-6 18:43

sobereva 发表于 2014-12-6 18:32
至于Fe，你用6-31G* 输出一下，一比较就能看出，一个S和两个P都是因为加了弥散函数所致
注意这种加弥散函 ...

确实，我没有在数据库找到6-31+g（d）Fe的基组，可能这个基组没用涵盖铁，我在重新测试一下，找到合理一些的铁全电子基组，还有，我这么做会不会导致计算的错误？

作者
Author: chenxin199261 时间: 2014-12-6 19:23

chenxin199261 发表于 2014-12-6 18:40
是的，我拿单个的碳原子进行测试了，基组是正常的。
我在仔细的check一下我的输入文件，谢谢你了
我用 ...

我定义基组时候重叠定义了

作者
Author: sobereva 时间: 2014-12-6 23:38

chenxin199261 发表于 2014-12-6 18:43
确实，我没有在数据库找到6-31+g（d）Fe的基组，可能这个基组没用涵盖铁，我在重新测试一下，找到合理一 ...

不至于计算出错，但是比起6-31G*只会白多费时间，结果也没改进。
如果要考虑弥散函数，应该用已有弥散函数的基组，或者把6-31G*指数最小的SP壳层的指数除以2.5作为一个新的SP壳层加上。

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