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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » 急急急！高斯生成不了wfn文件怎么办
81/1返回列表
查看: 4406  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +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
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭