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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » 急急急！高斯生成不了wfn文件怎么办
81/1返回列表
查看: 4408  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zyr3365754

木虫 (著名写手)

[求助] 急急急！高斯生成不了wfn文件怎么办

想用Multiwfn作电子密度差图，但是Multiwfn需要读入wfn文件，
所以在高斯中用output=wfn
命令想得到wfn文件，最后原子坐标部分完了以后也加了路径和名字
F:\GaussOutput\XX.wfn
但是计算时没有生成wfn文件，我使用的ecp。

#opt freq=noraman tpsstpss/genecp output=wfn
算的Ag原子和Si原子，其中Ag原子用的Lanl2Tz(f)

我在坛子上看到一个帖子说output=wfn，不能得到ecp的波函数信息，
我的问题是
1.我得不到wfn文件是因为ecp的原因吗？
2.如果不使用ecp就不能得到wfn文件，那么我想使用Multiwfn。该怎么解决？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-13 11:07:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+2 2012-09-13 17:32:36
zyr3365754: 金币+5, 感谢交流， 2012-09-14 10:12:06
没看到你的输入输出文件 所以无法判断你的问题所在
使用了ecp当然也可以得到wfn文件
Multiwfn可以不使用wfn文件而使用fchk文件
赞一下 回复此楼
2楼2012-09-13 12:29:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

【答案】应助回帖

★ ★ ★ ★ ★
zyr3365754: 金币+5, 谢谢 2012-09-14 10:12:22
其实你这种情况 最好改用wfx文件 详见Multiwfn的说明书
赞一下 回复此楼
3楼2012-09-13 12:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zyr3365754: 金币+10, 我明天把我的输入文件贴出来你看看。非常感谢 2012-09-14 10:13:00
不是因为ECP的原因。
没能生成wfn文件我估计是因为你在定义基组、赝势以及wfn文件路径的时候格式弄乱了。

以下是个例子以资参考
%chk=C:\gtest\CuCO+.chk
# b3lyp/genecp out=wfn

b3lyp/cc-pvdz+lanl2dz

1 1
cu                -2.92658731    0.85317459    0.00000000
C                 -0.71658731    0.85317459    0.00000000
O                  0.39881269    0.85317459    0.00000000

cu 0
lanl2dz
****
C O 0
cc-pvdz
****

cu 0
lanl2dz

c:\t.wfn
赞一下 回复此楼
4楼2012-09-13 23:30:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

书万里

木虫 (著名写手)

历史控+原著控
哇，全是大牛回的贴，楼主很幸运哦。
赞一下 回复此楼
喂人民服雾。
5楼2012-09-14 08:12:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyr3365754

木虫 (著名写手)
引用回帖:
4楼: Originally posted by sobereva at 2012-09-13 23:30:49
不是因为ECP的原因。
没能生成wfn文件我估计是因为你在定义基组、赝势以及wfn文件路径的时候格式弄乱了。

以下是个例子以资参考
%chk=C:\gtest\CuCO+.chk
# b3lyp/genecp out=wfn

b3lyp/cc-pvdz+lanl2dz
...

# tpsstpss/genecp nosymm geom=connectivity output=wfn

Title Card Required

0 1
Ag                 0.00000000    1.48897900   -0.40098300
Ag                 0.00000000   -1.48897900   -0.40098300
Si                -1.14673100    0.00000000    1.34615900
Si                 1.14673100    0.00000000    1.34615900

1 3 1.0 4 1.0
2 3 1.0 4 1.0
3 4 1.0
4

Si  0
6-311+g(d)
****
Ag     0
S   1   1.00
      2.9500000              1.0000000        
S   1   1.00
      2.1490000              1.0000000        
S   1   1.00
      0.6684000              1.0000000        
S   1   1.00
      0.0997000              1.0000000        
S   1   1.00
      0.0347000              1.0000000        
P   1   1.00
      6.5530000              1.0000000        
P   1   1.00
      1.5650000              1.0000000        
P   1   1.00
      0.5748000              1.0000000        
P   1   1.00
      0.0833000              1.0000000        
P   1   1.00
      0.0252000              1.0000000        
D   2   1.00
      3.3910000              0.1396938        
      1.5990000              0.4744421        
D   1   1.00
      0.6282000              1.0000000        
D   1   1.00
      0.2108000              1.0000000        
F   1   1.00
      1.6110000              1.0000000        
****

AG     0
AG-ECP     3     28
f potential
  5
0    568.7006237             -0.0587930        
1    162.3579066            -20.1145146        
2     51.1025755           -104.2733114        
2     16.9205822            -40.4539787        
2      6.1669596             -3.4420009        
s-f potential
  5
0     76.0974658              2.9861527        
1     15.3327359             35.1576460        
2     18.7715345            450.1809906        
2     13.3663294           -866.0248308        
2      9.8236948            523.1110176        
p-f potential
  5
0     56.3318043              4.9640671        
1     69.0609098             21.5028219        
2     19.2717998            546.0275453        
2     12.5770654           -600.3822556        
2      8.7956670            348.2949289        
d-f potential
  5
0     53.4641078              3.0467486        
1     40.1975457             23.3656705        
2     11.9086073            777.2540117        
2      9.7528183          -1238.8602423        
2      8.1788997            608.0677121        

****

Ag  0
lanl2dTz(f)

F:\SiAg.wfn

以上是我的输入文件，Ag我用的LanLTz(f)基组。Si用的6-311+g(d)。以上导入高斯计算，只有out文件输出，out文件最后提示是：
1|1|UNPC-UNK|SP|RTPSSTPSS|GenECP|Ag2Si2|PCUSER|15-Sep-2012|0||# tpsstp
ss/genecp nosymm geom=connectivity output=wfn||Title Card Required||0,
1|Ag,0,0.,1.488979,-0.400983|Ag,0,0.,-1.488979,-0.400983|Si,0,-1.14673
1,0.,1.346159|Si,0,1.146731,0.,1.346159||Version=IA32W-G03RevE.01|HF=-
870.3878491|RMSD=2.384e-005|Thermal=0.|Dipole=0.,0.,-1.241521|PG=C02V
[SGV(Ag2),SGV'(Si2)]||@

Writing a WFN file to ****.
我自己就搞不清楚Writing a WFN file to ****.是什么意思？？
恳请大侠帮助
赞一下 回复此楼
6楼2012-09-15 15:32:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianwenk

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
zyr3365754: 金币+5, ★★★很有帮助, 非常感谢 2012-09-16 12:17:23
引用回帖:
6楼: Originally posted by zyr3365754 at 2012-09-15 15:32:58
# tpsstpss/genecp nosymm geom=connectivity output=wfn

Title Card Required

0 1
Ag                 0.00000000    1.48897900   -0.40098300
Ag                 0.00000000   -1.48897900   -0.400 ...

我最近也出现了这个问题，我是因为wFN写入的路径与上面空行太多了，后来发现只能空一行，就是说“F:\SiAg.wfn”与上面只能空一行，你可以试一下
赞一下 回复此楼
7楼2012-09-15 16:11:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫

【答案】应助回帖

★ ★ ★ ★ ★
zyr3365754: 金币+5, ★★★★★最佳答案, 非常感谢，以后多指导呀！再次感谢 2012-09-16 12:17:11
****

Ag  0
lanl2dTz(f)

F:\SiAg.wfn
将以上六行中的前五行删掉，只保留最后一行
高斯的输入文件对空行数目、空行位置是很敏感的，不要乱空行。对星号也是。
由于你已经以自定义方式定义了赝势，就不要再用诸如这样
Ag  0
lanl2dTz(f)
来再次定义赝势，这就重复了，而且也不会生效。而且高斯自己并不支持lanl2dTz(f)这种写法。高斯支持的赝势在手册里都已写明。
赞一下 回复此楼
8楼2012-09-15 18:48:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zyr3365754 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 17:31 by 4FFAWE8HcgUD
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 17:30 by 4FFAWE8HcgUD
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 5/250 2026-08-29 16:16 by manplx
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:02 by jCd0dEvKHShX
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 13:51 by jCd0dEvKHShX
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
关闭
请填处理意见
关闭