Fe，C原子的6-31+g(d)基组形式的问题

chenxin199261

我在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给出的基组一样呢？





谢谢各位了!!!!

sobereva

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

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

chenxin199261

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

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

chenxin199261

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

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

chenxin199261

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

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

谢谢你了

sobereva

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

sobereva

先说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怀疑是你用了其它关键词额外加的。