24小时热门版块排行榜    

查看: 3919  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +14 于海海海海 2026-03-19 14/700 2026-03-19 19:40 by maocaozhuxi
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 274求调剂 +6 S.H1 2026-03-18 6/300 2026-03-19 09:34 by 花店相见
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
信息提示
请填处理意见