24小时热门版块排行榜    

查看: 895  |  回复: 5

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

[交流] 【求助】求助高斯中的Wachters-Hay基组 已有2人参与

从一篇文献上发现Wachters (augmented by two diffuse p and one diffuse d sets, 8s7p4d, full contraction sheme 62111111/3311111/3111)基组适合于第一过渡金属配合物的几何优化,作者对Wachters基组给出的参考文献如下:
JCP, 1970,52,1033-1036 和JCP, 1977,66,4377-4384。我想请教的是,它的基组应该如何表示了?麻烦您以Cu为例,以高斯格式给出,谢谢!

在EMSL基组库检索到的Wachters基组为:
****
Cu     0
S   6   1.00
337200.0000000              0.0002800        
  50072.9000000              0.0021700        
  11373.4000000              0.0112200        
   3239.8200000              0.0446210        
   1071.9700000              0.1383600        
    395.0990000              0.3125500        
S   2   1.00
    158.3990000              0.4106000        
     67.3591000              0.2145300        
S   1   1.00
     22.2983000              1.0000000        
S   1   1.00
      9.3935700              1.0000000        
S   1   1.00
      2.5784800              1.0000000        
S   1   1.00
      0.9640800              1.0000000        
S   1   1.00
      0.1133030              1.0000000        
S   1   1.00
      0.0407910              1.0000000        
P   3   1.00
   2245.2900000              0.0022600        
    532.1060000              0.0184100        
    172.1950000              0.0864600        
P   3   1.00
     65.3239000              0.2545000        
     27.0551000              0.4341900        
     11.7435000              0.3339000        
P   1   1.00
      4.6938200              1.0000000        
P   2   1.00
      1.9066700              0.5693600        
      0.7114450              0.2470200        
P   1   1.00
      0.2325980              1.0000000        
P   1   1.00
      0.0692990              1.0000000        
D   3   1.00
     48.5439000              0.0309700        
     13.5490000              0.1619600        
      4.5162800              0.3781200        
D   1   1.00
      1.4732900              1.0000000        
D   1   1.00
      0.4148750              1.0000000        
D   1   1.00
      0.1491000              1.0000000  

但似乎p轨道的contraction 方案似乎不对,所以请虫友多多指教。是不是我的理解有误。
回复此楼

» 猜你喜欢

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

前途光明,出路难觅!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fooo

金虫 (正式写手)


lihb734(金币+1):感谢您的参与!谢谢! 2010-05-23 21:01:49
aylayl08(金币+1):谢谢讨论 2010-05-24 08:14:15
你怎么看出来乎p轨道的contraction 方案不对了,请说清楚些。
另外我个人的感觉,这个基组似乎太老了,现在应该有更好的基组,为什么不用而非要用这个?
引用回帖:
Originally posted by lihb734 at 2010-05-23 19:04:39:
从一篇文献上发现Wachters (augmented by two diffuse p and one diffuse d sets, 8s7p4d, full contraction sheme 62111111/3311111/3111)基组适合于第一过渡金属配合物的几何优化,作者对Wachters基组给出的参考 ...

2楼2010-05-23 19:26:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

引用回帖:
Originally posted by fooo at 2010-05-23 19:26:46:
你怎么看出来乎p轨道的contraction 方案不对了,请说清楚些。
另外我个人的感觉,这个基组似乎太老了,现在应该有更好的基组,为什么不用而非要用这个?


你怎么看出来乎p轨道的contraction 方案不对了——作者的是33收缩,而EMSL给出的是32收缩,我是这么理解的,不一定对啊。
这个基组确实有点老了,但有人证明它要优于常用的LANL2DZ, SDD等含ECP的基组,您有什么好的基组推荐码?
前途光明,出路难觅!
3楼2010-05-23 21:01:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fooo

金虫 (正式写手)


aylayl08(金币+1):谢谢讨论 2010-05-24 08:14:30
P   3   1.00
   2245.2900000              0.0022600        
    532.1060000              0.0184100        
    172.1950000              0.0864600        
P   3   1.00
     65.3239000              0.2545000        
     27.0551000              0.4341900        
     11.7435000              0.3339000        
这两个不都是33收缩吗?怎么是32收缩呢?
这些年出了不少新基组,在EMSL上也有不少,这方面的文献检索起来比较麻烦,似乎发表在JCP上的比较多。
引用回帖:
Originally posted by lihb734 at 2010-05-23 21:01:33:

你怎么看出来乎p轨道的contraction 方案不对了——作者的是33收缩,而EMSL给出的是32收缩,我是这么理解的,不一定对啊。
这个基组确实有点老了,但有人证明它要优于常用的LANL2DZ, SDD等含ECP的基组,您有什么 ...

4楼2010-05-23 21:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

引用回帖:
Originally posted by fooo at 2010-05-23 21:07:47:
P   3   1.00
   2245.2900000              0.0022600        
    532.1060000              0.0184100        
    172.1950000              0.0864600        
P   3   1.00
     65.3239000           ...

P   3   1.00
   2245.2900000              0.0022600        
    532.1060000              0.0184100        
    172.1950000              0.0864600        
P   3   1.00
     65.3239000              0.2545000        
     27.0551000              0.4341900        
     11.7435000              0.3339000        
P   1   1.00
      4.6938200              1.0000000        
P   2   1.00
      1.9066700              0.5693600        
      0.7114450              0.2470200        
P   1   1.00
      0.2325980              1.0000000        
P   1   1.00
      0.0692990              1.0

不好意思,是3312,而文献中给出的是3311,谢谢!对不起!给您添麻烦了。
前途光明,出路难觅!
5楼2010-05-23 21:47:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fooo

金虫 (正式写手)


lihb734(金币+3):感谢您的指点!说实话,我现在还不知道什么是contraction了。就算看了原始文献也搞不明白。 2010-05-24 07:50:44
aylayl08(金币+1):谢谢讨论 2010-05-24 08:14:54
这可能是有问题,最好查一下原始文献。不过如果参数没错,也可能关系不大,因为你不采用收缩也可以。
引用回帖:
Originally posted by lihb734 at 2010-05-23 21:47:30:

P   3   1.00
   2245.2900000              0.0022600        
    532.1060000              0.0184100        
    172.1950000              0.0864600        
P   3   1.00
     65.3239000         ...

6楼2010-05-24 03:38:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lihb734 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见