24小时热门版块排行榜    

查看: 1663  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
信息提示
请填处理意见