应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » 想用multiwfn画电荷密度图，但是获取wfn文件失败
61/1返回列表
查看: 2078  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-10-23 21:55:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

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

实际上你也可以用.fch文件作为Multiwfn的输入，效果一样，生成.wfn文件不是必需的。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

brothers
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朵）

brothers
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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭