24小时热门版块排行榜    

查看: 705  |  回复: 4

lastzealot

新虫 (著名写手)

【乱世的奸雄】

[求助] 关于外部基组的问题

https://bse.pnl.gov/bse/portal网站上有很多基组,计算时常需要用到这些,但是不大清楚怎么输入进去。比如说Ce,找到Stuttgart RSC 1997 ECP基组的gaussian94版本,进去后如下:
!  Stuttgart RSC 1997 ECP  EMSL  Basis Set Exchange Library   9/20/14 9:00 AM
! ELEMENTS                      REFERENCES
! ---------                       ----------
!

****

CE     0
CE-ECP     5     28
h-ul potential
1
2      1.000000000            0.000000000      
s-ul potential
1
2     20.137829000          580.083457000      
p-ul potential
1
2     15.998482000          310.302833000      
d-ul potential
1
2     14.974187000          167.813944000      
f-ul potential
1
2     23.402455000          -49.390229000      
g-ul potential
1
2     16.570553000          -21.331879000  

有俩问题:
1. 这只是赝势部分,基组部分在哪里?用这个赝势,和其它基组的基组部分,可以吗?
2. 在input文件中输入这些基组时,直接copy在最下面就可以了吗?是否需要进行修改?p-ul中的ul是什么啊?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lastzealot

新虫 (著名写手)

【乱世的奸雄】

刚才可能表达的不够专业,各位虫友可能没大明白。
我的意思是input里面的外部基组应该有两个部分吧
Ce 0
...
****  (基组)

Ce 0
...       (赝势)
但是这两个部分的作用都分别是啥呢?
如果用网站上给出的EMSL Basis Set,是不是直接拷贝到input里面就可以了?还是要进行其他修改呢?
2楼2014-09-21 00:45:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lastzealot

新虫 (著名写手)

【乱世的奸雄】

EMSL Basis Set中有的基组只有赝势,那前面那部分(如下所示)去哪找啊?
Ce 0
...
****  (基组)
3楼2014-09-21 00:47:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴世杰

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
lastzealot: 金币+10, ★★★很有帮助, 谢谢,请问那个g-ul中的ul是指什么啊?在input文件中可以这么写吗? 2014-09-21 17:17:45
这是我在https://bse.pnl.gov/bse/portal网站上找到的
Ge     0
S   1   1.00
      1.0180000              1.0000000        
S   1   1.00
      0.3188000              1.0000000        
S   1   1.00
      0.1030000              1.0000000        
P   1   1.00
      1.8360000              1.0000000        
P   1   1.00
      0.2622000              1.0000000        
P   1   1.00
      0.0797600              1.0000000        
D   1   1.00
     44.8800000              1.0000000        
D   1   1.00
     12.2500000              1.0000000        
D   1   1.00
      3.8320000              1.0000000        
D   1   1.00
      1.1060000              1.0000000        
****
Ce     0
S   1   1.00
      0.9018000              1.0000000        
S   1   1.00
      0.3558000              1.0000000        
S   1   1.00
      0.2580000              1.0000000        
S   1   1.00
      0.0671900              1.0000000        
S   1   1.00
      0.0316900              1.0000000        
S   1   1.00
      0.0147100              1.0000000        
P   1   1.00
      1.8540000              1.0000000        
P   1   1.00
      0.2073000              1.0000000        
P   1   1.00
      0.1229000              1.0000000        
P   1   1.00
      0.0581400              1.0000000        
P   1   1.00
      0.0283300              1.0000000        
P   1   1.00
      0.0123400              1.0000000        
D   1   1.00
      3.6790000              1.0000000        
D   1   1.00
      0.5340000              1.0000000        
D   1   1.00
      0.2463000              1.0000000        
D   1   1.00
      0.1281000              1.0000000        
D   1   1.00
      0.0643300              1.0000000        
D   1   1.00
      0.0295600              1.0000000        
F   1   1.00
     54.7400000              1.0000000        
F   1   1.00
     18.2400000              1.0000000        
F   1   1.00
      6.9560000              1.0000000        
F   1   1.00
      2.7930000              1.0000000        
F   1   1.00
      1.0680000              1.0000000        
F   1   1.00
      0.3499000              1.0000000        
****

GE     0
GE-ECP     2     18
D POTENTIAL
  5
2      1.35459995            -1.40196705      
2      5.02229977            -9.45403576      
2     15.32830048           -23.62409592      
2     47.47309875           -62.44671631      
1    136.24209595           -12.96366024      
S-D POTENTIAL
  6
2      2.30509996           -27.65140533      
2      2.67100000           157.55535889      
2      3.49440002          -208.90888977      
2      4.66540003           146.75300598      
1     14.92539978            26.79065895      
0     44.04529953             2.94686103      
P-D POTENTIAL
  6
2      1.72140002           -55.52495575      
2      1.98930001           194.35308838      
2      2.55629992          -249.20588684      
2      3.29259992           154.47665405      
1      8.26930046            13.02904606      
0     10.14990044             5.37435579      
CE     0
CE-ECP     3     54
F POTENTIAL
  6
2      0.18850000            -0.31545401      
2      0.69250000            -4.03798485      
2      2.45420003           -20.49917410      
2      9.75559997           -81.90377040      
2     50.36529920          -294.86666900      
1    141.28579700           -30.30850030      
S-F POTENTIAL
  8
2      0.62040001            87.62731930      
2      0.75950003          -175.92492700      
2      1.04569995           173.75613400      
2      1.52219999          -118.99214900      
2      2.23819995            80.98940270      
2      8.08870029            -2.07819796      
1      6.39750004            36.81498720      
0     17.84670070             6.47752524      
P-F POTENTIAL
  8
2      0.48100001            73.16195680      
2      0.54939997           -98.24731440      
2      0.83929998            47.34008030      
2      1.63699997           -20.87126540      
2      2.63610005            44.86721040      
2      5.01919985            18.21543310      
1      8.89319992            37.83178710      
0     24.34210010             5.46868611      
D-F POTENTIAL
  8
2      0.08600000            -0.01027300      
2      0.84719998            -2.24907804      
2      0.17870000            -0.02766900      
2      1.88349998            98.89129640      
2      2.53220010          -233.11683700      
2      3.19779992           202.36198400      
1     10.85770040            45.45883560      
0     31.21439930             6.81029177      
第一部分是价电子的描述,第二部分是内部原子轨道势函数参数。
4楼2014-09-21 10:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+3 2014-09-21 15:00:12
lastzealot: 金币+20, ★★★★★最佳答案, 非常感谢您的帮助,为你特意增加悬赏。问一个弱弱的问题,在windows版本的gaussian中,可以自己制作.gbs文件,每次在input文件中@这个文件吗? 2014-09-21 17:19:09
见此贴最后一节
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://hi.baidu.com/sobereva/item/88107793764a87bdcc80e5dc
5楼2014-09-21 11:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lastzealot 的主题更新
信息提示
请填处理意见