24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1482  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试324 中药学 一志愿天中医 求调剂 +3 李的Lucky 2026-04-10 3/150 2026-04-15 23:54 by 陈皮皮
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 材料299专硕求调剂 +17 +21 2026-04-09 17/850 2026-04-15 17:20 by zs92450
[考研] 310求调剂 +16 666真好 2026-04-11 18/900 2026-04-15 13:28 by 黑科技矿业
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 机械还有还有名额吗?太难了 +8 笑笑袁 2026-04-10 8/400 2026-04-14 08:44 by screening
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 一志愿矿大,材料工程专硕314分,0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
信息提示
请填处理意见