24小时热门版块排行榜    

查看: 2612  |  回复: 5

brothers

银虫 (小有名气)

[交流] 想用multiwfn画电荷密度图,但是获取wfn文件失败 已有2人参与

在gaussian中运行以下gjf文件,往往写不出wfn文件,报错为2070错误,最后一行输出为writing a wfn file to.....
仔细核对了每一个参数,最后发现是基组的原因。
当使用较小的基组时可以正常写出wfn文件,但是使用cc-pvtz-pp则失败,请问这是为什么?怎么解决呢?
实在是不明白啊
ps:当不使用out=wfn关键词的时候是完全正常的。

#p B3LYP/genecp out=wfn nosymm

Au

1 1
Au 0.000000   0.000000    0.000000

Au     0
S   9   1.00
     36.0126000              0.0190290        
     22.5645000             -0.1678000        
     14.1620000              0.4966770        
      8.8548700             -0.2662560        
      5.0607200             -0.6998980        
      1.4740500              0.8731080        
      0.7049130              0.5025820        
      0.2177870              0.0321960        
      0.0974080             -0.0055030        
S   9   1.00
     36.0126000             -0.0067770        
     22.5645000              0.0600080        
     14.1620000             -0.1833280        
      8.8548700              0.1147190        
      5.0607200              0.2488220        
      1.4740500             -0.4376710        
      0.7049130             -0.3311080        
      0.2177870              0.2741810        
      0.0974080              0.6377380        
S   9   1.00
     36.0126000              0.0039090        
     22.5645000              0.0337600        
     14.1620000             -0.1660600        
      8.8548700             -0.0517890        
      5.0607200              0.7950230        
      1.4740500             -1.8059270        
      0.7049130              0.5612390        
      0.2177870              1.8784270        
      0.0974080             -0.9692290        
S   9   1.00
     36.0126000             -0.0119140        
     22.5645000              0.1330690        
     14.1620000             -0.4342870        
      8.8548700              0.1701570        
      5.0607200              1.0647560        
      1.4740500             -4.0055360        
      0.7049130              4.4028320        
      0.2177870             -1.4321630        
      0.0974080             -1.1432550        
S   1   1.00
      0.0412140              1.0000000        
P   8   1.00
     21.5954000             -0.0132180        
     13.5145000              0.0801960        
      6.2551400             -0.2966860        
      1.7181400              0.5289590        
      0.8694390              0.4788890        
      0.4267590              0.1576200        
      0.1842970              0.0138970        
      0.0763790              0.0008350        
P   8   1.00
     21.5954000              0.0030800        
     13.5145000             -0.0214060        
      6.2551400              0.0867760        
      1.7181400             -0.1897740        
      0.8694390             -0.1858350        
      0.4267590              0.0047800        
      0.1842970              0.3638360        
      0.0763790              0.5417640        
P   8   1.00
     21.5954000              0.0070860        
     13.5145000             -0.0474930        
      6.2551400              0.1918930        
      1.7181400             -0.4636370        
      0.8694390             -0.4351380        
      0.4267590              0.4944030        
      0.1842970              0.6751140        
      0.0763790              0.1315360        
P   8   1.00
     21.5954000              0.0129200        
     13.5145000             -0.0820290        
      6.2551400              0.3386190        
      1.7181400             -1.2549920        
      0.8694390              0.2258130        
      0.4267590              1.5758790        
      0.1842970             -0.8043540        
      0.0763790             -0.5073130        
P   1   1.00
      0.0312200              1.0000000        
D   7   1.00
     87.0208000              0.0000440        
     11.0668000              0.0142880        
      6.6604700             -0.0688760        
      1.9469600              0.2395590        
      0.9847060              0.3873230        
      0.4746550              0.3487560        
      0.2165580              0.1984920        
D   7   1.00
     87.0208000             -0.0000310        
     11.0668000             -0.0199590        
      6.6604700              0.0953200        
      1.9469600             -0.4213430        
      0.9847060             -0.5069330        
      0.4746550              0.2854660        
      0.2165580              0.5825650        
D   7   1.00
     87.0208000             -0.0000440        
     11.0668000              0.0350830        
      6.6604700             -0.1619420        
      1.9469600              1.0421550        
      0.9847060             -0.2648650        
      0.4746550             -1.1163050        
      0.2165580              0.5457620        
D   1   1.00
      0.0914720              1.0000000        
F   1   1.00
      1.4385000              1.0000000        
F   1   1.00
      0.4878000              1.0000000        
G   1   1.00
      1.2257000              1.0000000        
****

AU     0
AU-ECP     5     60
h-ul potential
  1
2      1.0000000              0.0000000        
s-ul potential
  2
2     13.5232180            426.6418670        
2      6.2643840             36.8006680        
p-ul potential
  4
2     11.4138670             87.0020910        
2     10.3292150            174.0043700        
2      5.7074240              8.8706100        
2      4.8281650             17.9024380        
d-ul potential
  4
2      7.4309630             49.8836550        
2      8.3219900             74.6845490        
2      4.6096420              6.4862270        
2      3.5115070              9.5468210        
f-ul potential
  2
2      3.0846390              8.7916400        
2      3.0247430             11.6584560        
g-ul potential
  2
2      3.9784420             -5.2343370        
2      4.0114910             -6.7381420        

E:\Multiwfn\Au00.wfn
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
gmy1990: 金币+2 2014-10-23 22:26:13
应确保E:\Multiwfn这个目录已经存在

实际上你也可以用.fch文件作为Multiwfn的输入,效果一样,生成.wfn文件不是必需的。

» 本帖已获得的红花(最新10朵)

2楼2014-10-23 22:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brothers

银虫 (小有名气)

引用回帖:
2楼: Originally posted by sobereva at 2014-10-23 22:09:52
应确保E:\Multiwfn这个目录已经存在

实际上你也可以用.fch文件作为Multiwfn的输入,效果一样,生成.wfn文件不是必需的。

真快!!谢谢你的回复,文件夹已经确认存在了。将cc-pvtz-pp基组直接替换成一个更小的基组,别的都不做修改,就成功了....
既然可以用fch来完成,那我以后不写wfn了吧,不过仍然想知道到底是为什么cc-pvtz-pp不可以呢?(我调了大半天才发现问题所在啊.....)
ps:正在使用您的multiwfn,真厉害!!已经第一次画出了自己的电荷密度差图(换了个小基组先试了一下),非常感谢~
3楼2014-10-23 22:21:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brothers

银虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by sobereva at 2014-10-23 22:09:52
应确保E:\Multiwfn这个目录已经存在

实际上你也可以用.fch文件作为Multiwfn的输入,效果一样,生成.wfn文件不是必需的。

用fch成功了,谢谢
4楼2014-10-23 22:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
wfn格式太老,最早的ampac只支持spdf函数。gh函数是最近几年别的aim程序自己加上的,不包含在gaussian中。因此如果出现了gh等高角动量的基函数,gaussian就会出错 。

» 本帖已获得的红花(最新10朵)

beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2014-10-24 00:23:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brothers

银虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by beefly at 2014-10-24 00:23:56
wfn格式太老,最早的ampac只支持spdf函数。gh函数是最近几年别的aim程序自己加上的,不包含在gaussian中。因此如果出现了gh等高角动量的基函数,gaussian就会出错 。

原来如此。wfn最高角动量只支持到f,在带有g角动量基函数的情况下若试图让Gaussian输出wfn文件就会报错。
非常感谢!

[ 发自手机版 http://muchong.com/3g ]
6楼2014-10-24 01:01:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 brothers 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见