24小时热门版块排行榜    

查看: 1381  |  回复: 5

huosuhong

铜虫 (初入文坛)

[求助] 重金属元素基组问题

各位虫友大家好,请问在做拓扑分析用gaussian生成wfn时重金属元素用哪个基组比较好,比如Pb元素,我采用了aug-cc-pvdz-pp基组,据老师说这个基组考虑了部分内层电子,但是在算Pb时提示基函数个数小于轨道数目导致错误,请问各位有没有好的建议,谢谢!
回复此楼
加油
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangf44

禁虫 (小有名气)

★ ★ ★ ★ ★
感谢参与,应助指数 +1
huosuhong: 金币+5 2012-09-14 08:31:25
本帖内容被屏蔽

2楼2012-09-13 22:25:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
huosuhong: 金币+5 2012-09-14 08:31:40
gmy1990: 金币+3 2012-09-14 18:40:42
对于主族元素,如果在赝势下进行波函数分析,(aug)-cc-pVDZ-PP是最为推荐的。它是极少几个对于主族是小核的赝势,能够将亚价层密度表现出来,这是很重要的。可参见本文的讨论http://hi.baidu.com/sobereva/item/ed2f57e743c1be11595dd8c8

至于错误提示,如LS所说,应当是写法有误。
3楼2012-09-13 23:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huosuhong

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by sobereva at 2012-09-13 23:23:46
对于主族元素,如果在赝势下进行波函数分析,(aug)-cc-pVDZ-PP是最为推荐的。它是极少几个对于主族是小核的赝势,能够将亚价层密度表现出来,这是很重要的。可参见本文的讨论http://hi.baidu.com/sobereva/item/ed2 ...

非常感谢大侠相助,我一直在使用您的multiwfn程序,非常敬佩和仰慕您的才学和品性。我现在将我的gjf文件和out文件贴出来,能否帮我看看到底是什么问题吗,非常感谢,我们老师说是基组问题,所以建议我向卢大侠您请教一下
Pb.out(6KB)
http://kuai.xunlei.com/d/GBSYLEIWSVBO?p=130497
Pb.gjf(2.12KB)
http://kuai.xunlei.com/d/GAPAWTSZHKIP?p=130497
加油
4楼2012-09-14 08:29:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huosuhong

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by wangf44 at 2012-09-13 22:25:39
这个基组是含赝势的基组,你要把赝势加上去啊!

谢谢虫友的帮助,请问怎么加上去啊,能否说的详细点,谢谢
加油
5楼2012-09-14 08:30:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★
gmy1990: 金币+2 2012-09-14 18:40:51
写成这样。赝势信息是从EMSL里aug-cc-pVDZ-PP的条目里取的。
# b3lyp/genecp  output=wfx nosymm

Title Card Required

0 1
Pb                 1.59400000    0.00000000    0.00000000
Pb                -1.59400000    0.00000000    0.00000000
H                  1.80941937    2.89198798   -0.00000000
H                 -1.80941937   -2.89198798   -0.00000000

H 0
6-311++g(d,p)
****
Pb     0
S   6   1.00
     37.4583000              0.0248240        
     23.3170000             -0.1810620        
     14.5759000              0.5007060        
      6.4337600             -0.9797370        
      1.6671900              0.9666300        
      0.8077320              0.4364870        
S   6   1.00
     37.4583000              0.0067850        
     23.3170000             -0.0601870        
     14.5759000              0.1889710        
      6.4337600             -0.4248380        
      1.6671900              0.6220350        
      0.8077320              0.3160880        
S   1   1.00
      0.2270840              1.0000000        
S   1   1.00
      0.0842520              1.0000000        
S   1   1.00
      0.0278000              1.0000000        
P   5   1.00
     10.6186000              0.1701740        
      7.1878100             -0.4241790        
      1.8943900              0.6934430        
      0.8167090              0.4737730        
      0.2036880              0.0290470        
P   5   1.00
     10.6186000             -0.0479540        
      7.1878100              0.1245560        
      1.8943900             -0.2597760        
      0.8167090             -0.1414800        
      0.2036880              0.5317690        
P   1   1.00
      0.0649730              1.0000000        
P   1   1.00
      0.0197000              1.0000000        
D   5   1.00
     13.7266000              0.0107870        
      6.9898300             -0.0666500        
      2.1735200              0.3342590        
      1.0281200              0.4826840        
      0.4558570              0.3004140        
D   1   1.00
      0.1822000              1.0000000        
D   1   1.00
      0.0609000              1.0000000        
****

PB     0
PB-ECP     5     60
h-ul potential
  1
2      1.0000000              0.0000000        
s-ul potential
  2
2     12.2963030            281.2854990        
2      8.6326340             62.5202170        
p-ul potential
  4
2     10.2417900             72.2768970        
2      8.9241760            144.5910830        
2      6.5813420              4.7586930        
2      6.2554030              9.9406210        
d-ul potential
  4
2      7.7543360             35.8485070        
2      7.7202810             53.7243420        
2      4.9702640             10.1152560        
2      4.5637890             14.8337310        
f-ul potential
  2
2      3.8875120             12.2098920        
2      3.8119630             16.1902910        
g-ul potential
  2
2      5.6915770             -9.0966650        
2      5.7155670            -11.5319960      

./pb.wfx
6楼2012-09-14 16:43:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huosuhong 的主题更新
信息提示
请填处理意见