24小时热门版块排行榜    

查看: 4398  |  回复: 7

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。该怎么解决?
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +3 yuleib84 2026-08-25 4/200 2026-08-25 13:31 by zjjcn2001
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 +3 ziyangfang 2026-08-25 3/150 2026-08-25 13:14 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 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见