24小时热门版块排行榜    

查看: 1010  |  回复: 2

cug_zhang

金虫 (正式写手)

[求助] 高斯计算td总出错 已有2人参与

输入文件如下:
%chk=1011.chk
%nprocs=4
%mem=3GB
# rb3lyp/genecp td=(nstates=100)

Title Card Required

0 1
Si                 0.00000000    0.00000000    1.54206500
Pd                 0.00000000    1.36878700   -0.23466200
Pd                 0.00000000   -1.36878700   -0.23466200

1
2
3

Si  0
6-311+g(d)
****
Pd     0
S   3   1.00
      2.7870000             -1.6102393        
      1.9650000              1.8489842        
      0.6243000              0.6037492        
S   4   1.00
      2.7870000              1.3540775        
      1.9650000             -1.6780848        
      0.6243000             -0.8559381        
      0.1496000              1.0200299        
S   1   1.00
      0.0436000              1.0000000        
P   3   1.00
      5.9990000             -0.1034910        
      1.4430000              0.7456952        
      0.5264000              0.3656494        
P   2   1.00
      0.7368000              0.0763285        
      0.0899000              0.9740065        
P   1   1.00
      0.0262000              1.0000000        
D   3   1.00
      6.0910000              0.0376146        
      1.7190000              0.5200479        
      0.6056000              0.5706071        
D   1   1.00
      0.1883000              1.0000000
****

PD     0
PD-ECP     3     28
f-ul potential
  5
0    598.3336444             -0.0563177        
1    162.4298290            -20.1288036        
2     51.5714771           -105.8197923        
2     16.4888260            -42.5733345        
2      5.8287656             -3.6165086        
s-ul potential
  5
0     73.3806304              3.0003651        
1     14.7550438             32.4350093        
2     17.8350204            459.0830383        
2     12.7111477           -868.0629029        
2      9.3292063            514.4726098        
p-ul potential
  5
0     55.6689877              4.9593099        
1     64.2337771             21.1711029        
2     17.6254952            605.0560092        
2     11.9058155           -726.9641846        
2      8.5100832            396.3274883        
d-ul potential
  5
0     49.9994728              3.0508745        
1     39.7477547             22.2506580        
2     11.4321366            674.8357698        
2      9.1790080          -1040.8554048        
2      7.5624429            505.9375147        

****

Pd  0
lanl2dz ecp

错误如下:
QPErr --- A syntax error was detected in the input line.
2
  '
Last state="NG"
TCursr= 1629 LCursr=    1
Error reading general basis specification.
Error termination via Lnk1e in /public/home/zhangshuai/g09/l301.exe at Sun Sep 27 17:52:40 2015.



输入文件修改:
%chk=1011.chk
%nprocs=4
%mem=3GB
# rb3lyp/genecp  td=(nstates=100)

1 Energy

0 1
Si                 0.00000000    0.00000000    1.54206500
Pd                 0.00000000    1.36878700   -0.23466200
Pd                 0.00000000   -1.36878700   -0.23466200

Si 0
6-311+g(d)
****
Pd 0
lanl2dz
****

Pd 0
lanl2dz

错误信息如下:
QPErr --- A syntax error was detected in the input line.
2
  '
Last state="NG"
TCursr= 1629 LCursr=    1
Error reading general basis specification.
Error termination via Lnk1e in /public/home/zhangshuai/g09/l301.exe at Sun Sep 27 17:52:40 2015.


同样一个输入文件:
%chk=1011.chk
%mem=3GB
%nprocshared=4
# b3lyp/genecp td=(nstates=100)

Title Card Required

0 2
Si                 0.00000000    1.26917800   -1.15457600
Si                 0.00000000   -1.26917800   -1.15457600
Nb                 0.00000000    0.00000000    0.78849100

Si 0
6-311+g(d)
****
Nb 0
lanl2dz
****

Nb 0
lanl2dz

这个文件就可以正常计算。


请高手帮忙看看错误发生在哪里?谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oinkmasta

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
td计算出的excited states基本都是open shell啊,rb3lyp应该不行。但它说是syntax error我也没想到。
16届本科在北美
2楼2015-09-27 23:44:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
cug_zhang: 金币+20, ★★★★★最佳答案, 答案很好,谢谢帮助! 2015-09-28 09:17:51
你的问题描述得不清楚,给人以很大误解,希望以后能加强。
首先我测试了你的第二个输入文件,没有问题,于是我想是不是你的输入法的关系,换了中文输入法输入lanl2dz,仍然没有问题。然后我再回过来看你的第一个出错信息和第二个出错信息,发现它们的内容和时间节点是完全一样的,也就是说你给的还是第一个问题的出错信息。综上所述,第二个输入文件是可以运行的。
回到第一个出入文件,确实有问题,问题出在输入语法上,数字2周围。回去找2,发现你的分子坐标下有1,2,3一组数,但是你的命令行内没有geom=connectivity,于是删除1,2,3,再次运行,一切正常。综上所述,geom=connectivity与对应的描述部分要么一起删除,要么一起保留。问题解决。
3楼2015-09-28 08:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cug_zhang 的主题更新
信息提示
请填处理意见