24小时热门版块排行榜    

查看: 1687  |  回复: 5

huosuhong

铜虫 (初入文坛)

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

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

» 收录本帖的淘帖专辑推荐

Linux,高斯等相关知识

» 猜你喜欢

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

加油
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 21:14 by MHf40yA0Uhyb
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 20:50 by MHf40yA0Uhyb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-16 18:39 by rRvxEj3zorCy
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 QUjhNVAcOSff 2026-09-13 11/550 2026-09-16 18:37 by rRvxEj3zorCy
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 16:51 by i3KBRD19ERkG
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:00 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:24 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:27 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:27 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 03:02 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
请填处理意见