24小时热门版块排行榜    

查看: 3917  |  回复: 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 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 化学工程321分求调剂 +9 大米饭! 2026-03-15 12/600 2026-03-16 21:59 by 大米饭!
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703化学调剂 ,六级已过,有科研经历 +7 曦熙兮 2026-03-15 7/350 2026-03-16 16:34 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 求调剂,一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见